1f6f390a9SLois Curfman McInnes 2f6f390a9SLois Curfman McInnes #if !defined(__DENSE_H) 3f6f390a9SLois Curfman McInnes #define __DENSE_H 4da33ede1SBarry Smith #include "src/mat/matimpl.h" 5da33ede1SBarry Smith 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 */ 13273d9f13SBarry Smith PetscTruth roworiented; /* if true, row oriented input (default) */ 1413f74950SBarry Smith PetscInt pad; /* padding */ 154ce68768SBarry Smith PetscBLASInt *pivots; /* pivots in LU factorization */ 16*284134d9SBarry Smith PetscBLASInt lda; /* Lapack leading dimension of data */ 17*284134d9SBarry Smith PetscBLASInt Nmax; /* indicates the second dimension of data */ 18273d9f13SBarry Smith PetscTruth user_alloc; /* true if the user provided the dense data */ 19f6f390a9SLois Curfman McInnes } Mat_SeqDense; 20f6f390a9SLois Curfman McInnes 21dfbe8321SBarry Smith EXTERN PetscErrorCode MatMult_SeqDense(Mat A,Vec,Vec); 22dfbe8321SBarry Smith EXTERN PetscErrorCode MatMultAdd_SeqDense(Mat A,Vec,Vec,Vec); 23dfbe8321SBarry Smith EXTERN PetscErrorCode MatMultTranspose_SeqDense(Mat A,Vec,Vec); 24dfbe8321SBarry Smith EXTERN PetscErrorCode MatMultTransposeAdd_SeqDense(Mat A,Vec,Vec,Vec); 2518476b09SLois Curfman McInnes 26f6f390a9SLois Curfman McInnes #endif 27