From: Wenyu Lei 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 quadrature_formula(2 * fe.degree); + const Simplex::QGauss quadrature_formula(fe.degree + 1); #else const QGauss quadrature_formula(2 * fe.degree); #endif @@ -379,7 +379,7 @@ namespace Step38 Solution(), difference_per_cell, #ifdef USE_SIMPLEX - Simplex::QGauss(2 * fe.degree + 1), + Simplex::QGauss(fe.degree + 1), #else QGauss(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