#include <OpenTissue/configuration.h>
#include <OpenTissue/core/math/math_random.h>
#include <OpenTissue/core/math/big/big_types.h>
#include <OpenTissue/core/math/big/big_cholesky.h>
#include <boost/numeric/ublas/triangular.hpp>
#include <boost/numeric/ublas/banded.hpp>
#include <boost/numeric/ublas/symmetric.hpp>
#include <cassert>
#include <limits>
#include <OpenTissue/utility/utility_push_boost_filter.h>
#include <boost/test/auto_unit_test.hpp>
#include <boost/test/unit_test_suite.hpp>
#include <boost/test/floating_point_comparison.hpp>
#include <boost/test/test_tools.hpp>
#include <OpenTissue/utility/utility_pop_boost_filter.h>
Defines | |
#define | BOOST_AUTO_TEST_MAIN |
Functions | |
template<class TYPE , class MATRIX > | |
ublas::triangular_adaptor < const MATRIX, TYPE > | triangular (const MATRIX &A, const TYPE &uplo=TYPE()) |
template<class MATRIX > | |
ublas::banded_adaptor< const MATRIX > | banded (const MATRIX &A, const size_t lower, const size_t upper) |
template<class TYPE , class MATRIX > | |
ublas::symmetric_adaptor < const MATRIX, TYPE > | symmetric (const MATRIX &A, const TYPE &uplo=TYPE()) |
template<class MATRIX > | |
void | fill_symm (MATRIX &L, const size_t bands=std::numeric_limits< size_t >::max()) |
template<typename matrix_type , typename vector_type > | |
void | test (matrix_type const &A, vector_type &x, vector_type const &b, vector_type const &y) |
BOOST_AUTO_TEST_SUITE (opentissue_math_big_cholesky) | |
BOOST_AUTO_TEST_CASE (Gunter_Winkler_Konstantin_Kutzkow_Test) | |
BOOST_AUTO_TEST_CASE (random_test_case) | |
BOOST_AUTO_TEST_SUITE_END () |
#define BOOST_AUTO_TEST_MAIN |
ublas::banded_adaptor<const MATRIX> banded | ( | const MATRIX & | A, | |
const size_t | lower, | |||
const size_t | upper | |||
) |
Make a immutable banded adaptor from a matrix
: A = banded(B, lower, upper);
BOOST_AUTO_TEST_CASE | ( | random_test_case | ) |
BOOST_AUTO_TEST_CASE | ( | Gunter_Winkler_Konstantin_Kutzkow_Test | ) |
BOOST_AUTO_TEST_SUITE | ( | opentissue_math_big_cholesky | ) |
BOOST_AUTO_TEST_SUITE_END | ( | ) |
void fill_symm | ( | MATRIX & | L, | |
const size_t | bands = std::numeric_limits<size_t>::max() | |||
) |
Fill lower triangular matrix L
ublas::symmetric_adaptor<const MATRIX, TYPE> symmetric | ( | const MATRIX & | A, | |
const TYPE & | uplo = TYPE() | |||
) |
Make a immutable symmetric adaptor from a matrix
: A = symmetric< lower >(B); A = symmetric(B, lower());
void test | ( | matrix_type const & | A, | |
vector_type & | x, | |||
vector_type const & | b, | |||
vector_type const & | y | |||
) |
ublas::triangular_adaptor<const MATRIX, TYPE> triangular | ( | const MATRIX & | A, | |
const TYPE & | uplo = TYPE() | |||
) |
Make a immutable triangular adaptor from a matrix
: A = triangular< lower >(B); A = triangular(B, lower());