From 27ed3f3a89764cd54f7f8565304f8705e5efb4a7 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Thu, 28 Mar 2019 09:09:11 -0600 Subject: [PATCH] Use the correct quadrature formula order. --- MCMC-Laplace/mcmc-laplace.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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 quadrature_formula(2); + const QGauss quadrature_formula(fe.degree+1); const unsigned int n_q_points = quadrature_formula.size(); FEValues fe_values(fe, -- 2.39.5