diff options
| author |  David Grote <grote1@llnl.gov>
| 2022-02-04 17:12:55 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2022-02-05 01:12:55 +0000 |
| commit | 5692067b01252146b8d0a90137a1bd952644c8ca (patch) | |
| tree | bfb704a14588a54207b883254264a8c474d65a14 /Source/Particles/Collision/MCCProcess.cpp | |
| parent | 9263e6a9832f7ec69a8cc54fda3f6f1de51caaf1 (diff) | |
| download | WarpX-5692067b01252146b8d0a90137a1bd952644c8ca.tar.gz WarpX-5692067b01252146b8d0a90137a1bd952644c8ca.tar.zst WarpX-5692067b01252146b8d0a90137a1bd952644c8ca.zip | |
Fixes removing implicit double to float conversions (#2687)
* Fixes removing implicit double to float conversions
* Minor fix in Evolve
* Removed comments
* Reverted amrex::second to original version.
This change will be made in a separate PR
* More fixes
* Performance: pow2 & pow3
* One more pow2
* Cleaning: auto where types are repetitive
std::pow(float, int): see comments and
https://stackoverflow.com/a/5627278/2719194
Diffstat (limited to 'Source/Particles/Collision/MCCProcess.cpp')
| -rw-r--r-- | Source/Particles/Collision/MCCProcess.cpp | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/Source/Particles/Collision/MCCProcess.cpp b/Source/Particles/Collision/MCCProcess.cpp index a059dd2f5..9f027d056 100644 --- a/Source/Particles/Collision/MCCProcess.cpp +++ b/Source/Particles/Collision/MCCProcess.cpp @@ -37,6 +37,7 @@ MCCProcess::MCCProcess ( void MCCProcess::init (const std::string& scattering_process, const amrex::Real energy) { + using namespace amrex::literals; m_exe_h.m_sigmas_data = m_sigmas_h.data(); // save energy grid parameters for easy use @@ -45,7 +46,7 @@ MCCProcess::init (const std::string& scattering_process, const amrex::Real energ m_exe_h.m_energy_hi = m_energies[m_grid_size-1]; m_exe_h.m_sigma_lo = m_sigmas_h[0]; m_exe_h.m_sigma_hi = m_sigmas_h[m_grid_size-1]; - m_exe_h.m_dE = (m_exe_h.m_energy_hi - m_exe_h.m_energy_lo)/(m_grid_size - 1.); + m_exe_h.m_dE = (m_exe_h.m_energy_hi - m_exe_h.m_energy_lo)/(m_grid_size - 1._rt); m_exe_h.m_energy_penalty = energy; m_exe_h.m_type = parseProcessType(scattering_process); @@ -99,7 +100,7 @@ MCCProcess::readCrossSectionFile ( std::ifstream infile(cross_section_file); if(!infile.is_open()) amrex::Abort("Failed to open cross-section data file"); - double energy, sigma; + amrex::Real energy, sigma; while (infile >> energy >> sigma) { energies.push_back(energy); sigmas.push_back(sigma); |