1*8a729477SBarry Smith 2*8a729477SBarry Smith #include "matimpl.h" 3*8a729477SBarry Smith #include "math.h" 4*8a729477SBarry Smith 5*8a729477SBarry Smith /* The i[] and j[] arrays start at 1 not zero to support Fortran 77 */ 6*8a729477SBarry Smith /* In Fortran j[i[k]+p-1] is the pth column in row k */ 7*8a729477SBarry Smith 8*8a729477SBarry Smith /* 9*8a729477SBarry Smith singlemalloc indicates that a, i and j where all obtained with 10*8a729477SBarry Smith one big malloc 11*8a729477SBarry Smith */ 12*8a729477SBarry Smith typedef struct { 13*8a729477SBarry Smith int sorted, roworiented, nonew, singlemalloc; 14*8a729477SBarry Smith int m,n,nz,mem,*diag, /* rows and columns */ 15*8a729477SBarry Smith *i,*imax, *ilen, /* j + i[k] - 1 is start of row k */ 16*8a729477SBarry Smith *j; /* ilen is actual lenght of row */ 17*8a729477SBarry Smith Scalar *a; 18*8a729477SBarry Smith } Matiaij; 19*8a729477SBarry Smith 20*8a729477SBarry Smith 21