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 */ 13d3042a70SBarry Smith PetscScalar *unplacedarray; /* if one called MatDensePlaceArray(), this is where it stashed the original */ 14ace3abfcSBarry Smith PetscBool roworiented; /* if true, row oriented input (default) */ 1513f74950SBarry Smith PetscInt pad; /* padding */ 164ce68768SBarry Smith PetscBLASInt *pivots; /* pivots in LU factorization */ 17a49dc2a2SStefano Zampini PetscBLASInt lfwork; /* length of work array in factorization */ 18a49dc2a2SStefano Zampini PetscScalar *fwork; /* work array in factorization */ 19284134d9SBarry Smith PetscBLASInt lda; /* Lapack leading dimension of data */ 20ace3abfcSBarry Smith PetscBool changelda; /* change lda on resize? Default unless user set lda */ 2121a2c019SBarry Smith PetscBLASInt Mmax,Nmax; /* indicates the largest dimensions of data possible */ 22ace3abfcSBarry Smith PetscBool user_alloc; /* true if the user provided the dense data */ 23d3042a70SBarry Smith PetscBool unplaced_user_alloc; 24abc3b08eSStefano Zampini Mat ptapwork; /* workspace (SeqDense matrix) for PtAP */ 25afea5027SHong Zhang 26afea5027SHong Zhang Mat_MatTransMatMult *atb; /* used by MatTransposeMatMult_SeqAIJ_SeqDense */ 27f6f390a9SLois Curfman McInnes } Mat_SeqDense; 28f6f390a9SLois Curfman McInnes 2909573ac7SBarry Smith extern PetscErrorCode MatMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 3009573ac7SBarry Smith extern PetscErrorCode MatMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 3175648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 3275648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 3375648e8dSHong Zhang extern PetscErrorCode MatTransposeMatMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 3469f65d41SStefano Zampini extern PetscErrorCode MatMatTransposeMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 3569f65d41SStefano Zampini extern PetscErrorCode MatMatTransposeMultSymbolic_SeqDense_SeqDense(Mat,Mat,PetscReal,Mat*); 3669f65d41SStefano Zampini extern PetscErrorCode MatMatTransposeMultNumeric_SeqDense_SeqDense(Mat,Mat,Mat); 3709573ac7SBarry Smith extern PetscErrorCode MatMatMultSymbolic_SeqAIJ_SeqDense(Mat,Mat,PetscReal,Mat*); 3809573ac7SBarry Smith extern PetscErrorCode MatMatMultNumeric_SeqAIJ_SeqDense(Mat,Mat,Mat); 3918476b09SLois Curfman McInnes 40afea5027SHong Zhang PETSC_INTERN PetscErrorCode MatDestroy_SeqDense_MatTransMatMult(Mat); 41afea5027SHong Zhang 42150d2497SBarry Smith PETSC_INTERN PetscErrorCode MatMatMult_SeqAIJ_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 43150d2497SBarry Smith PETSC_INTERN PetscErrorCode MatMatMult_SeqDense_SeqDense(Mat,Mat,MatReuse,PetscReal,Mat*); 4405709791SSatish Balay PETSC_EXTERN PetscErrorCode MatSeqDenseInvertFactors_Private(Mat); 45ec01deb9SMatthew Knepley 46*d528f656SJakub Kruzik PETSC_INTERN PetscErrorCode MatCreateMPIMatConcatenateSeqMat_SeqDense(MPI_Comm,Mat,PetscInt,MatReuse,Mat*); 47*d528f656SJakub Kruzik PETSC_INTERN PetscErrorCode MatCreateMPIMatConcatenateSeqMat_MPIDense(MPI_Comm,Mat,PetscInt,MatReuse,Mat*); 48*d528f656SJakub Kruzik 49f6f390a9SLois Curfman McInnes #endif 50