2: #include <../src/mat/impls/aij/seq/aij.h> 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;