Fix sign on H component of PML

5 years ago · 1e874cb0c0
parent ea5e298023
commit 1e874cb0c0
1 changed files with 3 additions and 3 deletions
--- a/opencl_fdtd/kernels/update_h.cl
+++ b/opencl_fdtd/kernels/update_h.cl
@ -96,9 +96,9 @@ ftype Hz_old = Hz[i];
 {%- endif %}
 // H update equations
-Hx[i] -= dt * (dEzy - dEyz + pHxi);
+Hx[i] -= dt * (dEzy - dEyz - pHxi);
-Hy[i] -= dt * (dExz - dEzx + pHyi);
+Hy[i] -= dt * (dExz - dEzx - pHyi);
-Hz[i] -= dt * (dEyx - dExy + pHzi);
+Hz[i] -= dt * (dEyx - dExy - pHzi);
 {% if do_poynting -%}
 // Average H across timesteps