Add diagnostic plots

guyer · Apr 13, 2021 · 58feb9f · 58feb9f
1 parent 5ef7823
commit 58feb9f
Showing 1 changed file with 15 additions and 0 deletions.
diff --git a/examples/diffusion/circleQuad.py b/examples/diffusion/circleQuad.py
@@ -159,6 +159,21 @@
 >>> print(phi.allclose(x, atol = 0.035)) # doctest: +GMSH
 1
 
+>>> from fipy.meshes.gmshMesh import gmshVersion
+>>> version = gmshVersion()
+
+>>> Viewer(vars=CellVariable(mesh=mesh,
+...                          name="error - Gmsh {}".format(version),
+...                          value=phi - x)).plot(filename="error_gmsh{}.png".format(version))
+
+>>> Viewer(vars=CellVariable(mesh=mesh,
+...                          name="volume - Gmsh {}".format(version),
+...                          value=mesh.cellVolumes)).plot(filename="volume_gmsh{}.png".format(version))
+
+>>> Viewer(vars=CellVariable(name = "non-orthogonality - Gmsh {}".format(version),
+...                          mesh = mesh,
+...                          value = mesh._nonOrthogonality)).plot(filename="nonorthogonality_gmsh{}.png".format(version))
+
 Display the results if run as a script.
 
 >>> from fipy import input