This is the mail archive of the
fortran@gcc.gnu.org
mailing list for the GNU Fortran project.
Re: legacy code with real*16
- From: Pierre-Matthieu anglade <pierre dot matthieu dot anglade at gmail dot com>
- To: James Mitchell Allmond <gte203h at prism dot gatech dot edu>
- Cc: fortran at gcc dot gnu dot org
- Date: Tue, 29 Nov 2005 08:47:01 +0100
- Subject: Re: legacy code with real*16
- Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=beta; d=gmail.com; h=received:message-id:date:from:to:subject:cc:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=qvbNHVRhGElwgxDM5W5pWjlHBNFXSdoE8a+3hOjkNBdBcO/u6vqAPZWIrPwmDHDbhFV9lQhsFxgV+Rhmb+7nYTZfxIWwtjJV/E46cWig37br1hk2YuW4AZfsC1F6T3luxEKOzCheHBj07wHuwr+LodsPyyuZD7HgZThA4NgyG18=
- References: <8A093AC4-7079-4FA8-809E-088C7FE768AE@prism.gatech.edu> <72f4de7342dcd370ca0ba4864ef7385a@nasa.gov>
If you keep using x86 hardware, you might not find some hardware
quadruple precision numbers. Still you may give a try to libs that
implement it in software.
I don't know if this solution might be suitable for your case but if
you translate to C (you can use f2c to do that) then use cephes
(http://www.moshier.net/) to implement extended precision you might
then get (slowly) the precision you want.
--
Pierre-Matthieu Anglade