source for python script - it is not equivalent to 1.ctl
#!/usr/bin/python
import meep
import getopt
import sys
m=meep.MPI(sys.argv)
if meep.my_rank()==0:
print 'Count processors:', meep.count_processors(), 'My rank is',
meep.my_rank()
class epsilon(meep.Callback):
def __init__(self,inner_eps=1.0,outer_eps=1.0):
meep.Callback.__init__(self)
self.inner_eps=inner_eps
self.outer_eps=outer_eps
def eps(self,vec):
if
abs(vec.y()-4)<=0.5:
return
self.inner_eps
else:
return
self.outer_eps
vl=meep.vol2d(16,8,10)
eps=epsilon(12)
meep.setCallback(eps.__disown__())
pml=meep.pml(1.0)
#st=meep.structure(vl,meep.EPS,pml)
st=meep.structure(vl,meep.EPS,pml, meep.identity(),2)
fl=meep.fields(st)
sr=meep.continuous_src_time(0.15)
r_pos=meep.vec(1.0,4.0)
fl.add_point_source(meep.Ey,sr,r_pos)
fl.output_hdf5(meep.Dielectric,vl.interior())
while fl.t<=10000:
fl.step()
print m.xxx.elapsed_time()
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
