Fix loop over particle coordinate levels for source domain rejection (#…

…2751)
openmc-dev · Nov 2, 2023 · 7aaccbb · 7aaccbb
1 parent 19f121a
commit 7aaccbb
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/source.cpp b/src/source.cpp
@@ -222,10 +222,10 @@ SourceSite IndependentSource::sample(uint64_t* seed) const
             found = contains(domain_ids_, model::materials[mat_index]->id());
           }
         } else {
-          for (const auto& coord : p.coord()) {
+          for (int i = 0; i < p.n_coord(); i++) {
             auto id = (domain_type_ == DomainType::CELL)
-                        ? model::cells[coord.cell]->id_
-                        : model::universes[coord.universe]->id_;
+                        ? model::cells[p.coord(i).cell]->id_
+                        : model::universes[p.coord(i).universe]->id_;
             if ((found = contains(domain_ids_, id)))
               break;
           }