From f01a55721c38dcfb9ea2caac9e3fa72318bbd9ed Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 13 Dec 2006 16:04:34 +0000
Subject: [PATCH] added check whether the derivatives of the directors are
orthogonal to the directors
[[Imported from SVN: r1089]]
---
src/rodassembler.cc | 5 +++++
1 file changed, 5 insertions(+)
diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index ba24ea45..aac12b13 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -140,6 +140,11 @@ getFirstDerivativesOfDirectors(const Quaternion<double>& q,
dd_dvj[2][j] *= 2;
}
+
+ // Check: The derivatives of the directors must be orthogonal to the directors
+ for (int i=0; i<3; i++)
+ for (int j=0; j<3; j++)
+ assert (std::abs(q.director(i) * dd_dvj[i][j]) < 1e-7);
}
--
GitLab