4: #include <petsc/private/matimpl.h>
  5: #include <../src/mat/impls/aij/seq/aij.h> /* Mat_MatTransMatMult is defined here */

  7: /*
  8:   MATSEQDENSE format - conventional dense Fortran storage (by columns)
  9: */

 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 */

 21:   Mat_MatTransMatMult *atb;       /* used by MatTransposeMatMult_SeqAIJ_SeqDense */
 22: } Mat_SeqDense;

 24: extern PetscErrorCode MatMult_SeqDense(Mat A,Vec,Vec);
 25: extern PetscErrorCode MatMultAdd_SeqDense(Mat A,Vec,Vec,Vec);
 26: extern PetscErrorCode MatMultTranspose_SeqDense(Mat A,Vec,Vec);
 27: extern PetscErrorCode MatMultTransposeAdd_SeqDense(Mat A,Vec,Vec,Vec);
 28: extern PetscErrorCode MatMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*);
 29: extern PetscErrorCode MatMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat);
 30: extern PetscErrorCode MatTransposeMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*);
 31: extern PetscErrorCode MatTransposeMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*);
 32: extern PetscErrorCode MatTransposeMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat);
 33: extern PetscErrorCode MatMatMultSymbolic_SeqAIJ_SeqDense(Mat,Mat,PetscReal,Mat*);
 34: extern PetscErrorCode MatMatMultNumeric_SeqAIJ_SeqDense(Mat,Mat,Mat);

 36: PETSC_INTERN PetscErrorCode MatDestroy_SeqDense_MatTransMatMult(Mat);

 38: PETSC_EXTERN PetscErrorCode MatMatMult_SeqAIJ_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*);
 39: PETSC_EXTERN PetscErrorCode MatMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*);

 41: #endif