1 /* 2 Intel Copyright (C) .... 3 */ 4 5 #include "sparse_solver.h" 6 #include <Eigen/PardisoSupport> 7 8 template<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 23 void 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
