From d6096bf940998a5b11b41aae47a7ba651ea11941 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 22 Apr 2009 08:43:07 +0000
Subject: [PATCH] norm of the derivative in a dimension-independent way
[[Imported from SVN: r4077]]
---
src/harmonicenergystiffness.hh | 16 +++++++++-------
1 file changed, 9 insertions(+), 7 deletions(-)
diff --git a/src/harmonicenergystiffness.hh b/src/harmonicenergystiffness.hh
index c0306a84..c2b2780c 100644
--- a/src/harmonicenergystiffness.hh
+++ b/src/harmonicenergystiffness.hh
@@ -52,7 +52,7 @@ energy(const Entity& element,
LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(element.type(), localSolution);
- int quadOrder = gridDim;
+ int quadOrder = 1;//gridDim;
const Dune::QuadratureRule<double, gridDim>& quad
= Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder);
@@ -87,13 +87,15 @@ energy(const Entity& element,
}
#if 0
- double derivativenorm0 = 0;
- double derivativenorm1 = 0;
- for (int i=0; i<4; i++) {
- derivativenorm0 += derivative[i][0]*derivative[i][0];
- derivativenorm1 += derivative[i][1]*derivative[i][1];
+ std::cout << "Derivative norm: ";
+
+ for (int i=0; i<gridDim; i++) {
+ double derivativenorm = 0;
+ for (int j=0; j<4; j++)
+ derivativenorm += derivative[j][i]*derivative[j][i];
+ std::cout << derivativenorm << ", ";
}
- std::cout << "Derivative norm: " << derivativenorm0 << ", " << derivativenorm1 << std::endl;
+ std::cout << std::endl;
#endif
for (int comp=0; comp<4; comp++) {
--
GitLab