use quadratic in simplex/step-38

diff --git a/tests/simplex/step-38.cc b/tests/simplex/step-38.cc
index 82f919a5c3a04a6be599e84ff03aa1881cdc4f7c..cda26794bb4d3979a26afc940afc6516120f7ae9 100644 (file)
--- 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 9bd139ac6b80e58b9702b99a685d0ce2ebc3f164..90c80ee663981980d1b1f9c94f6690c6c26abb4f 100644 (file)
--- 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