#pragma once

#include <../src/mat/impls/aij/seq/aij.h>

typedef struct {
  PetscInt     nz;
  PetscInt     m;     /* number of rows */
  PetscInt     rmax;  /* maximum number of columns in a row */
  PetscInt     ncols; /* number of columns in each row */
  PetscInt    *icols; /* columns of nonzeros, stored one column at a time */
  PetscScalar *acols; /* values of nonzeros, stored as icols */
  /* these are only needed for the parallel case */
  Vec          xwork, fwork;
  VecScatter   xscat; /* gathers the locally needed part of global vector */
  PetscScalar *array; /* array used to create xwork */
} Mat_AIJCRL;