Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py

    r5729 r5847  
    199199        self.CFL = CFL
    200200        self.set_timestepping_method(timestepping_method)
    201        
     201        self.set_beta(beta_w)
    202202        self.boundary_map = None  # Will be populated by set_boundary       
    203203       
     
    297297
    298298        return self.time + self.starttime
     299
     300    def set_beta(self,beta):
     301        """Set default beta for limiting
     302        """
     303
     304        self.beta = beta
     305        for name in self.quantities:
     306            #print 'setting beta for quantity ',name
     307            Q = self.quantities[name]
     308            Q.set_beta(beta)
     309
     310    def get_beta(self):
     311        """Get default beta for limiting
     312        """
     313
     314        return self.beta
    299315
    300316    def set_default_order(self, n):
     
    12061222        self.update_ghosts()
    12071223
     1224
     1225        # Update time
     1226        self.time += self.timestep
     1227
    12081228        # Update vertex and edge values
    12091229        self.distribute_to_vertices_and_edges()
     
    12121232        self.update_boundary()
    12131233
    1214         # Update time
    1215         self.time += self.timestep
    12161234
    12171235       
     
    12431261        self.update_ghosts()
    12441262
     1263        # Update time
     1264        self.time += self.timestep
     1265
    12451266        # Update vertex and edge values
    12461267        self.distribute_to_vertices_and_edges()
     
    12481269        # Update boundary values
    12491270        self.update_boundary()
    1250 
    1251         # Update time
    1252         self.time += self.timestep
    12531271
    12541272        #------------------------------------
     
    13091327        self.update_ghosts()
    13101328
     1329        # Update time
     1330        self.time += self.timestep
     1331
    13111332        # Update vertex and edge values
    13121333        self.distribute_to_vertices_and_edges()
     
    13151336        self.update_boundary()
    13161337
    1317         # Update time
    1318         self.time += self.timestep
    13191338
    13201339        #------------------------------------
     
    13391358        self.update_ghosts()
    13401359
     1360
     1361        # Set substep time
     1362        self.time = initial_time + self.timestep*0.5
     1363
    13411364        # Update vertex and edge values
    13421365        self.distribute_to_vertices_and_edges()
     
    13451368        self.update_boundary()
    13461369
    1347         # Set substep time
    1348         self.time = initial_time + self.timestep*0.5
    13491370
    13501371        #------------------------------------
     
    13691390        self.update_ghosts()
    13701391
     1392        # Set new time
     1393        self.time = initial_time + self.timestep       
     1394
    13711395        # Update vertex and edge values
    13721396        self.distribute_to_vertices_and_edges()
     
    13751399        self.update_boundary()
    13761400
    1377         # Set new time
    1378         self.time = initial_time + self.timestep       
    1379        
    1380 
     1401       
    13811402    def evolve_to_end(self, finaltime = 1.0):
    1382         """Iterate evolve all the way to the end
    1383         """
     1403        """Iterate evolve all the way to the end      """
    13841404
    13851405        for _ in self.evolve(yieldstep=None, finaltime=finaltime):
Note: See TracChangeset for help on using the changeset viewer.