Lines Matching refs:j
38 do j = 1, nshl ! May be worth eliminating zeros in A(prim) matrices
39 tmp=shp(:,j)*WdetJ
40 BDiagl(:,j,1,1) = BDiagl(:,j,1,1)
42 & shg(:,j,1) * A1(:,1,1) +
43 & shg(:,j,2) * A2(:,1,1) +
44 & shg(:,j,3) * A3(:,1,1)
46 BDiagl(:,j,1,2) = BDiagl(:,j,1,2)
48 & shg(:,j,1) * A1(:,1,2)
52 BDiagl(:,j,1,3) = BDiagl(:,j,1,3)
55 & +shg(:,j,2) * A2(:,1,3)
58 BDiagl(:,j,1,4) = BDiagl(:,j,1,4)
62 & shg(:,j,3) * A3(:,1,4)
64 BDiagl(:,j,1,5) = BDiagl(:,j,1,5)
66 & shg(:,j,1) * A1(:,1,5) +
67 & shg(:,j,2) * A2(:,1,5) +
68 & shg(:,j,3) * A3(:,1,5)
70 BDiagl(:,j,2,1) = BDiagl(:,j,2,1)
72 & shg(:,j,1) * A1(:,2,1) +
73 & shg(:,j,2) * A2(:,2,1) +
74 & shg(:,j,3) * A3(:,2,1)
76 BDiagl(:,j,2,2) = BDiagl(:,j,2,2)
78 & shg(:,j,1) * A1(:,2,2) +
79 & shg(:,j,2) * A2(:,2,2) +
80 & shg(:,j,3) * A3(:,2,2)
82 BDiagl(:,j,2,3) = BDiagl(:,j,2,3)
85 & +shg(:,j,2) * A2(:,2,3)
88 BDiagl(:,j,2,4) = BDiagl(:,j,2,4)
92 & shg(:,j,3) * A3(:,2,4)
94 BDiagl(:,j,2,5) = BDiagl(:,j,2,5)
96 & shg(:,j,1) * A1(:,2,5) +
97 & shg(:,j,2) * A2(:,2,5) +
98 & shg(:,j,3) * A3(:,2,5)
100 BDiagl(:,j,3,1) = BDiagl(:,j,3,1)
102 & shg(:,j,1) * A1(:,3,1) +
103 & shg(:,j,2) * A2(:,3,1) +
104 & shg(:,j,3) * A3(:,3,1)
106 BDiagl(:,j,3,2) = BDiagl(:,j,3,2)
108 & shg(:,j,1) * A1(:,3,2)
112 BDiagl(:,j,3,3) = BDiagl(:,j,3,3)
114 & shg(:,j,1) * A1(:,3,3) +
115 & shg(:,j,2) * A2(:,3,3) +
116 & shg(:,j,3) * A3(:,3,3)
118 BDiagl(:,j,3,4) = BDiagl(:,j,3,4)
122 & shg(:,j,3) * A3(:,3,4)
124 BDiagl(:,j,3,5) = BDiagl(:,j,3,5)
126 & shg(:,j,1) * A1(:,3,5) +
127 & shg(:,j,2) * A2(:,3,5) +
128 & shg(:,j,3) * A3(:,3,5)
130 BDiagl(:,j,4,1) = BDiagl(:,j,4,1)
132 & shg(:,j,1) * A1(:,4,1) +
133 & shg(:,j,2) * A2(:,4,1) +
134 & shg(:,j,3) * A3(:,4,1)
136 BDiagl(:,j,4,2) = BDiagl(:,j,4,2)
138 & shg(:,j,1) * A1(:,4,2)
142 BDiagl(:,j,4,3) = BDiagl(:,j,4,3)
145 & +shg(:,j,2) * A2(:,4,3)
148 BDiagl(:,j,4,4) = BDiagl(:,j,4,4)
150 & shg(:,j,1) * A1(:,4,4) +
151 & shg(:,j,2) * A2(:,4,4) +
152 & shg(:,j,3) * A3(:,4,4)
154 BDiagl(:,j,4,5) = BDiagl(:,j,4,5)
156 & shg(:,j,1) * A1(:,4,5) +
157 & shg(:,j,2) * A2(:,4,5) +
158 & shg(:,j,3) * A3(:,4,5)
160 BDiagl(:,j,5,1) = BDiagl(:,j,5,1)
162 & shg(:,j,1) * A1(:,5,1) +
163 & shg(:,j,2) * A2(:,5,1) +
164 & shg(:,j,3) * A3(:,5,1)
166 BDiagl(:,j,5,2) = BDiagl(:,j,5,2)
168 & shg(:,j,1) * A1(:,5,2) +
169 & shg(:,j,2) * A2(:,5,2) +
170 & shg(:,j,3) * A3(:,5,2)
172 BDiagl(:,j,5,3) = BDiagl(:,j,5,3)
174 & shg(:,j,1) * A1(:,5,3) +
175 & shg(:,j,2) * A2(:,5,3) +
176 & shg(:,j,3) * A3(:,5,3)
178 BDiagl(:,j,5,4) = BDiagl(:,j,5,4)
180 & shg(:,j,1) * A1(:,5,4) +
181 & shg(:,j,2) * A2(:,5,4) +
182 & shg(:,j,3) * A3(:,5,4)
184 BDiagl(:,j,5,5) = BDiagl(:,j,5,5)
186 & shg(:,j,1) * A1(:,5,5) +
187 & shg(:,j,2) * A2(:,5,5) +
188 & shg(:,j,3) * A3(:,5,5)
200 do j=1,nshl ! take advantage of zeros in A0(Prim)
201 BDiagl(:,j,2,2) = BDiagl(:,j,2,2) + tmp * A0(:,2,2)
202 BDiagl(:,j,3,3) = BDiagl(:,j,3,3) + tmp * A0(:,3,3)
203 BDiagl(:,j,4,4) = BDiagl(:,j,4,4) + tmp * A0(:,4,4)
204 BDiagl(:,j,1,5) = BDiagl(:,j,1,5) + tmp * A0(:,1,5)
205 BDiagl(:,j,2,5) = BDiagl(:,j,2,5) + tmp * A0(:,2,5)
206 BDiagl(:,j,3,5) = BDiagl(:,j,3,5) + tmp * A0(:,3,5)
207 BDiagl(:,j,4,5) = BDiagl(:,j,4,5) + tmp * A0(:,4,5)
208 BDiagl(:,j,1,1) = BDiagl(:,j,1,1) + tmp * A0(:,1,1)
209 BDiagl(:,j,2,1) = BDiagl(:,j,2,1) + tmp * A0(:,2,1)
210 BDiagl(:,j,3,1) = BDiagl(:,j,3,1) + tmp * A0(:,3,1)
211 BDiagl(:,j,4,1) = BDiagl(:,j,4,1) + tmp * A0(:,4,1)
212 BDiagl(:,j,5,1) = BDiagl(:,j,5,1) + tmp * A0(:,5,1)
213 BDiagl(:,j,5,2) = BDiagl(:,j,5,2) + tmp * A0(:,5,2)
214 BDiagl(:,j,5,3) = BDiagl(:,j,5,3) + tmp * A0(:,5,3)
215 BDiagl(:,j,5,4) = BDiagl(:,j,5,4) + tmp * A0(:,5,4)
216 BDiagl(:,j,5,5) = BDiagl(:,j,5,5) + tmp * A0(:,5,5)
226 do j = 1, nshl
228 tmp2 = (shp(:,j)*shp(:,j)) * tmp
230 BDiagl(:,j,2,2) = BDiagl(:,j,2,2) + tmp2 * A0(:,2,2)
231 BDiagl(:,j,3,3) = BDiagl(:,j,3,3) + tmp2 * A0(:,3,3)
232 BDiagl(:,j,4,4) = BDiagl(:,j,4,4) + tmp2 * A0(:,4,4)
233 BDiagl(:,j,1,5) = BDiagl(:,j,1,5) + tmp2 * A0(:,1,5)
234 BDiagl(:,j,2,5) = BDiagl(:,j,2,5) + tmp2 * A0(:,2,5)
235 BDiagl(:,j,3,5) = BDiagl(:,j,3,5) + tmp2 * A0(:,3,5)
236 BDiagl(:,j,4,5) = BDiagl(:,j,4,5) + tmp2 * A0(:,4,5)
237 BDiagl(:,j,1,1) = BDiagl(:,j,1,1) + tmp2 * A0(:,1,1)
238 BDiagl(:,j,2,1) = BDiagl(:,j,2,1) + tmp2 * A0(:,2,1)
239 BDiagl(:,j,3,1) = BDiagl(:,j,3,1) + tmp2 * A0(:,3,1)
240 BDiagl(:,j,4,1) = BDiagl(:,j,4,1) + tmp2 * A0(:,4,1)
241 BDiagl(:,j,5,1) = BDiagl(:,j,5,1) + tmp2 * A0(:,5,1)
242 BDiagl(:,j,5,2) = BDiagl(:,j,5,2) + tmp2 * A0(:,5,2)
243 BDiagl(:,j,5,3) = BDiagl(:,j,5,3) + tmp2 * A0(:,5,3)
244 BDiagl(:,j,5,4) = BDiagl(:,j,5,4) + tmp2 * A0(:,5,4)
245 BDiagl(:,j,5,5) = BDiagl(:,j,5,5) + tmp2 * A0(:,5,5)
258 do j = 1, nflow
260 Atau1(:,i,j) = Atau1(:,i,j) + A1(:,i,k)*PTau(:,k,j)
267 do j = 1, nflow
269 A1tauA0(:,i,j) =
270 & Atau1(:,i,1)*A0(:,1,j) +
271 & Atau1(:,i,2)*A0(:,2,j) +
272 & Atau1(:,i,3)*A0(:,3,j) +
273 & Atau1(:,i,4)*A0(:,4,j) +
274 & Atau1(:,i,5)*A0(:,5,j)
283 do j = 1, nflow
285 Atau2(:,i,j) = Atau2(:,i,j) + A2(:,i,k)*PTau(:,k,j)
293 do j = 1, nflow
295 A2tauA0(:,i,j) =
296 & Atau2(:,i,1)*A0(:,1,j) +
297 & Atau2(:,i,2)*A0(:,2,j) +
298 & Atau2(:,i,3)*A0(:,3,j) +
299 & Atau2(:,i,4)*A0(:,4,j) +
300 & Atau2(:,i,5)*A0(:,5,j)
309 do j = 1, nflow
311 Atau3(:,i,j) = Atau3(:,i,j) + A3(:,i,k)*PTau(:,k,j)
320 do j = 1, nflow
322 A3tauA0(:,i,j) =
323 & Atau3(:,i,1)*A0(:,1,j) +
324 & Atau3(:,i,2)*A0(:,2,j) +
325 & Atau3(:,i,3)*A0(:,3,j) +
326 & Atau3(:,i,4)*A0(:,4,j) +
327 & Atau3(:,i,5)*A0(:,5,j)
574 do j = 1, nshl
575 tmp = WdetJ * shg(:,j,1) * shg(:,j,1)
578 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
791 do j = 1, nshl
792 tmp = WdetJ * shg(:,j,2) * shg(:,j,2)
795 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
1008 do j = 1, nshl
1009 tmp = WdetJ * shg(:,j,3) * shg(:,j,3)
1012 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
1348 do j = 1, nshl
1349 tmp = WdetJ * shg(:,j,1) * shg(:,j,2)
1352 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
1691 do j = 1, nshl
1692 tmp = WdetJ * shg(:,j,1) * shg(:,j,3)
1695 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
2034 do j = 1, nshl
2035 tmp = WdetJ * shg(:,j,2) * shg(:,j,3)
2038 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)