Changeset 4095

Timestamp:

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

Author:

nick

Message:

add number of trigs in mesh

Location:

anuga_core/source/anuga

Files:

• geospatial_data/geospatial_data.py (modified) (1 diff)
• utilities/polygon.py (modified) (1 diff)

anuga_core/source/anuga/geospatial_data/geospatial_data.py

@@ -610 +610 @@
     def split(self, factor=0.5):
         """Returns two geospatial_data object, first is size of the 'factor'
-        smaller the original one and the second is the remainer. The two new 
+        smaller the original and the second is the remainer. The two new 
         object are disjoin set of each other. 
         

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)
+