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
Priority:  5 - Normal Operating System:  Linux
Status:  None Assigned to:  None
Open/Closed:  Open Release:  2.1

Discussion

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.wolfr ... icPi(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.wolfr ... icPi(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-pr ... es/gsl/index.html), so I don't see any error messages.

Best,
Dennis

Dennis Amelunxen <damelunx>

 

Attached Files

This item currently has no attached files.

 

Dependencies

This item does not depend on any other items.

No items depend on this one.

 

Mail Notification Carbon-Copy List

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

    History

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