On Tue, 2010-07-27 at 08:11 -0400, Jeff Squyres wrote:
> On Jul 26, 2010, at 11:06 PM, Hugo Gagnon wrote:
>
> > 8 integer, parameter :: dp = kind(1.d0)
> > 9 real(kind=dp) :: inside(5), outside(5)
>
> I'm not a fortran expert -- is kind(1.d0) really double precision? According to http://gcc.gnu.org/onlinedocs/gcc-3.4.6/g77/Kind-Notation.html, kind(2) is double precision (but that's for a different compiler, and I don't quite grok the ".d0" notation).
>
Urgh! Thank heavens gcc have moved away from that stupid idea.
kind=8 is normally double precision (and is with gfortran). kind(1.0d0)
is always double precision.
The d (as opposed to e) means DP.
|