Skip to content

Commit

Permalink
fix issues of tab
Browse files Browse the repository at this point in the history
  • Loading branch information
Eya Dammak committed Jan 3, 2025
1 parent 425082b commit 56b4b2d
Show file tree
Hide file tree
Showing 2 changed files with 12 additions and 12 deletions.
2 changes: 1 addition & 1 deletion CMakeLists.txt
Original file line number Diff line number Diff line change
Expand Up @@ -342,7 +342,7 @@ if(BUILD_TESTING)
WORKING_DIRECTORY ${CMAKE_RUNTIME_OUTPUT_DIRECTORY}
)

add_test(NAME laser_ionization.1Rank
add_test(NAME laser_ionization.1Rank
COMMAND bash ${HiPACE_SOURCE_DIR}/tests/laser_ionization.1Rank.sh
$<TARGET_FILE:HiPACE> ${HiPACE_SOURCE_DIR}
WORKING_DIRECTORY ${CMAKE_RUNTIME_OUTPUT_DIRECTORY}
Expand Down
22 changes: 11 additions & 11 deletions src/particles/plasma/PlasmaParticleContainer.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -501,10 +501,10 @@ LaserIonization (const int islice,
(PhysConstSI::ep0 * PhysConstSI::m_e) );
const amrex::Real E0 = Hipace::m_normalized_units ?
wp * PhysConstSI::m_e * PhysConstSI::c / PhysConstSI::q_e : 1;
const amrex::Real lambda0 = laser.GetLambda0();
const amrex::Real lambda0 = laser.GetLambda0();
const amrex::Real omega0 = 2.0 * MathConst::pi * phys_const.c / lambda0;
const bool linear_polarization = laser.LinearPolarization();
const bool linear_polarization = laser.LinearPolarization();
int * const ion_lev = soa_ion.GetIntData(PlasmaIdx::ion_lev).data();
const amrex::Real * const x_prev = soa_ion.GetRealData(PlasmaIdx::x_prev).data();
const amrex::Real * const y_prev = soa_ion.GetRealData(PlasmaIdx::y_prev).data();
Expand Down Expand Up @@ -559,8 +559,8 @@ LaserIonization (const int islice,

amrex::Real Ep = std::sqrt( amrex::abs(Et*Et) + amrex::abs(El*El) );
Ep *= phys_const.m_e * phys_const.c / phys_const.q_e;
Ep *= E0;
Ep *= E0;
// Compute probability of ionization p
const amrex::Real gammap = (1.0_rt + uxp[ip] * uxp[ip] * clightsq
+ uyp[ip] * uyp[ip] * clightsq
Expand All @@ -571,12 +571,12 @@ LaserIonization (const int islice,
std::pow(Ep, adk_power[ion_lev_loc]) *
std::exp( adk_exp_prefactor[ion_lev_loc]/Ep );

amrex::Real w_dtau_ac;
if (linear_polarization) {
w_dtau_ac = w_dtau_dc * std::sqrt(Ep * laser_adk_prefactor[ion_lev_loc]);
} else {
w_dtau_ac = w_dtau_dc;
}
amrex::Real w_dtau_ac;
if (linear_polarization) {
w_dtau_ac = w_dtau_dc * std::sqrt(Ep * laser_adk_prefactor[ion_lev_loc]);
} else {
w_dtau_ac = w_dtau_dc;
}

amrex::Real p = 1._rt - std::exp( - w_dtau_ac );

Expand Down

0 comments on commit 56b4b2d

Please sign in to comment.