Example #1
0
name = "oned-test"
outdir = "result/%s" % name
prepare(outdir)

string = String(31)
string.source = source.HardSource(source.sin_oft, (f,), string.shape[0]/2)

hdf5 = h5py.File("%s/%s.hdf5" % (outdir, name),"w")
hdf5.attrs["name"]  = name
hdf5.attrs["dt"]    = dt
hdf5.attrs["dx"]    = dx
hdf5.require_group("timeline")

for t in range(0,300):
    string.update_dfield()
    string.update_efield()
    string.update_abc()
    # string.update_source(t*dt)
    string.efield[15] = source.sin_oft(t*dt, f)
    string.update_bfield()
    string.update_hfield()

    # plot(string.efield, "/tmp/%.3d.png", t)
    
    hdf5.require_group("timeline/"+str(t))
    hdf5["timeline"][str(t)]["ex"] = string.efield
    hdf5["timeline"][str(t)]["hy"] = string.hfield
    print(t)

hdf5.close()
Example #2
0
from yafdtd.grid import Cube
from yafdtd.utils import *
from scipy.constants import c, epsilon_0, mu_0


name = "threed-test"
outdir = "result/%s" % name
prepare(outdir)

hdf5 = h5py.File("%s/%s.hdf5" % (outdir, name), "w")
hdf5.attrs["name"] = name
hdf5.require_group("timeline")

length = 31

cube = Cube((length,length,length))

for t in range(100):
    cube.update_dfield()
    cube.update_efield()
    cube.ezfield[15,15,:] = sin_oft(0.005*t)
    cube.update_bfield()
    cube.update_hfield()
    
    hdf5.require_group("timeline/"+str(t))
    # hdf5["timeline"][str(t)]["ez"] = cube.ezfield[:,:,15]
    hdf5["timeline"][str(t)]["ez"] = cube.ezfield
    print t

hdf5.close()