From mboxrd@z Thu Jan  1 00:00:00 1970
X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on polar.synack.me
X-Spam-Level: 
X-Spam-Status: No, score=-0.3 required=5.0 tests=BAYES_00,
	REPLYTO_WITHOUT_TO_CC autolearn=no autolearn_force=no version=3.4.4
Path: 
 eternal-september.org!reader01.eternal-september.org!reader02.eternal-september.org!news.eternal-september.org!mx02.eternal-september.org!feeder.eternal-september.org!aioe.org!.POSTED!not-for-mail
From: "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de>
Newsgroups: comp.lang.ada
Subject: Re: getting same output as gfortran, long_float
Date: Fri, 1 May 2015 08:52:43 +0200
Organization: cbb software GmbH
Message-ID: <1i8x3r1feyzkt$.j85il7e3wpv9.dlg@40tude.net>
References: <mhu65r$ntj$1@speranza.aioe.org>
 <1kxou0nloqg9c$.1x0itzgdrlosm$.dlg@40tude.net>
 <mhujuf$jbs$1@speranza.aioe.org>
Reply-To: mailbox@dmitry-kazakov.de
NNTP-Posting-Host: evoS9sCOdnHjo0GRLLMU1Q.user.speranza.aioe.org
Mime-Version: 1.0
Content-Type: text/plain; charset="us-ascii"
Content-Transfer-Encoding: 7bit
X-Complaints-To: abuse@aioe.org
User-Agent: 40tude_Dialog/2.0.15.1
X-Notice: Filtered by postfilter v. 0.8.2
Xref: news.eternal-september.org comp.lang.ada:25678
Date: 2015-05-01T08:52:43+02:00
List-Id: <comp.lang.ada>

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

-- 
Regards,
Dmitry A. Kazakov
http://www.dmitry-kazakov.de