1 2 #include <../src/mat/impls/aij/seq/aij.h> 3 4 typedef struct { 5 PetscInt nz; 6 PetscInt m; /* number of rows */ 7 PetscInt rmax; /* maximum number of columns in a row */ 8 PetscInt ncols; /* number of columns in each row */ 9 PetscInt *icols; /* columns of nonzeros, stored one column at a time */ 10 PetscScalar *acols; /* values of nonzeros, stored as icols */ 11 /* these are only needed for the parallel case */ 12 Vec xwork,fwork; 13 VecScatter xscat; /* gathers the locally needed part of global vector */ 14 PetscScalar *array; /* array used to create xwork */ 15 } Mat_AIJCRL; 16