From 5c1603ea2fe87b8fd2ea61ca0bafa220a3efccab Mon Sep 17 00:00:00 2001
From: Luca Fedeli <luca.fedeli@cea.fr>
Date: Sat, 14 Sep 2024 01:13:44 +0200
Subject: [PATCH] Replace "std::endl" with "\n" (except for error messages)
 (#5183)

* replace std::endl with \n

* revert to std::endl when flushing the buffer seems intentional

* Update Source/Diagnostics/BTDiagnostics.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Source/Diagnostics/FlushFormats/FlushFormatCatalyst.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Source/Diagnostics/ReducedDiags/ChargeOnEB.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Source/Diagnostics/ReducedDiags/ChargeOnEB.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Source/FieldSolver/ImplicitSolvers/SemiImplicitEM.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Source/FieldSolver/ImplicitSolvers/ThetaImplicitEM.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Source/Particles/PhysicalParticleContainer.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Tools/QedTablesUtils/Source/ArgParser/QedTablesArgParser.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Tools/QedTablesUtils/Source/QedTableGenerator.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Update Tools/QedTablesUtils/Source/QedTableReader.cpp

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

---------

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>
---
 Source/Diagnostics/BTDiagnostics.cpp          |  6 ++--
 .../FlushFormats/FlushFormatCatalyst.cpp      |  2 +-
 .../Diagnostics/ReducedDiags/BeamRelevant.cpp |  6 ++--
 .../Diagnostics/ReducedDiags/ChargeOnEB.cpp   |  3 +-
 .../ReducedDiags/ColliderRelevant.cpp         |  2 +-
 .../ReducedDiags/DifferentialLuminosity.cpp   |  2 +-
 .../Diagnostics/ReducedDiags/FieldEnergy.cpp  |  2 +-
 .../Diagnostics/ReducedDiags/FieldMaximum.cpp |  2 +-
 .../ReducedDiags/FieldMomentum.cpp            |  2 +-
 .../Diagnostics/ReducedDiags/FieldProbe.cpp   |  4 +--
 .../ReducedDiags/FieldReduction.cpp           |  2 +-
 .../ReducedDiags/LoadBalanceCosts.cpp         |  6 ++--
 .../ReducedDiags/LoadBalanceEfficiency.cpp    |  2 +-
 .../ReducedDiags/ParticleEnergy.cpp           |  2 +-
 .../ReducedDiags/ParticleExtrema.cpp          |  2 +-
 .../ReducedDiags/ParticleHistogram.cpp        |  2 +-
 .../ReducedDiags/ParticleMomentum.cpp         |  2 +-
 .../ReducedDiags/ParticleNumber.cpp           |  2 +-
 .../Diagnostics/ReducedDiags/ReducedDiags.cpp |  2 +-
 .../Diagnostics/ReducedDiags/RhoMaximum.cpp   |  2 +-
 Source/Diagnostics/WarpXOpenPMD.cpp           |  2 +-
 .../ImplicitSolvers/SemiImplicitEM.cpp        | 19 ++++++-------
 .../ImplicitSolvers/ThetaImplicitEM.cpp       | 21 +++++++-------
 Source/Initialization/WarpXInitData.cpp       |  2 +-
 Source/NonlinearSolvers/NewtonSolver.H        | 28 +++++++++----------
 Source/NonlinearSolvers/PicardSolver.H        | 14 +++++-----
 .../Particles/PhysicalParticleContainer.cpp   |  6 ++--
 .../Source/ArgParser/QedTablesArgParser.cpp   |  4 +--
 Tools/QedTablesUtils/Source/QedTableCommons.H |  2 +-
 .../Source/QedTableGenerator.cpp              |  2 +-
 .../QedTablesUtils/Source/QedTableReader.cpp  |  2 +-
 31 files changed, 77 insertions(+), 80 deletions(-)

diff --git a/Source/Diagnostics/BTDiagnostics.cpp b/Source/Diagnostics/BTDiagnostics.cpp
index 1cee9909226..6fdb605f8dc 100644
--- a/Source/Diagnostics/BTDiagnostics.cpp
+++ b/Source/Diagnostics/BTDiagnostics.cpp
@@ -170,18 +170,18 @@ void BTDiagnostics::DerivedInitData ()
         if (final_snapshot_fill_iteration > warpx.maxStep()) {
             warpx.updateMaxStep(final_snapshot_fill_iteration);
             amrex::Print()<<"max_step insufficient to fill all BTD snapshots. Automatically increased to: "
-                << final_snapshot_fill_iteration << std::endl;
+                << final_snapshot_fill_iteration << "\n";
 
         }
         if (final_snapshot_fill_time > warpx.stopTime()) {
             warpx.updateStopTime(final_snapshot_fill_time);
             amrex::Print()<<"stop_time insufficient to fill all BTD snapshots. Automatically increased to: "
-                << final_snapshot_fill_time << std::endl;
+                << final_snapshot_fill_time << "\n";
 
         }
         if (warpx.maxStep() == std::numeric_limits<int>::max() && warpx.stopTime() == std::numeric_limits<amrex::Real>::max()) {
             amrex::Print()<<"max_step unspecified and stop time unspecified.  Setting max step to "
-                <<final_snapshot_fill_iteration<< " to fill all BTD snapshots." << std::endl;
+                <<final_snapshot_fill_iteration<< " to fill all BTD snapshots.\n";
             warpx.updateMaxStep(final_snapshot_fill_iteration);
         }
 
diff --git a/Source/Diagnostics/FlushFormats/FlushFormatCatalyst.cpp b/Source/Diagnostics/FlushFormats/FlushFormatCatalyst.cpp
index 5b1ccfac1a5..425c13de6a4 100644
--- a/Source/Diagnostics/FlushFormats/FlushFormatCatalyst.cpp
+++ b/Source/Diagnostics/FlushFormats/FlushFormatCatalyst.cpp
@@ -204,7 +204,7 @@ FlushFormatCatalyst::~FlushFormatCatalyst() {
         amrex::Abort(message);
     } else {
         // Temporary, remove for final
-        std::cout << "Successfully finalized Catalyst" << std::endl;
+        std::cout << "Successfully finalized Catalyst\n";
     }
 
 }
diff --git a/Source/Diagnostics/ReducedDiags/BeamRelevant.cpp b/Source/Diagnostics/ReducedDiags/BeamRelevant.cpp
index 9b5a28a3516..af18ab508a9 100644
--- a/Source/Diagnostics/ReducedDiags/BeamRelevant.cpp
+++ b/Source/Diagnostics/ReducedDiags/BeamRelevant.cpp
@@ -109,7 +109,7 @@ BeamRelevant::BeamRelevant (const std::string& rd_name)
             ofs << "[" << c++ << "]alpha_y()";        ofs << m_sep;
             ofs << "[" << c++ << "]beta_x(m)";        ofs << m_sep;
             ofs << "[" << c++ << "]beta_y(m)";        ofs << m_sep;
-            ofs << "[" << c++ << "]charge(C)";        ofs << std::endl;
+            ofs << "[" << c++ << "]charge(C)";        ofs << "\n";
 #elif (defined WARPX_DIM_XZ)
             int c = 0;
             ofs << "#";
@@ -131,7 +131,7 @@ BeamRelevant::BeamRelevant (const std::string& rd_name)
             ofs << "[" << c++ << "]emittance_z(m)";   ofs << m_sep;
             ofs << "[" << c++ << "]alpha_x()";        ofs << m_sep;
             ofs << "[" << c++ << "]beta_x(m)";        ofs << m_sep;
-            ofs << "[" << c++ << "]charge(C)";        ofs << std::endl;
+            ofs << "[" << c++ << "]charge(C)";        ofs << "\n";
 #elif (defined WARPX_DIM_1D_Z)
             int c = 0;
             ofs << "#";
@@ -148,7 +148,7 @@ BeamRelevant::BeamRelevant (const std::string& rd_name)
             ofs << "[" << c++ << "]pz_rms(kg*m/s)";   ofs << m_sep;
             ofs << "[" << c++ << "]gamma_rms()";      ofs << m_sep;
             ofs << "[" << c++ << "]emittance_z(m)";   ofs << m_sep;
-            ofs << "[" << c++ << "]charge(C)";        ofs << std::endl;
+            ofs << "[" << c++ << "]charge(C)";        ofs << "\n";
 #endif
             // close file
             ofs.close();
diff --git a/Source/Diagnostics/ReducedDiags/ChargeOnEB.cpp b/Source/Diagnostics/ReducedDiags/ChargeOnEB.cpp
index 01a2b5d8b23..975fed6d74c 100644
--- a/Source/Diagnostics/ReducedDiags/ChargeOnEB.cpp
+++ b/Source/Diagnostics/ReducedDiags/ChargeOnEB.cpp
@@ -79,8 +79,7 @@ ChargeOnEB::ChargeOnEB (const std::string& rd_name)
             ofs << m_sep;
             ofs << "[" << c++ << "]time(s)";
             ofs << m_sep;
-            ofs << "[" << c++ << "]Charge (C)";
-            ofs << std::endl;
+            ofs << "[" << c++ << "]Charge (C)\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/ColliderRelevant.cpp b/Source/Diagnostics/ReducedDiags/ColliderRelevant.cpp
index 1b0c74cf7fa..fb683e25319 100644
--- a/Source/Diagnostics/ReducedDiags/ColliderRelevant.cpp
+++ b/Source/Diagnostics/ReducedDiags/ColliderRelevant.cpp
@@ -181,7 +181,7 @@ ColliderRelevant::ColliderRelevant (const std::string& rd_name)
                 const auto& el = m_headers_indices[name];
                 ofs << m_sep << "[" << el.idx + off << "]" << el.header;
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/DifferentialLuminosity.cpp b/Source/Diagnostics/ReducedDiags/DifferentialLuminosity.cpp
index 111eca0f14c..59a32cf0545 100644
--- a/Source/Diagnostics/ReducedDiags/DifferentialLuminosity.cpp
+++ b/Source/Diagnostics/ReducedDiags/DifferentialLuminosity.cpp
@@ -115,7 +115,7 @@ DifferentialLuminosity::DifferentialLuminosity (const std::string& rd_name)
                 const Real b = m_bin_min + m_bin_size*(Real(i)+0.5_rt);
                 ofs << "bin" << 1+i << "=" << b << "(eV)";
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/FieldEnergy.cpp b/Source/Diagnostics/ReducedDiags/FieldEnergy.cpp
index 9c0fac101a2..fbfdaf11017 100644
--- a/Source/Diagnostics/ReducedDiags/FieldEnergy.cpp
+++ b/Source/Diagnostics/ReducedDiags/FieldEnergy.cpp
@@ -67,7 +67,7 @@ FieldEnergy::FieldEnergy (const std::string& rd_name)
                 ofs << m_sep;
                 ofs << "[" << c++ << "]B_lev" + std::to_string(lev) + "(J)";
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/FieldMaximum.cpp b/Source/Diagnostics/ReducedDiags/FieldMaximum.cpp
index 24d9dda3ea6..4fbbd1ec82c 100644
--- a/Source/Diagnostics/ReducedDiags/FieldMaximum.cpp
+++ b/Source/Diagnostics/ReducedDiags/FieldMaximum.cpp
@@ -92,7 +92,7 @@ FieldMaximum::FieldMaximum (const std::string& rd_name)
                 ofs << m_sep;
                 ofs << "[" << c++ << "]max_|B|_lev" + std::to_string(lev) + "(T)";
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/FieldMomentum.cpp b/Source/Diagnostics/ReducedDiags/FieldMomentum.cpp
index 7eb16efecff..72bece9265b 100644
--- a/Source/Diagnostics/ReducedDiags/FieldMomentum.cpp
+++ b/Source/Diagnostics/ReducedDiags/FieldMomentum.cpp
@@ -87,7 +87,7 @@ FieldMomentum::FieldMomentum (const std::string& rd_name)
                 ofs << "momentum_z_lev" << lev << "(kg*m/s)";
             }
 
-            ofs << std::endl;
+            ofs << "\n";
             ofs.close();
         }
     }
diff --git a/Source/Diagnostics/ReducedDiags/FieldProbe.cpp b/Source/Diagnostics/ReducedDiags/FieldProbe.cpp
index 1113075f4cd..f498cea7566 100644
--- a/Source/Diagnostics/ReducedDiags/FieldProbe.cpp
+++ b/Source/Diagnostics/ReducedDiags/FieldProbe.cpp
@@ -225,7 +225,7 @@ FieldProbe::FieldProbe (const std::string& rd_name)
                 ofs << m_sep;
                 ofs << "[" << c++ << "]part_S_lev" + std::to_string(lev) + u_map[FieldProbePIdx::S];
             }
-            ofs << std::endl;
+            ofs << "\n";
 
             // close file
             ofs.close();
@@ -678,7 +678,7 @@ void FieldProbe::WriteToFile (int step) const
             ofs << m_sep;
             ofs << sorted_data[i * noutputs + k];
         }
-        ofs << std::endl;
+        ofs << "\n";
     } // end loop over data size
     // close file
     ofs.close();
diff --git a/Source/Diagnostics/ReducedDiags/FieldReduction.cpp b/Source/Diagnostics/ReducedDiags/FieldReduction.cpp
index 9d0521eb181..fbf67be4e4a 100644
--- a/Source/Diagnostics/ReducedDiags/FieldReduction.cpp
+++ b/Source/Diagnostics/ReducedDiags/FieldReduction.cpp
@@ -77,7 +77,7 @@ FieldReduction::FieldReduction (const std::string& rd_name)
             ofs << m_sep;
             ofs << "[" << c++ << "]" + amrex::getEnumNameString(m_reduction_type) + " of " + parser_string + " (SI units)";
 
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/LoadBalanceCosts.cpp b/Source/Diagnostics/ReducedDiags/LoadBalanceCosts.cpp
index 882c94b12eb..a11a8d35853 100644
--- a/Source/Diagnostics/ReducedDiags/LoadBalanceCosts.cpp
+++ b/Source/Diagnostics/ReducedDiags/LoadBalanceCosts.cpp
@@ -275,7 +275,7 @@ void LoadBalanceCosts::WriteToFile (int step) const
     // end loop over data size
 
     // end line
-    ofs << std::endl;
+    ofs << "\n";
 
     // close file
     ofs.close();
@@ -329,7 +329,7 @@ void LoadBalanceCosts::WriteToFile (int step) const
             ofstmp << m_sep;
             ofstmp << "[" << c++ << "]hostname_box_" + std::to_string(boxNumber) + "()";
         }
-        ofstmp << std::endl;
+        ofstmp << "\n";
 
         // open the data-containing file
         const std::string fileDataName = m_path + m_rd_name + "." + m_extension;
@@ -360,7 +360,7 @@ void LoadBalanceCosts::WriteToFile (int step) const
             {
                 ofstmp << m_sep << "NaN";
             }
-            ofstmp << std::endl;
+            ofstmp << "\n";
         }
 
         // close files
diff --git a/Source/Diagnostics/ReducedDiags/LoadBalanceEfficiency.cpp b/Source/Diagnostics/ReducedDiags/LoadBalanceEfficiency.cpp
index b6f3a58e670..6c1f46a6c54 100644
--- a/Source/Diagnostics/ReducedDiags/LoadBalanceEfficiency.cpp
+++ b/Source/Diagnostics/ReducedDiags/LoadBalanceEfficiency.cpp
@@ -50,7 +50,7 @@ LoadBalanceEfficiency::LoadBalanceEfficiency (const std::string& rd_name)
                 ofs << m_sep;
                 ofs << "[" << c++ << "]lev" + std::to_string(lev);
             }
-            ofs << std::endl;
+            ofs << "\n";
 
             // close file
             ofs.close();
diff --git a/Source/Diagnostics/ReducedDiags/ParticleEnergy.cpp b/Source/Diagnostics/ReducedDiags/ParticleEnergy.cpp
index 9446af5e20c..41b58ea8fe3 100644
--- a/Source/Diagnostics/ReducedDiags/ParticleEnergy.cpp
+++ b/Source/Diagnostics/ReducedDiags/ParticleEnergy.cpp
@@ -78,7 +78,7 @@ ParticleEnergy::ParticleEnergy (const std::string& rd_name)
                 ofs << m_sep;
                 ofs << "[" << c++ << "]" << species_names[i] + "_mean(J)";
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/ParticleExtrema.cpp b/Source/Diagnostics/ReducedDiags/ParticleExtrema.cpp
index 811fee7a9de..842dd3c9efd 100644
--- a/Source/Diagnostics/ReducedDiags/ParticleExtrema.cpp
+++ b/Source/Diagnostics/ReducedDiags/ParticleExtrema.cpp
@@ -140,7 +140,7 @@ ParticleExtrema::ParticleExtrema (const std::string& rd_name)
                     const auto& el = m_headers_indices[name];
                     ofs << m_sep << "[" << el.idx + off << "]" << el.header;
                 }
-                ofs << std::endl;
+                ofs << "\n";
                 // close file
                 ofs.close();
             }
diff --git a/Source/Diagnostics/ReducedDiags/ParticleHistogram.cpp b/Source/Diagnostics/ReducedDiags/ParticleHistogram.cpp
index 4aeb2af526f..8ba78b30cf2 100644
--- a/Source/Diagnostics/ReducedDiags/ParticleHistogram.cpp
+++ b/Source/Diagnostics/ReducedDiags/ParticleHistogram.cpp
@@ -146,7 +146,7 @@ ParticleHistogram::ParticleHistogram (const std::string& rd_name)
                 ofs << "bin" + std::to_string(1+i)
                              + "=" + std::to_string(b) + "()";
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/ParticleMomentum.cpp b/Source/Diagnostics/ReducedDiags/ParticleMomentum.cpp
index 47d1ec2cdd4..8cd4a9bc00c 100644
--- a/Source/Diagnostics/ReducedDiags/ParticleMomentum.cpp
+++ b/Source/Diagnostics/ReducedDiags/ParticleMomentum.cpp
@@ -107,7 +107,7 @@ ParticleMomentum::ParticleMomentum (const std::string& rd_name)
                 ofs << species_names[i] + "_mean_z(kg*m/s)";
             }
 
-            ofs << std::endl;
+            ofs << "\n";
             ofs.close();
         }
     }
diff --git a/Source/Diagnostics/ReducedDiags/ParticleNumber.cpp b/Source/Diagnostics/ReducedDiags/ParticleNumber.cpp
index bfdf1daffe3..83dc3d118f6 100644
--- a/Source/Diagnostics/ReducedDiags/ParticleNumber.cpp
+++ b/Source/Diagnostics/ReducedDiags/ParticleNumber.cpp
@@ -75,7 +75,7 @@ ParticleNumber::ParticleNumber (const std::string& rd_name)
                 ofs << m_sep;
                 ofs << "[" << c++ << "]" << species_names[i] + "_weight()";
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/ReducedDiags/ReducedDiags.cpp b/Source/Diagnostics/ReducedDiags/ReducedDiags.cpp
index e3835dde416..ec31d9de81c 100644
--- a/Source/Diagnostics/ReducedDiags/ReducedDiags.cpp
+++ b/Source/Diagnostics/ReducedDiags/ReducedDiags.cpp
@@ -121,7 +121,7 @@ void ReducedDiags::WriteToFile (int step) const
     // end loop over data size
 
     // end line
-    ofs << std::endl;
+    ofs << "\n";
 
     // close file
     ofs.close();
diff --git a/Source/Diagnostics/ReducedDiags/RhoMaximum.cpp b/Source/Diagnostics/ReducedDiags/RhoMaximum.cpp
index 9902415d9a9..66460613a56 100644
--- a/Source/Diagnostics/ReducedDiags/RhoMaximum.cpp
+++ b/Source/Diagnostics/ReducedDiags/RhoMaximum.cpp
@@ -113,7 +113,7 @@ RhoMaximum::RhoMaximum (const std::string& rd_name)
                                          + "_|rho|_lev" + std::to_string(lev) + "(C/m^3)";
                 }
             }
-            ofs << std::endl;
+            ofs << "\n";
             // close file
             ofs.close();
         }
diff --git a/Source/Diagnostics/WarpXOpenPMD.cpp b/Source/Diagnostics/WarpXOpenPMD.cpp
index 36827bd316a..4e1b6238adf 100644
--- a/Source/Diagnostics/WarpXOpenPMD.cpp
+++ b/Source/Diagnostics/WarpXOpenPMD.cpp
@@ -466,7 +466,7 @@ void WarpXOpenPMDPlot::CloseStep (bool isBTD, bool isLastBTDFlush)
             std::string const filename = GetFileName(filepath);
 
             std::ofstream pv_helper_file(m_dirPrefix + "/paraview.pmd");
-            pv_helper_file << filename << std::endl;
+            pv_helper_file << filename << "\n";
             pv_helper_file.close();
         }
     }
diff --git a/Source/FieldSolver/ImplicitSolvers/SemiImplicitEM.cpp b/Source/FieldSolver/ImplicitSolvers/SemiImplicitEM.cpp
index cfd18354878..33e2ced7b53 100644
--- a/Source/FieldSolver/ImplicitSolvers/SemiImplicitEM.cpp
+++ b/Source/FieldSolver/ImplicitSolvers/SemiImplicitEM.cpp
@@ -36,21 +36,20 @@ void SemiImplicitEM::Define ( WarpX*  a_WarpX )
 void SemiImplicitEM::PrintParameters () const
 {
     if (!m_WarpX->Verbose()) { return; }
-    amrex::Print() << std::endl;
-    amrex::Print() << "-----------------------------------------------------------" << std::endl;
-    amrex::Print() << "----------- SEMI IMPLICIT EM SOLVER PARAMETERS ------------" << std::endl;
-    amrex::Print() << "-----------------------------------------------------------" << std::endl;
-    amrex::Print() << "max particle iterations:    " << m_max_particle_iterations << std::endl;
-    amrex::Print() << "particle tolerance:         " << m_particle_tolerance << std::endl;
+    amrex::Print() << "\n";
+    amrex::Print() << "-----------------------------------------------------------\n";
+    amrex::Print() << "----------- SEMI IMPLICIT EM SOLVER PARAMETERS ------------\n";
+    amrex::Print() << "-----------------------------------------------------------\n";
+    amrex::Print() << "max particle iterations:    " << m_max_particle_iterations << "\n";
+    amrex::Print() << "particle tolerance:         " << m_particle_tolerance << "\n";
     if (m_nlsolver_type==NonlinearSolverType::Picard) {
-        amrex::Print() << "Nonlinear solver type:      Picard" << std::endl;
+        amrex::Print() << "Nonlinear solver type:      Picard\n";
     }
     else if (m_nlsolver_type==NonlinearSolverType::Newton) {
-        amrex::Print() << "Nonlinear solver type:      Newton" << std::endl;
+        amrex::Print() << "Nonlinear solver type:      Newton\n";
     }
     m_nlsolver->PrintParams();
-    amrex::Print() << "-----------------------------------------------------------" << std::endl;
-    amrex::Print() << std::endl;
+    amrex::Print() << "-----------------------------------------------------------\n\n";
 }
 
 void SemiImplicitEM::OneStep ( amrex::Real  a_time,
diff --git a/Source/FieldSolver/ImplicitSolvers/ThetaImplicitEM.cpp b/Source/FieldSolver/ImplicitSolvers/ThetaImplicitEM.cpp
index 4c86389797f..3d74ddfde69 100644
--- a/Source/FieldSolver/ImplicitSolvers/ThetaImplicitEM.cpp
+++ b/Source/FieldSolver/ImplicitSolvers/ThetaImplicitEM.cpp
@@ -56,22 +56,21 @@ void ThetaImplicitEM::Define ( WarpX* const  a_WarpX )
 void ThetaImplicitEM::PrintParameters () const
 {
     if (!m_WarpX->Verbose()) { return; }
-    amrex::Print() << std::endl;
-    amrex::Print() << "-----------------------------------------------------------" << std::endl;
-    amrex::Print() << "----------- THETA IMPLICIT EM SOLVER PARAMETERS -----------" << std::endl;
-    amrex::Print() << "-----------------------------------------------------------" << std::endl;
-    amrex::Print() << "Time-bias parameter theta:  " << m_theta << std::endl;
-    amrex::Print() << "max particle iterations:    " << m_max_particle_iterations << std::endl;
-    amrex::Print() << "particle tolerance:         " << m_particle_tolerance << std::endl;
+    amrex::Print() << "\n";
+    amrex::Print() << "-----------------------------------------------------------\n";
+    amrex::Print() << "----------- THETA IMPLICIT EM SOLVER PARAMETERS -----------\n";
+    amrex::Print() << "-----------------------------------------------------------\n";
+    amrex::Print() << "Time-bias parameter theta:  " << m_theta << "\n";
+    amrex::Print() << "max particle iterations:    " << m_max_particle_iterations << "\n";
+    amrex::Print() << "particle tolerance:         " << m_particle_tolerance << "\n";
     if (m_nlsolver_type==NonlinearSolverType::Picard) {
-        amrex::Print() << "Nonlinear solver type:      Picard" << std::endl;
+        amrex::Print() << "Nonlinear solver type:      Picard\n";
     }
     else if (m_nlsolver_type==NonlinearSolverType::Newton) {
-        amrex::Print() << "Nonlinear solver type:      Newton" << std::endl;
+        amrex::Print() << "Nonlinear solver type:      Newton\n";
     }
     m_nlsolver->PrintParams();
-    amrex::Print() << "-----------------------------------------------------------" << std::endl;
-    amrex::Print() << std::endl;
+    amrex::Print() << "-----------------------------------------------------------\n\n";
 }
 
 void ThetaImplicitEM::OneStep ( const amrex::Real  a_time,
diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp
index de763831d98..49b0d439c50 100644
--- a/Source/Initialization/WarpXInitData.cpp
+++ b/Source/Initialization/WarpXInitData.cpp
@@ -830,7 +830,7 @@ WarpX::computeMaxStepBoostAccelerator() {
         static_cast<int>(interaction_time_boost/dt[maxLevel()]);
     max_step = computed_max_step;
     Print()<<"max_step computed in computeMaxStepBoostAccelerator: "
-           <<max_step<<std::endl;
+           <<max_step<<"\n";
 }
 
 void
diff --git a/Source/NonlinearSolvers/NewtonSolver.H b/Source/NonlinearSolvers/NewtonSolver.H
index 93ad432208a..814ef2b8dc3 100644
--- a/Source/NonlinearSolvers/NewtonSolver.H
+++ b/Source/NonlinearSolvers/NewtonSolver.H
@@ -69,16 +69,16 @@ public:
 
     void PrintParams () const override
     {
-        amrex::Print()     << "Newton verbose:             " << (this->m_verbose?"true":"false") << std::endl;
-        amrex::Print()     << "Newton max iterations:      " << m_maxits << std::endl;
-        amrex::Print()     << "Newton relative tolerance:  " << m_rtol << std::endl;
-        amrex::Print()     << "Newton absolute tolerance:  " << m_atol << std::endl;
-        amrex::Print()     << "Newton require convergence: " << (m_require_convergence?"true":"false") << std::endl;
-        amrex::Print()     << "GMRES verbose:            " << m_gmres_verbose_int << std::endl;
-        amrex::Print()     << "GMRES restart length:     " << m_gmres_restart_length << std::endl;
-        amrex::Print()     << "GMRES max iterations:     " << m_gmres_maxits << std::endl;
-        amrex::Print()     << "GMRES relative tolerance: " << m_gmres_rtol << std::endl;
-        amrex::Print()     << "GMRES absolute tolerance: " << m_gmres_atol << std::endl;
+        amrex::Print()     << "Newton verbose:             " << (this->m_verbose?"true":"false") << "\n";
+        amrex::Print()     << "Newton max iterations:      " << m_maxits << "\n";
+        amrex::Print()     << "Newton relative tolerance:  " << m_rtol << "\n";
+        amrex::Print()     << "Newton absolute tolerance:  " << m_atol << "\n";
+        amrex::Print()     << "Newton require convergence: " << (m_require_convergence?"true":"false") << "\n";
+        amrex::Print()     << "GMRES verbose:            " << m_gmres_verbose_int << "\n";
+        amrex::Print()     << "GMRES restart length:     " << m_gmres_restart_length << "\n";
+        amrex::Print()     << "GMRES max iterations:     " << m_gmres_maxits << "\n";
+        amrex::Print()     << "GMRES relative tolerance: " << m_gmres_rtol << "\n";
+        amrex::Print()     << "GMRES absolute tolerance: " << m_gmres_atol << "\n";
     }
 
 private:
@@ -261,19 +261,19 @@ void NewtonSolver<Vec,Ops>::Solve ( Vec&         a_U,
 
         if (norm_abs < m_rtol) {
             amrex::Print() << "Newton: exiting at iteration = " << std::setw(3) << iter
-                           << ". Satisfied absolute tolerance " << m_atol << std::endl;
+                           << ". Satisfied absolute tolerance " << m_atol << "\n";
             break;
         }
 
         if (norm_rel < m_rtol) {
             amrex::Print() << "Newton: exiting at iteration = " << std::setw(3) << iter
-                           << ". Satisfied relative tolerance " << m_rtol << std::endl;
+                           << ". Satisfied relative tolerance " << m_rtol << "\n";
             break;
         }
 
         if (norm_abs > 100._rt*norm0) {
             amrex::Print() << "Newton: exiting at iteration = " << std::setw(3) << iter
-                           << ". SOLVER DIVERGED! relative tolerance = " << m_rtol << std::endl;
+                           << ". SOLVER DIVERGED! relative tolerance = " << m_rtol << "\n";
             std::stringstream convergenceMsg;
             convergenceMsg << "Newton: exiting at iteration " << std::setw(3) << iter <<
                               ". SOLVER DIVERGED! absolute norm = " << norm_abs <<
@@ -291,7 +291,7 @@ void NewtonSolver<Vec,Ops>::Solve ( Vec&         a_U,
         iter++;
         if (iter >= m_maxits) {
             amrex::Print() << "Newton: exiting at iter = " << std::setw(3) << iter
-                           << ". Maximum iteration reached: iter = " << m_maxits << std::endl;
+                           << ". Maximum iteration reached: iter = " << m_maxits << "\n";
             break;
         }
 
diff --git a/Source/NonlinearSolvers/PicardSolver.H b/Source/NonlinearSolvers/PicardSolver.H
index f05b9a106e6..4eed4d6c2e0 100644
--- a/Source/NonlinearSolvers/PicardSolver.H
+++ b/Source/NonlinearSolvers/PicardSolver.H
@@ -55,10 +55,10 @@ public:
 
     void PrintParams () const override
     {
-        amrex::Print() << "Picard max iterations:      " << m_maxits << std::endl;
-        amrex::Print() << "Picard relative tolerance:  " << m_rtol << std::endl;
-        amrex::Print() << "Picard absolute tolerance:  " << m_atol << std::endl;
-        amrex::Print() << "Picard require convergence: " << (m_require_convergence?"true":"false") << std::endl;
+        amrex::Print() << "Picard max iterations:      " << m_maxits << "\n";
+        amrex::Print() << "Picard relative tolerance:  " << m_rtol << "\n";
+        amrex::Print() << "Picard absolute tolerance:  " << m_atol << "\n";
+        amrex::Print() << "Picard require convergence: " << (m_require_convergence?"true":"false") << "\n";
     }
 
 private:
@@ -179,19 +179,19 @@ void PicardSolver<Vec,Ops>::Solve ( Vec&         a_U,
 
         if (norm_abs < m_atol) {
             amrex::Print() << "Picard: exiting at iter = " << std::setw(3) << iter
-                           << ". Satisfied absolute tolerance " << m_atol << std::endl;
+                           << ". Satisfied absolute tolerance " << m_atol << "\n";
             break;
         }
 
         if (norm_rel < m_rtol) {
             amrex::Print() << "Picard: exiting at iter = " << std::setw(3) << iter
-                           << ". Satisfied relative tolerance " << m_rtol << std::endl;
+                           << ". Satisfied relative tolerance " << m_rtol << "\n";
             break;
         }
 
         if (iter >= m_maxits) {
             amrex::Print() << "Picard: exiting at iter = " << std::setw(3) << iter
-                           << ". Maximum iteration reached: iter = " << m_maxits << std::endl;
+                           << ". Maximum iteration reached: iter = " << m_maxits << "\n";
             break;
         }
 
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
index d1a19f06993..0617b36a273 100644
--- a/Source/Particles/PhysicalParticleContainer.cpp
+++ b/Source/Particles/PhysicalParticleContainer.cpp
@@ -2976,10 +2976,10 @@ PhysicalParticleContainer::ImplicitPushXP (WarpXParIter& pti,
 #if !defined(AMREX_USE_GPU)
                 std::stringstream convergenceMsg;
                 convergenceMsg << "Picard solver for particle failed to converge after " <<
-                    iter << " iterations. " << std::endl;
+                    iter << " iterations.\n";
                 convergenceMsg << "Position step norm is " << step_norm <<
-                    " and the tolerance is " << particle_tolerance << std::endl;
-                convergenceMsg << " ux = " << ux[ip] << ", uy = " << uy[ip] << ", uz = " << uz[ip] << std::endl;
+                    " and the tolerance is " << particle_tolerance << "\n";
+                convergenceMsg << " ux = " << ux[ip] << ", uy = " << uy[ip] << ", uz = " << uz[ip] << "\n";
                 convergenceMsg << " xp = " << xp     << ", yp = " << yp     << ", zp = " << zp;
                 ablastr::warn_manager::WMRecordWarning("ImplicitPushXP", convergenceMsg.str());
 #endif
diff --git a/Tools/QedTablesUtils/Source/ArgParser/QedTablesArgParser.cpp b/Tools/QedTablesUtils/Source/ArgParser/QedTablesArgParser.cpp
index 41d27477dcf..a92f191037d 100644
--- a/Tools/QedTablesUtils/Source/ArgParser/QedTablesArgParser.cpp
+++ b/Tools/QedTablesUtils/Source/ArgParser/QedTablesArgParser.cpp
@@ -85,7 +85,7 @@ ArgParser::ParseArgs (const std::vector<Key>& keys, const int argc, char const*
 void
 ArgParser::PrintHelp (const vector<ArgParser::Key>& cmd_list)
 {
-    cout << "Command line options: " << endl;
+    cout << "Command line options:\n";
 
     for (const auto& el : cmd_list){
         const auto type = get<1>(el);
@@ -102,7 +102,7 @@ ArgParser::PrintHelp (const vector<ArgParser::Key>& cmd_list)
         cout << get<0>(el)
             << " " << stype
             << " " << get<2>(el)
-            << endl;
+            << "\n";
     }
 
 }
diff --git a/Tools/QedTablesUtils/Source/QedTableCommons.H b/Tools/QedTablesUtils/Source/QedTableCommons.H
index 40551b9e13c..2233513bc97 100644
--- a/Tools/QedTablesUtils/Source/QedTableCommons.H
+++ b/Tools/QedTablesUtils/Source/QedTableCommons.H
@@ -19,7 +19,7 @@ void AbortWithMessage(const std::string& msg)
 
 void SuccessExit()
 {
-    std::cout << "___________________________" << std::endl;
+    std::cout << "___________________________\n";
     exit(0);
 }
 
diff --git a/Tools/QedTablesUtils/Source/QedTableGenerator.cpp b/Tools/QedTablesUtils/Source/QedTableGenerator.cpp
index 1ea62b5c6ed..7bfa5787ec8 100644
--- a/Tools/QedTablesUtils/Source/QedTableGenerator.cpp
+++ b/Tools/QedTablesUtils/Source/QedTableGenerator.cpp
@@ -49,7 +49,7 @@ void GenerateTableQS (const ParsedArgs& args, const string& outfile_name);
 
 int main (int argc, char** argv)
 {
-    cout << "### QED Table Generator ###" << endl;
+    cout << "### QED Table Generator ###\n";
     const auto args_map = ParseArgs(line_commands, argc, argv);
 
     if (args_map.empty() || Contains(args_map, "-h")){
diff --git a/Tools/QedTablesUtils/Source/QedTableReader.cpp b/Tools/QedTablesUtils/Source/QedTableReader.cpp
index ba9d58775f2..27b284f34c2 100644
--- a/Tools/QedTablesUtils/Source/QedTableReader.cpp
+++ b/Tools/QedTablesUtils/Source/QedTableReader.cpp
@@ -57,7 +57,7 @@ class qs_photon_emission_table_wrapper :
 
 int main (int argc, char** argv)
 {
-    cout << "### QED Table Reader ###" << endl;
+    cout << "### QED Table Reader ###\n";
     const auto args_map = ParseArgs(line_commands, argc, argv);
 
     if (args_map.empty() || Contains(args_map, "-h")){