read_array.m

来自「用于计算矩阵的特征值,及矩阵的其他运算.可以用与稀疏矩阵」· M 代码 · 共 32 行

function A = read_array(filename);

fid = fopen(filename,'r');
dim = fread(fid,1,'uint');
for i = 1:dim
  d(i) = fread(fid,1,'uint');
end;
[B,count] = fread(fid,prod(d),'float');
B = reshape(B,fliplr(d));

dim = fread(fid,1,'uint');
if (~feof(fid)) %collection of arrays.
  A{1} = B;
  ct = 2;
  while (~feof(fid))
    for i = 1:dim
      d(i) = fread(fid,1,'uint');
    end;
    [B,count] = fread(fid,prod(d),'float');
    A{ct} = reshape(B,fliplr(d));
    ct = ct + 1;
    dim = fread(fid,1,'uint');
  end;
else %else a single array 
  A = B;
end;

fclose(fid);