Improved the bracketing method for SFD slope

MintonGroup · Feb 23, 2023 · 2423034 · 2423034
1 parent 90009c7
commit 2423034
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/fraggle/fraggle_util.f90 b/src/fraggle/fraggle_util.f90
@@ -182,8 +182,8 @@ module subroutine fraggle_util_set_mass_dist(self, param)
                   y1 = y1 - (i-1)**(-x1/3.0_DP)
                end do
                if (y0*y1 < 0.0_DP) exit
-               x1 = x1 * 1.1_DP
-               x0 = x0 * 1.1_DP
+               if (y0 > 0.0_DP) x0 = x0 * 1.6_DP
+               if (y1 < 0.0_DP) x1 = x1 * 1.6_DP
             end do
 
             ! Find the mass scaling factor with bisection