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

Precompute offsets in random ray source update (#3775)

Browse source
This commit is contained in:
Jonathan Shimwell 2026-02-06 16:25:14 +01:00 • committed by GitHub
parent f2c936cf5b
commit 6f625f3dea
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
1 changed files with 7 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/random_ray/flat_source_domain.cpp
View file

@ -112,19 +112,19 @@ void FlatSourceDomain::update_single_neutron_source(SourceRegionHandle& srh)
double density_mult = srh.density_mult();
if (material != MATERIAL_VOID) {
double inverse_k_eff = 1.0 / k_eff_;
int material_offset = material * negroups_;
int scatter_offset = material * negroups_ * negroups_;
for (int g_out = 0; g_out < negroups_; g_out++) {
double sigma_t = sigma_t_[material * negroups_ + g_out] * density_mult;
double sigma_t = sigma_t_[material_offset + g_out] * density_mult;
double scatter_source = 0.0;
double fission_source = 0.0;
for (int g_in = 0; g_in < negroups_; g_in++) {
double scalar_flux = srh.scalar_flux_old(g_in);
double sigma_s = sigma_s_[material * negroups_ * negroups_ +
g_out * negroups_ + g_in] *
density_mult;
double nu_sigma_f =
nu_sigma_f_[material * negroups_ + g_in] * density_mult;
double chi = chi_[material * negroups_ + g_out];
double sigma_s =
sigma_s_[scatter_offset + g_out * negroups_ + g_in] * density_mult;
double nu_sigma_f = nu_sigma_f_[material_offset + g_in] * density_mult;
double chi = chi_[material_offset + g_out];
scatter_source += sigma_s * scalar_flux;
if (settings::create_fission_neutrons) {