dgens = [] for i in range(N): if i == 0: truefalse = True else: truefalse = False conf = sanitize_config(init_offset=truefalse, tag=i+1, kRef=kRef, rad=rads[i], phys=phys[i], center=centers[i], ) dgen = generate_concentric_dict(conf) cmds = write_params_geo(conf) call(cmds) dgens.append(dgen) doms = merge_msh_bubbles(dgens) myd = Domains(doms) myd.write2dot('graph.dot') call(['dot', '-Teps', 'graph.dot'], stdout=open('graph.eps', 'wb')) dd = myd print(N) ##################################
bem.global_parameters.hmat.coarsening = False ################################################# kRef_rc = 0.1 * np.pi eps_rc = 2 N = 1 geoconf = sanitize_config(Nints=N, names=[ str(i) for i in range(N+1) ], phys=(eps_rc, 1, 1), kRef=kRef_rc, meshname='sphere', init_offset=True) dd = generate_disjoint_dict(geoconf) cmds = write_params_geo(geoconf) nlambdas = [5, 10, 20, 40, 80, 160] ################################################# kRef = geoconf['kRef'] eps, _, _ = geoconf['phys'][1] Ndom = len(geoconf['names']) ################################################# iincident = 1 def dir_data(x, normal, dom_ind, result): result[0] = -np.exp( 1j * kRef * x[iincident])
#!/usr/bin/env python # coding: utf8 from subprocess import call import domains as d dgens = [] N, tag = 3, 1 conf = d.sanitize_config(Nints=N, init_offset=True, tag=tag) dgen = d.generate_concentric_dict(conf) cmds = d.write_params_geo(conf, file_geo='geo/ellipse-concentric.script.geo') call(cmds) dgens.append(dgen) tag += N N = 3 conf = d.sanitize_config(Nints=N, init_offset=False, tag=tag) conf['center'] = (3, 0, 0) dgen = d.generate_concentric_dict(conf) cmds = d.write_params_geo(conf, file_geo='geo/sphere-concentric.script.geo') call(cmds) dgens.append(dgen) tag += N N = 3 conf = d.sanitize_config(Nints=N, init_offset=False, tag=tag) conf['center'] = (-3, 0, 0) dgen = d.generate_concentric_dict(conf) cmds = d.write_params_geo(conf, file_geo='geo/sphere-concentric.script.geo')
#!/usr/bin/env python # coding: utf8 from subprocess import call import domains as d dgens = [] N, tag = 3, 1 conf = d.sanitize_config(Nints=N, init_offset=True, tag=tag) conf['center'] = (-2, 0, 0) dgen = d.generate_concentric_dict(conf) cmds = d.write_params_geo(conf, file_geo='geo/ellipse-concentric.script.geo') call(cmds) dgens.append(dgen) tag += N N = 3 conf = d.sanitize_config(Nints=N, init_offset=False, tag=tag) conf['center'] = (2, 0, 0) dgen = d.generate_concentric_dict(conf) cmds = d.write_params_geo(conf, file_geo='geo/sphere-concentric.script.geo') call(cmds) dgens.append(dgen) dicts = d.merge_msh(dgens)