From 679c496289d5b82311f932cbcc879de5f1f213b4 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 7 Jun 2011 09:17:17 +0000
Subject: [PATCH] abort if there is a difference between the analytical Hessian
and the fd Hessian
[[Imported from SVN: r7386]]
---
test/localgeodesicfestiffnesstest.cc | 9 +++++++++
1 file changed, 9 insertions(+)
diff --git a/test/localgeodesicfestiffnesstest.cc b/test/localgeodesicfestiffnesstest.cc
index 73dabd80..08f0b7fc 100644
--- a/test/localgeodesicfestiffnesstest.cc
+++ b/test/localgeodesicfestiffnesstest.cc
@@ -179,6 +179,15 @@ void testHessian()
std::cout << "hessian:" << std::endl;
printmatrix(std::cout, hessian, "hessian", "--");
+ ///////////////////////////////////////////////////////////////////////////////////
+ // Abort if there is a difference
+ ///////////////////////////////////////////////////////////////////////////////////
+ Matrix<FieldMatrix<double,spaceDim,spaceDim> > difference = hessian;
+ difference -= fdHessian;
+
+ if (difference.infinity_norm() > 1e-4)
+ assert(false);
+
}
}
--
GitLab