From 56bb3d1f76583b2fd3e6dd3f0a012ab18df4c72c Mon Sep 17 00:00:00 2001 From: Denis Davydov Date: Thu, 3 Nov 2016 14:54:34 +0100 Subject: [PATCH] use fe_degree+2 quadrature points in matrix_free/mass_operator_0x --- tests/matrix_free/mass_operator_01.cc | 6 +++--- tests/matrix_free/mass_operator_02.cc | 6 +++--- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/tests/matrix_free/mass_operator_01.cc b/tests/matrix_free/mass_operator_01.cc index 2e95bd51a9..4cb4d31ba9 100644 --- a/tests/matrix_free/mass_operator_01.cc +++ b/tests/matrix_free/mass_operator_01.cc @@ -107,7 +107,7 @@ void test () MatrixFree mf_data; { - const QGauss<1> quad (fe_degree+1); + const QGauss<1> quad (fe_degree+2); typename MatrixFree::AdditionalData data; data.mpi_communicator = MPI_COMM_WORLD; data.tasks_parallel_scheme = @@ -116,7 +116,7 @@ void test () mf_data.reinit (dof, constraints, quad, data); } - MatrixFreeOperators::MassOperator mf; + MatrixFreeOperators::MassOperator mf; mf.initialize(mf_data); mf.compute_diagonal(); LinearAlgebra::distributed::Vector in, out, ref; @@ -147,7 +147,7 @@ void test () sparse_matrix.reinit (csp); } { - QGauss quadrature_formula(fe_degree+1); + QGauss quadrature_formula(fe_degree+2); FEValues fe_values (dof.get_fe(), quadrature_formula, update_values | update_gradients | diff --git a/tests/matrix_free/mass_operator_02.cc b/tests/matrix_free/mass_operator_02.cc index b7db2771e1..b082200bd5 100644 --- a/tests/matrix_free/mass_operator_02.cc +++ b/tests/matrix_free/mass_operator_02.cc @@ -71,7 +71,7 @@ void test () MatrixFree mf_data; { - const QGauss<1> quad (fe_degree+1); + const QGauss<1> quad (fe_degree+2); typename MatrixFree::AdditionalData data; data.mpi_communicator = MPI_COMM_WORLD; data.tasks_parallel_scheme = @@ -80,7 +80,7 @@ void test () mf_data.reinit (dof, constraints, quad, data); } - MatrixFreeOperators::MassOperator mf; + MatrixFreeOperators::MassOperator mf; mf.initialize(mf_data); mf.compute_diagonal(); const LinearAlgebra::distributed::Vector &diagonal @@ -114,7 +114,7 @@ void test () sparse_matrix.reinit (csp); } { - QGauss quadrature_formula(fe_degree+1); + QGauss quadrature_formula(fe_degree+2); FEValues fe_values (dof.get_fe(), quadrature_formula, update_values | update_gradients | -- 2.39.5