Intel Fortran compiler

ksasidhar — Tue, 25 Mar 2008 00:12:24 GMT

Hello,

I am using intel fortran 10.0.023 for my simulations. I am trying to carry out specific computation which works on few computers and does not work on others. The statement is as follows.

DO i=1,3
IF (position(i) == TWOPI) THEN
position(i) = TWOPI - 1E-6
END IF
END DO

Can you please let me know ifthere is anythingto be changed here?

Thanks in advance.

Regards,

Sasidhar

Re: Intel Fortran compiler

jimdempseyatthecove — Tue, 25 Mar 2008 04:25:19 GMT

Sasidhar,

Is position real(4) or real(8)?
is TWOPI real(4) or real(8)?

If TWOPI = 3.14159..., is 3.14159... a single precisionor double precision constant? (Match the literal type with the variable type as constants are not promoted)

1E-6 is single precision. Subtracting a 1E-6 from a 3E+0 will require 7 digits of precision which are not available in single precision.

If position and TWOPI are real(8) then use 1D-6 for the literal or perhapse you are looking for

position(i) = TWOPI- TINY(TWOPI)

This will produce the largest number less than TWOPI representable in the precision of what is used for TWOPI. (Hopefully position is the same precision as TWOPI).

Jim Dempsey

topic Intel Fortran compiler in Intel® Fortran Compiler

Intel Fortran compiler

Re: Intel Fortran compiler