Skip to content

Clean-up: moved current correction print statement inside PSATD solver condition #6721

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
oshapoval wants to merge 2 commits into
base: development
Choose a base branch
from
+56 −54
Open

Clean-up: moved current correction print statement inside PSATD solver condition #6721

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
040036d
Clean-up: moved current correction print statement inside PSATD solve…
oshapoval Mar 27, 2026
3007f11
Merge remote-tracking branch 'origin/development' into update_runtime…
oshapoval Apr 24, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
110 changes: 56 additions & 54 deletions Source/Initialization/WarpXInitData.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -647,61 +647,63 @@ WarpX::PrintMainPICparameters ()
else if (WarpX::electromagnetic_solver_id == ElectromagneticSolverAlgo::HybridPIC){
amrex::Print() << "Maxwell Solver: | Hybrid-PIC (Ohm's law) \n";
}
#ifdef WARPX_USE_FFT
// Print PSATD solver's configuration
if (WarpX::electromagnetic_solver_id == ElectromagneticSolverAlgo::PSATD){
amrex::Print() << "Maxwell Solver: | PSATD \n";
}
if ((m_v_galilean[0]!=0) or (m_v_galilean[1]!=0) or (m_v_galilean[2]!=0)) {
amrex::Print() << " | - Galilean \n" <<
" | - v_galilean = (" << m_v_galilean[0] << "," <<
m_v_galilean[1] << "," << m_v_galilean[2] << ")\n";
}
if ((m_v_comoving[0]!=0) or (m_v_comoving[1]!=0) or (m_v_comoving[2]!=0)) {
amrex::Print() << " | - comoving \n" <<
" | - v_comoving = (" << m_v_comoving[0] << "," <<
m_v_comoving[1] << "," << m_v_comoving[2] << ")\n";
}
if (WarpX::update_with_rho) {
amrex::Print() << " | - update with rho is ON \n";
}
if (current_correction) {
amrex::Print() << " | - current correction is ON \n";
}
if (WarpX::do_dive_cleaning) {
amrex::Print() << " | - div(E) cleaning is ON \n";
}
if (WarpX::do_divb_cleaning) {
amrex::Print() << " | - div(B) cleaning is ON \n";
}
if (m_JRhom == 1){
amrex::Print() << " | - PSATD-JRhom deposition is ON \n";
amrex::Print() << " | - m_JRhom_subintervals = "
<< WarpX::m_JRhom_subintervals << "\n";
if (time_dependency_J == TimeDependencyJ::Linear){
amrex::Print() << " | - time_dependency_J = linear \n";
}
else if (time_dependency_J == TimeDependencyJ::Constant){
amrex::Print() << " | - time_dependency_J = constant \n";
}
else if (time_dependency_J == TimeDependencyJ::Quadratic){
amrex::Print() << " | - time_dependency_J = quadratic \n";
}
if (time_dependency_rho == TimeDependencyRho::Linear){
amrex::Print() << " | - time_dependency_rho = linear \n";
}
else if (time_dependency_rho == TimeDependencyRho::Constant){
amrex::Print() << " | - time_dependency_rho = constant \n";
}
else if (time_dependency_rho == TimeDependencyRho::Quadratic){
amrex::Print() << " | - time_dependency_rho = quadratic \n";
}
}
if (fft_do_time_averaging){
amrex::Print()<<" | - time-averaged is ON \n";
}
#endif // WARPX_USE_FFT
#ifdef WARPX_USE_FFT
// All FFT/PSATD-related prints only when PSATD is selected at runtime
if (WarpX::electromagnetic_solver_id == ElectromagneticSolverAlgo::PSATD){
amrex::Print() << "Maxwell Solver: | PSATD \n";

if ((m_v_galilean[0]!=0) or (m_v_galilean[1]!=0) or (m_v_galilean[2]!=0)) {
amrex::Print() << " | - Galilean \n" <<
" | - v_galilean = (" << m_v_galilean[0] << "," <<
m_v_galilean[1] << "," << m_v_galilean[2] << ")\n";
}
if ((m_v_comoving[0]!=0) or (m_v_comoving[1]!=0) or (m_v_comoving[2]!=0)) {
amrex::Print() << " | - comoving \n" <<
" | - v_comoving = (" << m_v_comoving[0] << "," <<
m_v_comoving[1] << "," << m_v_comoving[2] << ")\n";
}
if (WarpX::update_with_rho) {
amrex::Print() << " | - update with rho is ON \n";
}
if (current_correction) {
amrex::Print() << " | - current correction is ON \n";
}
if (m_JRhom == 1){
amrex::Print() << " | - PSATD-JRhom deposition is ON \n";
amrex::Print() << " | - m_JRhom_subintervals = "
<< WarpX::m_JRhom_subintervals << "\n";
if (time_dependency_J == TimeDependencyJ::Linear){
amrex::Print() << " | - time_dependency_J = linear \n";
}
else if (time_dependency_J == TimeDependencyJ::Constant){
amrex::Print() << " | - time_dependency_J = constant \n";
}
else if (time_dependency_J == TimeDependencyJ::Quadratic){
amrex::Print() << " | - time_dependency_J = quadratic \n";
}
}

if (time_dependency_rho == TimeDependencyRho::Linear){
amrex::Print() << " | - time_dependency_rho = linear \n";
}
else if (time_dependency_rho == TimeDependencyRho::Constant){
amrex::Print() << " | - time_dependency_rho = constant \n";
}
else if (time_dependency_rho == TimeDependencyRho::Quadratic){
amrex::Print() << " | - time_dependency_rho = quadratic \n";
}

if (fft_do_time_averaging){
amrex::Print()<<" | - time-averaged is ON \n";
}
}
#endif // WARPX_USE_FFT
if (WarpX::do_dive_cleaning) {
amrex::Print() << " | - div(E) cleaning is ON \n";
}
if (WarpX::do_divb_cleaning) {
amrex::Print() << " | - div(B) cleaning is ON \n";
}
if (grid_type == GridType::Collocated){
amrex::Print() << " | - collocated grid \n";
}
Expand Down
Loading
⚡