diff options
author | 2019-07-08 10:24:40 -0700 | |
---|---|---|
committer | 2019-07-08 10:24:40 -0700 | |
commit | 20c785a67312b88e497a94e1a27a343918cdc9ce (patch) | |
tree | 96053ff7651e6d14376013b53fbf2ee71b297835 /Source/Particles/PhysicalParticleContainer.cpp | |
parent | af336d881e641b43df054d2701e6e0095fe48d8c (diff) | |
download | WarpX-20c785a67312b88e497a94e1a27a343918cdc9ce.tar.gz WarpX-20c785a67312b88e497a94e1a27a343918cdc9ce.tar.zst WarpX-20c785a67312b88e497a94e1a27a343918cdc9ce.zip |
Cleaning code for slices in boosted frame diag
Diffstat (limited to 'Source/Particles/PhysicalParticleContainer.cpp')
-rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index 7e7c9534e..27155e1c6 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -1897,7 +1897,7 @@ void PhysicalParticleContainer::GetParticleSlice(const int direction, const Real for (int lev = 0; lev < nlevs; ++lev) { - const Real* dx = Geom(lev).CellSize(); + const Real* dx = Geom(lev).CellSize(); const Real* plo = Geom(lev).ProbLo(); // first we touch each map entry in serial @@ -1951,7 +1951,7 @@ void PhysicalParticleContainer::GetParticleSlice(const int direction, const Real // timestep, skip it. if ( not (((zp_new[i] >= z_new) && (zp_old[i] <= z_old)) || ((zp_new[i] <= z_new) && (zp_old[i] >= z_old))) ) continue; - + // Lorentz transform particles to lab frame Real gamma_new_p = std::sqrt(1.0 + inv_c2*(uxp_new[i]*uxp_new[i] + uyp_new[i]*uyp_new[i] + uzp_new[i]*uzp_new[i])); Real t_new_p = WarpX::gamma_boost*t_boost - uzfrm*zp_new[i]*inv_c2; |