diff options
author | Gael Guennebaud <g.gael@free.fr> | 2009-07-14 23:00:53 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2009-07-14 23:00:53 +0200 |
commit | 7a9519a9be524a879b44a1912ae694cfe4ee4ef7 (patch) | |
tree | da5d64e9b9a8c2658f17f54997777b28fa8fd0d2 /test/adjoint.cpp | |
parent | f5d2317b12f3d6eea68db5fb48743ca1801d10d3 (diff) |
fix typo in blue norm
Diffstat (limited to 'test/adjoint.cpp')
-rw-r--r-- | test/adjoint.cpp | 32 |
1 files changed, 9 insertions, 23 deletions
diff --git a/test/adjoint.cpp b/test/adjoint.cpp index 14ee44a0f..403f16a45 100644 --- a/test/adjoint.cpp +++ b/test/adjoint.cpp @@ -114,29 +114,15 @@ template<typename MatrixType> void adjoint(const MatrixType& m) void test_adjoint() { -// for(int i = 0; i < g_repeat; i++) { -// CALL_SUBTEST( adjoint(Matrix<float, 1, 1>()) ); -// CALL_SUBTEST( adjoint(Matrix3d()) ); -// CALL_SUBTEST( adjoint(Matrix4f()) ); -// CALL_SUBTEST( adjoint(MatrixXcf(4, 4)) ); -// CALL_SUBTEST( adjoint(MatrixXi(8, 12)) ); -// CALL_SUBTEST( adjoint(MatrixXf(21, 21)) ); -// } - // test a large matrix only once -// CALL_SUBTEST( adjoint(Matrix<float, 100, 100>()) ); - for(int i = 0; i < g_repeat; i++) - { - std::cerr.precision(20); - int s = 1000000; - double y = 1.131242353467546478463457843445677435233e23 * ei_abs(ei_random<double>()); - VectorXf v = VectorXf::Ones(s) * y; -// Vector4f x(v.segment(0,s/4).blueNorm(), v.segment(s/4+1,s/4).blueNorm(), -// v.segment((s/2)+1,s/4).blueNorm(), v.segment(3*s/4+1,s - 3*s/4-1).blueNorm()); -// std::cerr << v.norm() << " == " << v.stableNorm() << " == " << v.blueNorm() << " == " << x.norm() << "\n"; - std::cerr << v.norm() << "\n" << v.stableNorm() << "\n" << v.blueNorm() << "\n" << ei_sqrt(double(s)) * y << "\n\n\n"; - -// VectorXd d = VectorXd::Ones(s) * y;//v.cast<double>(); -// std::cerr << d.norm() << "\n" << d.stableNorm() << "\n" << d.blueNorm() << "\n" << ei_sqrt(double(s)) * y << "\n\n\n"; + for(int i = 0; i < g_repeat; i++) { + CALL_SUBTEST( adjoint(Matrix<float, 1, 1>()) ); + CALL_SUBTEST( adjoint(Matrix3d()) ); + CALL_SUBTEST( adjoint(Matrix4f()) ); + CALL_SUBTEST( adjoint(MatrixXcf(4, 4)) ); + CALL_SUBTEST( adjoint(MatrixXi(8, 12)) ); + CALL_SUBTEST( adjoint(MatrixXf(21, 21)) ); } + // test a large matrix only once + CALL_SUBTEST( adjoint(Matrix<float, 100, 100>()) ); } |