Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025 #include "main.h"
00026
00027
00028
00029
00030 template<typename MatrixType> void selfadjoint(const MatrixType& m)
00031 {
00032 typedef typename MatrixType::Index Index;
00033 typedef typename MatrixType::Scalar Scalar;
00034 typedef typename NumTraits<Scalar>::Real RealScalar;
00035
00036 Index rows = m.rows();
00037 Index cols = m.cols();
00038
00039 MatrixType m1 = MatrixType::Random(rows, cols),
00040 m3(rows, cols);
00041
00042 m1.diagonal() = m1.diagonal().real().template cast<Scalar>();
00043
00044
00045 m3 = m1.template selfadjointView<Upper>();
00046 VERIFY_IS_APPROX(MatrixType(m3.template triangularView<Upper>()), MatrixType(m1.template triangularView<Upper>()));
00047 VERIFY_IS_APPROX(m3, m3.adjoint());
00048
00049
00050 m3 = m1.template selfadjointView<Lower>();
00051 VERIFY_IS_APPROX(MatrixType(m3.template triangularView<Lower>()), MatrixType(m1.template triangularView<Lower>()));
00052 VERIFY_IS_APPROX(m3, m3.adjoint());
00053 }
00054
00055 void bug_159()
00056 {
00057 Matrix3d m = Matrix3d::Random().selfadjointView<Lower>();
00058 }
00059
00060 void test_selfadjoint()
00061 {
00062 for(int i = 0; i < g_repeat ; i++)
00063 {
00064 int s = internal::random<int>(1,20); EIGEN_UNUSED_VARIABLE(s);
00065
00066 CALL_SUBTEST_1( selfadjoint(Matrix<float, 1, 1>()) );
00067 CALL_SUBTEST_2( selfadjoint(Matrix<float, 2, 2>()) );
00068 CALL_SUBTEST_3( selfadjoint(Matrix3cf()) );
00069 CALL_SUBTEST_4( selfadjoint(MatrixXcd(s,s)) );
00070 CALL_SUBTEST_5( selfadjoint(Matrix<float,Dynamic,Dynamic,RowMajor>(s, s)) );
00071 }
00072
00073 CALL_SUBTEST_1( bug_159() );
00074 }