ground/gcs/src/libs/eigen/test/pardiso_support.cpp

   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 }