From 27ba6f716900eceaf11a9c5b0916817fe1c1021c Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 22 Apr 2011 16:29:17 +0000
Subject: [PATCH] test method orthonormalFrame
[[Imported from SVN: r7201]]
---
test/unitvectortest.cc | 25 +++++++++++++++++++++----
1 file changed, 21 insertions(+), 4 deletions(-)
diff --git a/test/unitvectortest.cc b/test/unitvectortest.cc
index 0b1880fe..ded44a5f 100644
--- a/test/unitvectortest.cc
+++ b/test/unitvectortest.cc
@@ -82,6 +82,17 @@ FieldMatrix<double,worldDim,worldDim> getSecondDerivativeOfSecondArgumentFD(cons
return ret2;
}
+template <class TargetSpace, int worldDim>
+void testOrthonormalFrame(const TargetSpace& a)
+{
+ const size_t spaceDim = TargetSpace::dim;
+ FieldMatrix<double,spaceDim,worldDim> B = a.orthonormalFrame();
+
+ for (int i=0; i<spaceDim; i++)
+ for (int j=0; j<spaceDim; j++)
+ assert( std::fabs(B[i]*B[j] - (i==j)) < 1e-10 );
+}
+
template <class TargetSpace, int dim>
void testDerivativeOfSquaredDistance(const TargetSpace& a, const TargetSpace& b)
{
@@ -291,10 +302,13 @@ void testUnitVector2d()
// Set up elements of S^1
for (int i=0; i<nTestPoints; i++) {
+ Dune::array<double,2> w0 = {testPoints[i][0], testPoints[i][1]};
+ UnitVector<2> v0(w0);
+
+ testOrthonormalFrame<UnitVector<2>, 2>(v0);
+
for (int j=0; j<nTestPoints; j++) {
- Dune::array<double,2> w0 = {testPoints[i][0], testPoints[i][1]};
- UnitVector<2> v0(w0);
Dune::array<double,2> w1 = {testPoints[j][0], testPoints[j][1]};
UnitVector<2> v1(w1);
@@ -319,10 +333,13 @@ void testUnitVector3d()
// Set up elements of S^1
for (int i=0; i<nTestPoints; i++) {
+ Dune::array<double,3> w0 = {testPoints[i][0], testPoints[i][1], testPoints[i][2]};
+ UnitVector<3> uv0(w0);
+
+ testOrthonormalFrame<UnitVector<3>, 3>(uv0);
+
for (int j=0; j<nTestPoints; j++) {
- Dune::array<double,3> w0 = {testPoints[i][0], testPoints[i][1], testPoints[i][2]};
- UnitVector<3> uv0(w0);
Dune::array<double,3> w1 = {testPoints[j][0], testPoints[j][1], testPoints[j][2]};
UnitVector<3> uv1(w1);
--
GitLab