[math] Remez algorithm in Boost Bessel implementation

znailcode--- via Boost <[email protected]>
Newsgroups gmane.comp.lib.boost.devel
Message-ID <[email protected]>
Hi,

I am currently working on the Struve function implementation alongside the Bessel functions.

To approximate H₀(x) and H₁(x), I plan to use the Remez minimax approximation algorithm, similar to what is used in the current Bessel implementation:
https://github.com/boostorg/math/blob/develop/include/boost/math/special_functions/detail/bessel_j0.hpp (for example, in the region x ≤ 8).

I see that an implementation of the Remez algorithm already exists in Boost here:
https://github.com/boostorg/math/blob/develop/include_private/boost/math/tools/remez.hpp

I would like to ask which implementation of the Remez algorithm was used to derive the Bessel function approximations. Was it this Boost implementation, or a different one? Additionally, how were the reference values computed for the approximation process?

Currently, I am using mpmath in Python to generate reference values for the Struve function with the required precision.

Best regards,
Zaripov Nail
_______________________________________________
Boost mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://lists.boost.org/mailman3/lists/boost.lists.boost.org/
Archived at: https://lists.boost.org/archives/list/[email protected]/message/LF6HHQCG527SUXC7FAE2CSGUF4HQVPX5/
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.