pietc_s.f90

Go to the documentation of this file.

module pietc_s
use pkind, only: sp,spc
implicit none
logical ,parameter:: T=.true.
logical ,parameter:: F=.false.
real(sp),parameter:: u0=0_sp
real(sp),parameter:: u1=1_sp
real(sp),parameter:: mu1=-u1
real(sp),parameter:: u2=2_sp
real(sp),parameter:: mu2=-u2
real(sp),parameter:: u3=3_sp
real(sp),parameter:: mu3=-u3
real(sp),parameter:: u4=4_sp
real(sp),parameter:: mu4=-u4
real(sp),parameter:: u5=5_sp
real(sp),parameter:: mu5=-u5
real(sp),parameter:: u6=6_sp
real(sp),parameter:: mu6=-u6
real(sp),parameter:: o2=u1/u2
real(sp),parameter:: o3=u1/u3
real(sp),parameter:: o4=u1/u4
real(sp),parameter:: o5=u1/u5
real(sp),parameter:: o6=u1/u6
real(sp),parameter:: mo2=-o2
real(sp),parameter:: mo3=-o3
real(sp),parameter:: mo4=-o4
real(sp),parameter:: mo5=-o5
real(sp),parameter:: mo6=-06
real(sp),parameter:: pi =3.1415926535897932384626433832795028841971693993751058209749e0_sp
real(sp),parameter:: pi2=6.2831853071795864769252867665590057683943387987502116419498e0_sp
real(sp),parameter:: pih=1.5707963267948966192313216916397514420985846996875529104874e0_sp
real(sp),parameter:: rpi=1.7724538509055160272981674833411451827975494561223871282138e0_sp
real(sp),parameter:: r2 =1.4142135623730950488016887242096980785696718753769480731766e0_sp
real(sp),parameter:: r3 =1.7320508075688772935274463415058723669428052538103806280558e0_sp
real(sp),parameter:: r5 =2.2360679774997896964091736687312762354406183596115257242708e0_sp
real(sp),parameter:: or2=u1/r2
real(sp),parameter:: or3=u1/r3
real(sp),parameter:: or5=u1/r5
real(sp),parameter:: phi=1.6180339887498948482045868343656381177203091798057628621354e0_sp
real(sp),parameter:: euler=0.57721566490153286060651209008240243104215933593992359880e0_sp
! Degree to radians; radians to degrees:
real(sp),parameter:: dtor=pi/180
real(sp),parameter:: rtod=180/pi
! Sines of all main fractions of 90 degrees (down to ninths):
real(sp),parameter:: s10=.173648177666930348851716626769314796000375677184069387236241e0_sp
real(sp),parameter:: s11=.195090322016128267848284868477022240927691617751954807754502e0_sp
real(sp),parameter:: s13=.222520933956314404288902564496794759466355568764544955311987e0_sp
real(sp),parameter:: s15=.258819045102520762348898837624048328349068901319930513814003e0_sp
real(sp),parameter:: s18=.309016994374947424102293417182819058860154589902881431067724e0_sp
real(sp),parameter:: s20=.342020143325668733044099614682259580763083367514160628465048e0_sp
real(sp),parameter:: s22=.382683432365089771728459984030398866761344562485627041433800e0_sp
real(sp),parameter:: s26=.433883739117558120475768332848358754609990727787459876444547e0_sp
real(sp),parameter:: s30=o2
real(sp),parameter:: s34=.555570233019602224742830813948532874374937190754804045924153e0_sp
real(sp),parameter:: s36=.587785252292473129168705954639072768597652437643145991072272e0_sp
real(sp),parameter:: s39=.623489801858733530525004884004239810632274730896402105365549e0_sp
real(sp),parameter:: s40=.642787609686539326322643409907263432907559884205681790324977e0_sp
real(sp),parameter:: s45=or2
real(sp),parameter:: s50=.766044443118978035202392650555416673935832457080395245854045e0_sp
real(sp),parameter:: s51=.781831482468029808708444526674057750232334518708687528980634e0_sp
real(sp),parameter:: s54=.809016994374947424102293417182819058860154589902881431067724e0_sp
real(sp),parameter:: s56=.831469612302545237078788377617905756738560811987249963446124e0_sp
real(sp),parameter:: s60=r3*o2
real(sp),parameter:: s64=.900968867902419126236102319507445051165919162131857150053562e0_sp
real(sp),parameter:: s68=.923879532511286756128183189396788286822416625863642486115097e0_sp
real(sp),parameter:: s70=.939692620785908384054109277324731469936208134264464633090286e0_sp
real(sp),parameter:: s72=.951056516295153572116439333379382143405698634125750222447305e0_sp
real(sp),parameter:: s75=.965925826289068286749743199728897367633904839008404550402343e0_sp
real(sp),parameter:: s77=.974927912181823607018131682993931217232785800619997437648079e0_sp
real(sp),parameter:: s79=.980785280403230449126182236134239036973933730893336095002916e0_sp
real(sp),parameter:: s80=.984807753012208059366743024589523013670643251719842418790025e0_sp
!... and their minuses:
real(sp),parameter:: ms10=-s10
real(sp),parameter:: ms11=-s11
real(sp),parameter:: ms13=-s13
real(sp),parameter:: ms15=-s15
real(sp),parameter:: ms18=-s18
real(sp),parameter:: ms20=-s20
real(sp),parameter:: ms22=-s22
real(sp),parameter:: ms26=-s26
real(sp),parameter:: ms30=-s30
real(sp),parameter:: ms34=-s34
real(sp),parameter:: ms36=-s36
real(sp),parameter:: ms39=-s39
real(sp),parameter:: ms40=-s40
real(sp),parameter:: ms45=-s45
real(sp),parameter:: ms50=-s50
real(sp),parameter:: ms51=-s51
real(sp),parameter:: ms54=-s54
real(sp),parameter:: ms56=-s56
real(sp),parameter:: ms60=-s60
real(sp),parameter:: ms64=-s64
real(sp),parameter:: ms68=-s68
real(sp),parameter:: ms70=-s70
real(sp),parameter:: ms72=-s72
real(sp),parameter:: ms75=-s75
real(sp),parameter:: ms77=-s77
real(sp),parameter:: ms79=-s79
real(sp),parameter:: ms80=-s80

complex(spc),parameter::  c0=(u0,u0)
complex(spc),parameter::  c1=(u1,u0)
complex(spc),parameter::  mc1=-c1
complex(spc),parameter::  ci=(u0,u1)
complex(spc),parameter::  mci=-ci
complex(spc),parameter::  cipi=ci*pi
! Main fractional rotations, as unimodualr complex numbers:
complex(spc),parameter::  z000=c1
complex(spc),parameter::  z010=( s80,s10)
complex(spc),parameter::  z011=( s79,s11)
complex(spc),parameter::  z013=( s77,s13)
complex(spc),parameter::  z015=( s75,s15)
complex(spc),parameter::  z018=( s72,s18)
complex(spc),parameter::  z020=( s70,s20)
complex(spc),parameter::  z022=( s68,s22)
complex(spc),parameter::  z026=( s64,s26)
complex(spc),parameter::  z030=( s60,s30)
complex(spc),parameter::  z034=( s56,s34)
complex(spc),parameter::  z036=( s54,s36)
complex(spc),parameter::  z039=( s51,s39)
complex(spc),parameter::  z040=( s50,s40)
complex(spc),parameter::  z045=( s45,s45)
complex(spc),parameter::  z050=( s40,s50)
complex(spc),parameter::  z051=( s39,s51)
complex(spc),parameter::  z054=( s36,s54)
complex(spc),parameter::  z056=( s34,s56)
complex(spc),parameter::  z060=( s30,s60)
complex(spc),parameter::  z064=( s26,s64)
complex(spc),parameter::  z068=( s22,s68)
complex(spc),parameter::  z070=( s20,s70)
complex(spc),parameter::  z072=( s18,s72)
complex(spc),parameter::  z075=( s15,s75)
complex(spc),parameter::  z077=( s13,s77)
complex(spc),parameter::  z079=( s11,s79)
complex(spc),parameter::  z080=( s10,s80)
complex(spc),parameter::  z090=ci
complex(spc),parameter::  z100=(ms10,s80)
complex(spc),parameter::  z101=(ms11,s79)
complex(spc),parameter::  z103=(ms13,s77)
complex(spc),parameter::  z105=(ms15,s75)
complex(spc),parameter::  z108=(ms18,s72)
complex(spc),parameter::  z110=(ms20,s70)
complex(spc),parameter::  z112=(ms22,s68)
complex(spc),parameter::  z116=(ms26,s64)
complex(spc),parameter::  z120=(ms30,s60)
complex(spc),parameter::  z124=(ms34,s56)
complex(spc),parameter::  z126=(ms36,s54)
complex(spc),parameter::  z129=(ms39,s51)
complex(spc),parameter::  z130=(ms40,s50)
complex(spc),parameter::  z135=(ms45,s45)
complex(spc),parameter::  z140=(ms50,s40)
complex(spc),parameter::  z141=(ms51,s39)
complex(spc),parameter::  z144=(ms54,s36)
complex(spc),parameter::  z146=(ms56,s34)
complex(spc),parameter::  z150=(ms60,s30)
complex(spc),parameter::  z154=(ms64,s26)
complex(spc),parameter::  z158=(ms68,s22)
complex(spc),parameter::  z160=(ms70,s20)
complex(spc),parameter::  z162=(ms72,s18)
complex(spc),parameter::  z165=(ms75,s15)
complex(spc),parameter::  z167=(ms77,s13)
complex(spc),parameter::  z169=(ms79,s11)
complex(spc),parameter::  z170=(ms80,s10)
complex(spc),parameter::  z180=-z000
complex(spc),parameter::  z190=-z010
complex(spc),parameter::  z191=-z011
complex(spc),parameter::  z193=-z013
complex(spc),parameter::  z195=-z015
complex(spc),parameter::  z198=-z018
complex(spc),parameter::  z200=-z020
complex(spc),parameter::  z202=-z022
complex(spc),parameter::  z206=-z026
complex(spc),parameter::  z210=-z030
complex(spc),parameter::  z214=-z034
complex(spc),parameter::  z216=-z036
complex(spc),parameter::  z219=-z039
complex(spc),parameter::  z220=-z040
complex(spc),parameter::  z225=-z045
complex(spc),parameter::  z230=-z050
complex(spc),parameter::  z231=-z051
complex(spc),parameter::  z234=-z054
complex(spc),parameter::  z236=-z056
complex(spc),parameter::  z240=-z060
complex(spc),parameter::  z244=-z064
complex(spc),parameter::  z248=-z068
complex(spc),parameter::  z250=-z070
complex(spc),parameter::  z252=-z072
complex(spc),parameter::  z255=-z075
complex(spc),parameter::  z257=-z077
complex(spc),parameter::  z259=-z079
complex(spc),parameter::  z260=-z080
complex(spc),parameter::  z270=-z090
complex(spc),parameter::  z280=-z100
complex(spc),parameter::  z281=-z101
complex(spc),parameter::  z283=-z103
complex(spc),parameter::  z285=-z105
complex(spc),parameter::  z288=-z108
complex(spc),parameter::  z290=-z110
complex(spc),parameter::  z292=-z112
complex(spc),parameter::  z296=-z116
complex(spc),parameter::  z300=-z120
complex(spc),parameter::  z304=-z124
complex(spc),parameter::  z306=-z126
complex(spc),parameter::  z309=-z129
complex(spc),parameter::  z310=-z130
complex(spc),parameter::  z315=-z135
complex(spc),parameter::  z320=-z140
complex(spc),parameter::  z321=-z141
complex(spc),parameter::  z324=-z144
complex(spc),parameter::  z326=-z146
complex(spc),parameter::  z330=-z150
complex(spc),parameter::  z334=-z154
complex(spc),parameter::  z338=-z158
complex(spc),parameter::  z340=-z160
complex(spc),parameter::  z342=-z162
complex(spc),parameter::  z345=-z165
complex(spc),parameter::  z347=-z167
complex(spc),parameter::  z349=-z169
complex(spc),parameter::  z350=-z170
end module pietc_s