From: Andrew Pinski <pinskia@physics.uc.edu>
To: Steve Kargl <sgk@troutmask.apl.washington.edu>
Cc: gcc@gcc.gnu.org, fortran@gcc.gnu.org
Subject: Re: Someone broke complex arithmetic
Date: Tue, 19 Jul 2005 03:02:00 -0000 [thread overview]
Message-ID: <8543fda2eb87c81a7ab47c78a5c7b187@physics.uc.edu> (raw)
In-Reply-To: <20050719024516.GA88459@troutmask.apl.washington.edu>
On Jul 18, 2005, at 10:45 PM, Steve Kargl wrote:
> On Mon, Jul 18, 2005 at 07:31:27PM -0700, Steve Kargl wrote:
>>
>> Here's the output from a program that brought the problem
>> to my attention. It uses downward recursion to compute
>> spherical Bessel functions.
>>
>> NAG's F95 compiler
>> n x jn(x) jn(cmplx(x,0))
>> 0 2.2900E+01 -3.4445E-02 -3.4445E-02 0.0000E+00
>> 1 2.2900E+01 2.5337E-02 2.5337E-02 0.0000E+00
>> 2 2.2900E+01 3.7765E-02 3.7765E-02 0.0000E+00
>> 3 2.2900E+01 -1.7091E-02 -1.7091E-02 0.0000E+00
>> 4 2.2900E+01 -4.2989E-02 -4.2989E-02 0.0000E+00
>>
>
> gfortran in 4.0.1 gives the correct results shown above.
This is most likely the same problem as PR 22504.
Could you attach your program to that PR?
-- Pinski
next prev parent reply other threads:[~2005-07-19 3:02 UTC|newest]
Thread overview: 10+ messages / expand[flat|nested] mbox.gz Atom feed top
2005-07-19 2:31 Steve Kargl
2005-07-19 2:45 ` Steve Kargl
2005-07-19 3:02 ` Andrew Pinski [this message]
2005-07-19 3:11 ` Steve Kargl
[not found] <42D7F643.2080202@wanadoo.fr>
[not found] ` <200507160334.j6G3Yad26764@makai.watson.ibm.com>
[not found] ` <42DA8C72.2020207@wanadoo.fr>
[not found] ` <200507172119.j6HLJPd26676@makai.watson.ibm.com>
[not found] ` <42DAD899.2090403@wanadoo.fr>
[not found] ` <200507172332.j6HNWmd27626@makai.watson.ibm.com>
[not found] ` <42DB3313.2070707@wanadoo.fr>
[not found] ` <200507182101.j6IL1fd27934@makai.watson.ibm.com>
[not found] ` <42DC9232.10807@wanadoo.fr>
2005-07-19 7:37 ` FX Coudert
2005-07-19 17:17 ` Daniel Berlin
2005-07-19 17:31 ` Steve Kargl
2005-07-19 18:21 ` Richard Henderson
2005-07-19 18:32 ` Steve Kargl
2005-07-25 10:49 ` Gerald Pfeifer
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=8543fda2eb87c81a7ab47c78a5c7b187@physics.uc.edu \
--to=pinskia@physics.uc.edu \
--cc=fortran@gcc.gnu.org \
--cc=gcc@gcc.gnu.org \
--cc=sgk@troutmask.apl.washington.edu \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).