#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: print 'No visualiser' from norms import linf_norm yieldstep = 1
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() # Start the evolve computions import hotshot profiler = hotshot.Profile("hotshot." + str(numprocs) + "." + str(myid) + ".prof") s = '''for t in domain.evolve(yieldstep = 0.1, finaltime = 2.0): if myid == 0: domain.write_time() '''
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() # Start the evolve computions import hotshot profiler = hotshot.Profile("hotshot." + str(numprocs) + "." + str(myid) + ".prof") s = '''for t in domain.evolve(yieldstep = 0.1, finaltime = 2.0): if myid == 0: domain.write_time() '''
# 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: print 'No visualiser' from norms import linf_norm