Python Parallel_Domain.set_quantityの例

コード例 #1

ファイル: run_parallel_merimbula_test.py プロジェクト: ZijunZhang-Zoe/anuga_core_backup

#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

コード例 #2

ファイル: run_parallel_advection_prof.py プロジェクト: MattAndersonPE/anuga_core

                         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()
'''
    

コード例 #3

                         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()
'''

コード例 #4

ファイル: run_parallel_merimbula_test.py プロジェクト: MattAndersonPE/anuga_core

# 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