Intel® Fortran Compiler
Build applications that can scale for the future with optimized code designed for Intel® Xeon® and compatible processors.
29248 Discussions

double precision - loss of precision

rssmps
Beginner
375 Views
I've got aproblem with precision loss.
does anyone know if this isa 6.6B problem?
module common
double precision ::PI=3.141592653589793238D0, GC=32.174
contains
FUNCTION FCAP(N) RESULT(OUT)
USE GTAC
IMPLICIT NONE
DOUBLE PRECISION OUT
INTEGER N
OUT=RHO(N)*CP(N)*PI*D(N)*D(N)/4.0*(DX(N)+DY(N))
END FUNCTION
end module
when I call the fcap function, using the debuggerI noticed that
PI =3.14159250259399
3.141592653589793238D0
^^^
what happened here?
Has anyone else encountered this before?
Thanks!

0 Kudos
0 Replies
Reply