cp readslice.R ts_readslice.R
[wrf-fire-matlab.git] / util1_jan / read_s.m
blob6116c368f48016da05ef24ebc2ad26ac1bb56d80
1 function a=read_s(name,num1,num2)
2 % read sparse matrix from file created by call write_array subroutines with this name num1 num2
3 % the mode of reading may change; do not call the functions below directly
4 if ~exist('num1','var'),
5         num1=-1;
6 end
7 if ~exist('num2','var'),
8         num2=-1;
9 end
10 % read the pieces
11 filename=file_name([name,'_sz'],num1,num2);
12 nn=read_array(filename);
13 filename=file_name([name,'_ia'],num1,num2);
14 ia=read_array(filename);
15 filename=file_name([name,'_ja'],num1,num2);
16 ja=read_array(filename);
17 filename=file_name([name,'_da'],num1,num2);
18 da=read_array(filename);
19 % create the sparse magtrix
20 a=sparse(ia,ja,da,nn(1),nn(2));