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

changed particle init index to match history based method (1 vs. 0-based). Also added some print statements etc.

Browse source
This commit is contained in:
John Tramm 2019-11-26 20:47:36 +00:00
parent ff27ed2bbe
commit 5379b1c291
2 changed files with 5 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

5
src/simulation.cpp
View file

@ -421,6 +421,7 @@ void process_collision_events()
// Reset banked weight during collision
p->n_bank_ = 0;
p->n_bank_second_ = 0;
p->wgt_bank_ = 0.0;
for (int& v : p->n_delayed_bank_) v = 0;
@ -547,7 +548,7 @@ void transport()
// Initialize all histories
// TODO: Parallelize
for (int i = 0; i < n_particles; i++) {
initialize_history(particles + i, source_offset + i);
initialize_history(particles + i, source_offset + i + 1);
}
//std::cout << "Enqueing particles for XS Lookups..." << std::endl;
@ -568,7 +569,7 @@ void transport()
std::cout << "Collisions = " << collision_queue_length << std::endl;
*/
int max = std::max({calculate_fuel_xs_queue_length, calculate_nonfuel_xs_queue_length, advance_particle_queue_length, surface_crossing_queue_length, collision_queue_length});
check_energies();
//check_energies();
if (max == 0) {
break;
} else if (max == calculate_fuel_xs_queue_length) {