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']