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()
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()