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:
