print "nb=%d, nb_edge=%d surrnb=%d norm=%d" % \
(cnt_nb, cnt_nbedge, cnt_surrnb, cnt_norm)
if cnt_el or cnt_vc or cnt_ec or cnt_cc:
print "elen=%d, vc=%d ec=%d cc=%d" % \
(cnt_el, cnt_vc, cnt_ec, cnt_cc)
if cnt_se or cnt_xe or cnt_ye or cnt_be:
print "se=%d, xe=%d ye=%d be=%d" % \
(cnt_se, cnt_xe, cnt_ye, cnt_be)
if cnt_sv or cnt_xv or cnt_yv or cnt_bv:
print "sv=%d, xv=%d yv=%d bv=%d" % \
(cnt_sv, cnt_xv, cnt_yv, cnt_bv)
if __name__ == "__main__":
from anuga_cuda import generate_merimbula_domain
from anuga_cuda.compute_fluxes.compute_fluxes import spe_bubble_sort
domain1 = generate_merimbula_domain()
domain2 = generate_merimbula_domain()
#sort_domain(domain2)
#sort_domain_check(domain2)
#sort_domain(domain1)
print " **** rearrange_domain **** "
domain2=rearrange_domain(domain2, False)
print " **** check domain ****"
rearrange_domain_check(domain1, domain2)
#!/usr/bin/env python
import numpy
from pycuda import driver as drv
from anuga_cuda import generate_merimbula_domain, generate_cairns_domain, get_kernel_function_info
using_tsunami_domain = False
if using_tsunami_domain:
domain1 = generate_cairns_domain(False)
domain2 = generate_cairns_domain(True)
else:
domain1 = generate_merimbula_domain()
domain2 = generate_merimbula_domain(gpu=True)
domain2.equip_kernel_functions()
domain1.evolve(yieldstep = 50, finaltime = 500)
domain1.evolve(yieldstep = 50, finaltime = 500)
domain1.evolve(yieldstep = 50, finaltime = 500)
if using_tsunami_domain:
domain1.distribute_to_vertices_and_edges()
else:
from swb2_domain_ext import \
extrapolate_second_order_edge_sw as extrapol2_ext
# Shortcuts
Stage = domain1.quantities['stage']
Xmom = domain1.quantities['xmomentum']
Ymom = domain1.quantities['ymomentum']
