Re: isnan in fortran

[c.]

On 11 Apr 2012, at 14:09, John W. Eaton wrote:

> OK, that's a problem.  With the intrinsic ISNAN, the compiler chooses
> the right version based on the type of the argument.  But with Fortran
> 77 and our own function there is no overloading, so we would need to
> use different names.
> 
> jwe

In that case I would rather go back to the sed or cpp option
as (X.NE.X) works with both single and double precision.
c.