1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
import meshzoo
import pytest
import meshplex
def test_euler_characteristic():
points = [[0.0, 0.0], [1.0, 0.0], [0.0, 1.0]]
cells = [[0, 1, 2]]
mesh = meshplex.MeshTri(points, cells)
assert mesh.euler_characteristic == 1
assert mesh.genus == 0
points, cells = meshzoo.icosa_sphere(5)
mesh = meshplex.MeshTri(points, cells)
assert mesh.euler_characteristic == 2
assert mesh.genus == 0
points, cells = meshzoo.moebius(num_twists=1, nl=21, nw=6)
mesh = meshplex.MeshTri(points, cells)
assert mesh.euler_characteristic == 0
with pytest.raises(RuntimeError):
mesh.genus
|
