--- Comment #1 from pinskia at gcc dot gnu dot org 2009-10-26 15:10 ---
NaNs don't exist in Fortran's data types except for the Fortran 2003's IEEE
types.
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=41831
--- Comment #2 from kargl at gcc dot gnu dot org 2009-10-26 22:39 ---
> My conclusion is that there is a problem when the
> NAN is passed from the C to the fortran.
That's not my conclusion. :)
> double C2F(returnanan)(void)
> {
(snip)
>
> PROGRAM test
> DATA ONE/1.0D0
--- Comment #3 from michael dot baudin at scilab dot org 2009-10-27 13:32
---
You are right, this was a bug in the fortran source code.
I cannot add implicit none in this slatec/gamma routine, but
declaring the return type of the function fixes the bug.
Thank your for looking at this b