meep.vec(model.size_x / 2, model.size_y / 2,
model.size_z / 2 - model.pml_thickness * 0))
vol.center_origin()
## Define the Perfectly Matched Layers
#perfectly_matched_layers = meep.pml(model.pml_thickness, meep.Z) ## PML on both faces at Z axis
perfectly_matched_layers = meep.pml(model.pml_thickness) ## PML on all faces
if not sim_param['frequency_domain']:
meep.master_printf("== Time domain structure setup ==\n")
## Define each polarizability by redirecting the callback to the corresponding "where_material" function
## Define the frequency-independent epsilon for all materials (needed here, before defining s, or unstable)
model.double_vec = model.eps
meep.set_EPS_Callback(model.__disown__())
s = meep.structure(vol, meep.EPS, perfectly_matched_layers,
meep.identity())
## Add all the materials
model.build_polarizabilities(s)
## Add the source dependence
#srctype = meep.band_src_time(model.srcFreq/c, model.srcWidth/c, model.simtime*c/1.1)
srctype = meep.gaussian_src_time(model.srcFreq / c, model.srcWidth /
c) ## , 0, 1000e-12 ??
else:
meep.master_printf(
"== Frequency domain structure setup (for frequency of %g Hz) ==\n" %
sim_param['frequency'])
if (model.Kx != 0 or model.Ky != 0):
print "Warning: frequency-domain solver may be broken for nonperpendicular incidence"
## Initialize volume
vol = meep.vol3d(model.size_x, model.size_y, model.size_z, 1./model.resolution)
volume_except_pml = meep.volume(
meep.vec(-model.size_x/2, -model.size_y/2, -model.size_z/2+model.pml_thickness*0),
meep.vec(model.size_x/2, model.size_y/2, model.size_z/2-model.pml_thickness*0))
vol.center_origin()
## Define the Perfectly Matched Layers
perfectly_matched_layers = meep.pml(model.pml_thickness) ## PML on both faces at Z axis
if not sim_param['frequency_domain']:
meep.master_printf("== Time domain structure setup ==\n")
## Define each polarizability by redirecting the callback to the corresponding "where_material" function
## Define the frequency-independent epsilon for all materials (needed here, before defining s, or unstable)
model.double_vec = model.get_static_permittivity; meep.set_EPS_Callback(model.__disown__())
s = meep.structure(vol, meep.EPS, perfectly_matched_layers, meep.identity())
## Add all the materials
model.build_polarizabilities(s)
## Add the source dependence
#srctype = meep.band_src_time(model.srcFreq/c, model.srcWidth/c, model.simtime*c/1.1)
srctype = meep.gaussian_src_time(model.srcFreq/c, model.srcWidth/c) ## , 0, 1000e-12 ??
else:
meep.master_printf("== Frequency domain structure setup (for frequency of %g Hz) ==\n" % sim_param['frequency'])
if (model.Kx!=0 or model.Ky!=0): print "Warning: frequency-domain solver may be broken for nonperpendicular incidence"
## Frequency-domain simulation does not support dispersive materials yet. We must define each material by
## using the nondispersive permittivity and the nondispersive conductivity
## (both calculated from polarizabilities at given frequency)