Update update_grid.cc

artis-mcrt · May 1, 2024 · 8f528a9 · 8f528a9
1 parent 27a128e
commit 8f528a9
Showing 1 changed file with 9 additions and 10 deletions.
diff --git a/update_grid.cc b/update_grid.cc
@@ -1251,16 +1251,15 @@ void update_grid(FILE *estimators_file, const int nts, const int nts_prev, const
 #endif
 
     // generate a loop over all mgi
-    // if constexpr (DETAILED_BF_ESTIMATORS_ON) {
-    //  for (int mgi = 0; mgi < grid::get_npts_model(); mgi++) {
-    //    const double deltaV =
-    //        grid::get_modelcell_assocvolume_tmin(mgi) * pow(globals::timesteps[nts_prev].mid / globals::tmin, 3);
-    //
-    //    const double estimator_normfactor = 1 / deltaV / deltat / globals::nprocs;
-    //
-    //    radfield::normalise_bf_estimators(mgi, estimator_normfactor / H);
-    //  }
-    //}
+    if constexpr (DETAILED_BF_ESTIMATORS_ON) {
+      for (int mgi = 0; mgi < grid::get_npts_model(); mgi++) {
+        const double estimator_normfactor =
+            1 / (grid::get_modelcell_assocvolume_tmin(mgi) * pow(globals::timesteps[nts_prev].mid / globals::tmin, 3)) /
+            globals::timesteps[nts_prev].width / globals::nprocs;
+
+        radfield::normalise_bf_estimators(mgi, estimator_normfactor / H);
+      }
+    }
 
     for (int mgi = 0; mgi < grid::get_npts_model(); mgi++) {
       /// Check if this task should work on the current model grid cell.