Informe problema libc6 en cálculos matemáticos

[José F Zárate <al276005 en academ01 lag itesm mx>]

   Que tal... 

        Se ha reportado problemas de error de calculo con programas bajo la
libc6 (especmficamente con g77), que no se presentan en la lib5.

   Estos dmas a venido un hilo de discusisn en las listas Debian-spanish.

El problema es como sigue:
*************************************************************
Subject:      g77 logarithm problem
From:         Do-Hoon Kwon <kwon 22 en osu edu>
Date:         1998/01/29
Message-ID:   <34D0D4D5 1B28EC47 en osu edu>
Newsgroups:   comp.os.linux.misc,comp.lang.fortran
[More Headers]


Hello,
 I use g77 with glibc-2.0.6 on my RedHat 5.0 linux box at school.

[dhk@dhkpc tmp]$ g77 --version
GNU Fortran 0.5.21
Copyright (C) 1997 Free Software Foundation, Inc.

I'm having a weird problem calculating logarithm with complex argument.
My test program is to calculate log(1+i*1).

      complex xj
      xj=(0.0,1.0)
      write (*,*) log(1.0+xj*1.0)
      end

[dhk@dhkpc tmp]$ g77 test.f -o test.x
[dhk@dhkpc tmp]$ test.x
 (0.240239993,0.785398185)

Here is the output of MATLAB 5.1 on the same linux box.


[dhk@dhkpc cubekeep]$ matlab

                            < M A T L A B (R) >
                (c) Copyright 1984-97 The MathWorks, Inc.
                            All Rights Reserved
                             Version 5.1.0.421
                                May 25 1997

      
  To get started, type one of these commands: helpwin, helpdesk, or
demo.
  For information on all of the MathWorks products, type tour.

>> log(1+i)

ans =

   0.3466 + 0.7854i
  

Only the g77 on my linux box is giving me a wrong asnwer. Is this a
known problem? I tried clog() instead of log(), but no difference.
Any help or pointer is appreciated. Thank you in advance.

Do-Hoon Kwon
***********************************************************************

Otras pruebas seqalan lo siguiente:
***********************************************************************
	(Debian 1.3.1, stable)
 gcc             2.7.2.1-9
 g77             0.5.20-1
 libc5           5.4.33-6
	Respuesta: (0.346573591,0.785398185)

	(Debian pre2.0, unstable)
 gcc             2.7.2.3-3
 g77             0.5.21-1
 libc6           2.0.6-3
	Respuesta: (0.240239993,0.785398185)

-por E. Zanardi-
Ademas Enrique Zanardi comenta:
"
Parece que ese error en particular sslo aparece si haces logaritmo de
valores complex*8 (complejos en simple precisisn). Si usas double complex
da el resultado correcto. Tambiin hay por ahm una versisn "rarita" del
gcc+g77 que no da errores (la 2.7.2.3.f.1 o algo asm. La nuestra es la
2.7.2.3) pero no he conseguido encontrarla por ningzn lado. Si alguien me
pasa un URL se lo agradecerma.
"
***********************************************************************


Los que tengan en uso software que usa logaritmos estin pendientes de esto...

  Saludos!

PacoZ.