]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Adjusted settings for default FESeries::Legendre object. 11479/head
authorMarc Fehling <mafehling.git@gmail.com>
Wed, 6 Jan 2021 04:14:59 +0000 (21:14 -0700)
committerMarc Fehling <mafehling.git@gmail.com>
Wed, 6 Jan 2021 04:14:59 +0000 (21:14 -0700)
include/deal.II/numerics/smoothness_estimator.h
source/numerics/smoothness_estimator.cc

index 7aa7200d62a156bfff60a0704ce2fb8206d5a7f1..b0d11b6531ec89342d58280866b3a939f8740b90 100644 (file)
@@ -228,10 +228,10 @@ namespace SmoothnessEstimator
      * the default configuration for smoothness estimation purposes.
      *
      * For each finite element of the provided @p fe_collection, we use as many
-     * modes as its polynomial degree plus one, since we start with the first
-     * Legendre polynomial which is just a constant. Further for each element,
-     * we use a Gaussian quadrature designed to yield exact results for the
-     * highest order Legendre polynomial used.
+     * modes as its polynomial degree plus two. This includes the first Legendre
+     * polynomial which is just a constant. Further for each element, we use a
+     * Gaussian quadrature designed to yield exact results for the highest order
+     * Legendre polynomial used.
      */
     template <int dim, int spacedim>
     FESeries::Legendre<dim, spacedim>
index 82c651701c9f480e9f7e069bbdd057e0ccc22310..fa7c0c3c41fcac4af9db8a4e56e4e65036395330 100644 (file)
@@ -290,9 +290,13 @@ namespace SmoothnessEstimator
     default_fe_series(const hp::FECollection<dim, spacedim> &fe_collection)
     {
       // Default number of coefficients per direction.
+      //
+      // With a number of modes equal to the polynomial degree plus two for each
+      // finite element, the smoothness estimation algorithm tends to produce
+      // stable results.
       std::vector<unsigned int> n_coefficients_per_direction;
       for (unsigned int i = 0; i < fe_collection.size(); ++i)
-        n_coefficients_per_direction.push_back(fe_collection[i].degree + 1);
+        n_coefficients_per_direction.push_back(fe_collection[i].degree + 2);
 
       // Default quadrature collection.
       //

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.