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:

https://www.wolframalpha.com/input/?i=EllipticPi(0.1,1.3,0.25)

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

gsl_sf_ellint_Pcomp(1.3,0.5,-0.1)

But while

https://www.wolframalpha.com/input/?i=EllipticPi(1.1,1.3,0.25)

gives 4.72112..., the output of

gsl_sf_ellint_Pcomp(1.3,0.5,-1.1)

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.

Best,
Dennis

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-4448.
    Corresponding source code