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

   1 // This file is part of Eigen, a lightweight C++ template library
   2 // for linear algebra.
   3 //
   4 // Copyright (C) 2011 Gael Guennebaud <g.gael@free.fr>
   5 //
   6 // This Source Code Form is subject to the terms of the Mozilla
   7 // Public License v. 2.0. If a copy of the MPL was not distributed
   8 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
   9
  10 #define EIGEN_NO_DEBUG_SMALL_PRODUCT_BLOCKS
  11 #include "sparse_solver.h"
  12
  13 #include <Eigen/CholmodSupport>
  14
  15 template<typename T> void test_cholmod_T()
  16 {
  17   CholmodDecomposition<SparseMatrix<T>, Lower> g_chol_colmajor_lower; g_chol_colmajor_lower.setMode(CholmodSupernodalLLt);
  18   CholmodDecomposition<SparseMatrix<T>, Upper> g_chol_colmajor_upper; g_chol_colmajor_upper.setMode(CholmodSupernodalLLt);
  19   CholmodDecomposition<SparseMatrix<T>, Lower> g_llt_colmajor_lower;  g_llt_colmajor_lower.setMode(CholmodSimplicialLLt);
  20   CholmodDecomposition<SparseMatrix<T>, Upper> g_llt_colmajor_upper;  g_llt_colmajor_upper.setMode(CholmodSimplicialLLt);
  21   CholmodDecomposition<SparseMatrix<T>, Lower> g_ldlt_colmajor_lower; g_ldlt_colmajor_lower.setMode(CholmodLDLt);
  22   CholmodDecomposition<SparseMatrix<T>, Upper> g_ldlt_colmajor_upper; g_ldlt_colmajor_upper.setMode(CholmodLDLt);
  23
  24   CholmodSupernodalLLT<SparseMatrix<T>, Lower> chol_colmajor_lower;
  25   CholmodSupernodalLLT<SparseMatrix<T>, Upper> chol_colmajor_upper;
  26   CholmodSimplicialLLT<SparseMatrix<T>, Lower> llt_colmajor_lower;
  27   CholmodSimplicialLLT<SparseMatrix<T>, Upper> llt_colmajor_upper;
  28   CholmodSimplicialLDLT<SparseMatrix<T>, Lower> ldlt_colmajor_lower;
  29   CholmodSimplicialLDLT<SparseMatrix<T>, Upper> ldlt_colmajor_upper;
  30
  31   check_sparse_spd_solving(g_chol_colmajor_lower);
  32   check_sparse_spd_solving(g_chol_colmajor_upper);
  33   check_sparse_spd_solving(g_llt_colmajor_lower);
  34   check_sparse_spd_solving(g_llt_colmajor_upper);
  35   check_sparse_spd_solving(g_ldlt_colmajor_lower);
  36   check_sparse_spd_solving(g_ldlt_colmajor_upper);
  37
  38   check_sparse_spd_solving(chol_colmajor_lower);
  39   check_sparse_spd_solving(chol_colmajor_upper);
  40   check_sparse_spd_solving(llt_colmajor_lower);
  41   check_sparse_spd_solving(llt_colmajor_upper);
  42   check_sparse_spd_solving(ldlt_colmajor_lower);
  43   check_sparse_spd_solving(ldlt_colmajor_upper);
  44
  45   check_sparse_spd_determinant(chol_colmajor_lower);
  46   check_sparse_spd_determinant(chol_colmajor_upper);
  47   check_sparse_spd_determinant(llt_colmajor_lower);
  48   check_sparse_spd_determinant(llt_colmajor_upper);
  49   check_sparse_spd_determinant(ldlt_colmajor_lower);
  50   check_sparse_spd_determinant(ldlt_colmajor_upper);
  51 }
  52
  53 void test_cholmod_support()
  54 {
  55   CALL_SUBTEST_1(test_cholmod_T<double>());
  56   CALL_SUBTEST_2(test_cholmod_T<std::complex<double> >());
  57 }