Changeset 5847 for anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py

Timestamp:

Oct 18, 2008, 9:13:18 PM (15 years ago)

Author:

steve

Message:

Changed parallel_api so that global mesh only needs to
be constructed on processor 0

File:

• anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py (modified) (12 diffs)

anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py

@@ -199 +199 @@
         self.CFL = CFL
         self.set_timestepping_method(timestepping_method)
-        
+        self.set_beta(beta_w)
         self.boundary_map = None  # Will be populated by set_boundary        
         
@@ -297 +297 @@
 
         return self.time + self.starttime
+
+    def set_beta(self,beta):
+        """Set default beta for limiting
+        """
+
+        self.beta = beta
+        for name in self.quantities:
+            #print 'setting beta for quantity ',name
+            Q = self.quantities[name]
+            Q.set_beta(beta)
+
+    def get_beta(self):
+        """Get default beta for limiting
+        """
+
+        return self.beta
 
     def set_default_order(self, n):
@@ -1206 +1222 @@
         self.update_ghosts()
 
+
+        # Update time
+        self.time += self.timestep
+
         # Update vertex and edge values
         self.distribute_to_vertices_and_edges()
@@ -1212 +1232 @@
         self.update_boundary()
 
-        # Update time
-        self.time += self.timestep
 
         
@@ -1243 +1261 @@
         self.update_ghosts()
 
+        # Update time
+        self.time += self.timestep
+
         # Update vertex and edge values
         self.distribute_to_vertices_and_edges()
@@ -1248 +1269 @@
         # Update boundary values
         self.update_boundary()
-
-        # Update time
-        self.time += self.timestep
 
         #------------------------------------
@@ -1309 +1327 @@
         self.update_ghosts()
 
+        # Update time
+        self.time += self.timestep
+
         # Update vertex and edge values
         self.distribute_to_vertices_and_edges()
@@ -1315 +1336 @@
         self.update_boundary()
 
-        # Update time
-        self.time += self.timestep
 
         #------------------------------------
@@ -1339 +1358 @@
         self.update_ghosts()
 
+
+        # Set substep time
+        self.time = initial_time + self.timestep*0.5
+
         # Update vertex and edge values
         self.distribute_to_vertices_and_edges()
@@ -1345 +1368 @@
         self.update_boundary()
 
-        # Set substep time
-        self.time = initial_time + self.timestep*0.5
 
         #------------------------------------
@@ -1369 +1390 @@
         self.update_ghosts()
 
+        # Set new time
+        self.time = initial_time + self.timestep       
+
         # Update vertex and edge values
         self.distribute_to_vertices_and_edges()
@@ -1375 +1399 @@
         self.update_boundary()
 
-        # Set new time
-        self.time = initial_time + self.timestep       
-        
-
+        
     def evolve_to_end(self, finaltime = 1.0):
-        """Iterate evolve all the way to the end
-        """
+        """Iterate evolve all the way to the end      """
 
         for _ in self.evolve(yieldstep=None, finaltime=finaltime):