added external xtensor within PhotonInteraction for all shell photoelectric xs to utilize contiguous memory in calculate_xs

This commit is contained in:
Patrick Myers 2022-05-25 09:06:53 -05:00
parent 1be02f90bd
commit 687d43b0ab
3 changed files with 15 additions and 19 deletions

View file

@ -358,9 +358,9 @@ void sample_photon_reaction(Particle& p)
continue;
// Evaluation subshell photoionization cross section
double xs = std::exp(shell.cross_section(i_grid - i_start) +
f * (shell.cross_section(i_grid + 1 - i_start) -
shell.cross_section(i_grid - i_start)));
double xs = std::exp(element.cross_sections_(i_grid - i_start) +
f * (element.cross_sections_(i_grid + 1 - i_start) -
element.cross_sections_(i_grid - i_start)));
prob += xs;
if (prob > cutoff) {