[Pkg-octave-devel] Bug#418158: Bug#418158: Bug#418158: more info
* Rafael Laboissiere <rafael@debian.org> [2008-02-24 18:23]:
> * Ólafur Jens Sigurðsson <ojsbug@gmail.com> [2008-02-24 12:34]:
>
> > I just tried some values and the limit seems to be poisscdf(604,604),
> > after that the function fails, that is poisscdf(605,605) fails and
> > poisscdf(604,604) succeeds.
>
> It fails for smaller values than that, like poisscdf (599, 600). The
> culprit is the gammainc function, which ultimately calls the function
> defined in d9lgit.f.
Just for the record, the gamma_inc_P function in the gsl octave-forge
package works fine in the range for which gammainc from d9lgit.f fails. For
instance:
octave:1> gamma_inc_P (1e5,1e5)
ans = 0.50042
John, is there any chance of integrating GSL directly in Octave? At least
making the compilation in files like gammainc.cc conditional on the presence
of libgsl and, if not, falling back to Slatec.
--
Rafael
Reply to: