[Bug libfortran/24313] complex sqrt function does not return principal value

[harald dot vogt at desy dot de]

------- Comment #12 from harald dot vogt at desy dot de  2006-01-25 21:07 -------
(In reply to comment #9)
> Patch for libgfortran 
> 
> http://gcc.gnu.org/ml/gcc-patches/2005-10/msg00626.html
> 
It is a useful strategy to use libm but...

It seems that glibc has still problems with csqrt. See also
http://sources.redhat.com/bugzilla/show_bug.cgi?id=2181
http://sources.redhat.com/bugzilla/show_bug.cgi?id=2182 .
Therefore it may be useful to check in intrinsics/configure if libm is broken
in the architecture used, and if so use the code in intrinsics/c99_functions.c
.
A possible solution for csqrt is shown in
http://www-zeuthen.desy.de/~hvogt/gfortran/acinclude.m4.diff .
It is better to be independend from glibc coming with the OS. This enables to
install gcc4(gfortran) without patching libm if it has bugs.


-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=24313