diff --git a/test/averagedistanceassemblertest.cc b/test/averagedistanceassemblertest.cc index 682554942c34dbc1aa0920b549db04387883986a..6edde06a820b657cecb6e4a0334057475d26196d 100644 --- a/test/averagedistanceassemblertest.cc +++ b/test/averagedistanceassemblertest.cc @@ -37,10 +37,10 @@ void testPoint(const std::vector<TargetSpace>& corners, // test the 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.assembleHessianApproximation(argument, hessianApproximation); + //assembler.assembleHessianApproximation(argument, hessianApproximation); for (size_t i=0; i<hessian.N(); i++) for (size_t j=0; j<hessian.M(); j++) { @@ -48,6 +48,9 @@ void testPoint(const std::vector<TargetSpace>& corners, 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; diff -= hessianApproximation; std::cout << "Matrix inf diff: " << diff.infinity_norm() << std::endl;