Bug#121852: marked as done (g77-3.0 fails to compute a complex log of certain values)
Your message dated 07 Mar 2002 16:00:43 +0100
with message-id <yrywuwoh17o.fsf@llc2.psi.ch>
and subject line
has caused the attached Bug report to be marked as done.
This means that you claim that the problem has been dealt with.
If this is not the case it is now your responsibility to reopen the
Bug report if necessary, and/or fix the problem forthwith.
(NB: If you are a system administrator and have no idea what I am
talking about this indicates a serious mail system misconfiguration
somewhere. Please contact me immediately.)
Debian bug tracking system administrator
(administrator, Debian Bugs database)
--------------------------------------
Received: (at submit) by bugs.debian.org; 30 Nov 2001 15:44:34 +0000
>From guasch@itp.uni-karlsruhe.de Fri Nov 30 09:44:34 2001
Return-path: <guasch@itp.uni-karlsruhe.de>
Received: from (localhost) [129.13.104.206] (mail)
by master.debian.org with esmtp (Exim 3.12 1 (Debian))
id 169pq6-0002vJ-00; Fri, 30 Nov 2001 09:44:34 -0600
Received: from guasch by localhost with local (Exim 3.32 #1 (Debian))
id 169puX-000628-00; Fri, 30 Nov 2001 16:49:09 +0100
From: Jaume Guasch <guasch@itp.uni-karlsruhe.de>
To: Debian Bug Tracking System <submit@bugs.debian.org>
Subject: g77-3.0 fails to compute a complex log of certain values
X-Reportbug-Version: 1.36
X-Mailer: reportbug 1.36
Date: Fri, 30 Nov 2001 16:49:08 +0100
Message-Id: <E169puX-000628-00@localhost>
Sender: Jaume Guasch <guasch@jasi.physik.uni-karlsruhe.de>
X-BadReturnPath: guasch@jasi.physik.uni-karlsruhe.de rewritten as guasch@itp.uni-karlsruhe.de
using "From" header
Delivered-To: submit@bugs.debian.org
Package: g77-3.0
Version: 1:3.0.2-3
Severity: normal
g77-3.0 fails to compute a complex log of certain values
computing log(x), for certain complex values of x (e.g. x=(1.00390625,0.) )
the fortran program enters a loop (no error, no core ... ), does not
compute the value, and the program has to be killed (with Ctrl-C).
g77-2.95 can compute the value.
Below there is a script session (together with the program) to reproduce
the bug. g77-2.95 runs the program fine, g77-3.0 doesn't
Best regards,
Jaume
==================================================
Script started on Fri Nov 30 16:36:03 2001
jasi $ g77 --version
GNU Fortran 0.5.25 20010319 (prerelease)
Copyright (C) 1997 Free Software Foundation, Inc.
For more version information on components of the GNU Fortran
compilation system, especially useful when reporting bugs,
type the command `g77 --verbose'.
GNU Fortran comes with NO WARRANTY, to the extent permitted by law.
You may redistribute copies of GNU Fortran
under the terms of the GNU General Public License.
For more information about these matters, see the file named COPYING
or type the command `info -f g77 Copying'.
jasi $ gcc --version
2.95.4
jasi $ cat prova.F
program prova
implicit none
double complex cs,x1,x2
double precision precc,x0
x0=0.d0
precc = 0.01d0
print *,"precc,x0:",precc,x0
C This is fine
cs = exp(log(DCMPLX(1+precc,x0)))
write(6,FMT=*), "cs=exp(log(1+precc))=",dble(cs)
x0=0.d0
precc = 0.001d0
print *,"precc,x0:",precc,x0
C This is fine
cs = exp(log(DCMPLX(1+precc,x0)))
write(6,*), "cs=exp(log(1+precc))=",dble(cs)
C This value is not fine
precc=0.00390625d0
print *,"precc,x0:",precc,x0
x1=DCMPLX(1+precc,x0)
write(6,*) "x1=",x1
x2=log(x1)
write(6,*) "x2=log(x1)=",x2
cs=exp(x2)
write(6,*), "cs=exp(x2)=",dble(cs)
cs = exp(log(DCMPLX(1+precc,x0)))
write(6,*), "cs=exp(log(1+precc))=",dble(cs)
stop
end
jasi $ g77 -o prova295 prova.F
jasi $ ./prova295
precc,x0: 0.01 0.
cs=exp(log(1+precc))= 1.01
precc,x0: 0.001 0.
cs=exp(log(1+precc))= 1.001
precc,x0: 0.00390625 0.
x1= (1.00390625,0.)
x2=log(x1)= (0.00389864042,0.)
cs=exp(x2)= 1.00390625
cs=exp(log(1+precc))= 1.00390625
jasi $ g77-3.0 --version
GNU Fortran 0.5.26 20011023 (release)
Copyright (C) 2001 Free Software Foundation, Inc.
For more version information on components of the GNU Fortran
compilation system, especially useful when reporting bugs,
type the command `g77 --verbose'.
GNU Fortran comes with NO WARRANTY, to the extent permitted by law.
You may redistribute copies of GNU Fortran
under the terms of the GNU General Public License.
For more information about these matters, see the file named COPYING
or type the command `info -f g77 Copying'.
jasi $ gcc-2 3, .0 --version
3.0.2
jasi $ g77-3.0 -o prova3 prova.F
jasi $ ./prova3
precc,x0: 0.01 0.
cs=exp(log(1+precc))= 1.01
precc,x0: 0.001 0.
cs=exp(log(1+precc))= 1.001
precc,x0: 0.00390625 0.
x1= (1.00390625,0.)
Interrupt
jasi $ exit
Script done on Fri Nov 30 16:37:32 2001
-- System Information
Debian Release: testing/unstable
Architecture: i386
Kernel: Linux jasi 2.4.14 #1 Sun Nov 18 18:19:45 CET 2001 i686
Locale: LANG=C, LC_CTYPE=
Versions of packages g77-3.0 depends on:
ii gcc-3.0 1:3.0.2-3 The GNU C compiler.
ii gcc-3.0-base 1:3.0.2-3 The GNU Compiler Collection (base
ii libc6 2.2.4-5 GNU C Library: Shared libraries an
---------------------------------------
Received: (at 121852-done) by bugs.debian.org; 7 Mar 2002 15:00:47 +0000
>From jaume.guasch@psi.ch Thu Mar 07 09:00:47 2002
Return-path: <jaume.guasch@psi.ch>
Received: from llc2.psi.ch [129.129.16.72]
by master.debian.org with esmtp (Exim 3.12 1 (Debian))
id 16izNv-0005j2-00; Thu, 07 Mar 2002 09:00:47 -0600
Received: (from guasch@localhost)
by llc2.psi.ch (8.9.3/8.9.3) id QAA09001;
Thu, 7 Mar 2002 16:00:45 +0100
X-Authentication-Warning: llc2.psi.ch: guasch set sender to jaume.guasch@psi.ch using -f
Sender: guasch@llc2.psi.ch
To: 121852-done@bugs.debian.org
Subject:
References: <E169puX-000628-00@localhost>
<ib97kp0cnyr.fsf@itppc15.physik.uni-karlsruhe.de>
From: Jaume Guasch <jaume.guasch@psi.ch>
X-No-Archive: Yes
Date: 07 Mar 2002 16:00:43 +0100
In-Reply-To: <ib97kp0cnyr.fsf@itppc15.physik.uni-karlsruhe.de>
Message-ID: <yrywuwoh17o.fsf@llc2.psi.ch>
Lines: 0
User-Agent: Gnus/5.0808 (Gnus v5.8.8) Emacs/20.7
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Delivered-To: 121852-done@bugs.debian.org
Reply to: