diff options
Diffstat (limited to 'Source/Particles/PhysicalParticleContainer.cpp')
| -rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 24 |
1 files changed, 2 insertions, 22 deletions
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index 0efb6f61c..8845fb9e2 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -342,18 +342,6 @@ PhysicalParticleContainer::PhysicalParticleContainer (AmrCore* amr_core, int isp pp_boundary.query("reflect_all_velocities", flag); m_boundary_conditions.Set_reflect_all_velocities(flag); - // Get Galilean velocity - ParmParse pp_psatd("psatd"); - bool use_default_v_galilean = false; - pp_psatd.query("use_default_v_galilean", use_default_v_galilean); - if (use_default_v_galilean) { - m_v_galilean[2] = -std::sqrt(1._rt - 1._rt / (WarpX::gamma_boost * WarpX::gamma_boost)); - } else { - queryArrWithParser(pp_psatd, "v_galilean", m_v_galilean, 0, 3); - } - // Scale the Galilean velocity by the speed of light - for (int i=0; i<3; i++) m_v_galilean[i] *= PhysConst::c; - } PhysicalParticleContainer::PhysicalParticleContainer (AmrCore* amr_core) @@ -2223,7 +2211,7 @@ PhysicalParticleContainer::PushP (int lev, Real dt, const auto getExternalEB = GetExternalEBField(pti); - const auto& xyzmin = WarpX::GetInstance().LowerCornerWithGalilean(box,m_v_galilean,lev); + const std::array<amrex::Real,3>& xyzmin = WarpX::LowerCorner(box, lev, 0._rt); const Dim3 lo = lbound(box); @@ -2576,14 +2564,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, const auto getExternalEB = GetExternalEBField(pti, offset); // Lower corner of tile box physical domain (take into account Galilean shift) - Real cur_time = WarpX::GetInstance().gett_new(lev); - const auto& time_of_last_gal_shift = WarpX::GetInstance().time_of_last_gal_shift; - Real time_shift = (cur_time - time_of_last_gal_shift); - amrex::Array<amrex::Real,3> galilean_shift ={ - m_v_galilean[0]*time_shift, - m_v_galilean[1]*time_shift, - m_v_galilean[2]*time_shift }; - const std::array<Real, 3>& xyzmin = WarpX::LowerCorner(box, galilean_shift, gather_lev); + const std::array<amrex::Real, 3>& xyzmin = WarpX::LowerCorner(box, gather_lev, 0._rt); const Dim3 lo = lbound(box); @@ -2797,7 +2778,6 @@ PhysicalParticleContainer::getIonizationFunc (const WarpXParIter& pti, WARPX_PROFILE("PhysicalParticleContainer::getIonizationFunc()"); return IonizationFilterFunc(pti, lev, ngEB, Ex, Ey, Ez, Bx, By, Bz, - m_v_galilean, ionization_energies.dataPtr(), adk_prefactor.dataPtr(), adk_exp_prefactor.dataPtr(), |