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
|
import numpy as np
from helpers import assert_norm_equality, save
import dmsh
def generate_halfspace(show=False):
geo = dmsh.Intersection(
[
dmsh.HalfSpace(np.sqrt(0.5) * np.array([1.0, 1.0])),
dmsh.Circle([0.0, 0.0], 1.0),
]
)
X, cells = dmsh.generate(geo, 0.1, show=show, max_steps=100)
return X, cells
def test_halfspace(show=False):
X, cells = generate_halfspace(show)
ref_norms = [1.6399670188761661e02, 1.0011048291798387e01, 9.9959986881486440e-01]
assert_norm_equality(X.flatten(), ref_norms, 1.0e-6)
if __name__ == "__main__":
X, cells = generate_halfspace(show=True)
save("halfspace.png", X, cells)
|
