diff --git a/src/maxwell/maxwell_cpu.cpp b/src/maxwell/maxwell_cpu.cpp index ab8a541cd7719582280579bc934634fbebc4d69b..0d29adf73788cab02922e6b75f4ace7307a38f32 100644 --- a/src/maxwell/maxwell_cpu.cpp +++ b/src/maxwell/maxwell_cpu.cpp @@ -118,11 +118,11 @@ compute_curls(solver_state& state, const field& curr, field& next) { for (const auto& ed : state.eds) { - compute_field_curls(ed, curr.Ex, curr.Ey, curr.Ez, next.Hx, next.Hy, next.Hz); + compute_field_curls(ed, curr.Ex, curr.Ey, curr.Ez, -1.0, next.Hx, next.Hy, next.Hz); #ifndef ADD_SOURCES compute_field_curls(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez); #else - compute_field_curls_sources(ed, curr.Hx, curr.Hy, curr.Hz, next.Ex, next.Ey, next.Ez, state.Jx_src, state.Jy_src, state.Jz_src); + compute_field_curls_sources(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez, state.Jx_src, state.Jy_src, state.Jz_src); #endif } @@ -138,7 +138,7 @@ compute_curls_E(solver_state& state, const field& curr, field& next) { for (const auto& ed : state.eds) { - compute_field_curls(ed, curr.Ex, curr.Ey, curr.Ez, next.Hx, next.Hy, next.Hz); + compute_field_curls(ed, curr.Ex, curr.Ey, curr.Ez, -1.0, next.Hx, next.Hy, next.Hz); } } @@ -150,7 +150,7 @@ compute_curls_H(solver_state& state, const field& curr, field& next) #ifndef ADD_SOURCES compute_field_curls(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez); #else - compute_field_curls_sources(ed, curr.Hx, curr.Hy, curr.Hz, next.Ex, next.Ey, next.Ez, state.Jx_src, state.Jy_src, state.Jz_src); + compute_field_curls_sources(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez, state.Jx_src, state.Jy_src, state.Jz_src); #endif }