diff options
| author |  Luca Fedeli <luca.fedeli@cea.fr>
| 2021-01-12 20:42:59 +0100 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2021-01-12 11:42:59 -0800 |
| commit | cb8dc2b37bb31094a686a377975eccd75b551244 (patch) | |
| tree | 31692eb529b2227da0ca47d52ef0cf177cce2945 /Source/Particles/PhysicalParticleContainer.cpp | |
| parent | 3f20f88bddafd98d87f0a22fd173fd27dc6272a4 (diff) | |
| download | WarpX-cb8dc2b37bb31094a686a377975eccd75b551244.tar.gz WarpX-cb8dc2b37bb31094a686a377975eccd75b551244.tar.zst WarpX-cb8dc2b37bb31094a686a377975eccd75b551244.zip | |
Fix some other warnings (#1618)
* fix warnings
* fixed bug
Diffstat (limited to '')
| -rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 16 |
1 files changed, 9 insertions, 7 deletions
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index 417f76bef..9b8184823 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -316,19 +316,19 @@ PhysicalParticleContainer::AddGaussianBeam ( u.z *= PhysConst::c; if (do_symmetrize){ // Add four particles to the beam: - CheckAndAddParticle(x, y, z, { u.x, u.y, u.z}, weight/4., + CheckAndAddParticle(x, y, z, { u.x, u.y, u.z}, weight/4._rt, particle_x, particle_y, particle_z, particle_ux, particle_uy, particle_uz, particle_w); - CheckAndAddParticle(x, -y, z, { u.x, -u.y, u.z}, weight/4., + CheckAndAddParticle(x, -y, z, { u.x, -u.y, u.z}, weight/4._rt, particle_x, particle_y, particle_z, particle_ux, particle_uy, particle_uz, particle_w); - CheckAndAddParticle(-x, y, z, { -u.x, u.y, u.z}, weight/4., + CheckAndAddParticle(-x, y, z, { -u.x, u.y, u.z}, weight/4._rt, particle_x, particle_y, particle_z, particle_ux, particle_uy, particle_uz, particle_w); - CheckAndAddParticle(-x, -y, z, { -u.x, -u.y, u.z}, weight/4., + CheckAndAddParticle(-x, -y, z, { -u.x, -u.y, u.z}, weight/4._rt, particle_x, particle_y, particle_z, particle_ux, particle_uy, particle_uz, particle_w); @@ -1941,11 +1941,13 @@ PhysicalParticleContainer::InitIonizationModule () // For now, we assume l=0 and m=0. // The approximate expressions are used, // without Gamma function - Real wa = std::pow(PhysConst::alpha,3) * PhysConst::c / PhysConst::r_e; + constexpr auto a3 = PhysConst::alpha*PhysConst::alpha*PhysConst::alpha; + constexpr auto a4 = a3 * PhysConst::alpha; + Real wa = a3 * PhysConst::c / PhysConst::r_e; Real Ea = PhysConst::m_e * PhysConst::c*PhysConst::c /PhysConst::q_e * - std::pow(PhysConst::alpha,4)/PhysConst::r_e; + a4/PhysConst::r_e; Real UH = table_ionization_energies[0]; - Real l_eff = std::sqrt(UH/h_ionization_energies[0]) - 1.; + Real l_eff = std::sqrt(UH/h_ionization_energies[0]) - 1._rt; const Real dt = WarpX::GetInstance().getdt(0); |