jan/opencl_fdtd
2
2
Fork 1
Code Issues Pull Requests Releases Wiki Activity

call update_s during loop

Browse Source
This commit is contained in:
Jan Petykiewicz 2020-10-16 19:32:35 -07:00
parent 4e0bd8b3c6
commit 8cbb0e9864
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opencl_fdtd/simulation.py
View File

@ -26,7 +26,7 @@ class Simulation(object):
""" """
Constructs and holds the basic FDTD operations and related fields Constructs and holds the basic FDTD operations and related fields
After constructing this object, call the (update_E, update_H) members After constructing this object, call the (update_E, update_H, update_S) members
to perform FDTD updates on the stored (E, H, S) fields: to perform FDTD updates on the stored (E, H, S) fields:
pmls = [{'axis': a, 'polarity': p} for a in 'xyz' for p in 'np'] pmls = [{'axis': a, 'polarity': p} for a in 'xyz' for p in 'np']
@ -43,6 +43,8 @@ class Simulation(object):
# Perturb the field (i.e., add a soft current source) # Perturb the field (i.e., add a soft current source)
sim.E[ind] += numpy.sin(omega * t * sim.dt) sim.E[ind] += numpy.sin(omega * t * sim.dt)
event = sim.update_H([]) event = sim.update_H([])
if sim.update_S:
event = sim.update_S([event])
event.wait() event.wait()
with lzma.open('saved_simulation', 'wb') as f: with lzma.open('saved_simulation', 'wb') as f: