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
|