コード例 #1
0
ファイル: eg-dist2sph.py プロジェクト: zimoun/mtf
elif len(rads) != len(phys):
    raise ValueError('different length')
else:
    N = len(rads)


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'))
コード例 #2
0
from miesphere import mie_D4grid, mie_N4grid
from krylov import gmres

from my_relative_error import relative_error

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

#################################################
コード例 #3
0
ファイル: star.py プロジェクト: zimoun/mtf
#!/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')
コード例 #4
0
ファイル: add_graph_bug.py プロジェクト: zimoun/mtf
#!/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)