Actual source code: sbaij.h

  1: /* $Id: sbaij.h,v 1.21 2001/08/07 03:03:01 balay Exp $ */

 5:  #include src/mat/matimpl.h
 6:  #include src/mat/impls/baij/seq/baij.h

  8: /*  
  9:   MATSEQSBAIJ format - Block compressed row storage. The i[] and j[] 
 10:   arrays start at 0.
 11: */

 13: typedef struct {
 14:   SEQBAIJHEADER

 16:   int              *inew;        /* pointer to beginning of each row of reordered matrix */
 17:   int              *jnew;        /* column values: jnew + i[k] is start of row k */
 18:   MatScalar        *anew;        /* nonzero diagonal and superdiagonal elements of reordered matrix */
 19:   PetscScalar      *solves_work; /* work space used in MatSolves */
 20:   int              solves_work_n;/* size of solves_work */
 21:   int              *a2anew;        /* map used for symm permutation */
 22:   PetscTruth       permute;        /* if true, a non-trivial permutation is used for factorization */

 24:   /* carry MatFactorInfo from symbolic factor to numeric factor */
 25:   int              factor_levels;
 26:   PetscReal        factor_damping;
 27:   PetscReal        factor_shift;
 28:   PetscReal        factor_zeropivot;
 29: } Mat_SeqSBAIJ;

 31: EXTERN int MatICCFactorSymbolic_SeqSBAIJ(Mat,IS,MatFactorInfo*,Mat *);
 32: EXTERN int MatDuplicate_SeqSBAIJ(Mat,MatDuplicateOption,Mat*);
 33: EXTERN int MatMarkDiagonal_SeqSBAIJ(Mat);

 35: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering(Mat,Mat*);
 36: EXTERN int MatSolve_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec);
 37: EXTERN int MatSolveTranspose_SeqSBAIJ_1_NaturalOrdering(Mat,Vec,Vec);

 39: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering(Mat,Mat*);
 40: EXTERN int MatSolve_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec);
 41: EXTERN int MatSolveTranspose_SeqSBAIJ_2_NaturalOrdering(Mat,Vec,Vec);

 43: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering(Mat,Mat*);
 44: EXTERN int MatSolve_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec);
 45: EXTERN int MatSolveTranspose_SeqSBAIJ_3_NaturalOrdering(Mat,Vec,Vec);

 47: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering(Mat,Mat*);
 48: EXTERN int MatSolve_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec);
 49: EXTERN int MatSolveTranspose_SeqSBAIJ_4_NaturalOrdering(Mat,Vec,Vec);

 51: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering(Mat,Mat*);
 52: EXTERN int MatSolve_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec);
 53: EXTERN int MatSolveTranspose_SeqSBAIJ_5_NaturalOrdering(Mat,Vec,Vec);

 55: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering(Mat,Mat*);
 56: EXTERN int MatSolve_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec);
 57: EXTERN int MatSolveTranspose_SeqSBAIJ_6_NaturalOrdering(Mat,Vec,Vec);

 59: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering(Mat,Mat*);
 60: EXTERN int MatSolve_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec);
 61: EXTERN int MatSolveTranspose_SeqSBAIJ_7_NaturalOrdering(Mat,Vec,Vec);

 63: EXTERN int MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering(Mat,Mat*);
 64: EXTERN int MatSolve_SeqSBAIJ_N_NaturalOrdering(Mat,Vec,Vec);

 66: EXTERN int MatRelax_SeqSBAIJ(Mat,Vec,PetscReal,MatSORType,PetscReal,int,int,Vec);
 67: EXTERN int MatLoad_SeqSBAIJ(PetscViewer,const MatType,Mat*);

 69: #endif