bugGNU Scientific Library - Bugs: bug #47027, gsl_sf_coulomb_wave_FG_e returns...


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

bug #47027: gsl_sf_coulomb_wave_FG_e returns NaN but with success flag

Submitted by:  Patrick Alken <psa>
Submitted on:  Sun 31 Jan 2016 12:49:29 AM UTC  
Category: NoneSeverity: 3 - Normal
Operating System: Status: None
Assigned to: NoneOpen/Closed: Open

Sun 31 Jan 2016 12:49:29 AM UTC, original submission:

from jenssss =at= phys =dot= au =dot= dk

For certain parameters the gsl_sf_coulomb_wave_FG_e returns a NaN value
for Gp (the derivative of the irregular Coulomb function), but the
function does not throw an error.

The problem seems to be similar to bug #30885.

I used gsl 1.16, from the repositories of Ubuntu 15.10.

The following example program

#include <stdio.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_sf_coulomb.h>
main (void)
gsl_sf_result F,G,Fp,Gp;
double exp_F,exp_G;
double eta=340;
double L=0;
double x=48.524525790349422;

int iret = gsl_sf_coulomb_wave_FG_e

printf ("status = %s\n", gsl_strerror(iret));
printf ("F = %.18e\n"
" ± % .18e\n",
F.val, F.err);
printf ("Fp = %.18e\n"
" ± % .18e\n",
Fp.val, Fp.err);
printf ("G = %.18e\n"
" ± % .18e\n",
G.val, G.err);
printf ("Gp = %.18e\n"
" ± % .18e\n",
Gp.val, Gp.err);
return iret;

Gives the following output:

status = success
F = 2.593734965617792327e-309
± 2.593734965617787386e-312
Fp = 9.371127281856027958e-309
± 9.371127281859619815e-312
G = 5.343744177206742137e+307
± 5.343744177206741934e+304
Gp = -nan
± inf

Patrick Alken <psa>
Project Administrator


No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -unavailable- added by psa (Submitted the item)

    Please enter the title of George Orwell's famous dystopian book (it's a date):



    No Changes Have Been Made to This Item

    Back to the top

    Powered by Savane 3.1-cleanup1