]> https://gitweb.dealii.org/ - code-gallery.git/commitdiff
Use the correct quadrature formula order.
authorWolfgang Bangerth <bangerth@colostate.edu>
Thu, 28 Mar 2019 15:09:11 +0000 (09:09 -0600)
committerWolfgang Bangerth <bangerth@colostate.edu>
Thu, 28 Mar 2019 15:09:11 +0000 (09:09 -0600)
MCMC-Laplace/mcmc-laplace.cc

index 9d593737721b40cefd2807abd550d09a3ccf1807..a50ca6a98ff6f08e3cc506ce58f26ae587978eb9 100644 (file)
@@ -225,7 +225,7 @@ namespace ForwardSimulator
       cell_matrix.reinit(dofs_per_cell, dofs_per_cell);
       cell_rhs.reinit(dofs_per_cell);
 
-      const QGauss<dim>  quadrature_formula(2);
+      const QGauss<dim>  quadrature_formula(fe.degree+1);
       const unsigned int n_q_points = quadrature_formula.size();
 
       FEValues<dim> fe_values(fe,

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.