From 7aaccbbeee40e3cd3766de492a7ed4644c272dce Mon Sep 17 00:00:00 2001
From: Paul Romano <paul.k.romano@gmail.com>
Date: Wed, 1 Nov 2023 19:56:28 -0500
Subject: [PATCH] Fix loop over particle coordinate levels for source domain
 rejection (#2751)

---
 src/source.cpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/source.cpp b/src/source.cpp
index 3251ca0d0d..3ddd903bf4 100644
--- 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;
           }