Change indexing for prn sampling for URR ptables

This commit is contained in:
Paul Romano 2021-03-24 07:56:01 -05:00
parent de2cbe0c1e
commit f8241f567f

View file

@ -808,9 +808,7 @@ void Nuclide::calculate_urr_xs(int i_temp, Particle& p) const
// reuse random numbers for the same nuclide at different temperatures,
// therefore preserving correlation of temperature in probability tables.
p.stream_ = STREAM_URR_PTABLE;
//TODO: to maintain the same random number stream as the Fortran code this
//replaces, the seed is set with index_ + 1 instead of index_
double r = future_prn(static_cast<int64_t>(index_ + 1), *p.current_seed());
double r = future_prn(static_cast<int64_t>(index_), *p.current_seed());
p.stream_ = STREAM_TRACKING;
int i_low = 0;