Changeset 4131 for anuga_core/source/anuga/utilities

Timestamp:

Jan 5, 2007, 2:30:04 PM (18 years ago)

Author:

nick

Message:

add "number_mesh_triangles" to polygon, given interior regions and bounding poly returns approx number of triangles in mesh

also
editted screen_catcher in util.py

File:

• anuga_core/source/anuga/utilities/polygon.py (modified) (1 diff)

anuga_core/source/anuga/utilities/polygon.py

@@ -737 +737 @@
     FIXME: Add tests for this function
     """
-    import polygon_area
+    from anuga.utilities.polygon import polygon_area
     
     # TO DO check if any of the regions fall inside one another
     no_triangles = 0.0
     area = polygon_area(bounding_poly)
+#    print 'area of bounding_poly with res ',remainder_res,' is ', area/1000000.
     for i,j in interior_regions:
         this_area = polygon_area(i)
-        no_triangles += this_area/j
+        this_triangles = this_area/j
+        no_triangles += this_triangles
         area -= this_area
         #convert to square Kms
-        print 'area of ',j, this_area/1000000., area/1000000.
-    no_triangles += area/remainder_res
+        print 'area of',j, 'trigs', this_triangles, 'this area', this_area/1000000.
+    bound_triangles = area/remainder_res
+    no_triangles += bound_triangles
+    print 'area of',remainder_res,'bound triangles:', bound_triangles, 'bound area:', area/1000000.
     return int(no_triangles/0.7)
-