From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Thu, 28 Mar 2019 15:09:11 +0000 (-0600)
Subject: Use the correct quadrature formula order.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=27ed3f3a89764cd54f7f8565304f8705e5efb4a7;p=code-gallery.git

Use the correct quadrature formula order.
---

diff --git a/MCMC-Laplace/mcmc-laplace.cc b/MCMC-Laplace/mcmc-laplace.cc
index 9d59373..a50ca6a 100644
--- a/MCMC-Laplace/mcmc-laplace.cc
+++ b/MCMC-Laplace/mcmc-laplace.cc
@@ -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,