Re: getting same output as gfortran, long_float

[robin.vowels@gmail.com]

On Friday, May 1, 2015 at 4:52:52 PM UTC+10, Dmitry A. Kazakov wrote:
> On Thu, 30 Apr 2015 20:12:47 -0500, Nasser M. Abbasi wrote:
> 
> > I found that gfortran can do 128 bit floating point without
> > the use of the compiler switch -fdefault-real-8. Which will
> > map to similar thing as the above Ada construct:
> > 
> > -------------------
> > PROGRAM foo
> > IMPLICIT NONE
> > REAL(KIND = 16) :: x  !-- kind=16 tells it is double quad
> > x = 12.0D0 * 0.0001D0/(1.0D0 * (1.0D0 - 0.1D0)**4 )
> > PRINT *, x
> > END PROGRAM
> > ------------------
> 
> I didn't use FORTRAN for decades, but in good old FORTRAN-IV you declare
> specific lengths using T*n, e.g.
> 
> INTEGER*2
> REAL*4
> REAL*8
> 
> so, I would suggest, maybe naively,
> 
> REAL*16

All those statements are non-standard,
and some current compilers treat them as errors.