m = 3
n = 3
h = .6
d = 1
c1 = hn(h,n)
c2 = h^(n-1)
%c2 = h^(n)

%c1 = 0;
%c2 = 1;

%c1 = rand
%c2 = rand
global GG
good=1;
mf = 0;
while good
  %c1 = rand;
  %c2 = rand;
  %h = rand
  %d = rand
  %d=1, h=.1
  %[A,gam,fval,good] = gamcheck(m,h,c1,c2);
  %n = round(randinseg(3,4))
  %[A,fval,diamX] = gamcheck1(m,n,h)
  %[A,phi,fval,diamX] = gamcheck2(m,n,h)
  %Bound = sum(h.^((1:n)-1)) + h^(n-1)*diamX
  %[A,phi,fval] = gamcheck01(m,n,h,d)
  %[A,fval] = gamrmlp2(GG,.1)

  %[A,phi,fval] = pgamcheck(m,n,h,d)
  %if abs(fval) > 1.11+1e-12
  %  'bad'
  %  keyboard
  %end

  %[A,phi,fval] = pgamcheck(m,n,h,d);
  %af = abs(fval);
  %if af>mf
  %  mf = af;
  %  %if mf>.099
  %  %  keyboard
  %  %end
  %  mA = A;
  %  mphi = phi;
  %end
  %global GS1 GS2
  %fprintf('mf = %1.12f \n',mf);
   
  %c1 = rand
  %c2 = rand
  %h = rand
  %[A,gam,fval] = alfbetcheck(c1,c2,h);
  %g = (abs(fval) <= c1 + c2*h + 1e-10);
  %if ~g
  %  'BAAAAAAAAAAD'
  %  keyboard
  %end
  
  %[A,gam,fval] = gam3dfc(c1,c2,h,m);  
  %[A,gam,fval] = gam3ind(c1,c2,h,m);  
  %[A,gam,fval] = gam2ind(c1,c2,h);
  %if fval>mf
  %[A,gam,fval] = maxaff(c1,c2,m);
  %[A,gam,fval] = maxaff3(c1,c2,m);
  %[A,gam,fval] = gam2m(c1,c2,h,m);
  %if fval<mf
  %  mf = fval;
  %  %mf = mf/c2;
  %  mA = A;
  %  mgam = gam;
  %  cc = [c1 c2 h];
  %end
  %fprintf('mf = %1.15f \n',mf);

  if 0
    global cc;
    A0=pnormdim(rand(m,m,n-1),1);
    philp = optmasklp(A0);
    %global GG
    %fvalp = tryphi(A0,philp)
    fvalp = philp * cc
    philp = philp-min(philp);
    phe = extremizf1(philp,m,n);
    phe = round(phe-min(phe));
    %fvale = tryphi(A0,phe)
    fvale = phe * cc
    
    %mask = getgmask(phe,m,n);
    mask = getmask1(phe,m,n);
    phim = fillgmask(mask,m);
    %fvalm = tryphi(A0,phim)
    fvalm = phim * cc
    
    
    %keyboard
    %x = phe(:);
    %b = inconvhull(x,PHI1);
    %if ~b
    %  PHI1 = [PHI1 x];
    %  corners = size(PHI1,2)
    %end
    
    %fv = abs(fvalp-fvalm);
    %fv = totdif(phim,phe);
    ii = getphinds(m,n);

    fv = totdif(phim(ii),phe(ii));
    if fv > .01
      keyboard
    end
    %if fv>mf
    if 0
      mf = fv;
      mA = A0;
      mphi = philp;
      save BADDAT mf mA mphi
    end
    fprintf('mf = %1.15f \n',mf);
  end
  
  
  %function xx = randxx(m,nn)
  %for k=1:m
  %  xx{k} = ceil(m*rand(1,nn));
  %end

  if 0
    n = 3;
    m = 3;

    xx0 = {};
    xx1 = {};
    for k=1:m
      xx0{k} = [];
      xx1{k} = [k];
    end
    
    GGk = {};
    for k=1:m
      xxk = {};
      for l=1:m
        xxk{l} = [k];
      end
      
      GGk{k} = fullgamx(n,A0,xxk,mask);
      
    end

    mask = randmask(m,n+2);
    phi = fillgmask(mask,m);
    h = rand
    [A0,fval] = phifixed(phi,m,n,h,d);
    
    %k1 = 1;
    %k2 = 1;
    %for i=1:n-1
    %  h = strict(A0(:,:,i));
    %  k1 = k1 + k2*h;
    %  k2 = h;
    %end
    k1 = sum(h.^(0:n));
    k2 = h;
    
    
    %A0=pnormdim(rand(m,m,100),1);
    
    %GG0 = fullgamx(n-1,A0,xx0,mask);
    GG = fullgamx(n,A0,xx0,mask);    
    b0 = isbimodal(GG);
    b1 = all(abs(sum(GG'))<=k1+1e-10);
    
    s1 = abs(sum(max(GG)));
    s2 = abs(sum(min(GG)));
    
    %b2 = abs(sum(diag(GG)))<=k1+k2+1e-10;
    %b2 = max(s1,s2)<=k1+1e-10;
    b2 = max(s1,s2)<=k1+k2+1e-10;
    
    good = (b0 & b1 & b2)
    if ~good
      keyboard
    end
  end
  
  
  if 1
    m = 4;
    M = anybeta(m);
    k1 = rand;
    k2 = rand;
    h = rand;
    [fval,A,alf0,bet0] = abgcheck(m,h,k1,k2);
    
    bet1 = A * bet0;
    alf1 = A * alf0;

    
    g1 = (max(abs(M*bet1)) <= k2*h + 1e-8);
    g2 = (abs(sum(alf1)) <= k1 + 1e-8);
    
    good = (g1 & g2)
    if ~good
      keyboard
    end
  end
  
  
  if 0
    %cc = randn(m^n,1);
    %philp = optphi(cc,m,n);
    A0=pnormdim(rand(m,m,n-1),1);
    philp = optmasklp(A0);
    philp = philp-min(philp);
    extremizfr(philp,m,n);
    
    
    %if abs(fvalp-fvale) > 1e-8
    %if abs(fvalp-fvalm) > 1e-7
    %fv = abs(fvalp-fvalm);
    %if fv>mf
    %  mf = fv;
    %  mA = A0;
    %  mphi = philp;
    %  save BADDAT mf mA mphi
    %end
    %fprintf('mf = %1.15f \n',mf);
  end
    
  %B = d*(1+h+h^2)
  %good = good & ((fval <= B + 1e-8));
  %D = max(max(GG)-min(GG))
  %good = good & (
  %E = B - fval
  %D + E
  %keyboard
end

%global XXX
%save BADGAM c1 c2 h A gam fval XXX

%global XX FN
%save BADGAM8996 h A FN fval XX