]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Clarify comments about distort_random() in step-49. 9345/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Thu, 16 Jan 2020 01:21:49 +0000 (18:21 -0700)
committerWolfgang Bangerth <bangerth@colostate.edu>
Thu, 16 Jan 2020 01:21:49 +0000 (18:21 -0700)
examples/step-49/doc/intro.dox
examples/step-49/step-49.cc

index 54aee436cdbb5ac7fceb509e7d3fa240cf95663b..cd3d45138bb85999cea58c19c260135b77de9627 100644 (file)
@@ -209,9 +209,18 @@ demonstrated in <code>grid_7()</code> 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.)
 
 
 <h4>Merging Meshes</h4>
index df07d396f201ab6c289adac7411920513444633e..ed3cb3d033e3d46e45ae4bdd0ccc2057d86436c7 100644 (file)
@@ -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;

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.