diff options
| -rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.cpp | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.cpp index dd608fa18..a76311db3 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.cpp @@ -142,7 +142,8 @@ GalileanPsatdAlgorithmRZ::pushSpectralFields (SpectralFieldDataRZ & f) Complex const divE = kr*(Ep_old - Em_old) + I*kz*Ez_old; Complex const divJ = kr*(Jp - Jm) + I*kz*Jz; - rho_diff = T2*(X2 - X3)*PhysConst::ep0*divE + T_rho*X2*divJ; + auto const myeps0 = PhysConst::ep0; // temporary for NVCC + rho_diff = T2*(X2 - X3)*myeps0*divE + T_rho*X2*divJ; } // Update E (see WarpX online documentation: theory section) |