Skip to content
Snippets Groups Projects
Commit 4d5060f3 authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

comment out the call to method hessianApproximation, which doesn't exist anymore

[[Imported from SVN: r8033]]
parent 3fa30a78
No related branches found
No related tags found
No related merge requests found
...@@ -37,10 +37,10 @@ void testPoint(const std::vector<TargetSpace>& corners, ...@@ -37,10 +37,10 @@ void testPoint(const std::vector<TargetSpace>& corners,
// test the hessian // test the hessian
FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> hessian; FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> hessian;
FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> hessianApproximation; FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> hessianApproximation(0);
assembler.assembleHessian(argument, hessian); assembler.assembleHessian(argument, hessian);
assembler.assembleHessianApproximation(argument, hessianApproximation); //assembler.assembleHessianApproximation(argument, hessianApproximation);
for (size_t i=0; i<hessian.N(); i++) for (size_t i=0; i<hessian.N(); i++)
for (size_t j=0; j<hessian.M(); j++) { for (size_t j=0; j<hessian.M(); j++) {
...@@ -48,6 +48,9 @@ void testPoint(const std::vector<TargetSpace>& corners, ...@@ -48,6 +48,9 @@ void testPoint(const std::vector<TargetSpace>& corners,
assert(!std::isnan(hessianApproximation[i][j])); assert(!std::isnan(hessianApproximation[i][j]));
} }
std::cout << "WARNING: no approximation of the Hessian available, not testing" << std::endl;
return;
FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> diff = hessian; FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> diff = hessian;
diff -= hessianApproximation; diff -= hessianApproximation;
std::cout << "Matrix inf diff: " << diff.infinity_norm() << std::endl; std::cout << "Matrix inf diff: " << diff.infinity_norm() << std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment