@@ -948,6 +948,32 @@ def test_trirefine():
948948 assert_array_almost_equal (xyz_data [0 ], xyz_data [1 ])
949949
950950
951+ def test_cubic_trirefine_masked ():
952+ np .random .seed (19680801 )
953+ x = np .zeros (1000 )
954+ for i in range (10 ):
955+ x [500 + i ] = np .random .rand ()
956+ y = np .repeat (np .arange (100 ), 10 )
957+ z = np .random .rand (1000 )
958+ tri = mtri .Triangulation (x , y )
959+ refiner = mtri .UniformTriRefiner (tri )
960+ cinterp = mtri .CubicTriInterpolator (tri , z )
961+ refiner .refine_field (z , triinterpolator = cinterp , subdiv = 2 )
962+
963+
964+ def test_linear_trirefine_masked ():
965+ np .random .seed (19680801 )
966+ x = np .zeros (1000 )
967+ for i in range (10 ):
968+ x [500 + i ] = np .random .rand ()
969+ y = np .repeat (np .arange (100 ), 10 )
970+ z = np .random .rand (1000 )
971+ tri = mtri .Triangulation (x , y )
972+ refiner = mtri .UniformTriRefiner (tri )
973+ linterp = mtri .LinearTriInterpolator (tri , z )
974+ refiner .refine_field (z , triinterpolator = linterp , subdiv = 2 )
975+
976+
951977def meshgrid_triangles (n ):
952978 """
953979 Utility function.
0 commit comments