29 #ifndef __PASO_MKL_H__ 30 #define __PASO_MKL_H__ 36 #ifdef ESYS_INDEXTYPE_LONG 37 #define ES_PARDISO pardiso_64 38 #define ES_MKL_INT MKL_INT64 40 #define ES_PARDISO pardiso 41 #define ES_MKL_INT MKL_INT 45 #include <mkl_pardiso.h> 49 #define MKL_ERROR_NO 0 50 #define MKL_MTYPE_REAL_SYM -2 51 #define MKL_MTYPE_REAL_UNSYM 11 53 #define MKL_REORDERING_MINIMUM_DEGREE 0 54 #define MKL_REORDERING_NESTED_DISSECTION 2 55 #define MKL_REORDERING_NESTED_DISSECTION_OMP 3 56 #define MKL_PHASE_SYMBOLIC_FACTORIZATION 11 57 #define MKL_PHASE_FACTORIZATION 22 58 #define MKL_PHASE_SOLVE 33 59 #define MKL_PHASE_RELEASE_MEMORY -1 64 dim_t numRefinements,
bool verbose);
68 #endif // __PASO_MKL_H__ boost::shared_ptr< SparseMatrix > SparseMatrix_ptr
Definition: SparseMatrix.h:35
void MKL_free(SparseMatrix *A)
Definition: MKL.cpp:35
int index_t
type for array/matrix indices used both globally and on each rank
Definition: DataTypes.h:59
void MKL_solve(SparseMatrix_ptr A, double *out, double *in, index_t reordering, dim_t numRefinements, bool verbose)
Definition: MKL.cpp:66
index_t dim_t
Definition: DataTypes.h:64