123456789101112131415161718192021222324252627282930313233343536373839404142434445 |
- #ifdef EIGEN_DEFAULT_TO_ROW_MAJOR
- #undef EIGEN_DEFAULT_TO_ROW_MAJOR
- #endif
- #include "sparse_solver.h"
- #include <Eigen/SparseLU>
- #include <unsupported/Eigen/SparseExtra>
- template<typename T> void test_sparselu_T()
- {
- SparseLU<SparseMatrix<T, ColMajor> > sparselu_colamd;
- SparseLU<SparseMatrix<T, ColMajor>, AMDOrdering<int> > sparselu_amd;
- SparseLU<SparseMatrix<T, ColMajor, long int>, NaturalOrdering<long int> > sparselu_natural;
-
- check_sparse_square_solving(sparselu_colamd, 300, 100000, true);
- check_sparse_square_solving(sparselu_amd, 300, 10000, true);
- check_sparse_square_solving(sparselu_natural, 300, 2000, true);
-
- check_sparse_square_abs_determinant(sparselu_colamd);
- check_sparse_square_abs_determinant(sparselu_amd);
-
- check_sparse_square_determinant(sparselu_colamd);
- check_sparse_square_determinant(sparselu_amd);
- }
- EIGEN_DECLARE_TEST(sparselu)
- {
- CALL_SUBTEST_1(test_sparselu_T<float>());
- CALL_SUBTEST_2(test_sparselu_T<double>());
- CALL_SUBTEST_3(test_sparselu_T<std::complex<float> >());
- CALL_SUBTEST_4(test_sparselu_T<std::complex<double> >());
- }
|