Re: [Bug-gsl] GSL: suggestion, E_n

[Oleg V. Motygin]

> At Wed, 08 Oct 2008 14:21:55 +0100, Oleg V. Motygin wrote:
> > I think it would be great if you could improve the GSL routines for
> >  computations of exponential integrals so that they can accept 
> > complex arguments too. A good algorithm for this was published in
> >
> > ACM Transactions on Mathematical Software Donald E. Amos, 
> > Computation of Exponential Integrals of a Complex Argument, vol.16,
> >  no. 2, p.169--177, 1990, http://doi.acm.org/10.1145/78928.78933
> >
> > ACM Transactions on Mathematical Software Donald E. Amos, Algorithm
> >  683: A Portable FORTRAN Subroutine for Exponential Integrals of a
> >  Complex Argument, vol.16, no. 2, p.178--182, 1990, 
> > http://doi.acm.org/10.1145/78928.78934
>
> Thanks for the suggestion - I agree it could be a useful function to 
> have. Would you be able to work on this?

Dear Brian,

Thank you very much for the reply. I guess that I do not have sufficient
qualification for this work (especially in multi-architecture libraries
and C/C++ programming). Anyway, I would be happy to participate in some
way. I have used this algorithm in my numerical investigations with
Octave. For this I slightly adapted the fortran code supplying it with
machine constants for IEEE (IBM.PC) arithmetics. Besides, I wrote a C++
interface for the fortran code to use it with Octave and tested the
function against Mathematica's one. So, if it is of interest I can
send you the stuff with comments.

With best wishes,
Oleg