aboutsummaryrefslogtreecommitdiff
path: root/Source/Particles/Collision/BackgroundMCC/MCCProcess.H
diff options
context:
space:
mode:
Diffstat (limited to 'Source/Particles/Collision/BackgroundMCC/MCCProcess.H')
-rw-r--r--Source/Particles/Collision/BackgroundMCC/MCCProcess.H40
1 files changed, 20 insertions, 20 deletions
diff --git a/Source/Particles/Collision/BackgroundMCC/MCCProcess.H b/Source/Particles/Collision/BackgroundMCC/MCCProcess.H
index 68ab61aba..0b28adf75 100644
--- a/Source/Particles/Collision/BackgroundMCC/MCCProcess.H
+++ b/Source/Particles/Collision/BackgroundMCC/MCCProcess.H
@@ -27,7 +27,7 @@ public:
MCCProcess (
const std::string& scattering_process,
const std::string& cross_section_file,
- const amrex::Real energy
+ const amrex::ParticleReal energy
);
template <typename InputVector>
@@ -35,7 +35,7 @@ public:
const std::string& scattering_process,
const InputVector&& energies,
const InputVector&& sigmas,
- const amrex::Real energy
+ const amrex::ParticleReal energy
);
MCCProcess (MCCProcess const&) = delete;
@@ -55,14 +55,14 @@ public:
static
void readCrossSectionFile (
const std::string cross_section_file,
- amrex::Vector<amrex::Real>& energies,
- amrex::Gpu::HostVector<amrex::Real>& sigmas
+ amrex::Vector<amrex::ParticleReal>& energies,
+ amrex::Gpu::HostVector<amrex::ParticleReal>& sigmas
);
static
void sanityCheckEnergyGrid (
- const amrex::Vector<amrex::Real>& energies,
- amrex::Real dE
+ const amrex::Vector<amrex::ParticleReal>& energies,
+ amrex::ParticleReal dE
);
struct Executor {
@@ -74,7 +74,7 @@ public:
*
*/
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
- amrex::Real getCrossSection (amrex::Real E_coll) const
+ amrex::ParticleReal getCrossSection (amrex::ParticleReal E_coll) const
{
if (E_coll < m_energy_lo) {
return m_sigma_lo;
@@ -84,7 +84,7 @@ public:
using amrex::Math::floor;
using amrex::Math::ceil;
// calculate index of bounding energy pairs
- amrex::Real temp = (E_coll - m_energy_lo) / m_dE;
+ amrex::ParticleReal temp = (E_coll - m_energy_lo) / m_dE;
int idx_1 = static_cast<int>(floor(temp));
int idx_2 = static_cast<int>(ceil(temp));
@@ -94,9 +94,9 @@ public:
}
}
- amrex::Real* m_sigmas_data = nullptr;
- amrex::Real m_energy_lo, m_energy_hi, m_sigma_lo, m_sigma_hi, m_dE;
- amrex::Real m_energy_penalty;
+ amrex::ParticleReal* m_sigmas_data = nullptr;
+ amrex::ParticleReal m_energy_lo, m_energy_hi, m_sigma_lo, m_sigma_hi, m_dE;
+ amrex::ParticleReal m_energy_penalty;
MCCProcessType m_type;
};
@@ -108,15 +108,15 @@ public:
#endif
}
- amrex::Real getCrossSection (amrex::Real E_coll) const
+ amrex::ParticleReal getCrossSection (amrex::ParticleReal E_coll) const
{
return m_exe_h.getCrossSection(E_coll);
}
- amrex::Real getEnergyPenalty () const { return m_exe_h.m_energy_penalty; }
- amrex::Real getMinEnergyInput () const { return m_exe_h.m_energy_lo; }
- amrex::Real getMaxEnergyInput () const { return m_exe_h.m_energy_hi; }
- amrex::Real getEnergyInputStep () const { return m_exe_h.m_dE; }
+ amrex::ParticleReal getEnergyPenalty () const { return m_exe_h.m_energy_penalty; }
+ amrex::ParticleReal getMinEnergyInput () const { return m_exe_h.m_energy_lo; }
+ amrex::ParticleReal getMaxEnergyInput () const { return m_exe_h.m_energy_hi; }
+ amrex::ParticleReal getEnergyInputStep () const { return m_exe_h.m_dE; }
MCCProcessType type () const { return m_exe_h.m_type; }
@@ -125,15 +125,15 @@ private:
static
MCCProcessType parseProcessType(const std::string& process);
- void init (const std::string& scattering_process, const amrex::Real energy);
+ void init (const std::string& scattering_process, const amrex::ParticleReal energy);
- amrex::Vector<amrex::Real> m_energies;
+ amrex::Vector<amrex::ParticleReal> m_energies;
#ifdef AMREX_USE_GPU
- amrex::Gpu::DeviceVector<amrex::Real> m_sigmas_d;
+ amrex::Gpu::DeviceVector<amrex::ParticleReal> m_sigmas_d;
Executor m_exe_d;
#endif
- amrex::Gpu::HostVector<amrex::Real> m_sigmas_h;
+ amrex::Gpu::HostVector<amrex::ParticleReal> m_sigmas_h;
Executor m_exe_h;
int m_grid_size;