[suggestion] for gsl_sf_Ci_e (in specfunc/sinint.c)
Gerard Jungman
jungman@lanl.gov
Wed Dec 19 13:20:00 GMT 2001
Arvid Requate wrote:
>
> while using the special function routine gsl_sf_Ci_e, I noticed
> that given an argument of zero it aborts with an error, while
> mathematically it should simply return zero.
Ci(x) has a logarithmic singularity at x=0. There may
be varying definitions, but the standard one has this
behaviour, as does the implementation in GSL.
Thanks.
--
G. Jungman
More information about the Gsl-discuss
mailing list