bugGNU Scientific Library - Bugs: bug #53451, gsl_sf_ellint_Pcomp( k, n, mode )...


You are not allowed to post comments on this tracker with your current authentication level.

bug #53451: gsl_sf_ellint_Pcomp( k, n, mode ) returns NaN if mode < -1

Submitter:  Dennis Amelunxen <damelunx>
Submitted:  Sat 24 Mar 2018 10:00:17 AM UTC
Category:  None Severity:  3 - Normal
Operating System:  Linux Status:  None
Assigned to:  None Open/Closed:  Open
Release:  2.1

Sat 24 Mar 2018 10:00:17 AM UTC, original submission:  

Hi there,

I am using

gsl_sf_ellint_Pcomp( k, n, mode )

from the GSL, and only get NaNs when using mode<-1 and k>1.3.
But the integral is defined, compare for example with the result from WolframAlpha:


gives 1.43517..., which is also the output of


But while


gives 4.72112..., the output of


is NaN.

I'm using GSL in R (via the wrapper package https://cran.r-project.org/web/packages/gsl/index.html), so I don't see any error messages.


Dennis Amelunxen <damelunx>


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by damelunx (Submitted the item)
  • -email is unavailable- added by damelunx

    Follows 1 latest change.

    Date Changed by Updated Field Previous Value => Replaced by
    2018-03-24 damelunx Carbon-Copy- Added -email is unavailable-

    Back to the top

    Powered by Savane 3.13-57b4.
    Corresponding source code