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

Use FieldVector::dimension instead of FieldVector::size

[[Imported from SVN: r8026]]
parent 9d7410b6
No related branches found
No related tags found
No related merge requests found
......@@ -36,8 +36,8 @@ void testPoint(const std::vector<TargetSpace>& corners,
assembler.assembleGradient(argument, gradient);
// test the hessian
FieldMatrix<double, TargetSpace::TangentVector::size, TargetSpace::TangentVector::size> hessian;
FieldMatrix<double, TargetSpace::TangentVector::size, TargetSpace::TangentVector::size> hessianApproximation;
FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> hessian;
FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> hessianApproximation;
assembler.assembleHessian(argument, hessian);
assembler.assembleHessianApproximation(argument, hessianApproximation);
......@@ -48,7 +48,7 @@ void testPoint(const std::vector<TargetSpace>& corners,
assert(!std::isnan(hessianApproximation[i][j]));
}
FieldMatrix<double, TargetSpace::TangentVector::size, TargetSpace::TangentVector::size> diff = hessian;
FieldMatrix<double, TargetSpace::TangentVector::dimension, TargetSpace::TangentVector::dimension> diff = hessian;
diff -= hessianApproximation;
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