util1_jan/adiff.m

   1 function [d,varargout]=adiff(root1,tiles1,root2,tiles2,steps)
   2 % d=adiff(root1,tiles1,root2,tiles2,steps)
   3 n=length(steps);
   4 for i=1:n
   5         step=steps(i);
   6         a1=read_array_tiles(root1,tiles1,step);
   7         a2=read_array_tiles(root2,tiles2,step);
   8         d(i)=big(a1-a2)/max([big(a1),big(a2),realmin]);
   9         disp(d(i))
  10 end
  11 if nargout==3,
  12     varargout{1}=a1;
  13     varargout{2}=a2;
  14 end
  15 end
  16