vis = RealtimeVisualiser(domain)
vis.render_quantity_height("elevation", offset=0.001, dynamic=False)
vis.render_quantity_height("stage", dynamic=True)
vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
vis.start()
import time
time.sleep(2.0)
yieldstep = 0.05
finaltime = 2.0
#Check that the boundary value gets propagated to all elements
for t in domain.evolve(yieldstep=yieldstep, finaltime=finaltime):
if myid == 0:
domain.write_time()
#print_test_stats(domain, tri_full_flag)
if visualise:
vis.update()
if visualise: vis.evolveFinished()
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' % domain.communication_reduce_time
print 'Broadcast time %.2f seconds' % domain.communication_broadcast_time
if visualise: vis.join()
finalize()
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'\
%domain.communication_reduce_time
if visualise:
vis.join()
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'\
%domain.communication_reduce_time
if visualise: vis.join()
