From 48796666758cc2eaffc908e0545809f19fb2c9f2 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 22 Apr 2010 14:39:30 +0000
Subject: [PATCH] start a test for the Hessian of the harmonic energy
[[Imported from SVN: r5928]]
---
test/harmonicenergytest.cc | 11 +++++++++++
1 file changed, 11 insertions(+)
diff --git a/test/harmonicenergytest.cc b/test/harmonicenergytest.cc
index 402b86e6..7a15f469 100644
--- a/test/harmonicenergytest.cc
+++ b/test/harmonicenergytest.cc
@@ -85,4 +85,15 @@ int main(int argc, char** argv)
coefficients[2] = uv;
testEnergy<GridType>(grid, coefficients);
+
+ // //////////////////////////////////////////////////////////
+ // Test the approximation to the Hesse matrix
+ // //////////////////////////////////////////////////////////
+
+ HarmonicEnergyLocalStiffness<GridType::LeafGridView,TargetSpace> assembler;
+
+ assembler.assemble(*grid->leafbegin<0>(), coefficients);
+
+ std::cout << "Hessian: \n" << assembler.A[0][0] << std::endl;
+
}
--
GitLab