new prefactor for calculation of width_p for momentum in linear

Hi-PACE · Jan 17, 2025 · 56158d6 · 56158d6
1 parent 29c2a53
commit 56158d6
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/particles/plasma/PlasmaParticleContainerInit.cpp b/src/particles/plasma/PlasmaParticleContainerInit.cpp
@@ -431,11 +431,13 @@ InitIonizationModule (const amrex::Geometry& geom, const amrex::Real background_
     m_adk_prefactor.resize(ion_atomic_number);
     m_adk_exp_prefactor.resize(ion_atomic_number);
     m_laser_adk_prefactor.resize(ion_atomic_number);
+    m_laser_dp_prefactor.resize(ion_atomic_number);
 
     amrex::Gpu::PinnedVector<amrex::Real> h_adk_power(ion_atomic_number);
     amrex::Gpu::PinnedVector<amrex::Real> h_adk_prefactor(ion_atomic_number);
     amrex::Gpu::PinnedVector<amrex::Real> h_adk_exp_prefactor(ion_atomic_number);
     amrex::Gpu::PinnedVector<amrex::Real> h_laser_adk_prefactor(ion_atomic_number);
+    amrex::Gpu::PinnedVector<amrex::Real> h_laser_dp_prefactor(ion_atomic_number);
 
     for (int i=0; i<ion_atomic_number; ++i)
     {
@@ -448,6 +450,7 @@ InitIonizationModule (const amrex::Geometry& geom, const amrex::Real background_
             * std::pow(2*std::pow((Uion/UH),3./2.)*Ea,2*n_eff - 1);
         h_adk_exp_prefactor[i] = -2./3. * std::pow( Uion/UH,3./2.) * Ea;
         h_laser_adk_prefactor[i] = (3./MathConst::pi) * std::pow(Uion/UH, -3./2.) / Ea;
+	h_laser_dp_prefactor[i] = 3./2. * std::pow(Uion/UH, -3./2.) / Ea;
     }
 
     amrex::Gpu::copy(amrex::Gpu::hostToDevice,
@@ -458,4 +461,6 @@ InitIonizationModule (const amrex::Geometry& geom, const amrex::Real background_
         h_adk_exp_prefactor.begin(), h_adk_exp_prefactor.end(), m_adk_exp_prefactor.begin());
     amrex::Gpu::copy(amrex::Gpu::hostToDevice,
         h_laser_adk_prefactor.begin(), h_laser_adk_prefactor.end(), m_laser_adk_prefactor.begin());
+    amrex::Gpu::copy(amrex::Gpu::hostToDevice,
+         h_laser_dp_prefactor.begin(), h_laser_dp_prefactor.end(), m_laser_dp_prefactor.begin());
 }