aboutsummaryrefslogtreecommitdiffhomepage
path: root/test/adjoint.cpp
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2009-07-14 23:00:53 +0200
committerGravatar Gael Guennebaud <g.gael@free.fr>2009-07-14 23:00:53 +0200
commit7a9519a9be524a879b44a1912ae694cfe4ee4ef7 (patch)
treeda5d64e9b9a8c2658f17f54997777b28fa8fd0d2 /test/adjoint.cpp
parentf5d2317b12f3d6eea68db5fb48743ca1801d10d3 (diff)
fix typo in blue norm
Diffstat (limited to 'test/adjoint.cpp')
-rw-r--r--test/adjoint.cpp32
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>()) );
}