})
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
finaltime = 200
#yieldstep = 1000
#finaltime = 50000
for t in domain.evolve(yieldstep=yieldstep, finaltime=finaltime):
if myid == 0:
domain.write_time()
print_test_stats(domain, tri_full_flag)
if myid == 0:
print 'That took %.2f seconds' % (time.time() - t0)
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
finaltime = 200
#yieldstep = 1000
#finaltime = 50000
for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
if myid == 0:
domain.write_time()
print_test_stats(domain, tri_full_flag)
if myid == 0:
print 'That took %.2f seconds' %(time.time()-t0)
visualise = False
if visualise:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
vis.render_quantity_height("stage", zScale=5.0, dynamic=True)
vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
vis.start()
import time
time.sleep(2.0)
if myid == 0:
import time
t0 = time.time()
for t in domain.evolve(yieldstep=5.0, finaltime=30.0):
if myid == 0:
domain.write_time()
if visualise:
vis.update()
if visualise: vis.evolveFinished()
# Output some computation statistics
if myid == 0:
print 'That took %.2f seconds' % (time.time() - t0)
print 'Communication time %.2f seconds'\
%domain.communication_time
print 'Reduction Communication time %.2f seconds'\
if visualise:
from anuga.visualiser import RealtimeVisualiser
vis = RealtimeVisualiser(domain)
vis.render_quantity_height("stage", zScale = 5.0, dynamic=True)
vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
vis.start()
import time
time.sleep(2.0)
if myid == 0:
import time
t0 = time.time()
for t in domain.evolve(yieldstep = 5.0, finaltime = 30.0):
if myid == 0:
domain.write_time()
if visualise:
vis.update()
if visualise: vis.evolveFinished()
# Output some computation statistics
if myid == 0:
print 'That took %.2f seconds' %(time.time()-t0)
