Changeset 3021 for inundation/pyvolution/domain.py

Timestamp:

May 30, 2006, 4:17:18 PM (19 years ago)

Author:

steve

Message:

updated timestep using only full cells

File:

• inundation/pyvolution/domain.py (modified) (6 diffs)

inundation/pyvolution/domain.py

@@ -79 +79 @@
                       verbose=verbose)
 
-        if verbose: print 'Initialising Domain'        
-        from Numeric import zeros, Float, Int
+        if verbose: print 'Initialising Domain'
+        from Numeric import zeros, Float, Int, ones
         from quantity import Quantity, Conserved_quantity
 
@@ -117 +117 @@
         # List of other quantity names
         if ghost_recv_dict  is None:
-            self.ghost_recv_dict  = []
+            self.ghost_recv_dict  = {}
         else:
             self.ghost_recv_dict  = ghost_recv_dict
@@ -125 +125 @@
 
         # Setup Communication Buffers
-           
+
         if verbose: print 'Domain: Set up communication buffers (parallel)'
         self.nsys = len(self.conserved_quantities)
@@ -136 +136 @@
             buffer_shape = self.ghost_recv_dict[key][0].shape[0]
             self.ghost_recv_dict[key].append(zeros( (buffer_shape,self.nsys) ,Float))
+
+
+        # Setup cell full flag
+        # =1 for full
+        # =0 for ghost
+        N=self.number_of_elements
+        self.tri_full_flag = ones(N, Int)
+        for i in self.ghost_recv_dict.keys():
+            for id in self.ghost_recv_dict[i][0]:
+                self.tri_full_flag[id] = 0
+
 
         #Defaults
@@ -186 +197 @@
             self.set_quantity_vertices_dict(vertex_quantity_dict)
 
-           
+
         if verbose: print 'Domain: Done'
 
-            
+
 
 
@@ -730 +741 @@
                 self.number_of_steps = 0
                 self.number_of_first_order_steps = 0
-                
+
 
     def evolve_to_end(self, finaltime = 1.0):