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