[Pkg-octave-devel] Bug#418158: Bug#418158: Bug#418158: more info
rafael at debian.org
Fri Feb 29 22:46:52 UTC 2008
* Rafael Laboissiere <rafael at debian.org> [2008-02-24 18:23]:
> * Ólafur Jens Sigurðsson <ojsbug at 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
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.
More information about the Pkg-octave-devel