Adjusted settings for default FESeries::Legendre object.

diff --git a/include/deal.II/numerics/smoothness_estimator.h b/include/deal.II/numerics/smoothness_estimator.h
index 7aa7200d62a156bfff60a0704ce2fb8206d5a7f1..b0d11b6531ec89342d58280866b3a939f8740b90 100644 (file)
--- a/include/deal.II/numerics/smoothness_estimator.h
+++ b/include/deal.II/numerics/smoothness_estimator.h
@@ -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>

diff --git a/source/numerics/smoothness_estimator.cc b/source/numerics/smoothness_estimator.cc
index 82c651701c9f480e9f7e069bbdd057e0ccc22310..fa7c0c3c41fcac4af9db8a4e56e4e65036395330 100644 (file)
--- a/source/numerics/smoothness_estimator.cc
+++ b/source/numerics/smoothness_estimator.cc
@@ -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.
       //