From: bangerth Date: Fri, 5 Aug 2011 21:43:22 +0000 (+0000) Subject: Increase degree of quadrature formula since the integrand actually has X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a6d39a2a670f56b53b0bc2ed797766cf3c529b8a;p=dealii-svn.git Increase degree of quadrature formula since the integrand actually has higher polynomial degree. git-svn-id: https://svn.dealii.org/trunk@24025 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-47/step-47.cc b/deal.II/examples/step-47/step-47.cc index e2b2826546..df97516c7a 100644 --- a/deal.II/examples/step-47/step-47.cc +++ b/deal.II/examples/step-47/step-47.cc @@ -320,7 +320,7 @@ void LaplaceProblem::setup_system () template void LaplaceProblem::assemble_system () { - const QGauss quadrature_formula(2); + const QGauss quadrature_formula(3); FEValues plain_fe_values (fe_collection[0], quadrature_formula, @@ -716,7 +716,7 @@ LaplaceProblem::compute_quadrature (const Quadrature &plain_quadrature } } - Assert (xfem_quadrature.size() == 20, ExcInternalError()); + Assert (xfem_quadrature.size() == plain_quadrature.size() * 5, ExcInternalError()); return std::pair >(2, xfem_quadrature); } @@ -811,7 +811,7 @@ LaplaceProblem::compute_quadrature (const Quadrature &plain_quadrature xfem_quadrature.initialize(xfem_points, xfem_weights); } } - Assert (xfem_quadrature.size() == 8, ExcInternalError()); + Assert (xfem_quadrature.size() == plain_quadrature.size() * 2, ExcInternalError()); return std::pair >(3, xfem_quadrature); }