[Pkg-octave-devel] Bug#418158: Bug#418158: Bug#418158: Bug#418158: more info
John W. Eaton
jwe at bevo.che.wisc.edu
Mon Mar 3 18:14:59 UTC 2008
On 29-Feb-2008, Rafael Laboissiere wrote:
| * 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
| 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.
Yes, but maybe not until after 3.1 is released unless someone shows up
with a patch.
jwe
More information about the Pkg-octave-devel
mailing list