]> https://gitweb.dealii.org/ - dealii.git/commitdiff
use fe_degree+2 quadrature points in matrix_free/mass_operator_0x 3393/head
authorDenis Davydov <davydden@gmail.com>
Thu, 3 Nov 2016 13:54:34 +0000 (14:54 +0100)
committerDenis Davydov <davydden@gmail.com>
Thu, 3 Nov 2016 13:54:34 +0000 (14:54 +0100)
tests/matrix_free/mass_operator_01.cc
tests/matrix_free/mass_operator_02.cc

index 2e95bd51a95e85a78119ff1ebb8d33a407150e60..4cb4d31ba9d28adfd26fef53b10e38c6b2ff7794 100644 (file)
@@ -107,7 +107,7 @@ void test ()
 
   MatrixFree<dim,number> mf_data;
   {
-    const QGauss<1> quad (fe_degree+1);
+    const QGauss<1> quad (fe_degree+2);
     typename MatrixFree<dim,number>::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<dim,fe_degree, fe_degree+1, 1, number> mf;
+  MatrixFreeOperators::MassOperator<dim,fe_degree, fe_degree+2, 1, number> mf;
   mf.initialize(mf_data);
   mf.compute_diagonal();
   LinearAlgebra::distributed::Vector<number> in, out, ref;
@@ -147,7 +147,7 @@ void test ()
     sparse_matrix.reinit (csp);
   }
   {
-    QGauss<dim>  quadrature_formula(fe_degree+1);
+    QGauss<dim>  quadrature_formula(fe_degree+2);
 
     FEValues<dim> fe_values (dof.get_fe(), quadrature_formula,
                              update_values    |  update_gradients |
index b7db2771e1f1f3401e23d0eb6daf2725a2258155..b082200bd59aba44a30749d3aca7f8427c2dfedf 100644 (file)
@@ -71,7 +71,7 @@ void test ()
 
   MatrixFree<dim,number> mf_data;
   {
-    const QGauss<1> quad (fe_degree+1);
+    const QGauss<1> quad (fe_degree+2);
     typename MatrixFree<dim,number>::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<dim,fe_degree, fe_degree+1, 1, number> mf;
+  MatrixFreeOperators::MassOperator<dim,fe_degree, fe_degree+2, 1, number> mf;
   mf.initialize(mf_data);
   mf.compute_diagonal();
   const LinearAlgebra::distributed::Vector<double> &diagonal
@@ -114,7 +114,7 @@ void test ()
     sparse_matrix.reinit (csp);
   }
   {
-    QGauss<dim>  quadrature_formula(fe_degree+1);
+    QGauss<dim>  quadrature_formula(fe_degree+2);
 
     FEValues<dim> fe_values (dof.get_fe(), quadrature_formula,
                              update_values    |  update_gradients |

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.