Changeset 6410 for branches/numpy/anuga/utilities/test_polygon.py

Timestamp:

Feb 25, 2009, 9:37:22 AM (15 years ago)

Author:

rwilson

Message:

numpy changes.

File:

• branches/numpy/anuga/utilities/test_polygon.py (modified) (1 diff)

branches/numpy/anuga/utilities/test_polygon.py

@@ -1703 +1703 @@
         assert num.allclose(z, z_ref)
 
+    def test_inside_polygon_badtest_1(self):
+        '''Test failure found in a larger test in shallow_water.
+        (test_get_maximum_inundation in test_data_manager.py (line 10574))
+       
+        Data below, when fed into:
+            indices = inside_polygon(points, polygon)
+        returns indices == [], and it shouldn't.
+
+        However, it works fine here!?!?
+        '''
+
+        polygon = [[50, 1], [99, 1], [99, 49], [50, 49]]
+
+        points = [[  0.,  0.], [  0., 10.], [  0., 20.], [  0., 30.],
+                  [  0., 40.], [  0., 50.], [  5.,  0.], [  5., 10.],
+                  [  5., 20.], [  5., 30.], [  5., 40.], [  5., 50.],
+                  [ 10.,  0.], [ 10., 10.], [ 10., 20.], [ 10., 30.],
+                  [ 10., 40.], [ 10., 50.], [ 15.,  0.], [ 15., 10.],
+                  [ 15., 20.], [ 15., 30.], [ 15., 40.], [ 15., 50.],
+                  [ 20.,  0.], [ 20., 10.], [ 20., 20.], [ 20., 30.],
+                  [ 20., 40.], [ 20., 50.], [ 25.,  0.], [ 25., 10.],
+                  [ 25., 20.], [ 25., 30.], [ 25., 40.], [ 25., 50.],
+                  [ 30.,  0.], [ 30., 10.], [ 30., 20.], [ 30., 30.],
+                  [ 30., 40.], [ 30., 50.], [ 35.,  0.], [ 35., 10.],
+                  [ 35., 20.], [ 35., 30.], [ 35., 40.], [ 35., 50.],
+                  [ 40.,  0.], [ 40., 10.], [ 40., 20.], [ 40., 30.],
+                  [ 40., 40.], [ 40., 50.], [ 45.,  0.], [ 45., 10.],
+                  [ 45., 20.], [ 45., 30.], [ 45., 40.], [ 45., 50.],
+                  [ 50.,  0.], [ 50., 10.], [ 50., 20.], [ 50., 30.],
+                  [ 50., 40.], [ 50., 50.], [ 55.,  0.], [ 55., 10.],
+                  [ 55., 20.], [ 55., 30.], [ 55., 40.], [ 55., 50.],
+                  [ 60.,  0.], [ 60., 10.], [ 60., 20.], [ 60., 30.],
+                  [ 60., 40.], [ 60., 50.], [ 65.,  0.], [ 65., 10.],
+                  [ 65., 20.], [ 65., 30.], [ 65., 40.], [ 65., 50.],
+                  [ 70.,  0.], [ 70., 10.], [ 70., 20.], [ 70., 30.],
+                  [ 70., 40.], [ 70., 50.], [ 75.,  0.], [ 75., 10.],
+                  [ 75., 20.], [ 75., 30.], [ 75., 40.], [ 75., 50.],
+                  [ 80.,  0.], [ 80., 10.], [ 80., 20.], [ 80., 30.],
+                  [ 80., 40.], [ 80., 50.], [ 85.,  0.], [ 85., 10.],
+                  [ 85., 20.], [ 85., 30.], [ 85., 40.], [ 85., 50.],
+                  [ 90.,  0.], [ 90., 10.], [ 90., 20.], [ 90., 30.],
+                  [ 90., 40.], [ 90., 50.], [ 95.,  0.], [ 95., 10.],
+                  [ 95., 20.], [ 95., 30.], [ 95., 40.], [ 95., 50.],
+                  [100.,  0.], [100., 10.], [100., 20.], [100., 30.],
+                  [100., 40.], [100., 50.]]
+
+        points = num.array(points)
+
+        indices = inside_polygon(points, polygon)
+
+
 ################################################################################
 
 if __name__ == "__main__":
     suite = unittest.makeSuite(Test_Polygon,'test')
-    #suite = unittest.makeSuite(Test_Polygon,'test_inside_polygon_geospatial')
+    #suite = unittest.makeSuite(Test_Polygon,'test_inside_polygon_badtest_1')
     runner = unittest.TextTestRunner()
     runner.run(suite)