boostorg / boostorg/math

arm64_v8a: ibeta_inv exceeds 200 iterations for some arguments

未关闭
#442 8 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看
主要语言
C++
星标
374
派生
264
平均合并
1 天 23 小时
30 天内合并 PR
17

描述

Hi, when compiling for 64-bit ARM (arm64_v8a), ibeta_inv(alpha, beta, x) produces

boost::exception_detail::clone_impl >: Error in function boost::math::ibeta(long double, long double, long double): Root finding evaluation exceeded 200 iterations, giving up now.

when either alpha or beta are smaller than about 3e-3. For other target processors I've tried (arm_v7a, x86, x86_64) this works fine.

I've now hacked up my own custom root finder for these cases which converges in just a few iterations. My code's certainly not up to boost standards, but I'm happy to share it as a starting point if someone wants to use it to improve the ibeta_inv root finder in this regard.

Though I suspect that something else is going on here - the bug may be indicative of some deeper numerical issues on arm64_v8a?

贡献指南

这个仓库没有索引到贡献指南

评估

这个 Issue 还没有评估数据。

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。