1 2 #if !defined(__SBAIJ_H) 3 #define __SBAIJ_H 4 #include "src/mat/matimpl.h" 5 #include "src/mat/impls/baij/seq/baij.h" 6 7 /* 8 MATSEQSBAIJ format - Block compressed row storage. The i[] and j[] 9 arrays start at 0. 10 */ 11 12 typedef struct { 13 SEQBAIJHEADER 14 15 int *inew; /* pointer to beginning of each row of reordered matrix */ 16 int *jnew; /* column values: jnew + i[k] is start of row k */ 17 MatScalar *anew; /* nonzero diagonal and superdiagonal elements of reordered matrix */ 18 PetscScalar *solves_work; /* work space used in MatSolves */ 19 int solves_work_n;/* size of solves_work */ 20 int *a2anew; /* map used for symm permutation */ 21 PetscTruth permute; /* if true, a non-trivial permutation is used for factorization */ 22 23 /* carry MatFactorInfo from symbolic factor to numeric factor */ 24 int factor_levels; 25 PetscReal factor_damping; 26 PetscReal factor_shift; 27 PetscReal factor_zeropivot; 28 } Mat_SeqSBAIJ; 29 30 EXTERN PetscErrorCode MatICCFactorSymbolic_SeqSBAIJ(Mat,IS,MatFactorInfo*,Mat *); 31 EXTERN PetscErrorCode MatDuplicate_SeqSBAIJ(Mat,MatDuplicateOption,Mat*); 32 EXTERN PetscErrorCode MatMarkDiagonal_SeqSBAIJ(Mat); 33 34 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering(Mat,Mat*); 35 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec); 36 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec); 37 38 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering(Mat,Mat*); 39 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec); 40 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec); 41 42 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering(Mat,Mat*); 43 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec); 44 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec); 45 46 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering(Mat,Mat*); 47 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec); 48 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec); 49 50 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering(Mat,Mat*); 51 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec); 52 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec); 53 54 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering(Mat,Mat*); 55 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec); 56 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec); 57 58 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering(Mat,Mat*); 59 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec); 60 EXTERN PetscErrorCode MatSolveTranspose_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec); 61 62 EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering(Mat,Mat*); 63 EXTERN PetscErrorCode MatSolve_SeqSBAIJ_N_NaturalOrdering(Mat,Vec,Vec); 64 65 EXTERN PetscErrorCode MatRelax_SeqSBAIJ(Mat,Vec,PetscReal,MatSORType,PetscReal,int,int,Vec); 66 EXTERN PetscErrorCode MatLoad_SeqSBAIJ(PetscViewer,const MatType,Mat*); 67 68 #endif 69