1f6f390a9SLois Curfman McInnes 2f6f390a9SLois Curfman McInnes #if !defined(__DENSE_H) 3f6f390a9SLois Curfman McInnes #define __DENSE_H 4af0996ceSBarry Smith #include <petsc/private/matimpl.h> 5afea5027SHong 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 */ 20abc3b08eSStefano Zampini Mat ptapwork; /* workspace (SeqDense matrix) for PtAP */ 21afea5027SHong Zhang 22afea5027SHong Zhang Mat_MatTransMatMult *atb; /* used by MatTransposeMatMult_SeqAIJ_SeqDense */ 23f6f390a9SLois Curfman McInnes } Mat_SeqDense; 24f6f390a9SLois Curfman McInnes 2509573ac7SBarry Smith extern PetscErrorCode MatMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 2609573ac7SBarry Smith extern PetscErrorCode MatMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 2775648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 2875648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 2975648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 3009573ac7SBarry Smith extern PetscErrorCode MatMatMultSymbolic_SeqAIJ_SeqDense(Mat,Mat,PetscReal,Mat*); 3109573ac7SBarry Smith extern PetscErrorCode MatMatMultNumeric_SeqAIJ_SeqDense(Mat,Mat,Mat); 3218476b09SLois Curfman McInnes 33afea5027SHong Zhang PETSC_INTERN PetscErrorCode MatDestroy_SeqDense_MatTransMatMult(Mat); 34afea5027SHong Zhang 35*150d2497SBarry Smith PETSC_INTERN PetscErrorCode MatMatMult_SeqAIJ_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 36*150d2497SBarry Smith PETSC_INTERN PetscErrorCode MatMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 37ec01deb9SMatthew Knepley 38f6f390a9SLois Curfman McInnes #endif 39