]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Increase degree of quadrature formula since the integrand actually has
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 5 Aug 2011 21:43:22 +0000 (21:43 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 5 Aug 2011 21:43:22 +0000 (21:43 +0000)
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

index e2b282654660d14afad7f558ce5cf37d9d19e918..df97516c7af06e91a8959585a79e0e9ed8d1eb00 100644 (file)
@@ -320,7 +320,7 @@ void LaplaceProblem<dim>::setup_system ()
 template <int dim>
 void LaplaceProblem<dim>::assemble_system ()
 {
-  const QGauss<dim>  quadrature_formula(2);
+  const QGauss<dim>  quadrature_formula(3);
 
 
   FEValues<dim> plain_fe_values (fe_collection[0], quadrature_formula,
@@ -716,7 +716,7 @@ LaplaceProblem<dim>::compute_quadrature (const Quadrature<dim> &plain_quadrature
            }
        }
 
-      Assert (xfem_quadrature.size() == 20, ExcInternalError());
+      Assert (xfem_quadrature.size() == plain_quadrature.size() * 5, ExcInternalError());
       return std::pair<unsigned int, Quadrature<dim> >(2, xfem_quadrature);
     }
 
@@ -811,7 +811,7 @@ LaplaceProblem<dim>::compute_quadrature (const Quadrature<dim> &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<unsigned int, Quadrature<dim> >(3, xfem_quadrature);
     }
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.