1 2 #if !defined(__DENSE_H) 3 #define __DENSE_H 4 #include <petsc-private/matimpl.h> 5 6 7 /* 8 MATSEQDENSE format - conventional dense Fortran storage (by columns) 9 */ 10 11 typedef struct { 12 PetscScalar *v; /* matrix elements */ 13 PetscBool roworiented; /* if true, row oriented input (default) */ 14 PetscInt pad; /* padding */ 15 PetscBLASInt *pivots; /* pivots in LU factorization */ 16 PetscBLASInt lda; /* Lapack leading dimension of data */ 17 PetscBool changelda; /* change lda on resize? Default unless user set lda */ 18 PetscBLASInt Mmax,Nmax; /* indicates the largest dimensions of data possible */ 19 PetscBool user_alloc; /* true if the user provided the dense data */ 20 } Mat_SeqDense; 21 22 PETSC_INTERN PetscErrorCode MatMult_SeqDense(Mat A,Vec,Vec); 23 PETSC_INTERN PetscErrorCode MatMultAdd_SeqDense(Mat A,Vec,Vec,Vec); 24 PETSC_INTERN PetscErrorCode MatMultTranspose_SeqDense(Mat A,Vec,Vec); 25 PETSC_INTERN PetscErrorCode MatMultTransposeAdd_SeqDense(Mat A,Vec,Vec,Vec); 26 PETSC_INTERN PetscErrorCode MatMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 27 PETSC_INTERN PetscErrorCode MatMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 28 PETSC_INTERN PetscErrorCode MatTransposeMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 29 PETSC_INTERN PetscErrorCode MatTransposeMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 30 PETSC_INTERN PetscErrorCode MatTransposeMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 31 PETSC_INTERN PetscErrorCode MatMatMultSymbolic_SeqAIJ_SeqDense(Mat,Mat,PetscReal,Mat*); 32 PETSC_INTERN PetscErrorCode MatMatMultNumeric_SeqAIJ_SeqDense(Mat,Mat,Mat); 33 34 EXTERN_C_BEGIN 35 PETSC_INTERN PetscErrorCode MatMatMult_SeqAIJ_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 36 PETSC_INTERN PetscErrorCode MatMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 37 EXTERN_C_END 38 39 #endif 40