Lines Matching refs:j
60 integer i, j local
99 do j = 1, m
100 c(i,j) = a(i,j) * b(i,j)
116 integer m, n, i, j local
120 do j = 1, m
121 a(i,j) = 0.0e-0
135 integer m, n, i, j local
139 do j = 1, m
140 b(i,j) = a(i,j)
154 integer m, n, i, j local
158 do j = 1, m
159 a(i,j) = a(i,j) * s
173 integer m, n, i, j local
177 do j = 1, m
178 b(i,j) = a(i,j) * s
192 integer m, n, i, j local
196 do j = 1, m
197 b(i,j) = b(i,j) + a(i,j)
211 integer m, n, i, j local
215 do j = 1, m
216 b(i,j) = b(i,j) - a(i,j)
230 integer m, n, i, j local
235 do j = 1, m
236 flesDot1 = flesDot1 + a(i,j) * a(i,j)
250 integer m, n, i, j local
255 do j = 1, m
256 flesDot2 = flesDot2 + a(i,j) * b(i,j)
270 integer m, n, i, j local
275 do j= 1, m
276 y(i,j) = y(i,j) + a * x(i,j)
290 integer m, n, i, j local
295 do j = 1, m
296 y(i,j) = a * y(i,j) + x(i,j)
310 integer m, n, i, j local
314 do j = 1, m
315 if ( x(i,j) .ne. 0 ) x(i,j) = 1. / x(i,j)
338 integer i, j, m1 local
472 do j = m1, m, 8
482 tmp1 = tmp1 + x(i,j+0) * y(i)
483 tmp2 = tmp2 + x(i,j+1) * y(i)
484 tmp3 = tmp3 + x(i,j+2) * y(i)
485 tmp4 = tmp4 + x(i,j+3) * y(i)
486 tmp5 = tmp5 + x(i,j+4) * y(i)
487 tmp6 = tmp6 + x(i,j+5) * y(i)
488 tmp7 = tmp7 + x(i,j+6) * y(i)
489 tmp8 = tmp8 + x(i,j+7) * y(i)
491 c(j+0) = tmp1
492 c(j+1) = tmp2
493 c(j+2) = tmp3
494 c(j+3) = tmp4
495 c(j+4) = tmp5
496 c(j+5) = tmp6
497 c(j+6) = tmp7
498 c(j+7) = tmp8
520 integer i, j, m1 local
621 do j = m1, m, 8
622 tmp1 = c(j+0)
623 tmp2 = c(j+1)
624 tmp3 = c(j+2)
625 tmp4 = c(j+3)
626 tmp5 = c(j+4)
627 tmp6 = c(j+5)
628 tmp7 = c(j+6)
629 tmp8 = c(j+7)
632 1 + tmp1 * x(i,j+0) + tmp2 * x(i,j+1)
633 2 + tmp3 * x(i,j+2) + tmp4 * x(i,j+3)
634 3 + tmp5 * x(i,j+4) + tmp6 * x(i,j+5)
635 4 + tmp7 * x(i,j+6) + tmp8 * x(i,j+7)
658 integer i, j, m1 local
764 do j = m1, m, 8
765 tmp1 = c(j+0)
766 tmp2 = c(j+1)
767 tmp3 = c(j+2)
768 tmp4 = c(j+3)
769 tmp5 = c(j+4)
770 tmp6 = c(j+5)
771 tmp7 = c(j+6)
772 tmp8 = c(j+7)
775 1 + tmp1 * x(i,j+0) + tmp2 * x(i,j+1)
776 2 + tmp3 * x(i,j+2) + tmp4 * x(i,j+3)
777 3 + tmp5 * x(i,j+4) + tmp6 * x(i,j+5)
778 4 + tmp7 * x(i,j+6) + tmp8 * x(i,j+7)
801 integer i, j, m1 local
901 do j = m1, m, 8
902 tmp1 = c(j+0)
903 tmp2 = c(j+1)
904 tmp3 = c(j+2)
905 tmp4 = c(j+3)
906 tmp5 = c(j+4)
907 tmp6 = c(j+5)
908 tmp7 = c(j+6)
909 tmp8 = c(j+7)
912 1 - tmp1 * x(i,j+0) - tmp2 * x(i,j+1)
913 2 - tmp3 * x(i,j+2) - tmp4 * x(i,j+3)
914 3 - tmp5 * x(i,j+4) - tmp6 * x(i,j+5)
915 4 - tmp7 * x(i,j+6) - tmp8 * x(i,j+7)
936 integer i, j local
973 do j = 1, m
974 b(i,j) = a(i,j)
997 integer i, j local
1043 do j = 1, m
1044 c(j) = 0
1046 c(j) = c(j) + a(i,j) * b(i,j)
1069 integer i, j local
1105 do j = 1, m
1107 b(i,j) = b(i,j) + c(j) * a(i,j)
1141 do j = 1, nshl
1144 & + xGoC(1:npro,i0+1,j) * Ptemp(:,j,nPs)
1147 & + xGoC(1:npro,i0+2,j) * Ptemp(:,j,nPs)
1150 & + xGoC(1:npro,i0+3,j) * Ptemp(:,j,nPs)
1190 do j = 1, nshl
1191 j0 = (nsd) * (j - 1)
1194 & + xKebe(1:npro,i0+1,j0+1) * Ptemp(:,j,1)
1195 & + xKebe(1:npro,i0+1,j0+2) * Ptemp(:,j,2)
1196 & + xKebe(1:npro,i0+1,j0+3) * Ptemp(:,j,3)
1199 & + xKebe(1:npro,i0+2,j0+1) * Ptemp(:,j,1)
1200 & + xKebe(1:npro,i0+2,j0+2) * Ptemp(:,j,2)
1201 & + xKebe(1:npro,i0+2,j0+3) * Ptemp(:,j,3)
1203 & + xKebe(1:npro,i0+3,j0+1) * Ptemp(:,j,1)
1204 & + xKebe(1:npro,i0+3,j0+2) * Ptemp(:,j,2)
1205 & + xKebe(1:npro,i0+3,j0+3) * Ptemp(:,j,3)
1214 do j = 1, nshl
1217 & + xGoC(1:npro,i0+1,j) * Ptemp(:,j,nPs)
1219 & + xGoC(1:npro,i0+2,j) * Ptemp(:,j,nPs)
1221 & + xGoC(1:npro,i0+3,j) * Ptemp(:,j,nPs)
1258 do j = 1, nshl
1259 i0 = (nsd) * (j - 1)
1262 & - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1263 & - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1264 & - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1301 do j = 1, nshl
1302 i0 = (nsd) * (j - 1)
1305 & - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1306 & - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1307 & - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1318 do j = 1, nshl
1321 & + xGoC(1:npro,i0,j) * Ptemp(:,j,nPs)
1361 do j = 1, nshl
1362 j0 = (nsd) * (j - 1)
1365 & + xKebe(1:npro,i0+1,j0+1) * Ptemp(:,j,1)
1366 & + xKebe(1:npro,i0+1,j0+2) * Ptemp(:,j,2)
1367 & + xKebe(1:npro,i0+1,j0+3) * Ptemp(:,j,3)
1370 & + xKebe(1:npro,i0+2,j0+1) * Ptemp(:,j,1)
1371 & + xKebe(1:npro,i0+2,j0+2) * Ptemp(:,j,2)
1372 & + xKebe(1:npro,i0+2,j0+3) * Ptemp(:,j,3)
1374 & + xKebe(1:npro,i0+3,j0+1) * Ptemp(:,j,1)
1375 & + xKebe(1:npro,i0+3,j0+2) * Ptemp(:,j,2)
1376 & + xKebe(1:npro,i0+3,j0+3) * Ptemp(:,j,3)
1385 do j = 1, nshl
1388 & + xGoC(1:npro,i0+1,j) * Ptemp(:,j,nPs)
1390 & + xGoC(1:npro,i0+2,j) * Ptemp(:,j,nPs)
1392 & + xGoC(1:npro,i0+3,j) * Ptemp(:,j,nPs)
1400 do j = 1, nshl
1401 i0 = (nsd) * (j - 1)
1404 & - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1405 & - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1406 & - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1417 do j = 1, nshl
1420 & + xGoC(1:npro,i0,j) * Ptemp(:,j,nPs)
1480 do j = 1, nshl
1483 & + xTe(1:npro,i,j) * Ptemp(:,j,nPs)