On 07/ 7/10 01:46 PM, Chao Ding wrote:
If GSL_IEEE does NOT work in Windows, the initiation of the underflow
in my problem is that, I am computing a beta-binomial probability,
which is
gsl_sf_choose(n, k) * gsl_sf_beta(s+ k, f + n-k) /gsl_sf_beta(s, f);
When s, f get larger, the problem of underflow happens. Is there any
way to get rid of the underflow problem without GSL_IEEE_MODE setting?
Thanks,
Chao
Does re-writing in terms of gsl_sf_lnbeta help?
Rodney
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl