jan/meanas
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Reduce number of allocations during maxwell curls

Browse Source
This commit is contained in:
Jan Petykiewicz 2019-07-09 20:19:35 -07:00
parent 9d1d8fe869
commit 557e748356
1 changed files with 16 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
fdfd_tools/fdtd.py
View File

@ -26,10 +26,14 @@ def curl_h(dxes: dx_lists_t = None) -> functional_matrix:
return f - numpy.roll(f, 1, axis=ax) return f - numpy.roll(f, 1, axis=ax)
def ch_fun(h: field_t) -> field_t: def ch_fun(h: field_t) -> field_t:
e = [dh(h[2], 1) - dh(h[1], 2), output = numpy.empty_like(h)
dh(h[0], 2) - dh(h[2], 0), output[0] = dh(h[2], 1)
dh(h[1], 0) - dh(h[0], 1)] output[1] = dh(h[0], 2)
return e output[2] = dh(h[1], 0)
output[0] -= dh(h[1], 2)
output[1] -= dh(h[2], 0)
output[2] -= dh(h[0], 1)
return output
return ch_fun return ch_fun
@ -51,10 +55,14 @@ def curl_e(dxes: dx_lists_t = None) -> functional_matrix:
return numpy.roll(f, -1, axis=ax) - f return numpy.roll(f, -1, axis=ax) - f
def ce_fun(e: field_t) -> field_t: def ce_fun(e: field_t) -> field_t:
h = [de(e[2], 1) - de(e[1], 2), output = numpy.empty_like(e)
de(e[0], 2) - de(e[2], 0), output[0] = de(e[2], 1)
de(e[1], 0) - de(e[0], 1)] output[1] = de(e[0], 2)
return h output[2] = de(e[1], 0)
output[0] -= de(e[1], 2)
output[1] -= de(e[2], 0)
output[2] -= de(e[0], 1)
return output
return ce_fun return ce_fun