velocity=[0.1, 0.0]) # Make a notes of which triangles are full and which are ghost tri_full_flag = build_full_flag(domain, ghost_recv_dict) #domain.initialise_visualiser(rect=rect) #Boundaries T = Transmissive_boundary(domain) #R = Reflective_boundary(domain) domain.set_boundary({ 'outflow': T, 'inflow': T, 'inner': T, 'exterior': T, 'open': T, 'ghost': None }) domain.set_quantity('stage', quantities['stage']) #--------- # Evolution t0 = time.time() try: domain.initialise_visualiser(rect=rect) #domain.visualiser.coloring['stage'] = True #domain.visualiser.scale_z['stage'] = 0.2 except:
########################################### # 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 t0 = time.time()
# 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 t0 = time.time()
domain = Parallel_Domain(points, vertices, boundary, full_send_dict = full_send_dict, ghost_recv_dict = ghost_recv_dict, velocity = [0.1,0.0]) # Make a notes of which triangles are full and which are ghost tri_full_flag = build_full_flag(domain, ghost_recv_dict) #domain.initialise_visualiser(rect=rect) #Boundaries T = Transmissive_boundary(domain) #R = Reflective_boundary(domain) domain.set_boundary( {'outflow': T, 'inflow': T, 'inner':T, 'exterior': T, 'open':T, 'ghost':None} ) domain.set_quantity('stage', quantities['stage']) #--------- # Evolution t0 = time.time() try: domain.initialise_visualiser(rect=rect) #domain.visualiser.coloring['stage'] = True #domain.visualiser.scale_z['stage'] = 0.2 except: