Last change
on this file since 3868 was
2228,
checked in by lguez, 9 years ago
|
Correcting a problem from revision 2218. The type double precision
with option "-fdefault-real-8" of gfortran is promoted to 16-byte
precision and there is no specific procedure in arth with this
precision. Could not add a specific procedure in arth with double
precision because, with ifort, the option "-real-size 64" does not
promote the double precision, so that would make two identical
specific procedures in arth.
In module nrtype, replaced double precision by a parameterized real
kind so that the effective precision does not depend on a compiler
option.
In coefpoly, fxhyp, fyhyp and invert_zoom_x, use the parameterized
real kind defined in nrtype, instead of double precision.
Also, in module nrtype, removed unused derived types sprs2_sp and
sprs2_dp.
|
File size:
779 bytes
|
Line | |
---|
1 | MODULE nrtype |
---|
2 | |
---|
3 | implicit none |
---|
4 | |
---|
5 | integer, parameter:: k8 = selected_real_kind(13) |
---|
6 | |
---|
7 | ! Frequently used mathematical constants (with precision to spare): |
---|
8 | |
---|
9 | REAL, PARAMETER :: PI=3.141592653589793238462643383279502884197 |
---|
10 | REAL, PARAMETER :: PIO2=1.57079632679489661923132169163975144209858 |
---|
11 | REAL, PARAMETER :: TWOPI=6.283185307179586476925286766559005768394 |
---|
12 | REAL, PARAMETER :: SQRT2=1.41421356237309504880168872420969807856967 |
---|
13 | REAL, PARAMETER :: EULER=0.5772156649015328606065120900824024310422 |
---|
14 | |
---|
15 | REAL(K8), PARAMETER:: & |
---|
16 | PI_D = 3.141592653589793238462643383279502884197_k8 |
---|
17 | REAL(K8), PARAMETER:: & |
---|
18 | PIO2_D=1.57079632679489661923132169163975144209858_k8 |
---|
19 | REAL(K8), PARAMETER:: & |
---|
20 | TWOPI_D=6.283185307179586476925286766559005768394_k8 |
---|
21 | |
---|
22 | END MODULE nrtype |
---|
Note: See
TracBrowser
for help on using the repository browser.