From a6d39a2a670f56b53b0bc2ed797766cf3c529b8a Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 5 Aug 2011 21:43:22 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-47/step-47.cc | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) 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); } -- 2.39.5