prog

src/math.ml

@@ -90,7 +90,7 @@ let rec bisect_increasing ?(factor = 1e-10) ?(strict = true) ?(relax = 0.0) ?low
         (* go left *)
         (*Format.printf "left\n";*)
         bisect_increasing ~relax ~strict ~factor ?low
-          (bmin -. (relax *. width), m +. (relax *. width))
+          (max 0.0 (bmin -. (relax *. width)), m +. (relax *. width))
           f_to_evaluate
       else
         (* go right *)