femwind/read_array.m

   1 function a=read_array(name)
   2 % read matrix a from a given file name.txt
   3 filename=[name,'.txt'];
   4 if ~exist(filename,'file'),
   5     a=[];
   6     return
   7 end
   8 X=load(filename);
   9 magicnumber=X(1,:);
  10 if magicnumber~=456
  11     error(['cannot load matrix in file ',filename,' bad format'])
  12 end
  13 [m,n]=size(X);
  14 if n~=1,
  15     error('file must contain column vector')
  16 end
  17 nd=X(2);
  18 n1=X(3);
  19 n2=X(4);
  20 n3=X(5);
  21 fprintf('reading matrix size %g %g %g from file %s\n',n1,n2,n3,filename)
  22 a=reshape(X(6:m),n1,n2,n3);