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

Ensuring Particle::cell_instance_ is set to zero if not a distribcell.

Browse source
This commit is contained in:
Patrick Shriwise 2021-08-02 17:38:27 -05:00
parent 822858dfbb
commit 404dc1ecc9
1 changed files with 2 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

2
src/geometry.cpp
View file

@ -77,6 +77,7 @@ int cell_instance_at_level(const Particle& p, int level) {
// quick exit if this cell doesn't have distribcell instances
if (c.distribcell_index_ == C_NONE) return C_NONE;
// compute the cell's instance
int instance = 0;
for (int i = 0; i < level; i++) {
const auto& c_i {*model::cells[p.coord(i).cell]};
@ -159,6 +160,7 @@ find_cell_inner(Particle& p, const NeighborList* neighbor_list)
if (c.type_ == Fill::MATERIAL) {
// Found a material cell which means this is the lowest coord level.
p.cell_instance() = 0;
// Find the distribcell instance number.
if (c.distribcell_index_ >= 0) {
p.cell_instance() = cell_instance_at_level(p, p.n_coord() - 1);