At Thu, 02 Jun 2011 17:57:45 -0600, Curtis Hash wrote: > > Reply-To: [email protected] > From: Curtis Hash <[email protected]> > To: [email protected] > Subject: [Bug-gsl] alternative implementation of > gsl_ran_negative_binomial_pdf > Date: Thu, 02 Jun 2011 17:57:45 -0600 > > [1 <text/plain; ISO-8859-1 (7bit)>] > We believe there may be a more numerically stable way to write > gsl_ran_negative_binomial_pdf(). > > Test code:
Thanks for the bug report! I've applied your patch, much appreciated. One minor change - I rewrote log(1-p) as log1p(-p) to reduce cancellation error for p ~ 1. -- Brian Gough _______________________________________________ Bug-gsl mailing list [email protected] https://lists.gnu.org/mailman/listinfo/bug-gsl
