Update src/secondary_thermal.cpp

Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
This commit is contained in:
Olaf Schumann 2022-01-17 08:49:46 +01:00
parent 39e57619df
commit dc2e386a29
2 changed files with 4 additions and 6 deletions

View file

@ -4,6 +4,8 @@
#include "openmc/random_lcg.h"
#include "openmc/search.h"
#include <gsl/gsl-lite.hpp>
#include "xtensor/xview.hpp"
#include <cmath> // for log, exp
@ -38,11 +40,7 @@ void CoherentElasticAE::sample(
const auto& energies {xs_.bragg_edges()};
if (E_in < energies.front()) {
// XS is zero for E < energies[0] -> no scatter
mu = 1.0;
return;
}
Expects(E_in >= energies.front());
const int i = lower_bound_index(energies.begin(), energies.end(), E_in);

View file

@ -1,2 +1,2 @@
k-combined:
1.406055E+00 9.581397E-02
1.406055E+00 9.581396E-02