Allow for machine-precision differences in determining libration region

BlackHolePerturbationToolkit · Sep 15, 2023 · ef219d6 · ef219d6
1 parent 3906498
commit ef219d6
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/Kernel/TeukolskyMode.m b/Kernel/TeukolskyMode.m
@@ -117,6 +117,7 @@
   If[e != 0,
     rmin = p/(1+e);
     rmax = p/(1-e);
+    Module[{eps=2/10^Precision[{p,e}]}, rmin = (1-eps)rmin; rmax = (1+eps)rmax];
     R = TeukolskyRadial[s, l, m, a, \[Omega], Method->{"NumericalIntegration","Domain"-> {"In"->{rmin,rmax}, "Up"->{rmin,rmax}}},
         "Amplitudes" -> <|"In"-> R["In"]["UnscaledAmplitudes"], "Up"-> R["Up"]["UnscaledAmplitudes"]|>,
         "RenormalizedAngularMomentum"-> R["In"]["RenormalizedAngularMomentum"], "Eigenvalue" -> R["In"]["Eigenvalue"]];