remove flux-outter, not necessary

mcdc/kernel.py

@@ -2342,11 +2342,6 @@ def iqmc_prepare_source(mcdc):
     """
     flux_scatter = mcdc["technique"]["iqmc"]["score"]["flux"]
     flux_fission = mcdc["technique"]["iqmc"]["score"]["flux"]
-    if (
-        mcdc["setting"]["mode_eigenvalue"]
-        and mcdc["technique"]["iqmc"]["eigenmode_solver"] == "power_iteration"
-    ):
-        flux_fission = mcdc["technique"]["iqmc"]["score"]["flux-outter"]
     mesh = mcdc["technique"]["iqmc"]["mesh"]
     Nt = len(mesh["t"]) - 1
     Nx = len(mesh["x"]) - 1

mcdc/loop.py

@@ -529,7 +529,6 @@ def power_iteration(mcdc):
     # maximum number of iterations
     maxit = mcdc["technique"]["iqmc"]["maxitt"]
     score_bin = mcdc["technique"]["iqmc"]["score"]
-    score_bin["flux-outter"] = score_bin["flux"].copy()
     k_old = mcdc["k_eff"]
     solver = mcdc["technique"]["iqmc"]["fixed_source_solver"]
 
@@ -562,7 +561,6 @@ def power_iteration(mcdc):
         mcdc["technique"]["iqmc"]["res_outter"] = abs(mcdc["k_eff"] - k_old) / k_old
         k_old = mcdc["k_eff"]
         # store outter iteration values
-        score_bin["flux-outter"] = score_bin["flux"].copy()
         score_bin["effective-fission-outter"] = score_bin["effective-fission"].copy()
         fission_source_old = score_bin["fission-source"].copy()
         mcdc["technique"]["iqmc"]["itt_outter"] += 1

mcdc/type_.py

@@ -558,7 +558,6 @@ def make_type_technique(N_particle, G, card):
         scores_struct += [(name, float64, shape)]
     # TODO: make outter flux/effective fission size zero if not eigenmode
     # (causes problems with numba)
-    scores_struct += [("flux-outter", float64, (Ng, Nt, Nx, Ny, Nz))]
     scores_struct += [("effective-fission-outter", float64, (Ng, Nt, Nx, Ny, Nz))]
     scores = np.dtype(scores_struct)
     iqmc_list += [("score", scores)]