1f6f390a9SLois Curfman McInnes 2f6f390a9SLois Curfman McInnes #if !defined(__DENSE_H) 3f6f390a9SLois Curfman McInnes #define __DENSE_H 4b45d2f2cSJed Brown #include <petsc-private/matimpl.h> 5*afea5027SHong Zhang #include <../src/mat/impls/aij/seq/aij.h> /* Mat_MatTransMatMult is defined here */ 6f6f390a9SLois Curfman McInnes 7f6f390a9SLois Curfman McInnes /* 8f6f390a9SLois Curfman McInnes MATSEQDENSE format - conventional dense Fortran storage (by columns) 9f6f390a9SLois Curfman McInnes */ 10f6f390a9SLois Curfman McInnes 11f6f390a9SLois Curfman McInnes typedef struct { 12ea709b57SSatish Balay PetscScalar *v; /* matrix elements */ 13ace3abfcSBarry Smith PetscBool roworiented; /* if true, row oriented input (default) */ 1413f74950SBarry Smith PetscInt pad; /* padding */ 154ce68768SBarry Smith PetscBLASInt *pivots; /* pivots in LU factorization */ 16284134d9SBarry Smith PetscBLASInt lda; /* Lapack leading dimension of data */ 17ace3abfcSBarry Smith PetscBool changelda; /* change lda on resize? Default unless user set lda */ 1821a2c019SBarry Smith PetscBLASInt Mmax,Nmax; /* indicates the largest dimensions of data possible */ 19ace3abfcSBarry Smith PetscBool user_alloc; /* true if the user provided the dense data */ 20*afea5027SHong Zhang 21*afea5027SHong Zhang Mat_MatTransMatMult *atb; /* used by MatTransposeMatMult_SeqAIJ_SeqDense */ 22f6f390a9SLois Curfman McInnes } Mat_SeqDense; 23f6f390a9SLois Curfman McInnes 2409573ac7SBarry Smith extern PetscErrorCode MatMult_SeqDense(Mat A,Vec,Vec); 2509573ac7SBarry Smith extern PetscErrorCode MatMultAdd_SeqDense(Mat A,Vec,Vec,Vec); 2609573ac7SBarry Smith extern PetscErrorCode MatMultTranspose_SeqDense(Mat A,Vec,Vec); 2709573ac7SBarry Smith extern PetscErrorCode MatMultTransposeAdd_SeqDense(Mat A,Vec,Vec,Vec); 2809573ac7SBarry Smith extern PetscErrorCode MatMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 2909573ac7SBarry Smith extern PetscErrorCode MatMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 3075648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 3175648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 3275648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 3309573ac7SBarry Smith extern PetscErrorCode MatMatMultSymbolic_SeqAIJ_SeqDense(Mat,Mat,PetscReal,Mat*); 3409573ac7SBarry Smith extern PetscErrorCode MatMatMultNumeric_SeqAIJ_SeqDense(Mat,Mat,Mat); 3518476b09SLois Curfman McInnes 36*afea5027SHong Zhang PETSC_INTERN PetscErrorCode MatDestroy_SeqDense_MatTransMatMult(Mat); 37*afea5027SHong Zhang 388cc058d9SJed Brown PETSC_EXTERN PetscErrorCode MatMatMult_SeqAIJ_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 398cc058d9SJed Brown PETSC_EXTERN PetscErrorCode MatMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 40ec01deb9SMatthew Knepley 41f6f390a9SLois Curfman McInnes #endif 42