1 
2 #if !defined(__MPISBAIJ_H)
3 #define __MPISBAIJ_H
4 #include <../src/mat/impls/baij/seq/baij.h>
5 #include <../src/mat/impls/sbaij/seq/sbaij.h>
6 #include <../src/mat/impls/baij/mpi/mpibaij.h>
7 
8 typedef struct {
9   MPIBAIJHEADER;
10   Vec        slvec0,slvec1;               /* parallel vectors */
11   Vec        slvec0b,slvec1a,slvec1b;     /* seq vectors: local partition of slvec0 and slvec1 */
12   VecScatter sMvctx;                      /* scatter context for vector used for reducing communication */
13 
14   Vec diag;                               /* used in MatSOR_MPISBAIJ() with Eisenstat */
15   Vec bb1,xx1;
16 
17   /* these are used in MatSetValues() as tmp space before passing to the stasher */
18   PetscInt  n_loc,*in_loc;                /* nloc is length of in_loc and v_loc */
19   MatScalar *v_loc;
20 } Mat_MPISBAIJ;
21 
22 PETSC_INTERN PetscErrorCode MatLoad_MPISBAIJ(Mat,PetscViewer);
23 PETSC_INTERN PetscErrorCode MatMult_MPISBAIJ_Hermitian(Mat,Vec,Vec);
24 
25 PETSC_INTERN PetscErrorCode MatSetUpMultiply_MPISBAIJ(Mat);
26 PETSC_INTERN PetscErrorCode MatDisAssemble_MPISBAIJ(Mat);
27 PETSC_INTERN PetscErrorCode MatIncreaseOverlap_MPISBAIJ(Mat,PetscInt,IS[],PetscInt);
28 PETSC_INTERN PetscErrorCode MatGetRowMaxAbs_MPISBAIJ(Mat,Vec,PetscInt[]);
29 PETSC_INTERN PetscErrorCode MatSOR_MPISBAIJ(Mat,Vec,PetscReal,MatSORType,PetscReal,PetscInt,PetscInt,Vec);
30 
31 #endif
32