use MeshTri.init_circle in exx; fixes #476

kinnala · Sep 16, 2020 · 9eae8b9 · 9eae8b9
1 parent e292734
commit 9eae8b9
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 7 deletions.
diff --git a/docs/examples/ex12.py b/docs/examples/ex12.py
@@ -25,10 +25,9 @@
 
 import numpy as np
 
-m = from_file(Path(__file__).with_name("disk.json"))
-m.scale(1/np.linalg.norm(m.p, axis=0).max())  # unit radius
+m = MeshTri.init_circle(4)
 
-basis = InteriorBasis(m, ElementQuad2())
+basis = InteriorBasis(m, ElementTriP2())
 
 A = asm(laplace, basis)
 b = asm(unit_load, basis)

diff --git a/docs/examples/ex18.py b/docs/examples/ex18.py
@@ -51,8 +51,7 @@
 from scipy.sparse import bmat
 
 
-mesh = from_file(Path(__file__).with_name("disk.json")).to_meshtri()
-mesh.scale(1/np.linalg.norm(mesh.p, axis=0).max())  # unit radius
+mesh = MeshTri.init_circle(4)
 
 element = {'u': ElementVectorH1(ElementTriP2()),
            'p': ElementTriP1()}

diff --git a/docs/examples/ex20.py b/docs/examples/ex20.py
@@ -35,8 +35,7 @@
 import numpy as np
 
 
-mesh = from_file(Path(__file__).with_name("disk.json")).to_meshtri()
-mesh.scale(1/np.linalg.norm(mesh.p, axis=0).max())  # unit radius
+mesh = MeshTri.init_circle(4)
 element = ElementTriMorley()
 mapping = MappingAffine(mesh)
 ib = InteriorBasis(mesh, element, mapping, 2)