pietc.f90 File Reference

Double-precision constants. More...

Go to the source code of this file.

Modules
module	pietc
	Some of the commonly used constants (pi etc) mainly for double-precision subroutines.

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

Detailed Description

Double-precision constants.

Author

R. J. Purser

Date

2014

Definition in file pietc.f90.