boostorg / boostorg/math

arm64_v8a: ibeta_inv exceeds 200 iterations for some arguments

Đang mở
#442 8 bình luận 0 reaction 0 người được giao Xem trên GitHub
Ngôn ngữ chính
C++
Star
374
Fork
264
Merge trung bình
1 ngày 23 giờ
Pull request đã merge (30 ngày)
17

Mô tả

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?

Hướng dẫn đóng góp

Chưa lập chỉ mục được hướng dẫn đóng góp cho kho mã nguồn này

Đánh giá

Issue này chưa được đánh giá.

Nhận issue mới trong hộp thư của bạn

Bản tóm tắt ngắn những issue GitHub phù hợp với người mới.