Calling boost::math::ibeta_inv(1e50, 10, 0.1) produces:
libc++abi.dylib: terminating with uncaught exception of type boost::wrapexcept<std::domain_error>: Error in function boost::math::ibeta<long double>(long double, long double, long double): Parameter x outside the range [0,1] in the incomplete beta function (got x=835.800004005432128906).
I suspect the inverse iteration for ibeta is thrown off by the large alpha.
Calling
boost::math::ibeta_inv(1e50, 10, 0.1)produces:libc++abi.dylib: terminating with uncaught exception of type boost::wrapexcept<std::domain_error>: Error in function boost::math::ibeta<long double>(long double, long double, long double): Parameter x outside the range [0,1] in the incomplete beta function (got x=835.800004005432128906).I suspect the inverse iteration for ibeta is thrown off by the large alpha.