From af449f6ce13fed8f7daeee0aa5f579bfa3d9e00a Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 8 Mar 2012 10:42:20 +0000
Subject: [PATCH] A test for the OrthogonalMatrix class
[[Imported from SVN: r8557]]
---
test/Makefile.am | 3 ++
test/orthogonalmatrixtest.cc | 89 ++++++++++++++++++++++++++++++++++++
2 files changed, 92 insertions(+)
create mode 100644 test/orthogonalmatrixtest.cc
diff --git a/test/Makefile.am b/test/Makefile.am
index bd6d731c..829b0598 100644
--- a/test/Makefile.am
+++ b/test/Makefile.am
@@ -14,6 +14,7 @@ check_PROGRAMS = averagedistanceassemblertest \
localgeodesicfestiffnesstest \
localgfetestfunctiontest \
nestednesstest \
+ orthogonalmatrixtest \
rodassemblertest \
rotationtest \
svdtest \
@@ -41,6 +42,8 @@ cosseratenergytest_SOURCES = cosseratenergytest.cc
averagedistanceassemblertest_SOURCES = averagedistanceassemblertest.cc
+orthogonalmatrixtest_SOURCES = orthogonalmatrixtest.cc
+
rodassemblertest_SOURCES = rodassemblertest.cc
targetspacetest_SOURCES = targetspacetest.cc
diff --git a/test/orthogonalmatrixtest.cc b/test/orthogonalmatrixtest.cc
new file mode 100644
index 00000000..11b28d12
--- /dev/null
+++ b/test/orthogonalmatrixtest.cc
@@ -0,0 +1,89 @@
+#include <config.h>
+
+#include <dune/gfe/orthogonalmatrix.hh>
+
+#include "valuefactory.hh"
+
+
+/** \file
+ \brief Unit tests for the OrthogonalMatrix class
+*/
+
+using namespace Dune;
+
+double eps = 1e-6;
+
+/** \brief Test orthogonal projection onto the tangent space at q
+ */
+template <class T, int N>
+void testProjectionOntoTangentSpace(const OrthogonalMatrix<T,N>& p)
+{
+ std::vector<FieldMatrix<T,N,N> > testVectors;
+ ValueFactory<FieldMatrix<T,N,N> >::get(testVectors);
+
+ // Test each element in the list
+ for (size_t i=0; i<testVectors.size(); i++) {
+
+ // get projection
+ FieldMatrix<T,N,N> projected = p.projectOntoTangentSpace(testVectors[i]);
+
+ // Test whether the tangent vector is really tangent.
+ // The tangent space at q consist of all matrices of the form q * (a skew matrix).
+ // Hence we left-multiply by q^T and check whether the result is skew-symmetric.
+ FieldMatrix<T,N,N> skewPart(0);
+ for (int j=0; j<N; j++)
+ for (int k=0; k<N; k++)
+ for (int l=0; l<N; l++)
+ skewPart[j][k] += p.data()[l][j] * projected[l][k];
+
+ std::cout << "skew part\n" << skewPart << std::endl;
+ // is it skew?
+ for (int j=0; j<N; j++)
+ for (int k=0; k<=j; k++)
+ if ( std::abs(skewPart[j][k] + skewPart[k][j]) > eps ) {
+
+ std::cout << "matrix is not skew-symmetric\n" << skewPart << std::endl;
+ abort();
+
+ }
+ }
+
+
+}
+
+template <class T, int N>
+void test()
+{
+ std::cout << "Testing class " << className<OrthogonalMatrix<T,N> >() << std::endl;
+
+ // Use the ValueFactory for general matrices.
+ // I am too lazy to program a factory for orthogonal matrices
+ std::vector<FieldMatrix<T,N,N> > testPoints;
+ ValueFactory<FieldMatrix<T,N,N> >::get(testPoints);
+
+ int nTestPoints = testPoints.size();
+
+ // Test each element in the list
+ for (int i=0; i<nTestPoints; i++) {
+
+ // I cannot use matrices with determinant zero for testing
+ if (testPoints[i].determinant() > eps)
+ testProjectionOntoTangentSpace(OrthogonalMatrix<T,N>(testPoints[i]));
+
+ }
+
+}
+
+
+int main() try
+{
+ test<double,1>();
+ test<double,2>();
+ test<double,3>();
+
+} catch (Exception e) {
+
+ std::cout << e << std::endl;
+
+}
+
--
GitLab