function [AA1,C] = smartcollapse2(AA,C)
N = size(AA,1);
T = size(AA,3);
[h0,i0,j0] = SC(AA);
for t=1:T
  %[h0,i0,j0] = SC(AA(:,:,t));
  %pi0 = AA(:,i0,t);
  %pj0 = AA(:,j0,t);
  %[x,k] = max(abs(pi0-pj0));
  %a = pi0(k);
  a = h0;
  %b = pj0(k);
  b = 0;
  %A = [[a,b];[1-a,1-b]];
  AA1(:,:,t) = [[a,b];[1-a,1-b]];
end
%C = pnormdim(rand(2,1),1);
C(1) = 0;
C(2) = 1;
return


function [h,i0,j0] = SC(AA)
% the "stricture coefficient" of AA
[n,n,T] = size(AA);
h = 0;
for t=1:T
  A = AA(:,:,t);
  for i=1:n
    for j=i+1:n
      x = A(:,i);
      y = A(:,j);
      hxy = max(abs(x-y));
      if (hxy>h)
        i0=i;j0=j;h=hxy;
      end
    end
  end
end
return