From 340c0d8d2d016c16523b2200f4ac1a9cc24e846f Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 9 Feb 2011 17:02:49 +0000 Subject: [PATCH] Write a results section. git-svn-id: https://svn.dealii.org/trunk@23318 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-1/doc/results.dox | 31 +++++++++++++++++++++++++ 1 file changed, 31 insertions(+) diff --git a/deal.II/examples/step-1/doc/results.dox b/deal.II/examples/step-1/doc/results.dox index 74e6f11cef..7528d60d70 100644 --- a/deal.II/examples/step-1/doc/results.dox +++ b/deal.II/examples/step-1/doc/results.dox @@ -18,3 +18,34 @@ like this: The left one, well, is not very exciting. The right one is — at least — unconventional. +While the second mesh is entirely artificial and made-up, and +certainly not very practical in applications, to everyone's surprise it +has found its way into the literature: see the paper by M. Mu +titled "PDE.MART: A network-based problem-solving environment", ACM +Trans. Math. Software, vol. 31, pp. 508-531, 2005. Apparently it is +good for some things at least. + + +

Possible extensions

+ +This program obviously does not have a whole lot of functionality, but +in particular the second_grid function has a bunch of +places where you can play with it. For example, you could modify the +criterion by which we decide which cells to refine. An example would +be to change the condition to this: +@code + for (; cell!=endc; ++cell) + if (cell->center()[1] > 0) + cell->set_refine_flag (); +@endcode +This would refine all cells for which the $y$-coordinate of the cell's +center is greater than zero (the TriaAccessor::center +function that we call by dereferencing the cell iterator +returns a Point<2> object; subscripting [0] would give +the $x$-coordinate, subscripting [1] the +$y$-coordinate). By looking at the functions that TriaAccessor +provides, you can also use more complicated criteria for refinement. + +Another possibility would be to generate meshes of entirely different +geometries altogether. For this, take a look at the GridGenerator +namespace to see what it provides. -- 2.39.5