aboutsummaryrefslogtreecommitdiffhomepage
path: root/unsupported/test/mpreal_support.cpp
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2016-11-05 11:28:53 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2016-11-05 11:28:53 +0100
commit55b4fd1d401440d51bf5bc663385a48e4981ff66 (patch)
tree8d1d2689f2b2972d82c39e02cdc6e18a72cab764 /unsupported/test/mpreal_support.cpp
parenta354c3ca59878f9df5af45f86ac01824c6e551c5 (diff)
Extend mpreal unit test to check LLT with complexes.
Diffstat (limited to 'unsupported/test/mpreal_support.cpp')
-rw-r--r--unsupported/test/mpreal_support.cpp8
1 files changed, 8 insertions, 0 deletions
diff --git a/unsupported/test/mpreal_support.cpp b/unsupported/test/mpreal_support.cpp
index ffa5691eb..685e7ea45 100644
--- a/unsupported/test/mpreal_support.cpp
+++ b/unsupported/test/mpreal_support.cpp
@@ -12,6 +12,7 @@ void test_mpreal_support()
// set precision to 256 bits (double has only 53 bits)
mpreal::set_default_prec(256);
typedef Matrix<mpreal,Eigen::Dynamic,Eigen::Dynamic> MatrixXmp;
+ typedef Matrix<std::complex<mpreal>,Eigen::Dynamic,Eigen::Dynamic> MatrixXcmp;
std::cerr << "epsilon = " << NumTraits<mpreal>::epsilon() << "\n";
std::cerr << "dummy_precision = " << NumTraits<mpreal>::dummy_precision() << "\n";
@@ -25,6 +26,10 @@ void test_mpreal_support()
MatrixXmp B = MatrixXmp::Random(s,s);
MatrixXmp S = A.adjoint() * A;
MatrixXmp X;
+ MatrixXcmp Ac = MatrixXcmp::Random(s,s);
+ MatrixXcmp Bc = MatrixXcmp::Random(s,s);
+ MatrixXcmp Sc = Ac.adjoint() * Ac;
+ MatrixXcmp Xc;
// Basic stuffs
VERIFY_IS_APPROX(A.real(), A);
@@ -37,6 +42,9 @@ void test_mpreal_support()
// Cholesky
X = S.selfadjointView<Lower>().llt().solve(B);
VERIFY_IS_APPROX((S.selfadjointView<Lower>()*X).eval(),B);
+
+ Xc = Sc.selfadjointView<Lower>().llt().solve(Bc);
+ VERIFY_IS_APPROX((Sc.selfadjointView<Lower>()*Xc).eval(),Bc);
// partial LU
X = A.lu().solve(B);