Lines Matching refs:j
37 do j = 1, nshl ! May be worth eliminating zeros in A(prim) matrices
38 tmp=shp(:,j)*WdetJ
39 BDiagl(:,j,1,1) = BDiagl(:,j,1,1)
41 & shg(:,j,1) * A1(:,1,1) +
42 & shg(:,j,2) * A2(:,1,1) +
43 & shg(:,j,3) * A3(:,1,1)
45 BDiagl(:,j,1,2) = BDiagl(:,j,1,2)
47 & shg(:,j,1) * A1(:,1,2)
51 BDiagl(:,j,1,3) = BDiagl(:,j,1,3)
54 & +shg(:,j,2) * A2(:,1,3)
57 BDiagl(:,j,1,4) = BDiagl(:,j,1,4)
61 & shg(:,j,3) * A3(:,1,4)
63 BDiagl(:,j,1,5) = BDiagl(:,j,1,5)
65 & shg(:,j,1) * A1(:,1,5) +
66 & shg(:,j,2) * A2(:,1,5) +
67 & shg(:,j,3) * A3(:,1,5)
69 BDiagl(:,j,2,1) = BDiagl(:,j,2,1)
71 & shg(:,j,1) * A1(:,2,1) +
72 & shg(:,j,2) * A2(:,2,1) +
73 & shg(:,j,3) * A3(:,2,1)
75 BDiagl(:,j,2,2) = BDiagl(:,j,2,2)
77 & shg(:,j,1) * A1(:,2,2) +
78 & shg(:,j,2) * A2(:,2,2) +
79 & shg(:,j,3) * A3(:,2,2)
81 BDiagl(:,j,2,3) = BDiagl(:,j,2,3)
84 & +shg(:,j,2) * A2(:,2,3)
87 BDiagl(:,j,2,4) = BDiagl(:,j,2,4)
91 & shg(:,j,3) * A3(:,2,4)
93 BDiagl(:,j,2,5) = BDiagl(:,j,2,5)
95 & shg(:,j,1) * A1(:,2,5) +
96 & shg(:,j,2) * A2(:,2,5) +
97 & shg(:,j,3) * A3(:,2,5)
99 BDiagl(:,j,3,1) = BDiagl(:,j,3,1)
101 & shg(:,j,1) * A1(:,3,1) +
102 & shg(:,j,2) * A2(:,3,1) +
103 & shg(:,j,3) * A3(:,3,1)
105 BDiagl(:,j,3,2) = BDiagl(:,j,3,2)
107 & shg(:,j,1) * A1(:,3,2)
111 BDiagl(:,j,3,3) = BDiagl(:,j,3,3)
113 & shg(:,j,1) * A1(:,3,3) +
114 & shg(:,j,2) * A2(:,3,3) +
115 & shg(:,j,3) * A3(:,3,3)
117 BDiagl(:,j,3,4) = BDiagl(:,j,3,4)
121 & shg(:,j,3) * A3(:,3,4)
123 BDiagl(:,j,3,5) = BDiagl(:,j,3,5)
125 & shg(:,j,1) * A1(:,3,5) +
126 & shg(:,j,2) * A2(:,3,5) +
127 & shg(:,j,3) * A3(:,3,5)
129 BDiagl(:,j,4,1) = BDiagl(:,j,4,1)
131 & shg(:,j,1) * A1(:,4,1) +
132 & shg(:,j,2) * A2(:,4,1) +
133 & shg(:,j,3) * A3(:,4,1)
135 BDiagl(:,j,4,2) = BDiagl(:,j,4,2)
137 & shg(:,j,1) * A1(:,4,2)
141 BDiagl(:,j,4,3) = BDiagl(:,j,4,3)
144 & +shg(:,j,2) * A2(:,4,3)
147 BDiagl(:,j,4,4) = BDiagl(:,j,4,4)
149 & shg(:,j,1) * A1(:,4,4) +
150 & shg(:,j,2) * A2(:,4,4) +
151 & shg(:,j,3) * A3(:,4,4)
153 BDiagl(:,j,4,5) = BDiagl(:,j,4,5)
155 & shg(:,j,1) * A1(:,4,5) +
156 & shg(:,j,2) * A2(:,4,5) +
157 & shg(:,j,3) * A3(:,4,5)
159 BDiagl(:,j,5,1) = BDiagl(:,j,5,1)
161 & shg(:,j,1) * A1(:,5,1) +
162 & shg(:,j,2) * A2(:,5,1) +
163 & shg(:,j,3) * A3(:,5,1)
165 BDiagl(:,j,5,2) = BDiagl(:,j,5,2)
167 & shg(:,j,1) * A1(:,5,2) +
168 & shg(:,j,2) * A2(:,5,2) +
169 & shg(:,j,3) * A3(:,5,2)
171 BDiagl(:,j,5,3) = BDiagl(:,j,5,3)
173 & shg(:,j,1) * A1(:,5,3) +
174 & shg(:,j,2) * A2(:,5,3) +
175 & shg(:,j,3) * A3(:,5,3)
177 BDiagl(:,j,5,4) = BDiagl(:,j,5,4)
179 & shg(:,j,1) * A1(:,5,4) +
180 & shg(:,j,2) * A2(:,5,4) +
181 & shg(:,j,3) * A3(:,5,4)
183 BDiagl(:,j,5,5) = BDiagl(:,j,5,5)
185 & shg(:,j,1) * A1(:,5,5) +
186 & shg(:,j,2) * A2(:,5,5) +
187 & shg(:,j,3) * A3(:,5,5)
199 do j=1,nshl ! take advantage of zeros in A0(Prim)
200 BDiagl(:,j,2,2) = BDiagl(:,j,2,2) + tmp * A0(:,2,2)
201 BDiagl(:,j,3,3) = BDiagl(:,j,3,3) + tmp * A0(:,3,3)
202 BDiagl(:,j,4,4) = BDiagl(:,j,4,4) + tmp * A0(:,4,4)
203 BDiagl(:,j,1,5) = BDiagl(:,j,1,5) + tmp * A0(:,1,5)
204 BDiagl(:,j,2,5) = BDiagl(:,j,2,5) + tmp * A0(:,2,5)
205 BDiagl(:,j,3,5) = BDiagl(:,j,3,5) + tmp * A0(:,3,5)
206 BDiagl(:,j,4,5) = BDiagl(:,j,4,5) + tmp * A0(:,4,5)
207 BDiagl(:,j,1,1) = BDiagl(:,j,1,1) + tmp * A0(:,1,1)
208 BDiagl(:,j,2,1) = BDiagl(:,j,2,1) + tmp * A0(:,2,1)
209 BDiagl(:,j,3,1) = BDiagl(:,j,3,1) + tmp * A0(:,3,1)
210 BDiagl(:,j,4,1) = BDiagl(:,j,4,1) + tmp * A0(:,4,1)
211 BDiagl(:,j,5,1) = BDiagl(:,j,5,1) + tmp * A0(:,5,1)
212 BDiagl(:,j,5,2) = BDiagl(:,j,5,2) + tmp * A0(:,5,2)
213 BDiagl(:,j,5,3) = BDiagl(:,j,5,3) + tmp * A0(:,5,3)
214 BDiagl(:,j,5,4) = BDiagl(:,j,5,4) + tmp * A0(:,5,4)
215 BDiagl(:,j,5,5) = BDiagl(:,j,5,5) + tmp * A0(:,5,5)
225 do j = 1, nshl
227 tmp2 = (shp(:,j)*shp(:,j)) * tmp
229 BDiagl(:,j,2,2) = BDiagl(:,j,2,2) + tmp2 * A0(:,2,2)
230 BDiagl(:,j,3,3) = BDiagl(:,j,3,3) + tmp2 * A0(:,3,3)
231 BDiagl(:,j,4,4) = BDiagl(:,j,4,4) + tmp2 * A0(:,4,4)
232 BDiagl(:,j,1,5) = BDiagl(:,j,1,5) + tmp2 * A0(:,1,5)
233 BDiagl(:,j,2,5) = BDiagl(:,j,2,5) + tmp2 * A0(:,2,5)
234 BDiagl(:,j,3,5) = BDiagl(:,j,3,5) + tmp2 * A0(:,3,5)
235 BDiagl(:,j,4,5) = BDiagl(:,j,4,5) + tmp2 * A0(:,4,5)
236 BDiagl(:,j,1,1) = BDiagl(:,j,1,1) + tmp2 * A0(:,1,1)
237 BDiagl(:,j,2,1) = BDiagl(:,j,2,1) + tmp2 * A0(:,2,1)
238 BDiagl(:,j,3,1) = BDiagl(:,j,3,1) + tmp2 * A0(:,3,1)
239 BDiagl(:,j,4,1) = BDiagl(:,j,4,1) + tmp2 * A0(:,4,1)
240 BDiagl(:,j,5,1) = BDiagl(:,j,5,1) + tmp2 * A0(:,5,1)
241 BDiagl(:,j,5,2) = BDiagl(:,j,5,2) + tmp2 * A0(:,5,2)
242 BDiagl(:,j,5,3) = BDiagl(:,j,5,3) + tmp2 * A0(:,5,3)
243 BDiagl(:,j,5,4) = BDiagl(:,j,5,4) + tmp2 * A0(:,5,4)
244 BDiagl(:,j,5,5) = BDiagl(:,j,5,5) + tmp2 * A0(:,5,5)
265 do j = 1, nflow
267 A1tauA0(:,i,j) =
268 & Atau(:,i,1)*A0(:,1,j) +
269 & Atau(:,i,2)*A0(:,2,j) +
270 & Atau(:,i,3)*A0(:,3,j) +
271 & Atau(:,i,4)*A0(:,4,j) +
272 & Atau(:,i,5)*A0(:,5,j)
289 do j = 1, nflow
291 A2tauA0(:,i,j) =
292 & Atau(:,i,1)*A0(:,1,j) +
293 & Atau(:,i,2)*A0(:,2,j) +
294 & Atau(:,i,3)*A0(:,3,j) +
295 & Atau(:,i,4)*A0(:,4,j) +
296 & Atau(:,i,5)*A0(:,5,j)
310 do j = 1, nflow
312 A3tauA0(:,i,j) =
313 & Atau(:,i,1)*A0(:,1,j) +
314 & Atau(:,i,2)*A0(:,2,j) +
315 & Atau(:,i,3)*A0(:,3,j) +
316 & Atau(:,i,4)*A0(:,4,j) +
317 & Atau(:,i,5)*A0(:,5,j)
611 do j = 1, nshl
612 tmp = WdetJ * shg(:,j,1) * shg(:,j,1)
615 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
872 do j = 1, nshl
873 tmp = WdetJ * shg(:,j,2) * shg(:,j,2)
876 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
1133 do j = 1, nshl
1134 tmp = WdetJ * shg(:,j,3) * shg(:,j,3)
1137 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
1616 do j = 1, nshl
1617 tmp = WdetJ * shg(:,j,1) * shg(:,j,2)
1620 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
2099 do j = 1, nshl
2100 tmp = WdetJ * shg(:,j,1) * shg(:,j,3)
2103 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
2582 do j = 1, nshl
2583 tmp = WdetJ * shg(:,j,2) * shg(:,j,3)
2586 BDiagl(:,j,i,k) = BDiagl(:,j,i,k) + tmp * A0(:,i,k)
2633 do j = 1, nshl ! May be worth eliminating zeros in A(prim) matrices
2634 tmp=shp(:,j)*WdetJ
2635 BDiagtl(:,j) = BDiagtl(:,j)
2637 & shg(:,j,1) * A1t(:) +
2638 & shg(:,j,2) * A2t(:) +
2639 & shg(:,j,3) * A3t(:)
2652 do j=1,nshl ! take advantage of zeros in A0t(Prim)
2653 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)
2663 do j = 1, nshl
2665 tmp2 = (shp(:,j)*shp(:,j)) * tmp
2667 BDiagtl(:,j) = BDiagtl(:,j) + tmp2 * A0t(:)
2759 do j = 1, nshl
2760 tmp = WdetJ * shg(:,j,1) * shg(:,j,1)
2762 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)
2781 do j = 1, nshl
2782 tmp = WdetJ * shg(:,j,2) * shg(:,j,2)
2784 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)
2804 do j = 1, nshl
2805 tmp = WdetJ * shg(:,j,3) * shg(:,j,3)
2807 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)
2823 do j = 1, nshl
2824 tmp = WdetJ * shg(:,j,1) * shg(:,j,2)
2826 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)
2838 do j = 1, nshl
2839 tmp = WdetJ * shg(:,j,1) * shg(:,j,3)
2841 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)
2854 do j = 1, nshl
2855 tmp = WdetJ * shg(:,j,2) * shg(:,j,3)
2857 BDiagtl(:,j) = BDiagtl(:,j) + tmp * A0t(:)