Changeset 6001

Timestamp:

Nov 26, 2008, 3:29:51 PM (16 years ago)

Author:

ole

Message:

Make copy of the polygon in polygon_area to avoid destroying it

File:

• anuga_core/source/anuga/utilities/polygon.py (modified) (2 diffs)

anuga_core/source/anuga/utilities/polygon.py

@@ -482 +482 @@
 
 
-def polygon_area(polygon):
+def polygon_area(input_polygon):
     """ Determin area of arbitrary polygon
     Reference
@@ -489 +489 @@
     
     # Move polygon to origin (0,0) to avoid rounding errors
-    polygon = ensure_numeric(polygon)
-    min_x = min(polygon[:,0])
-    min_y = min(polygon[:,1])    
-    polygon -= [min_x, min_y]
+    # This makes a copy of the polygon to avoid destroying it
+    input_polygon = ensure_numeric(input_polygon)
+    min_x = min(input_polygon[:,0])
+    min_y = min(input_polygon[:,1])    
+    polygon = input_polygon - [min_x, min_y]
 
     # Compute area