This is the mail archive of the fortran@gcc.gnu.org mailing list for the GNU Fortran project.

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]
Other format:	[Raw text]

Re: Someone broke complex arithmetic

From: Andrew Pinski <pinskia at physics dot uc dot edu>
To: Steve Kargl <sgk at troutmask dot apl dot washington dot edu>
Cc: gcc at gcc dot gnu dot org, fortran at gcc dot gnu dot org
Date: Mon, 18 Jul 2005 23:02:12 -0400
Subject: Re: Someone broke complex arithmetic
References: <20050719023127.GA88291@troutmask.apl.washington.edu> <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

Follow-Ups:
- Re: Someone broke complex arithmetic
  - From: Steve Kargl

References:
- Someone broke complex arithmetic
  - From: Steve Kargl
- Re: Someone broke complex arithmetic
  - From: Steve Kargl

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]