search

JuliaIntervals / IntervalRootFinding.jl

Library for finding the roots of a function using interval arithmetic
https://juliaintervals.github.io/IntervalRootFinding.jl/
Other
129 stars 26 forks source link

Support for Bessel functions (or SpecialFunctions.jl in general) #175

Open merlinpelz opened 3 years ago

merlinpelz commented 3 years ago

It seems there is no ::Interval support for besselk and besseli of SpecialFunctions.jl since IntervalRootFinding.roots() gives a StackOverflowError when applied to (a combination of) these functions. This issue may be relevant for both IntervalArithmetic and IntervalRootFinding.

dpsanders commented 3 years ago

Yes, unfortunately implementing interval versions of special functions is Non-trivial.

However, many such functions are implemented in arb, so I think we should just wrap those.

merlinpelz commented 3 years ago

That would be highly appreciated - since the performance of IntervalRootFinding is so strong for "normal" functions and Bessel functions appear often in boundary value problems.

wsshin commented 2 years ago

I am suffering from this exact issue as well. It will be very nice to be able to find the roots of special functions by IntervalRootFinding.jl.