大参数值的gsl_cdf_beta_Pinv收敛问题

Question

我正在使用gsl 2.5版来计算Beta分布（documentation）的某些值的反cdf。但是，对于较大的参数值，我似乎遇到收敛问题。例如：

#include <math.h>
#include <stdlib.h>
#include <gsl/gsl_cdf.h>

int main(void){

   gsl_cdf_beta_Pinv(0.999810, 4.0, 499997.0);
   return(0);

   }

返回以下错误：

ERROR: inverse failed to converge

在几个地方提到了这个问题，但是我似乎找不到一致的解决方案。该问题已在2015年here中指出，但目前尚不清楚在最新版本的gsl中是否已解决此问题（我想不是）。我可以尝试按照电子邮件线程中的建议更改源（我不确定应该在源中的确切位置进行更改），但是我很好奇是否存在对此问题的可接受解决方案（或者是否存在其他库） C可以进行相同的计算而不会出现收敛问题。