Date: 2022-03-31T04:46:59+00:00

Actual source code: mpisell.h

      4: #endif
      5: #include <../src/mat/impls/sell/seq/sell.h>

      7: typedef struct {
      8:   Mat A,B;                             /* local submatrices: A (diag part),
      9:                                            B (off-diag part) */
     10:   PetscMPIInt size;                     /* size of communicator */
     11:   PetscMPIInt rank;                     /* rank of proc in communicator */

     13:   /* The following variables are used for matrix assembly */
     14:   PetscBool   donotstash;               /* PETSC_TRUE if off processor entries dropped */
     15:   MPI_Request *send_waits;              /* array of send requests */
     16:   MPI_Request *recv_waits;              /* array of receive requests */
     17:   PetscInt    nsends,nrecvs;           /* numbers of sends and receives */
     18:   PetscScalar *svalues,*rvalues;       /* sending and receiving data */
     19:   PetscInt    rmax;                     /* maximum message length */
     20: #if defined(PETSC_USE_CTABLE)
     21:   PetscTable colmap;
     22: #else
     23:   PetscInt *colmap;                     /* local col number of off-diag col */
     24: #endif
     25:   PetscInt *garray;                     /* global index of all off-processor columns */

     27:   /* The following variables are used for matrix-vector products */
     28:   Vec        lvec;                 /* local vector */
     29:   Vec        diag;
     30:   VecScatter Mvctx;                /* scatter context for vector */
     31:   PetscBool  roworiented;          /* if true, row-oriented input, default true */

     33:   /* The following variables are for MatGetRow() */
     34:   PetscInt    *rowindices;         /* column indices for row */
     35:   PetscScalar *rowvalues;          /* nonzero values in row */
     36:   PetscBool   getrowactive;        /* indicates MatGetRow(), not restored */

     38:   PetscInt *ld;                    /* number of entries per row left of diagona block */
     39: } Mat_MPISELL;

     41: PETSC_EXTERN PetscErrorCode MatCreate_MPISELL(Mat);
     42: PETSC_INTERN PetscErrorCode MatSetUpMultiply_MPISELL(Mat);

     44: PETSC_INTERN PetscErrorCode MatDisAssemble_MPISELL(Mat);
     45: PETSC_INTERN PetscErrorCode MatDuplicate_MPISELL(Mat,MatDuplicateOption,Mat*);

     47: PETSC_INTERN PetscErrorCode MatDestroy_MPISELL_PtAP(Mat);
     48: PETSC_INTERN PetscErrorCode MatDestroy_MPISELL(Mat);

     50: PETSC_INTERN PetscErrorCode MatSetValues_MPISELL(Mat,PetscInt,const PetscInt[],PetscInt,const PetscInt[],const PetscScalar [],InsertMode);
     51: PETSC_INTERN PetscErrorCode MatSetOption_MPISELL(Mat,MatOption,PetscBool);
     52: PETSC_INTERN PetscErrorCode MatGetSeqNonzeroStructure_MPISELL(Mat,Mat*);

     54: PETSC_INTERN PetscErrorCode MatSetFromOptions_MPISELL(PetscOptionItems*,Mat);
     55: PETSC_INTERN PetscErrorCode MatMPISELLSetPreallocation_MPISELL(Mat,PetscInt,const PetscInt[],PetscInt,const PetscInt[]);

     57: PETSC_INTERN PetscErrorCode MatConvert_MPISELL_MPIAIJ(Mat,MatType,MatReuse,Mat*);
     58: PETSC_INTERN PetscErrorCode MatConvert_MPIAIJ_MPISELL(Mat,MatType,MatReuse,Mat*);

     60: PETSC_INTERN PetscErrorCode MatSOR_MPISELL(Mat,Vec,PetscReal,MatSORType,PetscReal,PetscInt,PetscInt,Vec);

     62: PETSC_INTERN PetscErrorCode MatCreateColmap_MPISELL_Private(Mat);
     63: PETSC_INTERN PetscErrorCode MatDiagonalScaleLocal_MPISELL(Mat,Vec);