crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Make sure data gets initialized for particle restarts

Browse source
This commit is contained in:
Paul Romano 2020-06-17 22:22:01 -05:00
parent 2e29ffacfe
commit c8b99514b7
3 changed files with 16 additions and 13 deletions
Download patch file Download diff file Expand all files Collapse all files

22
src/simulation.cpp
View file

@ -69,19 +69,11 @@ int openmc_simulation_init()
// Skip if simulation has already been initialized
if (simulation::initialized) return 0;
// Determine limits for incident neutron/photon data
// Initialize nuclear data (energy limits, log grid)
if (settings::run_CE) {
determine_energy_limits();
initialize_data();
}
// Set up logarithmic grid for nuclides
for (auto& nuc : data::nuclides) {
nuc->init_grid();
}
int neutron = static_cast<int>(Particle::Type::neutron);
simulation::log_spacing = std::log(data::energy_max[neutron] /
data::energy_min[neutron]) / settings::n_log_bins;
// Determine how much work each process should do
calculate_work();
@ -569,7 +561,7 @@ void calculate_work()
}
}
void determine_energy_limits()
void initialize_data()
{
// Determine minimum/maximum energy for incident neutron/photon data
data::energy_max = {INFTY, INFTY};
@ -629,6 +621,14 @@ void determine_energy_limits()
}
}
}
// Set up logarithmic grid for nuclides
for (auto& nuc : data::nuclides) {
nuc->init_grid();
}
int neutron = static_cast<int>(Particle::Type::neutron);
simulation::log_spacing = std::log(data::energy_max[neutron] /
data::energy_min[neutron]) / settings::n_log_bins;
}
#ifdef OPENMC_MPI