########################################### # Create advection domain with direction (1,-1) # Initial condition is zero by default domain = Parallel_Domain(points, vertices, boundary, full_send_dict, ghost_recv_dict, velocity=[1.0, 0.0]) # Boundaries T = Transmissive_boundary(domain) domain.default_order = 2 domain.set_boundary( {'left': T, 'right': T, 'bottom': T, 'top': T, \ 'ghost': None} ) # Ensure that the domain definitions make sense domain.check_integrity() # Set the inititial conditions domain.set_quantity('stage', Set_Stage(0.2, 0.4, 1.0)) # Let processor 0 output some timing information if myid == 0: import time
"no vertices = ", len(vertices), "no boundaries = ", len(boundary) ########################################### # Start the computations on each subpartion ########################################### # Create advection domain with direction (1,-1) # Initial condition is zero by default domain = Parallel_Domain(points, vertices, boundary, full_send_dict, ghost_recv_dict, velocity=[1.0, 0.0]) # Boundaries T = Transmissive_boundary(domain) domain.default_order = 2 domain.set_boundary( {'left': T, 'right': T, 'bottom': T, 'top': T, \ 'ghost': None} ) # Ensure that the domain definitions make sense domain.check_integrity() # Set the inititial conditions domain.set_quantity('stage', Set_Stage(0.2,0.4,1.0)) # Let processor 0 output some timing information if myid == 0: import time