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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
import numpy as np
import meshzoo
from .helpers import is_near_equal, signed_simplex_volumes
def test_up():
x_range = np.linspace(0.0, 1.0, 11)
y_range = np.linspace(0.0, 1.0, 11)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="up")
assert len(points) == 121
assert is_near_equal(np.sum(points, axis=0), [60.5, 60.5])
assert len(cells) == 200
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
def test_up2():
x_range = np.linspace(0.0, 1.0, 3)
y_range = np.linspace(0.0, 1.0, 2)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="up")
assert len(points) == 6
assert is_near_equal(np.sum(points, axis=0), [3.0, 3.0])
assert len(cells) == 4
assert list(cells[0]) == [0, 1, 4]
assert list(cells[2]) == [0, 4, 3]
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
def test_zigzag():
x_range = np.linspace(0.0, 1.0, 11)
y_range = np.linspace(0.0, 1.0, 11)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="zigzag")
assert len(points) == 121
assert is_near_equal(np.sum(points, axis=0), [60.5, 60.5])
assert len(cells) == 200
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
def test_zigzag2():
x_range = np.linspace(0.0, 1.0, 2)
y_range = np.linspace(0.0, 1.0, 2)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="zigzag")
assert len(points) == 4
assert is_near_equal(np.sum(points, axis=0), [2.0, 2.0])
assert len(cells) == 2
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
def test_zigzag3():
x_range = np.linspace(0.0, 1.0, 3)
y_range = np.linspace(0.0, 1.0, 2)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="zigzag")
assert len(points) == 6
assert is_near_equal(np.sum(points, axis=0), [3.0, 3.0])
assert len(cells) == 4
assert list(cells[0]) == [0, 1, 4]
assert list(cells[2]) == [1, 2, 4]
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
def test_down():
x_range = np.linspace(0.0, 1.0, 5)
y_range = np.linspace(0.0, 1.0, 4)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="down")
assert len(points) == 20
assert len(cells) == 24
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
# meshzoo.show2d(points, cells)
def test_center():
x_range = np.linspace(0.0, 1.0, 11)
y_range = np.linspace(0.0, 1.0, 9)
points, cells = meshzoo.rectangle_tri(x_range, y_range, variant="center")
meshzoo.show2d(points, cells)
assert len(points) == 99
assert len(cells) == 160
assert np.all(signed_simplex_volumes(points, cells) > 0.0)
def test_quad1():
x_range = np.linspace(0.0, 1.0, 11)
y_range = np.linspace(0.0, 1.0, 11)
points, cells = meshzoo.rectangle_quad(x_range, y_range)
assert len(points) == 121
assert is_near_equal(np.sum(points, axis=0), [60.5, 60.5])
assert len(cells) == 100
# meshzoo.save2d("rectangle-quad.svg", points, cells)
def test_quad2():
x_range = np.linspace(0.0, 1.0, 3)
y_range = np.linspace(0.0, 1.0, 2)
points, cells = meshzoo.rectangle_quad(x_range, y_range)
assert len(points) == 6
assert is_near_equal(np.sum(points, axis=0), [3.0, 3.0])
assert np.all(cells == [[0, 1, 4, 3], [1, 2, 5, 4]])
def test_quad3():
x_range = np.linspace(0.0, 1.0, 3)
y_range = np.linspace(0.0, 1.0, 3)
points, cells = meshzoo.rectangle_quad(x_range, y_range)
assert len(points) == 9
assert np.all(cells == [[0, 1, 4, 3], [3, 4, 7, 6], [1, 2, 5, 4], [4, 5, 8, 7]])
def test_quad8():
x_range = np.linspace(0.0, 1.0, 3)
y_range = np.linspace(0.0, 1.0, 2)
points, cells = meshzoo.rectangle_quad(x_range, y_range, cell_type="quad8")
assert len(points) == 13
print(cells.tolist())
assert np.all(cells == [[0, 1, 4, 3, 6, 11, 8, 10], [1, 2, 5, 4, 7, 12, 9, 11]])
def test_quad9():
x_range = np.linspace(0.0, 1.0, 4)
y_range = np.linspace(0.0, 1.0, 3)
points, cells = meshzoo.rectangle_quad(x_range, y_range, cell_type="quad9")
assert len(points) == 35
print(cells.tolist())
assert np.all(
cells
== [
[0, 1, 5, 4, 12, 22, 15, 21, 29],
[4, 5, 9, 8, 15, 26, 18, 25, 32],
[1, 2, 6, 5, 13, 23, 16, 22, 30],
[5, 6, 10, 9, 16, 27, 19, 26, 33],
[2, 3, 7, 6, 14, 24, 17, 23, 31],
[6, 7, 11, 10, 17, 28, 20, 27, 34],
]
)
|
