pardiso_support.cpp revision c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd
1/* 2 Intel Copyright (C) .... 3*/ 4 5#include "sparse_solver.h" 6#include <Eigen/PardisoSupport> 7 8template<typename T> void test_pardiso_T() 9{ 10 PardisoLLT < SparseMatrix<T, RowMajor>, Lower> pardiso_llt_lower; 11 PardisoLLT < SparseMatrix<T, RowMajor>, Upper> pardiso_llt_upper; 12 PardisoLDLT < SparseMatrix<T, RowMajor>, Lower> pardiso_ldlt_lower; 13 PardisoLDLT < SparseMatrix<T, RowMajor>, Upper> pardiso_ldlt_upper; 14 PardisoLU < SparseMatrix<T, RowMajor> > pardiso_lu; 15 16 check_sparse_spd_solving(pardiso_llt_lower); 17 check_sparse_spd_solving(pardiso_llt_upper); 18 check_sparse_spd_solving(pardiso_ldlt_lower); 19 check_sparse_spd_solving(pardiso_ldlt_upper); 20 check_sparse_square_solving(pardiso_lu); 21} 22 23void test_pardiso_support() 24{ 25 CALL_SUBTEST_1(test_pardiso_T<float>()); 26 CALL_SUBTEST_2(test_pardiso_T<double>()); 27 CALL_SUBTEST_3(test_pardiso_T< std::complex<float> >()); 28 CALL_SUBTEST_4(test_pardiso_T< std::complex<double> >()); 29} 30