aboutsummaryrefslogtreecommitdiff
path: root/Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp')
-rw-r--r--Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp12
1 files changed, 6 insertions, 6 deletions
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp b/Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp
index ffa9fbfc8..fe9aa979a 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp
+++ b/Source/FieldSolver/FiniteDifferenceSolver/MacroscopicEvolveE.cpp
@@ -184,8 +184,8 @@ void FiniteDifferenceSolver::MacroscopicEvolveECartesian (
// Interpolated permittivity, epsilon, to Ex position on the grid
amrex::Real const epsilon_interp = ablastr::coarsen::sample::Interp(eps_arr, epsilon_stag,
Ex_stag, macro_cr, i, j, k, scomp);
- amrex::Real alpha = T_MacroAlgo::alpha( sigma_interp, epsilon_interp, dt);
- amrex::Real beta = T_MacroAlgo::beta( sigma_interp, epsilon_interp, dt);
+ const amrex::Real alpha = T_MacroAlgo::alpha( sigma_interp, epsilon_interp, dt);
+ const amrex::Real beta = T_MacroAlgo::beta( sigma_interp, epsilon_interp, dt);
Ex(i, j, k) = alpha * Ex(i, j, k)
+ beta * ( - T_Algo::DownwardDz(Hy, coefs_z, n_coefs_z, i, j, k,0)
+ T_Algo::DownwardDy(Hz, coefs_y, n_coefs_y, i, j, k,0)
@@ -208,8 +208,8 @@ void FiniteDifferenceSolver::MacroscopicEvolveECartesian (
// Interpolated permittivity, epsilon, to Ey position on the grid
amrex::Real const epsilon_interp = ablastr::coarsen::sample::Interp(eps_arr, epsilon_stag,
Ey_stag, macro_cr, i, j, k, scomp);
- amrex::Real alpha = T_MacroAlgo::alpha( sigma_interp, epsilon_interp, dt);
- amrex::Real beta = T_MacroAlgo::beta( sigma_interp, epsilon_interp, dt);
+ const amrex::Real alpha = T_MacroAlgo::alpha( sigma_interp, epsilon_interp, dt);
+ const amrex::Real beta = T_MacroAlgo::beta( sigma_interp, epsilon_interp, dt);
Ey(i, j, k) = alpha * Ey(i, j, k)
+ beta * ( - T_Algo::DownwardDx(Hz, coefs_x, n_coefs_x, i, j, k,0)
@@ -228,8 +228,8 @@ void FiniteDifferenceSolver::MacroscopicEvolveECartesian (
// Interpolated permittivity, epsilon, to Ez position on the grid
amrex::Real const epsilon_interp = ablastr::coarsen::sample::Interp(eps_arr, epsilon_stag,
Ez_stag, macro_cr, i, j, k, scomp);
- amrex::Real alpha = T_MacroAlgo::alpha( sigma_interp, epsilon_interp, dt);
- amrex::Real beta = T_MacroAlgo::beta( sigma_interp, epsilon_interp, dt);
+ const amrex::Real alpha = T_MacroAlgo::alpha( sigma_interp, epsilon_interp, dt);
+ const amrex::Real beta = T_MacroAlgo::beta( sigma_interp, epsilon_interp, dt);
Ez(i, j, k) = alpha * Ez(i, j, k)
+ beta * ( - T_Algo::DownwardDy(Hx, coefs_y, n_coefs_y, i, j, k,0)