Lines Matching refs:ay
21 DO ay=0,la-ax
22 sc(coset(ax,ay,az),1) = rap(1)*sc(coset(ax-1,ay,az),1) +&
23 f2 * kvec(1)*ss(coset(ax-1,ay,az),1)
24 ss(coset(ax,ay,az),1) = rap(1)*ss(coset(ax-1,ay,az),1) +&
25 f2 * kvec(1)*sc(coset(ax-1,ay,az),1)
94 DO ay=2,la
95 s(coset(0,ay,az),1,l) = rap(2)*s(coset(0,ay-1,az),1,l) +&
96 f2*REAL(ay-1,dp)*s(coset(0,ay-2,az),1,l)
97 IF ( ly1 > 0 ) s(coset(0,ay,az),1,l) = s(coset(0,ay,az),1,l) + &
98 f2y*s(coset(0,ay-1,az),1,ly1)
100 DO ay=0,la-1
101 IF ( lx1 > 0 ) s(coset(1,ay,az),1,l) = s(coset(1,ay,az),1,l) + &
102 f2x*s(coset(0,ay,az),1,lx1)
105 DO ay=0,la-ax
106 s(coset(ax,ay,az),1,l) = rap(1)*s(coset(ax-1,ay,az),1,l) +&
107 f3*s(coset(ax-2,ay,az),1,l)
108 IF ( lx1 > 0 ) s(coset(ax,ay,az),1,l) = s(coset(ax,ay,az),1,l) + &
109 f2x*s(coset(ax-1,ay,az),1,lx1)
121 DO ay=0,la-ax
122 s(coset(ax,ay,az),2,l) = s(coset(ax+1,ay,az),1,l) -&
123 rab(1)*s(coset(ax,ay,az),1,l)
124 s(coset(ax,ay,az),4,l) = s(coset(ax,ay,az+1),1,l) -&
125 rab(3)*s(coset(ax,ay,az),1,l)
130 DO ay=0,la_max-ax
132 s(coset(ax,ay,az),2,l) = rbp(1)*s(coset(ax,ay,az),1,l)
134 s(coset(ax,ay,az),2,l) = rbp(1)*s(coset(ax,ay,az),1,l) +&
135 fx*s(coset(ax-1,ay,az),1,l)
137 IF (lx1 > 0) s(coset(ax,ay,az),2,l) = s(coset(ax,ay,az),2,l) +&
138 f2x*s(coset(ax,ay,az),1,lx1)
139 IF (ay == 0) THEN
140 s(coset(ax,ay,az),3,l) = rbp(2)*s(coset(ax,ay,az),1,l)
142 s(coset(ax,ay,az),3,l) = rbp(2)*s(coset(ax,ay,az),1,l) +&
143 fy*s(coset(ax,ay-1,az),1,l)
145 IF (ly1 > 0) s(coset(ax,ay,az),3,l) = s(coset(ax,ay,az),3,l) +&
146 f2y*s(coset(ax,ay,az),1,ly1)
148 s(coset(ax,ay,az),4,l) = rbp(3)*s(coset(ax,ay,az),1,l)
150 s(coset(ax,ay,az),4,l) = rbp(3)*s(coset(ax,ay,az),1,l) +&
151 fz*s(coset(ax,ay,az-1),1,l)
153 IF (lz1 > 0) s(coset(ax,ay,az),4,l) = s(coset(ax,ay,az),4,l) +&
154 f2z*s(coset(ax,ay,az),1,lz1)
160 DO ay=0,la-ax
161 s(coset(ax,ay,az),coset(0,0,lb),l) =&
162 rab(3)*s(coset(ax,ay,az),coset(0,0,lb-1),l)
165 s(coset(ax,ay,az),coset(bx,by,bz),l) =&
166 rab(1)*s(coset(ax,ay,az),coset(bx-1,by,bz),l)
173 DO ay=0,la_max-ax
175 s(coset(ax,ay,az),coset(0,0,lb),l) =&
176 rbp(3)*s(coset(ax,ay,az),coset(0,0,lb-1),l) +&
177 f3*s(coset(ax,ay,az),coset(0,0,lb-2),l)
179 IF (lz1 > 0) s(coset(ax,ay,az),coset(0,0,lb),l) =&
180 f2z*s(coset(ax,ay,az),coset(0,0,lb-1),lz1)
181 IF (ay == 0) THEN
182 IF (ly1 > 0) s(coset(ax,ay,az),coset(0,1,bz),l) =&
183 f2y*s(coset(ax,ay,az),coset(0,0,bz),ly1)
185 s(coset(ax,ay,az),coset(0,by,bz),l) =&
186 f3*s(coset(ax,ay,az),coset(0,by-2,bz),l)
187 IF (ly1 > 0) s(coset(ax,ay,az),coset(0,by,bz),l) =&
188 f2y*s(coset(ax,ay,az),coset(0,by-1,bz),ly1)
190 s(coset(ax,ay,az),coset(0,1,bz),l) =&
191 fy*s(coset(ax,ay-1,az),coset(0,0,bz),l)
195 IF (lx1 > 0) s(coset(ax,ay,az),coset(1,by,bz),l) =&
196 f2x*s(coset(ax,ay,az),coset(0,by,bz),lx1)
200 s(coset(ax,ay,az),coset(bx,by,bz),l) =&
201 f3*s(coset(ax,ay,az),coset(bx-2,by,bz),l)
202 IF (lx1 > 0) s(coset(ax,ay,az),coset(bx,by,bz),l) =&
203 f2x*s(coset(ax,ay,az),coset(bx-1,by,bz),lx1)
207 IF (lx1 > 0) s(coset(ax,ay,az),coset(1,by,bz),l) =&
208 f2x*s(coset(ax,ay,az),coset(0,by,bz),lx1)
212 s(coset(ax,ay,az),coset(bx,by,bz),l) =&
213 f3*s(coset(ax,ay,az),coset(bx-2,by,bz),l)
214 IF (lx1 > 0) s(coset(ax,ay,az),coset(bx,by,bz),l) =&
215 f2x*s(coset(ax,ay,az),coset(bx-1,by,bz),lx1)