Added images

Author: JordiManyer

assets/amr/result.png

@@ -295,14 +295,17 @@ writevtk(new_grid,"half_cylinder_quadratic")
 # 1. Standard non-periodic mesh:
 model = cartesian_model((0,1,0,1),(3,3))
 plot_node_numbering(model)
+# ![](../assets/geometry/nodes_nonperiodic.png)
 
 # 2. Mesh with periodicity in x-direction:
 model = cartesian_model((0,1,0,1),(3,3),isperiodic=(true,false))
 plot_node_numbering(model)
+# ![](../assets/geometry/nodes_halfperiodic.png)
 
 # 3. Mesh with periodicity in both directions:
 model = cartesian_model((0,1,0,1),(3,3),isperiodic=(true,true))
 plot_node_numbering(model)
+# ![](../assets/geometry/nodes_fullperiodic.png)
 
 # Notice how the vertex numbers (displayed at node positions) show the topological
 # connectivity, while the nodes remain at their physical positions.
@@ -361,7 +364,7 @@ cell_vertex_ids = Table(cell_vertex_ids)
 
 # Get coordinates for the vertices:
 vertex_coords = node_coords[vertex_to_node]
-polys = map(get_polytope,reffes)
+polytopes = map(get_polytope,reffes)
 
 # #### Step 4: Create the Model
 #
@@ -385,3 +388,4 @@ mobius = UnstructuredDiscreteModel(grid,topo,labels)
 
 # Visualize the vertex numbering:
 plot_node_numbering(mobius)
+# ![](../assets/geometry/mobius.png)

assets/geometry/mobius.png

@@ -137,8 +137,8 @@ function amr_step(model,u_exact;order=1)
   η = estimate(ηh,uh)
 
   "Mark cells for refinement using Dörfler marking
-  This strategy marks cells containing a fixed fraction (0.8) of the total error"
-  m = DorflerMarking(0.8)
+  This strategy marks cells containing a fixed fraction (0.9) of the total error"
+  m = DorflerMarking(0.9)
   I = Adaptivity.mark(m,η)
 
   "Refine the mesh using newest vertex bisection"
@@ -179,11 +179,13 @@ for i in 1:nsteps
   )
 
   println("Error: $error, Error η: $(sum(η))")
-  @test (i < 3) || (error < last_error)
   last_error = error
   model = fmodel
 end
 
+# The final mesh gives the following result:
+# ![](../assets/amr/result.png)
+#
 # ## Conclusion
 #
 # In this tutorial, we have demonstrated how to: