In the current version, LinA loops when trying to run the exact solver with low precisions. To reproduce, consider the following piece of code
julia> f(x) = 3.0374999999999974e-7 * x ^ 3 + -0.0016874999999999985 * x ^ 2 + 3.0937499999999982 * x ^ 1 + 3750.0
julia> LinA.Linearize(f, 0.0, 5000.0, LinA.Relative(1e+1), LinA.ExactLin(); bounding = LinA.Under())
I am proposing a workaround in a separate branch.