[SciPy-user] exponential integral behaviour near -20

David Cournapeau david@ar.media.kyoto-u.ac...
Mon Dec 15 23:43:50 CST 2008


David Cournapeau wrote:
> Robert Kern wrote:
>   
>> Hmm. I dunno. Perhaps gfortran's runtime changed a branch cut for
>> CDLOG() between 4.2 and 4.3. That would explain the +/- pi*I in the
>> output.
>>     
>
> Actually, on my machine, the first goes through the code path with
> CDLOG, and the other one with the code path using CDEXP. I used a simple
> program fortran program:

Grr, thunderbird screwed the indentation, I attached the file instead,

David
-------------- next part --------------
A non-text attachment was scrubbed...
Name: yop.f
Type: text/x-fortran
Size: 1426 bytes
Desc: not available
Url : http://projects.scipy.org/pipermail/scipy-user/attachments/20081216/55b2f0e9/attachment.bin 


More information about the SciPy-user mailing list