diff options
Diffstat (limited to 'Source/FieldSolver/FiniteDifferenceSolver/HybridPICSolveE.cpp')
| -rw-r--r-- | Source/FieldSolver/FiniteDifferenceSolver/HybridPICSolveE.cpp | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/HybridPICSolveE.cpp b/Source/FieldSolver/FiniteDifferenceSolver/HybridPICSolveE.cpp index 8d74adb45..f4f236e61 100644 --- a/Source/FieldSolver/FiniteDifferenceSolver/HybridPICSolveE.cpp +++ b/Source/FieldSolver/FiniteDifferenceSolver/HybridPICSolveE.cpp @@ -196,7 +196,7 @@ void FiniteDifferenceSolver::HybridPICSolveE ( std::unique_ptr<amrex::MultiFab> const& Pefield, std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& edge_lengths, int lev, HybridPICModel const* hybrid_model, - DtType a_dt_type ) + const bool include_resistivity_term ) { // Select algorithm (The choice of algorithm is a runtime option, // but we compile code for each algorithm, using templates) @@ -205,14 +205,14 @@ void FiniteDifferenceSolver::HybridPICSolveE ( HybridPICSolveECylindrical <CylindricalYeeAlgorithm> ( Efield, Jfield, Jifield, Bfield, rhofield, Pefield, - edge_lengths, lev, hybrid_model, a_dt_type + edge_lengths, lev, hybrid_model, include_resistivity_term ); #else HybridPICSolveECartesian <CartesianYeeAlgorithm> ( Efield, Jfield, Jifield, Bfield, rhofield, Pefield, - edge_lengths, lev, hybrid_model, a_dt_type + edge_lengths, lev, hybrid_model, include_resistivity_term ); #endif @@ -233,14 +233,14 @@ void FiniteDifferenceSolver::HybridPICSolveECylindrical ( std::unique_ptr<amrex::MultiFab> const& Pefield, std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& edge_lengths, int lev, HybridPICModel const* hybrid_model, - DtType a_dt_type ) + const bool include_resistivity_term ) { #ifndef AMREX_USE_EB amrex::ignore_unused(edge_lengths); #endif amrex::ignore_unused( Efield, Jfield, Jifield, Bfield, rhofield, Pefield, edge_lengths, - lev, hybrid_model, a_dt_type + lev, hybrid_model, include_resistivity_term ); amrex::Abort(Utils::TextMsg::Err( "currently hybrid E-solve does not work for RZ")); @@ -258,7 +258,7 @@ void FiniteDifferenceSolver::HybridPICSolveECartesian ( std::unique_ptr<amrex::MultiFab> const& Pefield, std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& edge_lengths, int lev, HybridPICModel const* hybrid_model, - DtType a_dt_type ) + const bool include_resistivity_term ) { #ifndef AMREX_USE_EB amrex::ignore_unused(edge_lengths); @@ -435,7 +435,7 @@ void FiniteDifferenceSolver::HybridPICSolveECartesian ( Ex(i, j, k) = (enE_x - grad_Pe) / rho_val; // Add resistivity only if E field value is used to update B - if (a_dt_type != DtType::Full) Ex(i, j, k) += eta(rho_val) * Jx(i, j, k); + if (include_resistivity_term) Ex(i, j, k) += eta(rho_val) * Jx(i, j, k); }, // Ey calculation @@ -465,7 +465,7 @@ void FiniteDifferenceSolver::HybridPICSolveECartesian ( Ey(i, j, k) = (enE_y - grad_Pe) / rho_val; // Add resistivity only if E field value is used to update B - if (a_dt_type != DtType::Full) Ey(i, j, k) += eta(rho_val) * Jy(i, j, k); + if (include_resistivity_term) Ey(i, j, k) += eta(rho_val) * Jy(i, j, k); }, // Ez calculation @@ -490,7 +490,7 @@ void FiniteDifferenceSolver::HybridPICSolveECartesian ( Ez(i, j, k) = (enE_z - grad_Pe) / rho_val; // Add resistivity only if E field value is used to update B - if (a_dt_type != DtType::Full) Ez(i, j, k) += eta(rho_val) * Jz(i, j, k); + if (include_resistivity_term) Ez(i, j, k) += eta(rho_val) * Jz(i, j, k); } ); |