Actual source code: sbaij.h
1: /* $Id: sbaij.h,v 1.21 2001/08/07 03:03:01 balay Exp $ */
3: #include src/mat/matimpl.h
8: /*
9: MATSEQSBAIJ format - Block compressed row storage. The i[] and j[]
10: arrays start at 0.
11: */
13: typedef struct {
14: PetscTruth sorted; /* if true, rows are sorted by increasing columns */
15: PetscTruth roworiented; /* if true, row-oriented input, default */
16: int nonew; /* 1 don't add new nonzeros, -1 generate error on new */
17: PetscTruth singlemalloc; /* if true a, i, and j have been obtained with
18: one big malloc */
19: int bs,bs2; /* block size, square of block size */
20: int mbs,nbs; /* rows/bs or columns/bs */
21: int s_nz,s_maxnz; /* total diagonal and superdiagonal nonzero blocks,
22: total allocated diagonal and superdiagonal nonzero blocks*/
23: int *diag; /* pointers to diagonal elements */
24: int *i; /* pointer to beginning of each row */
25: int *inew; /* pointer to beginning of each row of reordered matrix */
26: int *imax; /* maximum space allocated for each row */
27: int *ilen; /* actual length of each row */
28: int *j; /* column values: j + i[k] is start of row k of reordered matrix */
29: int *jnew; /* column values: jnew + i[k] is start of row k */
30: MatScalar *a; /* nonzero diagonal and superdiagonal elements */
31: MatScalar *anew; /* nonzero diagonal and superdiagonal elements of reordered matrix */
32: IS row,icol; /* index sets, used for reorderings */
33: PetscScalar *solve_work; /* work space used in MatSolve */
34: PetscScalar *solves_work; /* work space used in MatSolves */
35: int solves_work_n;/* size of solves_work */
36: int reallocs; /* number of mallocs done during MatSetValues()
37: as more values are set then were preallocated */
38: PetscScalar *mult_work; /* work array for matrix vector product*/
39: PetscScalar *saved_values;
41: PetscTruth keepzeroedrows; /* if true, MatZeroRows() will not change nonzero structure */
42: int *a2anew; /* map used for symm permutation */
43: PetscTruth permute; /* if true, a non-trivial permutation is used for factorization */
44: PetscTruth pivotinblocks; /* pivot inside factorization of each diagonal block */
46: /* carry MatFactorInfo from symbolic factor to numeric factor */
47: int factor_levels;
48: PetscReal factor_damping;
49: PetscReal factor_shift;
50: PetscReal factor_zeropivot;
52: int *xtoy,*xtoyB; /* map nonzero pattern of X into Y's, used by MatAXPY() */
53: Mat XtoY; /* used by MatAXPY() */
54: } Mat_SeqSBAIJ;
56: EXTERN int MatICCFactorSymbolic_SeqSBAIJ(Mat,IS,MatFactorInfo*,Mat *);
57: EXTERN int MatDuplicate_SeqSBAIJ(Mat,MatDuplicateOption,Mat*);
58: EXTERN int MatMarkDiagonal_SeqSBAIJ(Mat);
60: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering(Mat,Mat*);
61: EXTERN int MatSolve_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec);
62: EXTERN int MatSolveTranspose_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec);
64: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering(Mat,Mat*);
65: EXTERN int MatSolve_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec);
66: EXTERN int MatSolveTranspose_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec);
68: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering(Mat,Mat*);
69: EXTERN int MatSolve_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec);
70: EXTERN int MatSolveTranspose_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec);
72: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering(Mat,Mat*);
73: EXTERN int MatSolve_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec);
74: EXTERN int MatSolveTranspose_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec);
76: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering(Mat,Mat*);
77: EXTERN int MatSolve_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec);
78: EXTERN int MatSolveTranspose_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec);
80: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering(Mat,Mat*);
81: EXTERN int MatSolve_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec);
82: EXTERN int MatSolveTranspose_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec);
84: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering(Mat,Mat*);
85: EXTERN int MatSolve_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec);
86: EXTERN int MatSolveTranspose_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec);
88: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering(Mat,Mat*);
89: EXTERN int MatSolve_SeqSBAIJ_N_NaturalOrdering(Mat,Vec,Vec);
91: EXTERN int MatRelax_SeqSBAIJ(Mat,Vec,PetscReal,MatSORType,PetscReal,int,int,Vec);
92: EXTERN int MatLoad_SeqSBAIJ(PetscViewer,MatType,Mat*);
94: #endif