1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
c
c-----------------------------------------------------------------------
c subroutine: const
c computes the multipliers for the various modules
c-----------------------------------------------------------------------
c
subroutine const(co3,co8,co16,co9,cdc,cdd)
double precision dtheta,dtwopi,dsq32,dsq2
double precision dcos1,dcos2,dcos3,dcos4
double precision dsin1,dsin2,dsin3,dsin4
dimension co3(3),co8(8),co16(18),co9(11),cdc(9),cdd(6)
dtwopi=8.0d0*datan(1.0d0)
dsq32=dsqrt(0.75d0)
dsq2=dsqrt(0.5d0)
c
c multipliers for the three point module
c
co3(1)=1.0
co3(2)=-1.5
co3(3)=-dsq32
c
c multipliers for the five point module
c
dtheta=dtwopi/5.0d0
dcos1=dcos(dtheta)
dcos2=dcos(2.0d0*dtheta)
dsin1=dsin(dtheta)
dsin2=dsin(2.0d0*dtheta)
cdd(1)=1.0
cdd(2)=-1.25
cdd(3)=-dsin1-dsin2
cdd(4)=0.5*(dcos1-dcos2)
cdd(5)=dsin1-dsin2
cdd(6)=dsin2
c
c
c multipliers for the seven point module
c
dtheta=dtwopi/7.0d0
dcos1=dcos(dtheta)
dcos2=dcos(2.0d0*dtheta)
dcos3=dcos(3.0d0*dtheta)
dsin1=dsin(dtheta)
dsin2=dsin(2.0d0*dtheta)
dsin3=dsin(3.0d0*dtheta)
cdc(1)=1.0
cdc(2)=-7.0d0/6.0d0
cdc(3)=-(dsin1+dsin2-dsin3)/3.0d0
cdc(4)=(dcos1+dcos2-2.0d0*dcos3)/3.0d0
cdc(5)=(2.0d0*dcos1-dcos2-dcos3)/3.0d0
cdc(6)=-(2.0d0*dsin1-dsin2+dsin3)/3.0d0
cdc(7)=-(dsin1+dsin2+2.0d0*dsin3)/3.0d0
cdc(8)=(dcos1-2.0d0*dcos2+dcos3)/3.0d0
cdc(9)=-(dsin1-2.0d0*dsin2-dsin3)/3.0d0
c
c multipliers for the eight point module
c
co8(1)=1.0
co8(2)=1.0
co8(3)=1.0
co8(4)=-1.0
co8(5)=1.0
co8(6)=-dsq2
co8(7)=-1.0
co8(8)=dsq2
c
c multipliers for the nine point module
c
dtheta=dtwopi/9.0d0
dcos1=dcos(dtheta)
dcos2=dcos(2.0d0*dtheta)
dcos4=dcos(4.0d0*dtheta)
dsin1=dsin(dtheta)
dsin2=dsin(2.0d0*dtheta)
dsin4=dsin(4.0d0*dtheta)
co9(1)=1.0
co9(2)=-1.5
co9(3)=-dsq32
co9(4)=0.5
co9(5)=(2.0d0*dcos1-dcos2-dcos4)/3.0d0
co9(6)=(dcos1-2.0d0*dcos2+dcos4)/3.0d0
co9(7)=(dcos1+dcos2-2.0d0*dcos4)/3.0d0
co9(8)=-(2.0d0*dsin1+dsin2-dsin4)/3.0d0
co9(9)=-(dsin1+2.0d0*dsin2+dsin4)/3.0d0
co9(10)=-(dsin1-dsin2-2.0d0*dsin4)/3.0d0
co9(11)=-dsq32
c
c multipliers for the sixteen point module
c
dtheta=dtwopi/16.0d0
dcos1=dcos(dtheta)
dcos3=dcos(3.0d0*dtheta)
dsin1=dsin(dtheta)
dsin3=dsin(3.0d0*dtheta)
co16(1)=1.0
co16(2)=1.0
co16(3)=1.0
co16(4)=-1.0
co16(5)=1.0
co16(6)=-dsq2
co16(7)=-1.0
co16(8)=dsq2
co16(9)=1.0
co16(10)=-(dsin1-dsin3)
co16(11)=-dsq2
co16(12)=-co16(10)
co16(13)=-1.0
co16(14)=-(dsin1+dsin3)
co16(15)=dsq2
co16(16)=-co16(14)
co16(17)=-dsin3
co16(18)=dcos3
return
end
|
