Changeset 4460

Timestamp:

May 18, 2007, 4:07:12 PM (17 years ago)

Author:

ole

Message:

Faster code for get_disconnected_triangles

File:

• anuga_core/source/anuga/abstract_2d_finite_volumes/general_mesh.py (modified) (2 diffs)

anuga_core/source/anuga/abstract_2d_finite_volumes/general_mesh.py

@@ -1 +1 @@
 
 from Numeric import concatenate, reshape, take, allclose
-from Numeric import array, zeros, Int, Float, sqrt, sum
+from Numeric import array, zeros, Int, Float, sqrt, sum, arange
 
 from anuga.coordinate_transforms.geo_reference import Geo_reference
@@ -359 +359 @@
         M = len(self) # Number of triangles
         K = 3*M       # Total number of unique vertices
-        T = reshape(array(range(K)).astype(Int), (M,3))
-
+        
+        #T = reshape(array(range(K)).astype(Int), (M,3))
+        T = reshape(arange(K).astype(Int), (M,3))  # Faster
+        
         return T