@@ -134,79 +134,83 @@ def test_no_modify():
134134
135135def test_trifinder ():
136136 # Test points within triangles of masked triangulation.
137- x ,y = np .meshgrid (np .arange (4 ), np .arange (4 ))
137+ x , y = np .meshgrid (np .arange (4 ), np .arange (4 ))
138138 x = x .ravel ()
139139 y = y .ravel ()
140- triangles = [[0 ,1 ,4 ], [1 ,5 ,4 ], [1 ,2 ,5 ], [2 ,6 ,5 ], [2 ,3 ,6 ], [3 ,7 ,6 ], [4 ,5 ,8 ],
141- [5 ,9 ,8 ], [5 ,6 ,9 ], [6 ,10 ,9 ], [6 ,7 ,10 ], [7 ,11 ,10 ], [8 ,9 ,12 ],
142- [9 ,13 ,12 ], [9 ,10 ,13 ], [10 ,14 ,13 ], [10 ,11 ,14 ], [11 ,15 ,14 ]]
140+ triangles = [[0 , 1 , 4 ], [1 , 5 , 4 ], [1 , 2 , 5 ], [2 , 6 , 5 ], [2 , 3 , 6 ],
141+ [3 , 7 , 6 ], [4 , 5 , 8 ], [5 , 9 , 8 ], [5 , 6 , 9 ], [6 , 10 , 9 ],
142+ [6 , 7 , 10 ], [7 , 11 , 10 ], [8 , 9 , 12 ], [9 , 13 , 12 ], [9 , 10 , 13 ],
143+ [10 , 14 , 13 ], [10 , 11 , 14 ], [11 , 15 , 14 ]]
143144 mask = np .zeros (len (triangles ))
144145 mask [8 :10 ] = 1
145146 triang = mtri .Triangulation (x , y , triangles , mask )
146147 trifinder = triang .get_trifinder ()
147148
148149 xs = [0.25 , 1.25 , 2.25 , 3.25 ]
149150 ys = [0.25 , 1.25 , 2.25 , 3.25 ]
150- xs ,ys = np .meshgrid (xs ,ys )
151+ xs , ys = np .meshgrid (xs , ys )
151152 xs = xs .ravel ()
152153 ys = ys .ravel ()
153154 tris = trifinder (xs , ys )
154- assert_array_equal (tris , [0 ,2 ,4 ,- 1 ,6 ,- 1 ,10 ,- 1 ,12 ,14 ,16 ,- 1 ,- 1 ,- 1 ,- 1 ,- 1 ])
155+ assert_array_equal (tris , [0 , 2 , 4 , - 1 , 6 , - 1 , 10 , - 1 ,
156+ 12 , 14 , 16 , - 1 , - 1 , - 1 , - 1 , - 1 ])
155157 tris = trifinder (xs - 0.5 , ys - 0.5 )
156- assert_array_equal (tris , [- 1 ,- 1 ,- 1 ,- 1 ,- 1 ,1 ,3 ,5 ,- 1 ,7 ,- 1 ,11 ,- 1 ,13 ,15 ,17 ])
158+ assert_array_equal (tris , [- 1 , - 1 , - 1 , - 1 , - 1 , 1 , 3 , 5 ,
159+ - 1 , 7 , - 1 , 11 , - 1 , 13 , 15 , 17 ])
157160
158161 # Test points exactly on boundary edges of masked triangulation.
159- xs = [0.5 , 1.5 , 2.5 , 0.5 , 1.5 , 2.5 , 1.5 , 1.5 , 0.0 , 1.0 , 2.0 , 3.0 ]
160- ys = [0.0 , 0.0 , 0.0 , 3.0 , 3.0 , 3.0 , 1.0 , 2.0 , 1.5 , 1.5 , 1.5 , 1.5 ]
162+ xs = [0.5 , 1.5 , 2.5 , 0.5 , 1.5 , 2.5 , 1.5 , 1.5 , 0.0 , 1.0 , 2.0 , 3.0 ]
163+ ys = [0.0 , 0.0 , 0.0 , 3.0 , 3.0 , 3.0 , 1.0 , 2.0 , 1.5 , 1.5 , 1.5 , 1.5 ]
161164 tris = trifinder (xs , ys )
162- assert_array_equal (tris , [0 ,2 , 4 , 13 ,15 ,17 ,3 , 14 ,6 , 7 , 10 ,11 ])
165+ assert_array_equal (tris , [0 , 2 , 4 , 13 , 15 , 17 , 3 , 14 , 6 , 7 , 10 , 11 ])
163166
164167 # Test points exactly on boundary corners of masked triangulation.
165168 xs = [0.0 , 3.0 ]
166169 ys = [0.0 , 3.0 ]
167170 tris = trifinder (xs , ys )
168- assert_array_equal (tris , [0 ,17 ])
171+ assert_array_equal (tris , [0 , 17 ])
169172
170173 # Test triangles with horizontal colinear points. These are not valid
171174 # triangulations, but we try to deal with the simplest violations.
172175 delta = 0.0 # If +ve, triangulation is OK, if -ve triangulation invalid,
173176 # if zero have colinear points but should pass tests anyway.
174177 x = [1.5 , 0 , 1 , 2 , 3 , 1.5 , 1.5 ]
175178 y = [- 1 , 0 , 0 , 0 , 0 , delta , 1 ]
176- triangles = [[0 ,2 , 1 ], [0 ,3 , 2 ], [0 ,4 , 3 ], [1 ,2 , 5 ], [2 ,3 , 5 ], [ 3 , 4 , 5 ], [ 1 , 5 , 6 ],
177- [4 , 6 , 5 ]]
179+ triangles = [[0 , 2 , 1 ], [0 , 3 , 2 ], [0 , 4 , 3 ], [1 , 2 , 5 ], [2 , 3 , 5 ],
180+ [3 , 4 , 5 ], [ 1 , 5 , 6 ], [ 4 , 6 , 5 ]]
178181 triang = mtri .Triangulation (x , y , triangles )
179182 trifinder = triang .get_trifinder ()
180183
181184 xs = [- 0.1 , 0.4 , 0.9 , 1.4 , 1.9 , 2.4 , 2.9 ]
182- ys = [- 0.1 ,0.1 ]
183- xs ,ys = np .meshgrid (xs , ys )
185+ ys = [- 0.1 , 0.1 ]
186+ xs , ys = np .meshgrid (xs , ys )
184187 tris = trifinder (xs , ys )
185- assert_array_equal (tris , [[- 1 ,0 ,0 ,1 ,1 ,2 ,- 1 ],[- 1 ,6 ,6 ,6 ,7 ,7 ,- 1 ]])
188+ assert_array_equal (tris , [[- 1 , 0 , 0 , 1 , 1 , 2 , - 1 ],
189+ [- 1 , 6 , 6 , 6 , 7 , 7 , - 1 ]])
186190
187191 # Test triangles with vertical colinear points. These are not valid
188192 # triangulations, but we try to deal with the simplest violations.
189193 delta = 0.0 # If +ve, triangulation is OK, if -ve triangulation invalid,
190194 # if zero have colinear points but should pass tests anyway.
191195 x = [- 1 , - delta , 0 , 0 , 0 , 0 , 1 ]
192196 y = [1.5 , 1.5 , 0 , 1 , 2 , 3 , 1.5 ]
193- triangles = [[0 ,1 , 2 ], [0 ,1 , 5 ], [1 ,2 , 3 ], [1 ,3 , 4 ], [1 ,4 , 5 ], [ 2 , 6 , 3 ], [ 3 , 6 , 4 ],
194- [4 , 6 , 5 ]]
197+ triangles = [[0 , 1 , 2 ], [0 , 1 , 5 ], [1 , 2 , 3 ], [1 , 3 , 4 ], [1 , 4 , 5 ],
198+ [2 , 6 , 3 ], [ 3 , 6 , 4 ], [ 4 , 6 , 5 ]]
195199 triang = mtri .Triangulation (x , y , triangles )
196200 trifinder = triang .get_trifinder ()
197201
198- xs = [- 0.1 ,0.1 ]
202+ xs = [- 0.1 , 0.1 ]
199203 ys = [- 0.1 , 0.4 , 0.9 , 1.4 , 1.9 , 2.4 , 2.9 ]
200- xs ,ys = np .meshgrid (xs , ys )
204+ xs , ys = np .meshgrid (xs , ys )
201205 tris = trifinder (xs , ys )
202- assert_array_equal (tris , [[- 1 ,- 1 ], [0 ,5 ], [0 ,5 ], [0 ,6 ], [1 ,6 ], [1 ,7 ],
203- [- 1 ,- 1 ]])
206+ assert_array_equal (tris , [[- 1 , - 1 ], [0 , 5 ], [0 , 5 ], [0 , 6 ], [1 , 6 ], [1 , 7 ],
207+ [- 1 , - 1 ]])
204208
205209 # Test that changing triangulation by setting a mask causes the trifinder
206210 # to be reinitialised.
207211 x = [0 , 1 , 0 , 1 ]
208212 y = [0 , 0 , 1 , 1 ]
209- triangles = [[0 ,1 , 2 ], [1 ,3 , 2 ]]
213+ triangles = [[0 , 1 , 2 ], [1 , 3 , 2 ]]
210214 triang = mtri .Triangulation (x , y , triangles )
211215 trifinder = triang .get_trifinder ()
212216
@@ -222,21 +226,22 @@ def test_trifinder():
222226
223227def test_triinterp ():
224228 # Test points within triangles of masked triangulation.
225- x ,y = np .meshgrid (np .arange (4 ), np .arange (4 ))
229+ x , y = np .meshgrid (np .arange (4 ), np .arange (4 ))
226230 x = x .ravel ()
227231 y = y .ravel ()
228232 z = 1.23 * x - 4.79 * y
229- triangles = [[0 ,1 ,4 ], [1 ,5 ,4 ], [1 ,2 ,5 ], [2 ,6 ,5 ], [2 ,3 ,6 ], [3 ,7 ,6 ], [4 ,5 ,8 ],
230- [5 ,9 ,8 ], [5 ,6 ,9 ], [6 ,10 ,9 ], [6 ,7 ,10 ], [7 ,11 ,10 ], [8 ,9 ,12 ],
231- [9 ,13 ,12 ], [9 ,10 ,13 ], [10 ,14 ,13 ], [10 ,11 ,14 ], [11 ,15 ,14 ]]
233+ triangles = [[0 , 1 , 4 ], [1 , 5 , 4 ], [1 , 2 , 5 ], [2 , 6 , 5 ], [2 , 3 , 6 ],
234+ [3 , 7 , 6 ], [4 , 5 , 8 ], [5 , 9 , 8 ], [5 , 6 , 9 ], [6 , 10 , 9 ],
235+ [6 , 7 , 10 ], [7 , 11 , 10 ], [8 , 9 , 12 ], [9 , 13 , 12 ], [9 , 10 , 13 ],
236+ [10 , 14 , 13 ], [10 , 11 , 14 ], [11 , 15 , 14 ]]
232237 mask = np .zeros (len (triangles ))
233238 mask [8 :10 ] = 1
234239 triang = mtri .Triangulation (x , y , triangles , mask )
235240 linear_interp = mtri .LinearTriInterpolator (triang , z )
236241
237242 xs = np .linspace (0.25 , 2.75 , 6 )
238243 ys = [0.25 , 0.75 , 2.25 , 2.75 ]
239- xs ,ys = np .meshgrid (xs ,ys )
244+ xs , ys = np .meshgrid (xs , ys )
240245 xs = xs .ravel ()
241246 ys = ys .ravel ()
242247 zs = linear_interp (xs , ys )
@@ -245,7 +250,7 @@ def test_triinterp():
245250 # Test points outside triangulation.
246251 xs = [- 0.25 , 1.25 , 1.75 , 3.25 ]
247252 ys = xs
248- xs , ys = np .meshgrid (xs ,ys )
253+ xs , ys = np .meshgrid (xs , ys )
249254 xs = xs .ravel ()
250255 ys = ys .ravel ()
251256 zs = linear_interp (xs , ys )
0 commit comments