|
From: | Oleg V. Motygin |
Subject: | Re: [Bug-gsl] GSL: suggestion, E_n |
Date: | Sun, 12 Oct 2008 13:48:32 +0100 |
User-agent: | Thunderbird 2.0.0.6 (X11/20070728) |
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 forcomputations of exponential integrals so that they can accept complex arguments too. A good algorithm for this was published inACM 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 aComplex Argument, vol.16, no. 2, p.178--182, 1990, http://doi.acm.org/10.1145/78928.78934Thanks 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
[Prev in Thread] | Current Thread | [Next in Thread] |