From af40fb3dad0065689127508b7b4899fc78ed3a9a Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 26 Oct 2011 15:43:43 +0000
Subject: [PATCH] abort when the test fails, instead of just printing loads of
numbers
[[Imported from SVN: r8057]]
---
test/cosseratenergytest.cc | 21 +++++++++++++++++----
1 file changed, 17 insertions(+), 4 deletions(-)
diff --git a/test/cosseratenergytest.cc b/test/cosseratenergytest.cc
index 9c6ce855..44e6bad0 100644
--- a/test/cosseratenergytest.cc
+++ b/test/cosseratenergytest.cc
@@ -129,6 +129,13 @@ void testEnergy(const GridType* grid, const std::vector<TargetSpace>& coefficien
CosseratEnergyLocalStiffness<typename GridType::LeafGridView,LocalFiniteElement,3> assembler(materialParameters,
NULL,
&zeroFunction);
+
+ // compute reference energy
+ double referenceEnergy = assembler.energy(*grid->template leafbegin<0>(),
+ feCache.get(grid->template leafbegin<0>()->type()),
+ coefficients);
+
+ // rotate the entire configuration
std::vector<TargetSpace> rotatedCoefficients(coefficients.size());
std::vector<Rotation<3> > testRotations;
@@ -151,9 +158,13 @@ void testEnergy(const GridType* grid, const std::vector<TargetSpace>& coefficien
rotatedCoefficients[j].q = testRotations[i].mult(coefficients[j].q);
}
- std::cout << "energy: " << assembler.energy(*grid->template leafbegin<0>(),
- feCache.get(grid->template leafbegin<0>()->type()),
- rotatedCoefficients) << std::endl;
+ double energy = assembler.energy(*grid->template leafbegin<0>(),
+ feCache.get(grid->template leafbegin<0>()->type()),
+ rotatedCoefficients);
+
+ assert(std::fabs(energy-referenceEnergy) < 1e-4);
+
+ //std::cout << "energy: " << energy << std::endl;
}
@@ -163,6 +174,8 @@ void testEnergy(const GridType* grid, const std::vector<TargetSpace>& coefficien
template <int domainDim>
void testFrameInvariance()
{
+ std::cout << " --- Testing frame invariance of the Cosserat energy, domain dimension: " << domainDim << " ---" << std::endl;
+
// ////////////////////////////////////////////////////////
// Make a test grid consisting of a single simplex
// ////////////////////////////////////////////////////////
@@ -215,7 +228,7 @@ void testFrameInvariance()
int main(int argc, char** argv) try
{
const int domainDim = 2;
- std::cout << " --- Testing Rotation<3>, domain dimension: " << domainDim << " ---" << std::endl;
+ std::cout << " --- Testing derivative of rotation matrix, domain dimension: " << domainDim << " ---" << std::endl;
std::vector<Rotation<3,double> > testPoints;
--
GitLab