Skip to content
Snippets Groups Projects
Commit d3a4a99d authored by Marco D'Antonio's avatar Marco D'Antonio
Browse files

Hotfix new curls computation

parent edee8b89
Branches
No related tags found
No related merge requests found
...@@ -118,11 +118,11 @@ compute_curls(solver_state& state, const field& curr, field& next) ...@@ -118,11 +118,11 @@ compute_curls(solver_state& state, const field& curr, field& next)
{ {
for (const auto& ed : state.eds) 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 #ifndef ADD_SOURCES
compute_field_curls(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez); compute_field_curls(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez);
#else #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 #endif
} }
...@@ -138,7 +138,7 @@ compute_curls_E(solver_state& state, const field& curr, field& next) ...@@ -138,7 +138,7 @@ compute_curls_E(solver_state& state, const field& curr, field& next)
{ {
for (const auto& ed : state.eds) 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) ...@@ -150,7 +150,7 @@ compute_curls_H(solver_state& state, const field& curr, field& next)
#ifndef ADD_SOURCES #ifndef ADD_SOURCES
compute_field_curls(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez); compute_field_curls(ed, curr.Hx, curr.Hy, curr.Hz, 1.0, next.Ex, next.Ey, next.Ez);
#else #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 #endif
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment