Phillip Helbig wrote:
>
> > real( double_k) :: pi = 3.1415926_double_k
> >
> > The "_double_k" on the literal is necessary
> > to get a double_k value. Otherwise, it's
> > a default real value, which may not be double_k.
>
> If it weren't for the restriction on initialisation expressions, one
> could write dble(3.1415926).
>
> Of course, if you are using such precision, you probably want to include
> more places.
as the 'double precision' type depends on the machine, the dble intrinsic should
never been used. in place of it, we have to use 'real(3.1415926,double_k)'.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|