Ignore:
Timestamp:
Nov 7, 2006, 1:13:03 PM (18 years ago)
Author:
ole
Message:

Parallel domains now store only full triangles in sww files.
Still need to remove ghost nodes.

Location:
anuga_core/source/anuga_parallel
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga_parallel/parallel_api.py

    r3926 r3928  
    171171
    172172    for p in range(numprocs):
     173        N = len(submesh['ghost_nodes'][p])               
    173174        M = len(submesh['ghost_triangles'][p])
    174         print 'There are %d ghost triangles on proc %d' %(M, p)
    175         N = len(submesh['ghost_nodes'][p])
    176         print 'There are %d ghost nodes on proc %d' %(N, p)
     175        print 'There are %d ghost nodes and %d ghost triangles on proc %d'\
     176              %(N, M, p)
    177177
    178178
  • anuga_core/source/anuga_parallel/parallel_shallow_water.py

    r3926 r3928  
    3434                 full_send_dict=None,
    3535                 ghost_recv_dict=None,
    36                  number_of_full_nodes=0,
    37                  number_of_full_triangles=0):
     36                 number_of_full_nodes=None,
     37                 number_of_full_triangles=None):
    3838
    3939        Domain.__init__(self,
     
    4848                        number_of_full_triangles=number_of_full_triangles)
    4949
    50         N = self.number_of_elements
     50        N = len(self) # number_of_triangles
    5151
    5252#        self.processor = pypar.rank()
  • anuga_core/source/anuga_parallel/test_parallel_sw_runup.py

    r3926 r3928  
    9696        local_interpolation_points.append(i)
    9797
    98 # Hack
     98# Hack before we excluded ghosts.
    9999if numprocs == 2:
    100100    if myid == 0:
    101101        del local_interpolation_points[0]               
    102102        #local_interpolation_points = [1,2,3]
    103 
    104103if numprocs == 3:
    105104    if myid == 1:
    106105        del local_interpolation_points[0]
    107 
    108 
    109106if numprocs == 4:
    110107    if myid == 0:
     
    122119
    123120for t in domain.evolve(yieldstep = 0.1, finaltime = 5.0):
    124     #domain.write_time()
    125 
     121    domain.write_time()
    126122   
    127123    # Record time series at known points
Note: See TracChangeset for help on using the changeset viewer.