diff options
Diffstat (limited to 'eigen/test/pastix_support.cpp')
| -rw-r--r-- | eigen/test/pastix_support.cpp | 44 |
1 files changed, 44 insertions, 0 deletions
diff --git a/eigen/test/pastix_support.cpp b/eigen/test/pastix_support.cpp new file mode 100644 index 0000000..14da094 --- /dev/null +++ b/eigen/test/pastix_support.cpp @@ -0,0 +1,44 @@ +// This file is part of Eigen, a lightweight C++ template library +// for linear algebra. +// +// Copyright (C) 2008-2012 Gael Guennebaud <gael.guennebaud@inria.fr> +// Copyright (C) 2012 Désiré Nuentsa-Wakam <desire.nuentsa_wakam@inria.fr> +// +// This Source Code Form is subject to the terms of the Mozilla +// Public License v. 2.0. If a copy of the MPL was not distributed +// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. +#include "sparse_solver.h" +#include <Eigen/PaStiXSupport> +#include <unsupported/Eigen/SparseExtra> + + +template<typename T> void test_pastix_T() +{ + PastixLLT< SparseMatrix<T, ColMajor>, Eigen::Lower > pastix_llt_lower; + PastixLDLT< SparseMatrix<T, ColMajor>, Eigen::Lower > pastix_ldlt_lower; + PastixLLT< SparseMatrix<T, ColMajor>, Eigen::Upper > pastix_llt_upper; + PastixLDLT< SparseMatrix<T, ColMajor>, Eigen::Upper > pastix_ldlt_upper; + PastixLU< SparseMatrix<T, ColMajor> > pastix_lu; + + check_sparse_spd_solving(pastix_llt_lower); + check_sparse_spd_solving(pastix_ldlt_lower); + check_sparse_spd_solving(pastix_llt_upper); + check_sparse_spd_solving(pastix_ldlt_upper); + check_sparse_square_solving(pastix_lu); +} + +// There is no support for selfadjoint matrices with PaStiX. +// Complex symmetric matrices should pass though +template<typename T> void test_pastix_T_LU() +{ + PastixLU< SparseMatrix<T, ColMajor> > pastix_lu; + check_sparse_square_solving(pastix_lu); +} + +void test_pastix_support() +{ + CALL_SUBTEST_1(test_pastix_T<float>()); + CALL_SUBTEST_2(test_pastix_T<double>()); + CALL_SUBTEST_3( (test_pastix_T_LU<std::complex<float> >()) ); + CALL_SUBTEST_4(test_pastix_T_LU<std::complex<double> >()); +} |