Electrostatic lab frame (#1531)

* Implemented lab frame electrostatic model

* Small clean up of ComputeSpaceChargeField

* Added ElectrostaticSphereLabFrame regression test

1 files changed, 3 insertions, 3 deletions

diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
index 51a615234..cd74d0926 100644
--- a/Source/Particles/PhysicalParticleContainer.cpp
+++ b/Source/Particles/PhysicalParticleContainer.cpp
@@ -1107,7 +1107,7 @@ PhysicalParticleContainer::Evolve (int lev,
                 //
                 // Current Deposition (only needed for electromagnetic solver)
                 //
-                if (!WarpX::do_electrostatic) {
+                if (WarpX::do_electrostatic == ElectrostaticSolverAlgo::None) {
                     int* AMREX_RESTRICT ion_lev;
                     if (do_field_ionization){
                         ion_lev = pti.GetiAttribs(particle_icomps["ionization_level"]).dataPtr();
@@ -1124,13 +1124,13 @@ PhysicalParticleContainer::Evolve (int lev,
                                        np_current, np-np_current, thread_num,
                                        lev, lev-1, dt);
                     }
-                } // end of "if !do_electrostatic"
+                } // end of "if do_electrostatic == ElectrostaticSolverAlgo::None"
             } // end of "if do_not_push"
 
             if (rho) {
                 // Deposit charge after particle push, in component 1 of MultiFab rho.
                 // (Skipped for electrostatic solver, as this may lead to out-of-bounds)
-                if (!WarpX::do_electrostatic) {
+                if (WarpX::do_electrostatic == ElectrostaticSolverAlgo::None) {
                     int* AMREX_RESTRICT ion_lev;
                     if (do_field_ionization){
                         ion_lev = pti.GetiAttribs(particle_icomps["ionization_level"]).dataPtr();