crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Fix for erroneously non-zero tally results of photon threshold reactions (#3242)

Browse source 
Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
This commit is contained in:
Adam Nelson 2025-01-08 13:58:59 -06:00 • committed by GitHub
parent 4492f9db10
commit 8ba66f9fe3
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
1 changed files with 10 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/photon.cpp
View file

@ -328,15 +328,19 @@ PhotonInteraction::PhotonInteraction(hid_t group)
}
// Take logarithm of energies and cross sections since they are log-log
// interpolated
// interpolated. Note that cross section libraries converted from ACE files
// represent zero as exp(-500) to avoid log-log interpolation errors. For
// values below exp(-499) we store the log as -900, for which exp(-900)
// evaluates to zero.
double limit = std::exp(-499.0);
energy_ = xt::log(energy_);
coherent_ = xt::where(coherent_ > 0.0, xt::log(coherent_), -500.0);
incoherent_ = xt::where(incoherent_ > 0.0, xt::log(incoherent_), -500.0);
coherent_ = xt::where(coherent_ > limit, xt::log(coherent_), -900.0);
incoherent_ = xt::where(incoherent_ > limit, xt::log(incoherent_), -900.0);
photoelectric_total_ = xt::where(
photoelectric_total_ > 0.0, xt::log(photoelectric_total_), -500.0);
photoelectric_total_ > limit, xt::log(photoelectric_total_), -900.0);
pair_production_total_ = xt::where(
pair_production_total_ > 0.0, xt::log(pair_production_total_), -500.0);
heating_ = xt::where(heating_ > 0.0, xt::log(heating_), -500.0);
pair_production_total_ > limit, xt::log(pair_production_total_), -900.0);
heating_ = xt::where(heating_ > limit, xt::log(heating_), -900.0);
}
PhotonInteraction::~PhotonInteraction()