From: Wenyu Lei <wlei@sissa.it>
Date: Mon, 18 Jan 2021 18:55:14 +0000 (+0100)
Subject: use quadratic in simplex/step-38
X-Git-Tag: v9.3.0-rc1~594^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F11579%2Fhead;p=dealii.git

use quadratic in simplex/step-38
---

diff --git a/tests/simplex/step-38.cc b/tests/simplex/step-38.cc
index 82f919a5c3..cda26794bb 100644
--- a/tests/simplex/step-38.cc
+++ b/tests/simplex/step-38.cc
@@ -63,7 +63,7 @@ namespace Step38
   class LaplaceBeltramiProblem
   {
   public:
-    LaplaceBeltramiProblem(const unsigned degree = 1);
+    LaplaceBeltramiProblem(const unsigned degree = 2);
     void
     run();
 
@@ -274,7 +274,7 @@ namespace Step38
     system_matrix = 0;
     system_rhs    = 0;
 #ifdef USE_SIMPLEX
-    const Simplex::QGauss<dim> quadrature_formula(2 * fe.degree);
+    const Simplex::QGauss<dim> quadrature_formula(fe.degree + 1);
 #else
     const QGauss<dim> quadrature_formula(2 * fe.degree);
 #endif
@@ -379,7 +379,7 @@ namespace Step38
                                       Solution<spacedim>(),
                                       difference_per_cell,
 #ifdef USE_SIMPLEX
-                                      Simplex::QGauss<dim>(2 * fe.degree + 1),
+                                      Simplex::QGauss<dim>(fe.degree + 1),
 #else
                                       QGauss<dim>(2 * fe.degree +
                                                   1), // This also works on
@@ -412,7 +412,7 @@ main()
     {
       using namespace Step38;
 
-      LaplaceBeltramiProblem<3> laplace_beltrami(1);
+      LaplaceBeltramiProblem<3> laplace_beltrami;
       laplace_beltrami.run();
     }
   catch (std::exception &exc)
diff --git a/tests/simplex/step-38.with_simplex_support=on.output b/tests/simplex/step-38.with_simplex_support=on.output
index 9bd139ac6b..90c80ee663 100644
--- a/tests/simplex/step-38.with_simplex_support=on.output
+++ b/tests/simplex/step-38.with_simplex_support=on.output
@@ -1,3 +1,3 @@
 Surface mesh has 10240 cells.
-Surface mesh has 5185 degrees of freedom.
-H1 error = 0.384717
+Surface mesh has 20609 degrees of freedom.
+H1 error = 0.00932553