[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