Changeset 3872

Timestamp:

Oct 26, 2006, 11:56:05 AM (18 years ago)

Author:

nick

Message:

added a commented out segment of code that could be used to test for
interior polygons that lie outside the bounding polygon.

Location:

anuga_core/source/anuga/pmesh

Files:

• mesh_interface.py (modified) (2 diffs)
• test_mesh_interface.py (modified) (1 diff)

anuga_core/source/anuga/pmesh/mesh_interface.py

@@ -37 +37 @@
 
     Boundary tags is a dictionary of symbolic tags. For every tag there
-    is a list of indices referring to segments associated with that tag
+    is a list of indices referring to segments associated with that tag.
+    If a segment is omitted it will be assigned the default tag ''.
 
     maximum_triangle_area is the maximal area per triangle
@@ -157 +158 @@
                       %(str(interior_polygon), str(bounding_polygon))
                 raise PolygonError, msg
+    
+# the following segment of code could be used to Test that all the 
+# interior polygons are inside the bounding_poly... however it might need 
+# to be change a bit   
+#
+#count = 0
+#for i in range(len(interior_regions)):
+#    region = interior_regions[i]
+#    interior_polygon = region[0]
+#    if len(inside_polygon(interior_polygon, bounding_polygon,
+#                   closed = True, verbose = False)) <> len(interior_polygon):
+#        print 'WARNING: interior polygon %d is outside bounding polygon' %(i)
+#        count += 1
+
+#if count == 0:
+#    print 'interior regions OK'
+#else:
+#    print 'check out your interior polygons'
+#    print 'check %s in production directory' %figname
+#    import sys; sys.exit()
+    
 
     if interior_holes is not None:        

anuga_core/source/anuga/pmesh/test_mesh_interface.py

@@ -647 +647 @@
         interior_regions = [[poly_onslow, 50000], [poly_region, 50000], [poly_coast,100000], [poly_thevenard, 100000]]
 
-#        boundary_tags = {'walls':[0,1],'bom':[2]}
-
+        boundary_tags = {'walls':[0,1],'bom':[2]}
+        
         try:
             m = create_mesh_from_regions(poly_all,