1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
|
from helpers import assert_norm_equality
import dmsh
def generate_pacman(show=False):
geo = dmsh.Difference(
dmsh.Circle([0.0, 0.0], 1.0),
dmsh.Polygon([[0.0, 0.0], [1.5, 0.4], [1.5, -0.4]]),
)
X, cells = dmsh.generate(geo, 0.1, show=show, tol=1.0e-5, max_steps=100)
return X, cells
def test_pacman(show=False):
X, cells = generate_pacman(show)
ref_norms = [3.0173012692535394e02, 1.3565685453257570e01, 9.9999999999884770e-01]
assert_norm_equality(X.flatten(), ref_norms, 1.0e-3)
if __name__ == "__main__":
test_pacman(show=False)
X, cells = generate_pacman(show=True)
# from helpers import save
# save("pacman.png", X, cells)
|