Changeset 4095 for anuga_core/source/anuga/utilities/polygon.py

Timestamp:

Dec 19, 2006, 3:30:43 PM (17 years ago)

Author:

nick

Message:

add number of trigs in mesh

File:

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

anuga_core/source/anuga/utilities/polygon.py

@@ -732 +732 @@
     return point
 
+def number_mesh_triangles(interior_regions, bounding_poly, remainder_res):
+    """Calcalutes the approximate number of triangles inside the bounding polygon
+    and the other interior regions
+    FIXME: Add tests for this function
+    """
+    import polygon_area
+    
+    # TO DO check if any of the regions fall inside one another
+    no_triangles = 0.0
+    area = polygon_area(bounding_poly)
+    for i,j in interior_regions:
+        this_area = polygon_area(i)
+        no_triangles += this_area/j
+        area -= this_area
+        #convert to square Kms
+        print 'area of ',j, this_area/1000000., area/1000000.
+    no_triangles += area/remainder_res
+    return int(no_triangles/0.7)
+