Changeset 2295

Timestamp:

Jan 27, 2006, 9:50:59 AM (19 years ago)

Author:

duncan

Message:

adding an actual test to test_mesh_interface

File:

• inundation/pmesh/test_mesh_interface.py (modified) (2 diffs)

inundation/pmesh/test_mesh_interface.py

@@ -43 +43 @@
                                      10000000,
                                      interior_regions=interior_regions,
-                                     poly_geo_reference=geo_ref_poly)
+                                     poly_geo_reference=geo_ref_poly,
+                                     mesh_geo_reference=mesh_geo)
 
         # Test the mesh instance
@@ -64 +65 @@
                         'FAILED!')
 
+        # Assuming the order of the region points is known.
+        # (This isn't true, if you consider create_mesh_from_regions
+        # a black box)
+        poly_point = m.getRegions()[0]
+        
+        #print "poly_point", poly_point
+        #print "polygon_absolute",polygon_absolute
+         
+        # poly_point values are relative to the mesh geo-ref
+        # make them absolute
+        self.failUnless(inside_polygon([poly_point.x+x,poly_point.y+y],
+                                       polygon_absolute, closed = False),
+                        'FAILED!')
+        
+        # Assuming the order of the region points is known.
+        # (This isn't true, if you consider create_mesh_from_regions
+        # a black box)
+        poly_point = m.getRegions()[1]
+        
+        #print "poly_point", poly_point
+        #print "polygon_absolute",polygon_absolute
+         
+        # poly_point values are relative to the mesh geo-ref
+        # make them absolute
+        self.failUnless(inside_polygon([poly_point.x+x,poly_point.y+y],
+                                       inner1_polygon_absolute,
+                                       closed = False),
+                        'FAILED!')
+        
+        # Assuming the order of the region points is known.
+        # (This isn't true, if you consider create_mesh_from_regions
+        # a black box)
+        poly_point = m.getRegions()[2]
+        
+        #print "poly_point", poly_point
+        #print "polygon_absolute",polygon_absolute
+         
+        # poly_point values are relative to the mesh geo-ref
+        # make them absolute
+        self.failUnless(inside_polygon([poly_point.x+x,poly_point.y+y],
+                                       inner2_polygon_absolute,
+                                       closed = False),
+                        'FAILED!')
 #-------------------------------------------------------------
 if __name__ == "__main__":