From 333eb12b5874169e8427072d8527605eee46c9d4 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 15 Jan 2020 18:21:49 -0700 Subject: [PATCH] Clarify comments about distort_random() in step-49. --- examples/step-49/doc/intro.dox | 15 ++++++++++++--- examples/step-49/step-49.cc | 6 ++++-- 2 files changed, 16 insertions(+), 5 deletions(-) diff --git a/examples/step-49/doc/intro.dox b/examples/step-49/doc/intro.dox index 54aee436cd..cd3d45138b 100644 --- a/examples/step-49/doc/intro.dox +++ b/examples/step-49/doc/intro.dox @@ -209,9 +209,18 @@ demonstrated in grid_7() and the result is as follows: This function is primarily intended to negate some of the superconvergence effects one gets when studying convergence on regular meshes, as well as to suppress some optimizations in deal.II that can exploit the fact that cells -are similar in shape. In practice, it is of course always better to -work with a sequence of unstructured meshes (see possible extensions at the -end of the this section). +are similar in shape. (Superconvergence refers to the fact that if a mesh +has certain symmetries -- for example, if the edges running into a vertex +are symmetric to this vertex, and if this is so for all vertices of a cell +-- that the solution is then often convergent with a higher order than one +would have expected from the usual error analysis. In the end, this +is a result of the fact that if one were to make a Taylor expansion of the +error, the symmetry leads to the fact that the expected next term of the +expansion happens to be zero, and the error order is determined by the +*second next* term. A distorted mesh does not have these symmetries and +consequently the error reflects what one will see when solving the equation +on *any* kind of mesh, rather than showing something that is only reflective +of a particular situation.)

Merging Meshes

diff --git a/examples/step-49/step-49.cc b/examples/step-49/step-49.cc index df07d396f2..ed3cb3d033 100644 --- a/examples/step-49/step-49.cc +++ b/examples/step-49/step-49.cc @@ -296,10 +296,12 @@ void grid_6() // In this last example, we create a mesh and then distort its (interior) // vertices by a random perturbation. This is not something you want to do for -// production computations, but it is a useful tool for testing +// production computations (because results are generally better on meshes +// with "nicely shaped" cells than on the deformed cells produced by +// GridTools::distort_random()), but it is a useful tool for testing // discretizations and codes to make sure they don't work just by accident // because the mesh happens to be uniformly structured and supporting -// super-convergence properties. +// superconvergence properties. void grid_7() { Triangulation<2> triangulation; -- 2.39.5