17 logical ,
parameter:: 
t=.true. 
    18 logical ,
parameter:: 
f=.false. 
    19 real(dp),
parameter:: 
u0=0_dp 
    20 real(dp),
parameter:: 
u1=1_dp 
    22 real(dp),
parameter:: 
u2=2_dp 
    24 real(dp),
parameter:: 
u3=3_dp 
    26 real(dp),
parameter:: 
u4=4_dp 
    28 real(dp),
parameter:: 
u5=5_dp 
    30 real(dp),
parameter:: 
u6=6_dp 
    42 real(dp),
parameter:: 
pi =3.1415926535897932384626433832795028841971693993751058209749e0_dp 
    43 real(dp),
parameter:: 
pi2=6.2831853071795864769252867665590057683943387987502116419498e0_dp 
    44 real(dp),
parameter:: 
pih=1.5707963267948966192313216916397514420985846996875529104874e0_dp 
    45 real(dp),
parameter:: 
rpi=1.7724538509055160272981674833411451827975494561223871282138e0_dp 
    46 real(dp),
parameter:: 
r2 =1.4142135623730950488016887242096980785696718753769480731766e0_dp 
    47 real(dp),
parameter:: 
r3 =1.7320508075688772935274463415058723669428052538103806280558e0_dp 
    48 real(dp),
parameter:: 
r5 =2.2360679774997896964091736687312762354406183596115257242708e0_dp 
    52 real(dp),
parameter:: 
phi=1.6180339887498948482045868343656381177203091798057628621354e0_dp 
    53 real(dp),
parameter:: 
euler=0.57721566490153286060651209008240243104215933593992359880e0_dp 
    57 real(dp),
parameter:: 
s10=.173648177666930348851716626769314796000375677184069387236241e0_dp 
    58 real(dp),
parameter:: 
s11=.195090322016128267848284868477022240927691617751954807754502e0_dp 
    59 real(dp),
parameter:: 
s13=.222520933956314404288902564496794759466355568764544955311987e0_dp 
    60 real(dp),
parameter:: 
s15=.258819045102520762348898837624048328349068901319930513814003e0_dp 
    61 real(dp),
parameter:: 
s18=.309016994374947424102293417182819058860154589902881431067724e0_dp 
    62 real(dp),
parameter:: 
s20=.342020143325668733044099614682259580763083367514160628465048e0_dp 
    63 real(dp),
parameter:: 
s22=.382683432365089771728459984030398866761344562485627041433800e0_dp 
    64 real(dp),
parameter:: 
s26=.433883739117558120475768332848358754609990727787459876444547e0_dp 
    66 real(dp),
parameter:: 
s34=.555570233019602224742830813948532874374937190754804045924153e0_dp 
    67 real(dp),
parameter:: 
s36=.587785252292473129168705954639072768597652437643145991072272e0_dp 
    68 real(dp),
parameter:: 
s39=.623489801858733530525004884004239810632274730896402105365549e0_dp 
    69 real(dp),
parameter:: 
s40=.642787609686539326322643409907263432907559884205681790324977e0_dp 
    71 real(dp),
parameter:: 
s50=.766044443118978035202392650555416673935832457080395245854045e0_dp 
    72 real(dp),
parameter:: 
s51=.781831482468029808708444526674057750232334518708687528980634e0_dp 
    73 real(dp),
parameter:: 
s54=.809016994374947424102293417182819058860154589902881431067724e0_dp 
    74 real(dp),
parameter:: 
s56=.831469612302545237078788377617905756738560811987249963446124e0_dp 
    76 real(dp),
parameter:: 
s64=.900968867902419126236102319507445051165919162131857150053562e0_dp 
    77 real(dp),
parameter:: 
s68=.923879532511286756128183189396788286822416625863642486115097e0_dp 
    78 real(dp),
parameter:: 
s70=.939692620785908384054109277324731469936208134264464633090286e0_dp 
    79 real(dp),
parameter:: 
s72=.951056516295153572116439333379382143405698634125750222447305e0_dp 
    80 real(dp),
parameter:: 
s75=.965925826289068286749743199728897367633904839008404550402343e0_dp 
    81 real(dp),
parameter:: 
s77=.974927912181823607018131682993931217232785800619997437648079e0_dp 
    82 real(dp),
parameter:: 
s79=.980785280403230449126182236134239036973933730893336095002916e0_dp 
    83 real(dp),
parameter:: 
s80=.984807753012208059366743024589523013670643251719842418790025e0_dp 
 complex(dpc), parameter z288
exp(2*pi*i*4/5) 
 
real(dp), parameter u4
four 
 
complex(dpc), parameter z347
exp(2*pi*i*27/28) 
 
real(dp), parameter ms60
minus-sine(60 deg) 
 
complex(dpc), parameter z162
exp(2*pi*i*9/20) 
 
real(dp), parameter s70
sine(70 deg) 
 
complex(dpc), parameter z051
exp(2*pi*i/7) 
 
integer, parameter dp
Double precision real kind. 
 
complex(dpc), parameter z220
exp(2*pi*i*11/18) 
 
complex(dpc), parameter z225
exp(2*pi*i*5/8) 
 
complex(dpc), parameter z020
exp(2*pi*i/18) 
 
complex(dpc), parameter z234
exp(2*pi*i*13/20) 
 
complex(dpc), parameter z244
exp(2*pi*i*19/28) 
 
complex(dpc), parameter z108
exp(2*pi*i*3/10) 
 
complex(dpc), parameter z270
exp(2*pi*i*3/4) 
 
real(dp), parameter s56
sine(90*5/8 deg) 
 
real(dp), parameter s15
sine(15 deg) 
 
real(dp), parameter u6
six 
 
Standard integer, real, and complex single and double precision kinds. 
 
complex(dpc), parameter z105
exp(2*pi*i*7/24) 
 
complex(dpc), parameter z000
exp(2*pi*i*0) 
 
real(dp), parameter ms40
minus-sine(40 deg) 
 
complex(dpc), parameter z030
exp(2*pi*i/12) 
 
complex(dpc), parameter z257
exp(2*pi*i*5/7) 
 
complex(dpc), parameter z210
exp(2*pi*i*7/12) 
 
real(dp), parameter s11
sine(90/8 deg) 
 
complex(dpc), parameter z338
exp(2*pi*i*15/16) 
 
real(dp), parameter mu1
minus-one 
 
complex(dpc), parameter z026
exp(2*pi*i/14) 
 
real(dp), parameter pih
pi*half 
 
complex(dpc), parameter z240
exp(2*pi*i*2/3) 
 
complex(dpc), parameter z206
exp(2*pi*i*4/7) 
 
complex(dpc), parameter z190
exp(2*pi*i*19/36) 
 
complex(dpc), parameter z018
exp(2*pi*i/20) 
 
real(dp), parameter mu5
minus-five 
 
complex(dpc), parameter z345
exp(2*pi*i*23/24) 
 
logical, parameter f
for pain-relief in logical ops 
 
complex(dpc), parameter z015
exp(2*pi*i/24) 
 
complex(dpc), parameter z334
exp(2*pi*i*13/14) 
 
real(dp), parameter s75
sine(75 deg) 
 
complex(dpc), parameter z090
exp(2*pi*i/4) 
 
real(dp), parameter o5
fifth 
 
real(dp), parameter s18
sine(18 deg) 
 
real(dp), parameter mu3
minus-three 
 
real(dp), parameter u0
zero 
 
complex(dpc), parameter z280
exp(2*pi*i*7/9) 
 
real(dp), parameter ms75
minus-sine(75 deg) 
 
real(dp), parameter ms15
minus-sine(15 deg) 
 
complex(dpc), parameter ci
complex i=sqrt(-1) 
 
real(dp), parameter u2
two 
 
real(dp), parameter s54
sine(54 deg) 
 
complex(dpc), parameter z195
exp(2*pi*i*13/24) 
 
real(dp), parameter s50
sine(50 deg) 
 
complex(dpc), parameter z248
exp(2*pi*i*11/16) 
 
complex(dpc), parameter z259
exp(2*pi*i*23/32) 
 
complex(dpc), parameter cipi
complex i*pi 
 
complex(dpc), parameter z169
exp(2*pi*i*15/32) 
 
complex(dpc), parameter z296
exp(2*pi*i*23/28) 
 
real(dp), parameter r2
Square root of 2. 
 
complex(dpc), parameter z250
exp(2*pi*i*25/36) 
 
real(dp), parameter or2
1. 
 
complex(dpc), parameter z260
exp(2*pi*i*13/18) 
 
complex(dpc), parameter z200
exp(2*pi*i*5/9) 
 
real(dp), parameter s34
sine(90*3/8 deg) 
 
complex(dpc), parameter z070
exp(2*pi*i*7/36) 
 
complex(dpc), parameter z231
exp(2*pi*i*9/14) 
 
complex(dpc), parameter z252
exp(2*pi*i*7/10) 
 
complex(dpc), parameter c1
complex one 
 
complex(dpc), parameter z167
exp(2*pi*i*13/28) 
 
real(dp), parameter s13
sine(90/7 deg) 
 
complex(dpc), parameter z290
exp(2*pi*i*29/36) 
 
complex(dpc), parameter z080
exp(2*pi*i*2/9) 
 
real(dp), parameter o6
sixth 
 
complex(dpc), parameter z130
exp(2*pi*i*13/36) 
 
complex(dpc), parameter z034
exp(2*pi*i*3/32) 
 
complex(dpc), parameter z144
exp(2*pi*i*2/5) 
 
real(dp), parameter ms36
minus-sine(36 deg) 
 
real(dp), parameter u5
five 
 
complex(dpc), parameter z230
exp(2*pi*i*23/36) 
 
real(dp), parameter mu2
minus-two 
 
complex(dpc), parameter z214
exp(2*pi*i*19/32) 
 
complex(dpc), parameter z193
exp(2*pi*i*15/28) 
 
real(dp), parameter phi
Golden number. 
 
complex(dpc), parameter z315
exp(2*pi*i*7/8) 
 
real(dp), parameter or3
1. 
 
complex(dpc), parameter z054
exp(2*pi*i*3/20) 
 
real(dp), parameter mu4
minus-four 
 
complex(dpc), parameter z056
exp(2*pi*i*5/32) 
 
real(dp), parameter pi
Pi. 
 
real(dp), parameter mu6
minus-six 
 
complex(dpc), parameter c0
complex zero 
 
complex(dpc), parameter z036
exp(2*pi*i/10) 
 
complex(dpc), parameter z160
exp(2*pi*i*4/9) 
 
complex(dpc), parameter z255
exp(2*pi*i*17/24) 
 
complex(dpc), parameter z120
exp(2*pi*i/3) 
 
real(dp), parameter s77
sine(90*6/7 deg) 
 
real(dp), parameter s80
sine(80 deg) 
 
real(dp), parameter ms45
minus-sine(45 deg) 
 
real(dp), parameter u3
three 
 
complex(dpc), parameter z324
exp(2*pi*i*9/10) 
 
complex(dpc), parameter z202
exp(2*pi*i*9/16) 
 
complex(dpc), parameter z079
exp(2*pi*i*7/32) 
 
complex(dpc), parameter z010
exp(2*pi*i/36) 
 
real(dp), parameter r5
Square root of 5. 
 
real(dp), parameter ms18
minus-sine(18 deg) 
 
complex(dpc), parameter z326
exp(2*pi*i*29/32) 
 
real(dp), parameter s20
sine(20 deg) 
 
complex(dpc), parameter z100
exp(2*pi*i*5/18) 
 
logical, parameter t
for pain-relief in logical ops 
 
real(dp), parameter ms50
minus-sine(50 deg) 
 
complex(dpc), parameter z141
exp(2*pi*i*11/28) 
 
complex(dpc), parameter z140
exp(2*pi*i*7/18) 
 
real(dp), parameter mo3
minus-third 
 
real(dp), parameter ms34
minus-sine(90*3/8 deg) 
 
real(dp), parameter s60
sine(60 deg) 
 
complex(dpc), parameter z022
exp(2*pi*i/16) 
 
complex(dpc), parameter mci
complex minus-i 
 
real(dp), parameter or5
1. 
 
complex(dpc), parameter z304
exp(2*pi*i*27/32) 
 
complex(dpc), parameter z340
exp(2*pi*i*17/18) 
 
real(dp), parameter dtor
Degrees to radians. 
 
real(dp), parameter s45
sine(45 deg) 
 
real(dp), parameter mo6
minus-sixth 
 
real(dp), parameter pi2
Pi*2. 
 
complex(dpc), parameter z060
exp(2*pi*i/6) 
 
real(dp), parameter ms80
minus-sine(80 deg) 
 
complex(dpc), parameter z112
exp(2*pi*i*5/16) 
 
complex(dpc), parameter z064
exp(2*pi*i*5/28) 
 
real(dp), parameter s22
sine(22.5 deg) 
 
real(dp), parameter mo4
minus-fourth 
 
complex(dpc), parameter z154
exp(2*pi*i*3/7) 
 
complex(dpc), parameter z330
exp(2*pi*i*11/12) 
 
real(dp), parameter s64
sine(90*5/7 deg) 
 
real(dp), parameter ms77
minus-sine(90*6/7 deg) 
 
complex(dpc), parameter z219
exp(2*pi*i*17/28) 
 
real(dp), parameter s79
sine(90*7/8 deg) 
 
real(dp), parameter ms26
minus-sine(90*2/7 deg) 
 
complex(dpc), parameter z124
exp(2*pi*i*11/32) 
 
complex(dpc), parameter z236
exp(2*pi*i*21/32) 
 
real(dp), parameter rtod
radians to degrees 
 
complex(dpc), parameter z300
exp(2*pi*i*5/6) 
 
real(dp), parameter ms70
minus-sine(70 deg) 
 
real(dp), parameter ms54
minus-sine(54 deg) 
 
complex(dpc), parameter z040
exp(2*pi*i/9) 
 
real(dp), parameter euler
Euler-Mascheroni constant. 
 
complex(dpc), parameter z072
exp(2*pi*i/5) 
 
complex(dpc), parameter z129
exp(2*pi*i*5/14) 
 
real(dp), parameter u1
one 
 
complex(dpc), parameter z321
exp(2*pi*i*25/28) 
 
complex(dpc), parameter z191
exp(2*pi*i*17/32) 
 
real(dp), parameter ms13
minus-sine(90/7 deg) 
 
complex(dpc), parameter z039
exp(2*pi*i*3/28) 
 
complex(dpc), parameter z116
exp(2*pi*i*9/28) 
 
complex(dpc), parameter z068
exp(2*pi*i*3/16) 
 
complex(dpc), parameter z342
exp(2*pi*i*19/20) 
 
complex(dpc), parameter z077
exp(2*pi*i*3/14) 
 
complex(dpc), parameter z158
exp(2*pi*i*7/16) 
 
real(dp), parameter r3
Square root of 3. 
 
real(dp), parameter s26
sine(90*2/7 deg) 
 
real(dp), parameter mo5
minus-fifth 
 
Some of the commonly used constants (pi etc) mainly for double-precision subroutines. 
 
real(dp), parameter s68
sine(67.5 deg) 
 
complex(dpc), parameter z150
exp(2*pi*i*5/12) 
 
complex(dpc), parameter z180
exp(2*pi*i/2) 
 
complex(dpc), parameter z075
exp(2*pi*i*5/24) 
 
complex(dpc), parameter z013
exp(2*pi*i/28) 
 
real(dp), parameter o3
third 
 
real(dp), parameter s72
sine(72 deg) 
 
real(dp), parameter ms22
minus-sine(22.5 deg) 
 
complex(dpc), parameter z309
exp(2*pi*i*6/7) 
 
complex(dpc), parameter z045
exp(2*pi*i/8) 
 
real(dp), parameter s36
sine(36 deg) 
 
real(dp), parameter ms11
minus-sine(90/8 deg) 
 
integer, parameter dpc
Double precision real kind. 
 
real(dp), parameter ms30
minus-sine(30 deg) 
 
complex(dpc), parameter z011
exp(2*pi*i/32) 
 
complex(dpc), parameter z101
exp(2*pi*i*9/32) 
 
real(dp), parameter s40
sine(40 deg) 
 
complex(dpc), parameter z165
exp(2*pi*i*11/24) 
 
real(dp), parameter o4
fourth 
 
complex(dpc), parameter z285
exp(2*pi*i*19/24) 
 
complex(dpc), parameter z146
exp(2*pi*i*13/32) 
 
complex(dpc), parameter z103
exp(2*pi*i*2/7) 
 
complex(dpc), parameter z349
exp(2*pi*i*31/32) 
 
real(dp), parameter ms68
minus-sine(67.5 deg) 
 
complex(dpc), parameter z310
exp(2*pi*i*31/36) 
 
complex(dpc), parameter z283
exp(2*pi*i*11/14) 
 
real(dp), parameter ms39
minus-sine(90*3/7 deg) 
 
complex(dpc), parameter z306
exp(2*pi*i*17/20) 
 
complex(dpc), parameter z281
exp(2*pi*i*25/32) 
 
real(dp), parameter ms64
minus-sine(90*5/7 deg) 
 
real(dp), parameter rpi
square root of pi 
 
complex(dpc), parameter z110
exp(2*pi*i*11/36) 
 
complex(dpc), parameter z050
exp(2*pi*i*5/36) 
 
complex(dpc), parameter z135
exp(2*pi*i*3/8) 
 
real(dp), parameter ms20
minus-sine(20 deg) 
 
complex(dpc), parameter mc1
complex minus-one 
 
real(dp), parameter ms10
minus-sine(10 deg) 
 
real(dp), parameter s30
sine(30 deg) 
 
complex(dpc), parameter z292
exp(2*pi*i*13/16) 
 
complex(dpc), parameter z198
exp(2*pi*i*11/20) 
 
complex(dpc), parameter z320
exp(2*pi*i*8/9) 
 
real(dp), parameter ms51
minus-sine(90*4/7 deg) 
 
real(dp), parameter o2
half 
 
real(dp), parameter ms56
minus-sine(90*5/8 deg) 
 
complex(dpc), parameter z170
exp(2*pi*i*17/36) 
 
real(dp), parameter mo2
minus-half 
 
complex(dpc), parameter z350
exp(2*pi*i*35/36) 
 
real(dp), parameter ms72
minus-sine(72 deg) 
 
real(dp), parameter s51
sine(90*4/7 deg) 
 
complex(dpc), parameter z216
exp(2*pi*i*3/5) 
 
real(dp), parameter s39
sine(90*3/7 deg) 
 
real(dp), parameter ms79
minus-sine(90*7/8 deg) 
 
complex(dpc), parameter z126
exp(2*pi*i*7/20) 
 
real(dp), parameter s10
sine(10 deg)