library_dirs = ['./code/MoM/'],
libraries = ['MoM','pthread'],
headers = ['<iostream>','<complex>','"Z_EJ_Z_HJ.h"'],
compiler = 'gcc',
extra_compile_args = ['-O3', '-pthread', '-w'])
return Z_CFIE_J
if __name__=="__main__":
path = './geo'
targetName = 'sphere'
f = 2.12e9
write_geo(path, targetName, 'lc', c/f/9.1)
write_geo(path, targetName, 'lx', 0.1)
write_geo(path, targetName, 'ly', 0.1)
write_geo(path, targetName, 'lz', 0.1)
executeGmsh(path, targetName, 0)
z_offset = 0.0
targetDimensions_scaling_factor = 1.0
languageForMeshConstruction = "Python"
meshFormat = 'GMSH'
meshFileTermination = '.msh'
target_mesh = MeshClass(path, targetName, targetDimensions_scaling_factor, z_offset, languageForMeshConstruction, meshFormat, meshFileTermination)
target_mesh.constructFromGmshFile()
N_RWG = target_mesh.N_RWG
w = 2. * pi * f
eps_r = 1.
mu_r = 1.
TDS_APPROX = 0
Z_s = 0.0
list_of_test_edges_numbers = arange(N_RWG).astype('i')
# and now the scaling tests
parameter_values = [0.1, 0.11, 0.12]
#parameter_values = [0.3, 0.4, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.2]
res = resClass(parameter_values)
for index in range(len(res.parameter_values)):
fileName = os.path.join(params_simu.pathToTarget,
params_simu.targetName)
write_geo(params_simu.pathToTarget, params_simu.targetName, 'lc',
c / params_simu.f / 10.)
write_geo(params_simu.pathToTarget, params_simu.targetName, 'lx',
res.parameter_values[index])
write_geo(params_simu.pathToTarget, params_simu.targetName, 'ly',
res.parameter_values[index])
write_geo(params_simu.pathToTarget, params_simu.targetName, 'lz',
res.parameter_values[index])
executeGmsh(params_simu.pathToTarget, params_simu.targetName, 0)
targetDimensions_scaling_factor = 1.0
target_mesh = MeshClass(params_simu.pathToTarget,
params_simu.targetName,
targetDimensions_scaling_factor,
params_simu.z_offset,
params_simu.languageForMeshConstruction)
N_RWG = target_mesh.E
sys.stdout.write("average edge length = %.5s" %
str(target_mesh.average_edge_length /
(c / params_simu.f)) + " lambdas \n")
## cubes computation
target_mesh.cubes_data_computation(a)
target_mesh.write_cubes()
# target_MLFMA construction