Changeset 2866 for inundation/pyvolution/mesh.py

Timestamp:

May 15, 2006, 2:26:21 PM (19 years ago)

Author:

ole

Message:

Implemented Domain(meshfile) variant widely, updated documentation and added more statistics

File:

• inundation/pyvolution/mesh.py (modified) (5 diffs)

inundation/pyvolution/mesh.py

@@ -60 +60 @@
 
 
-    def __init__(self, coordinates, triangles, boundary = None,
-                 tagged_elements = None, geo_reference = None,
-                 use_inscribed_circle = False):
+    def __init__(self, coordinates, triangles,
+                 boundary=None,
+                 tagged_elements=None,
+                 geo_reference=None,
+                 use_inscribed_circle=False,
+                 verbose=False):
         """
         Build triangles from x,y coordinates (sequence of 2-tuples or
@@ -73 +76 @@
         from Numeric import array, zeros, Int, Float, maximum, sqrt, sum
 
-        
-        General_mesh.__init__(self, coordinates, triangles, geo_reference)
+        General_mesh.__init__(self, coordinates, triangles,
+                              geo_reference, verbose=verbose)
+
+        if verbose: print 'Initialising mesh'         
 
         N = self.number_of_elements
@@ -94 +99 @@
 
         #Initialise each triangle
+        if verbose: print 'Mesh: Computing centroids and radii'        
         for i in range(N):
-            #if i % (N/10) == 0: print '(%d/%d)' %(i, N)
+            if verbose and i % ((N+10)/10) == 0: print '(%d/%d)' %(i, N)
 
             x0 = V[i, 0]; y0 = V[i, 1]
@@ -144 +150 @@
 
         #Build neighbour structure
+        if verbose: print 'Mesh: Building neigbour structure'                
         self.build_neighbour_structure()
 
         #Build surrogate neighbour structure
+        if verbose: print 'Mesh: Building surrogate neigbour structure'
         self.build_surrogate_neighbour_structure()
 
         #Build boundary dictionary mapping (id, edge) to symbolic tags
+        if verbose: print 'Mesh: Building boundary dictionary'
         self.build_boundary_dictionary(boundary)
 
         #Build tagged element  dictionary mapping (tag) to array of elements
+        if verbose: print 'Mesh: Building tagged elements dictionary'        
         self.build_tagged_elements_dictionary(tagged_elements)
 
@@ -159 +169 @@
 
         #FIXME check integrity?
+        if verbose: print 'Mesh: Done'                
 
     def __repr__(self):