View previous topic :: View next topic |
Author |
Message |
Wolfgang Hoeppner
Joined: 05 Dec 2006 Posts: 8 Location: Sibiu / Romania
|
Posted: Sat Jul 27, 2013 5:13 pm Post subject: Quad(ruple precision; REAL*16) support |
|
|
Will there be, in one of the next releases of FTN95, a support of quadruple precision arithmetic, at least for fundamental arithmetics (+ - * /), but preferably also for SQRT, SIN, COS, &c.?
Best wishes
Wolfgang. |
|
Back to top |
|
 |
PaulLaidler Site Admin
Joined: 21 Feb 2005 Posts: 8210 Location: Salford, UK
|
Posted: Sat Jul 27, 2013 8:35 pm Post subject: |
|
|
At the moment there are no plans to support REAL*16 arithmetic.
One approach might be to use the GMP (GNU) library. It would be helpful to know if anyone has experience of accessing this library from Fortran. |
|
Back to top |
|
 |
PaulLaidler Site Admin
Joined: 21 Feb 2005 Posts: 8210 Location: Salford, UK
|
Posted: Sun Jul 28, 2013 8:56 am Post subject: |
|
|
As a quick test I have added the following PRINT statement to a working 64 bit ClearWin+/gFortran program...
print*,sqrt(2.0_16)
and the output is
Code: | 1.41421356237309504880168872420969798E+0000 |
so it looks like you could work with 64 bit gFortran and ClearWin+.
p.s. It think that this value is correct up to 32 decimal places. |
|
Back to top |
|
 |
DanRRight
Joined: 10 Mar 2008 Posts: 2923 Location: South Pole, Antarctica
|
Posted: Mon Jul 29, 2013 6:02 pm Post subject: |
|
|
What order of exponents real*16 allows? |
|
Back to top |
|
 |
davidb
Joined: 17 Jul 2009 Posts: 560 Location: UK
|
Posted: Mon Jul 29, 2013 8:35 pm Post subject: |
|
|
Approximately 1.0E-4932 to 1.0E4932. _________________ Programmer in: Fortran 77/95/2003/2008, C, C++ (& OpenMP), java, Python, Perl |
|
Back to top |
|
 |
|