From babf2647fb613d9f8d3957b4d07dca48ff18ac66 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Fri, 11 May 2012 07:35:52 +0000 Subject: [PATCH] use non signaling version of NaN --- code/prec.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/prec.f90 b/code/prec.f90 index bdf3e5c7e..465a27391 100644 --- a/code/prec.f90 +++ b/code/prec.f90 @@ -41,9 +41,9 @@ module prec ! from http://www.hpc.unimelb.edu.au/doc/f90lrm/dfum_035.html ! copy can be found in documentation/Code/Fortran #ifdef LEGACY_COMPILER - real(pReal), parameter, public :: DAMASK_NaN = Z'7FF0000000000001' + real(pReal), parameter, public :: DAMASK_NaN = Z'7FF8000000000000' ! quiet NaN #else - real(pReal), parameter, public :: DAMASK_NaN = real(Z'7FF0000000000001', pReal) + real(pReal), parameter, public :: DAMASK_NaN = real(Z'7FF8000000000000', pReal) #endif type, public :: p_vec @@ -71,7 +71,7 @@ subroutine prec_init write(6,'(a,i3)') ' Bytes for pReal: ',pReal write(6,'(a,i3)') ' Bytes for pInt: ',pInt write(6,'(a,i3)') ' Bytes for pLongInt: ',pLongInt - write(6,'(a,e10.3)') ' NaN: ',DAMASK_NAN + write(6,'(a,e10.3)') ' NaN: ',DAMASK_NaN write(6,'(a,l3)') ' NaN /= NaN: ',DAMASK_NaN/=DAMASK_NaN if (DAMASK_NaN == DAMASK_NaN) call quit(9000) write(6,*)