curmin = 10
global m n
m=3;
save tmpdat A curmin

global HH PP Ainf Binf delta Lmf Lpd
global delta niter
delta = 1/m
rho = 1-1/(m-1) + .01
%rho = .5+(1-delta)/2;
niter = 1

global domask phi01
domask = 0; % mask doesn't yield opt
phi01 = 0; % restrict phi to [0,1]

global myopt
myopt = optimset('Diagnostics','off');
myopt = optimset('Display','off');
warning off


for n=2:6
  if domask
    Lmf = getLmf(m,n);
  elseif ~phi01    
    [Ainf,Binf] = lipab(m,n);
  end
  curmin = 10
  for t=1:20
    [A,p] = minmaxprdev(m,n,rho);
    if p<curmin
      curmin = p;
      Amin = A;
      fprintf('m=%d n=%d iter = %d  curmin = %1.9f \n',m,n,t,curmin);
      save tmpdat Amin curmin CTRX
    end
    if curmin < 1e-7
      break
    end
  end
  CTRX{n}=Amin;
end