check if nuclide is fissionable in eigenvalue tally

ilhamv · Nov 10, 2024 · c3c3e3b · c3c3e3b
1 parent bd510fe
commit c3c3e3b
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/mcdc/kernel.py b/mcdc/kernel.py
@@ -2126,6 +2126,8 @@ def eigenvalue_tally(P_arr, distance, mcdc):
             for i in range(material["N_nuclide"]):
                 ID_nuclide = material["nuclide_IDs"][i]
                 nuclide = mcdc["nuclides"][ID_nuclide]
+                if not nuclide["fissionable"]:
+                    continue
                 for j in range(J):
                     nu_d = get_nu_group(NU_FISSION_DELAYED, nuclide, E, j)
                     decay = nuclide["ce_decay"][j]