149b5e25fSSatish Balay 249b5e25fSSatish Balay #if !defined(__SBAIJ_H) 349b5e25fSSatish Balay #define __SBAIJ_H 4c833c1fbSBarry Smith #include "src/mat/matimpl.h" 56c6c5352SBarry Smith #include "src/mat/impls/baij/seq/baij.h" 649b5e25fSSatish Balay 749b5e25fSSatish Balay /* 849b5e25fSSatish Balay MATSEQSBAIJ format - Block compressed row storage. The i[] and j[] 949b5e25fSSatish Balay arrays start at 0. 1049b5e25fSSatish Balay */ 1149b5e25fSSatish Balay 1249b5e25fSSatish Balay typedef struct { 13421e10b8SBarry Smith SEQAIJHEADER(PetscScalar); 14e6b907acSBarry Smith SEQBAIJHEADER; 1513f74950SBarry Smith PetscInt *inew; /* pointer to beginning of each row of reordered matrix */ 1613f74950SBarry Smith PetscInt *jnew; /* column values: jnew + i[k] is start of row k */ 17a6175056SHong Zhang MatScalar *anew; /* nonzero diagonal and superdiagonal elements of reordered matrix */ 18d59c15a7SBarry Smith PetscScalar *solves_work; /* work space used in MatSolves */ 1913f74950SBarry Smith PetscInt solves_work_n;/* size of solves_work */ 20e2ee2a47SBarry Smith PetscScalar *relax_work; 2113f74950SBarry Smith PetscInt *a2anew; /* map used for symm permutation */ 22d36791b2SHong Zhang PetscTruth permute; /* if true, a non-trivial permutation is used for factorization */ 23941593c8SHong Zhang PetscTruth ignore_ltriangular; /* if true, ignore the lower triangular values inserted by users */ 24f5edf698SHong Zhang PetscTruth getrow_utriangular; /* if true, MatGetRow_SeqSBAIJ() is enabled to get the upper part of the row */ 2549b5e25fSSatish Balay } Mat_SeqSBAIJ; 2649b5e25fSSatish Balay 27f349cc81SSatish Balay EXTERN_C_BEGIN 28ab93d7beSBarry Smith EXTERN PetscErrorCode MatSeqSBAIJSetPreallocation_SeqSBAIJ(Mat,PetscInt,PetscInt,PetscInt*); 29f349cc81SSatish Balay EXTERN_C_END 30dfbe8321SBarry Smith EXTERN PetscErrorCode MatICCFactorSymbolic_SeqSBAIJ(Mat,IS,MatFactorInfo*,Mat *); 31dfbe8321SBarry Smith EXTERN PetscErrorCode MatDuplicate_SeqSBAIJ(Mat,MatDuplicateOption,Mat*); 32dfbe8321SBarry Smith EXTERN PetscErrorCode MatMarkDiagonal_SeqSBAIJ(Mat); 3349b5e25fSSatish Balay 34af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 35dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec); 364078b935SHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec); 374078b935SHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec); 38266135f8SHong Zhang 39af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 40dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec); 41*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec); 42*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec); 43266135f8SHong Zhang 44af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 45dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec); 46*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec); 47*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec); 48266135f8SHong Zhang 49af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 50dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec); 51*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec); 52*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec); 53266135f8SHong Zhang 54af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 55dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec); 56*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec); 57*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec); 58266135f8SHong Zhang 59af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 60dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec); 61*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec); 62*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec); 63266135f8SHong Zhang 64af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 65dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec); 66*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec); 67*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec); 6849b5e25fSSatish Balay 69af281ebdSHong Zhang EXTERN PetscErrorCode MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering(Mat,MatFactorInfo*,Mat*); 70dfbe8321SBarry Smith EXTERN PetscErrorCode MatSolve_SeqSBAIJ_N_NaturalOrdering(Mat,Vec,Vec); 71*759322afSHong Zhang EXTERN PetscErrorCode MatForwardSolve_SeqSBAIJ_N_NaturalOrdering(Mat,Vec,Vec); 72*759322afSHong Zhang EXTERN PetscErrorCode MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering(Mat,Vec,Vec); 73671cb588SHong Zhang 7413f74950SBarry Smith EXTERN PetscErrorCode MatRelax_SeqSBAIJ(Mat,Vec,PetscReal,MatSORType,PetscReal,PetscInt,PetscInt,Vec); 75f69a0ea3SMatthew Knepley EXTERN PetscErrorCode MatLoad_SeqSBAIJ(PetscViewer, MatType,Mat*); 76d06b337dSHong Zhang 7749b5e25fSSatish Balay #endif 78