diff options
| author |  Roelof Groenewald <40245517+roelof-groenewald@users.noreply.github.com>
| 2023-06-27 00:27:05 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2023-06-27 00:27:05 -0700 |
| commit | 9665008751b9cd7b97fc0e8005c95bc859a50bec (patch) | |
| tree | 4c82ef57f21c03f655ab4051b71abdea3f4fa79f /Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp | |
| parent | 902c79367dd6b503dcff5706f9ee6153cb7b1319 (diff) | |
| download | WarpX-9665008751b9cd7b97fc0e8005c95bc859a50bec.tar.gz WarpX-9665008751b9cd7b97fc0e8005c95bc859a50bec.tar.zst WarpX-9665008751b9cd7b97fc0e8005c95bc859a50bec.zip | |
refactored `HybridPICModel::HybridPICSolveE` to more readable (#4044)
Diffstat (limited to 'Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp')
| -rw-r--r-- | Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp b/Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp index 7160fa025..b50262064 100644 --- a/Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp +++ b/Source/FieldSolver/WarpXPushFieldsHybridPIC.cpp @@ -88,8 +88,8 @@ void WarpX::HybridPICEvolveFields () { m_hybrid_pic_model->CalculateCurrentAmpere(Bfield_fp, m_edge_lengths); m_hybrid_pic_model->HybridPICSolveE( - Efield_fp, current_fp, Bfield_fp, rho_fp, m_edge_lengths, - DtType::FirstHalf + Efield_fp, current_fp_temp, Bfield_fp, rho_fp_temp, m_edge_lengths, + true ); FillBoundaryE(guard_cells.ng_FieldSolver, WarpX::sync_nodal_points); EvolveB(0.5 / sub_steps * dt[0], DtType::FirstHalf); @@ -116,8 +116,8 @@ void WarpX::HybridPICEvolveFields () { m_hybrid_pic_model->CalculateCurrentAmpere(Bfield_fp, m_edge_lengths); m_hybrid_pic_model->HybridPICSolveE( - Efield_fp, current_fp, Bfield_fp, rho_fp, m_edge_lengths, - DtType::SecondHalf + Efield_fp, current_fp, Bfield_fp, rho_fp_temp, m_edge_lengths, + true ); FillBoundaryE(guard_cells.ng_FieldSolver, WarpX::sync_nodal_points); EvolveB(0.5 / sub_steps * dt[0], DtType::SecondHalf); @@ -148,8 +148,8 @@ void WarpX::HybridPICEvolveFields () // Update the E field to t=n+1 using the extrapolated J_i^n+1 value m_hybrid_pic_model->CalculateCurrentAmpere(Bfield_fp, m_edge_lengths); m_hybrid_pic_model->HybridPICSolveE( - Efield_fp, current_fp, Bfield_fp, rho_fp, m_edge_lengths, - DtType::Full + Efield_fp, current_fp_temp, Bfield_fp, rho_fp, m_edge_lengths, + false ); FillBoundaryE(guard_cells.ng_FieldSolver, WarpX::sync_nodal_points); |