normalize k by total_weight instead of n_particles

GuySten · GuySten · commit f228f27c634c · 2026-01-11T12:44:12.000+02:00
diff --git a/src/eigenvalue.cpp b/src/eigenvalue.cpp
@@ -75,9 +75,8 @@ void calculate_generation_keff()
 #endif
 
   // Normalize single batch estimate of k
-  // TODO: This should be normalized by total_weight, not by n_particles
   if (settings::solver_type != SolverType::RANDOM_RAY) {
-    keff_reduced /= settings::n_particles;
+    keff_reduced /= settings::total_weight;
   }
 
   simulation::k_generation.push_back(keff_reduced);

Original file line number	Diff line number	Diff line change
`@@ -75,9 +75,8 @@ void calculate_generation_keff()`
`75`	`75`	`#endif`
`76`	`76`
`77`	`77`	`// Normalize single batch estimate of k`
`78`		`- // TODO: This should be normalized by total_weight, not by n_particles`
`79`	`78`	`if (settings::solver_type != SolverType::RANDOM_RAY) {`
`80`		`- keff_reduced /= settings::n_particles;`
	`79`	`+ keff_reduced /= settings::total_weight;`
`81`	`80`	`}`
`82`	`81`
`83`	`82`	`simulation::k_generation.push_back(keff_reduced);`