Lines Matching refs:intp
32 do intp = 1, ngaussf ! Loop over quadrature points
44 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
45 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
46 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
47 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
48 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
49 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
50 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
51 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
52 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
83 wght=Qwtf(intp)
91 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
92 & + shgl(2,n,intp) * dxidx(:,2,1)
93 & + shgl(3,n,intp) * dxidx(:,3,1))
94 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
95 & + shgl(2,n,intp) * dxidx(:,2,2)
96 & + shgl(3,n,intp) * dxidx(:,3,2))
97 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
98 & + shgl(2,n,intp) * dxidx(:,2,3)
99 & + shgl(3,n,intp) * dxidx(:,3,3))
116 u1 = u1 + shp(i,intp)*yl(:,i,2)
117 u2 = u2 + shp(i,intp)*yl(:,i,3)
118 u3 = u3 + shp(i,intp)*yl(:,i,4)
222 do intp = 1, ngauss
235 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
236 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
237 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
238 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
239 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
240 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
241 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
242 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
243 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
273 wght=Qwt(lcsyst,intp) ! may be different now
279 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
280 & + shgl(2,n,intp) * dxidx(:,2,1)
281 & + shgl(3,n,intp) * dxidx(:,3,1))
282 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
283 & + shgl(2,n,intp) * dxidx(:,2,2)
284 & + shgl(3,n,intp) * dxidx(:,3,2))
285 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
286 & + shgl(2,n,intp) * dxidx(:,2,3)
287 & + shgl(3,n,intp) * dxidx(:,3,3))
304 fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1) ! hat{u1}
305 fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2) ! hat{u2}
306 fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3) ! hat{u3}
308 fresli(:,4) = fresli(:,4) + shp(i,intp)*hfresl(:,i,1)*
310 fresli(:,5) = fresli(:,5) + shp(i,intp)*hfresl(:,i,2)*
312 fresli(:,6) = fresli(:,6) + shp(i,intp)*hfresl(:,i,3)*
314 fresli(:,7) = fresli(:,7) + shp(i,intp)*hfresl(:,i,1)*
316 fresli(:,8) = fresli(:,8) + shp(i,intp)*hfresl(:,i,1)*
318 fresli(:,9) = fresli(:,9) + shp(i,intp)*hfresl(:,i,2)*
321 fresli(:,10) = fresli(:,10) + shp(i,intp)*hfresl(:,i,10) ! hat{S11}
322 fresli(:,11) = fresli(:,11) + shp(i,intp)*hfresl(:,i,11) ! hat{S22}
323 fresli(:,12) = fresli(:,12) + shp(i,intp)*hfresl(:,i,12) ! hat{S33}
324 fresli(:,13) = fresli(:,13) + shp(i,intp)*hfresl(:,i,13) ! hat{S12}
325 fresli(:,14) = fresli(:,14) + shp(i,intp)*hfresl(:,i,14) ! hat{S13}
326 fresli(:,15) = fresli(:,15) + shp(i,intp)*hfresl(:,i,15) ! hat{S23}
360 fresli(:,16) =fresli(:,16)+shp(i,intp)*SS11(:,i)! |hat{S}|*hat{S11}
361 fresli(:,17) =fresli(:,17)+shp(i,intp)*SS22(:,i)! |hat{S}|*hat{S22}
362 fresli(:,18) =fresli(:,18)+shp(i,intp)*SS33(:,i)! |hat{S}|*hat{S33}
363 fresli(:,19) =fresli(:,19)+shp(i,intp)*SS12(:,i)! |hat{S}|*hat{S12}
364 fresli(:,20) =fresli(:,20)+shp(i,intp)*SS13(:,i)! |hat{S}|*hat{S13}
365 fresli(:,21) =fresli(:,21)+shp(i,intp)*SS23(:,i)! |hat{S}|*hat{S23}
432 do intp = 1, ngaussf ! Loop over quadrature points
444 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
445 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
446 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
447 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
448 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
449 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
450 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
451 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
452 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
483 wght=Qwtf(intp)
491 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
492 & + shgl(2,n,intp) * dxidx(:,2,1)
493 & + shgl(3,n,intp) * dxidx(:,3,1))
494 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
495 & + shgl(2,n,intp) * dxidx(:,2,2)
496 & + shgl(3,n,intp) * dxidx(:,3,2))
497 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
498 & + shgl(2,n,intp) * dxidx(:,2,3)
499 & + shgl(3,n,intp) * dxidx(:,3,3))
516 u1 = u1 + shp(i,intp)*yl(:,i,2)
517 u2 = u2 + shp(i,intp)*yl(:,i,3)
518 u3 = u3 + shp(i,intp)*yl(:,i,4)
635 do intp = 1, ngaussf
648 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
649 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
650 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
651 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
652 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
653 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
654 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
655 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
656 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
687 wght=Qwtf(intp)
692 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
693 & + shgl(2,n,intp) * dxidx(:,2,1)
694 & + shgl(3,n,intp) * dxidx(:,3,1))
695 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
696 & + shgl(2,n,intp) * dxidx(:,2,2)
697 & + shgl(3,n,intp) * dxidx(:,3,2))
698 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
699 & + shgl(2,n,intp) * dxidx(:,2,3)
700 & + shgl(3,n,intp) * dxidx(:,3,3))
717 fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1) ! hat{u1}
718 fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2) ! hat{u2}
719 fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3) ! hat{u3}
721 fresli(:,4) = fresli(:,4) + shp(i,intp)*hfresl(:,i,4) ! hat{u1*u1}
722 fresli(:,5) = fresli(:,5) + shp(i,intp)*hfresl(:,i,5) ! hat{u2*u2}
723 fresli(:,6) = fresli(:,6) + shp(i,intp)*hfresl(:,i,6) ! hat{u3*u3}
724 fresli(:,7) = fresli(:,7) + shp(i,intp)*hfresl(:,i,7) ! hat{u1*u2}
725 fresli(:,8) = fresli(:,8) + shp(i,intp)*hfresl(:,i,8) ! hat{u1*u3}
726 fresli(:,9) = fresli(:,9) + shp(i,intp)*hfresl(:,i,9) ! hat{u2*u3}
728 fresli(:,10) = fresli(:,10) + shp(i,intp)*hfresl(:,i,10) ! hat{S11}
729 fresli(:,11) = fresli(:,11) + shp(i,intp)*hfresl(:,i,11) ! hat{S22}
730 fresli(:,12) = fresli(:,12) + shp(i,intp)*hfresl(:,i,12) ! hat{S33}
731 fresli(:,13) = fresli(:,13) + shp(i,intp)*hfresl(:,i,13) ! hat{S12}
732 fresli(:,14) = fresli(:,14) + shp(i,intp)*hfresl(:,i,14) ! hat{S13}
733 fresli(:,15) = fresli(:,15) + shp(i,intp)*hfresl(:,i,15) ! hat{S23}
735 fresli(:,16) = fresli(:,16) +shp(i,intp)*hfresl(:,i,17)! hat{S11*|S|}
736 fresli(:,17) = fresli(:,17) +shp(i,intp)*hfresl(:,i,18)! hat{S22*|S|}
737 fresli(:,18) = fresli(:,18) +shp(i,intp)*hfresl(:,i,19)! hat{S33*|S|}
738 fresli(:,19) = fresli(:,19) +shp(i,intp)*hfresl(:,i,20)! hat{S12*|S|}
739 fresli(:,20) = fresli(:,20) +shp(i,intp)*hfresl(:,i,21)! hat{S13*|S|}
740 fresli(:,21) = fresli(:,21) +shp(i,intp)*hfresl(:,i,22)! hat{S23*|S|}
811 do intp = 1, ngaussf ! Loop over quadrature points
823 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
824 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
825 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
826 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
827 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
828 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
829 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
830 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
831 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
862 wght=Qwtf(intp)
870 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
871 & + shgl(2,n,intp) * dxidx(:,2,1)
872 & + shgl(3,n,intp) * dxidx(:,3,1))
873 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
874 & + shgl(2,n,intp) * dxidx(:,2,2)
875 & + shgl(3,n,intp) * dxidx(:,3,2))
876 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
877 & + shgl(2,n,intp) * dxidx(:,2,3)
878 & + shgl(3,n,intp) * dxidx(:,3,3))
895 u1 = u1 + shp(i,intp)*yl(:,i,2)
896 u2 = u2 + shp(i,intp)*yl(:,i,3)
897 u3 = u3 + shp(i,intp)*yl(:,i,4)
926 fresli(:,i,1) = WdetJ * u1 * shp(i,intp) !G * u1 * WdetJ
927 fresli(:,i,2) = WdetJ * u2 * shp(i,intp) !G * u2 * WdetJ
928 fresli(:,i,3) = WdetJ * u3 * shp(i,intp) !G * u3 * WdetJ
930 fresli(:,i,4) = WdetJ * u1 * u1 * shp(i,intp) ! G*u1*u1*WdetJ
931 fresli(:,i,5) = WdetJ * u2 * u2 * shp(i,intp) ! G*u2*u2*WdetJ
932 fresli(:,i,6) = WdetJ * u3 * u3 * shp(i,intp) ! G*u3*u3*WdetJ
933 fresli(:,i,7) = WdetJ * u1 * u2 * shp(i,intp) ! G*u1*u2*WdetJ
934 fresli(:,i,8) = WdetJ * u1 * u3 * shp(i,intp) ! G*u1*u3*WdetJ
935 fresli(:,i,9) = WdetJ * u2 * u3 * shp(i,intp) ! G*u2*u3*WdetJ
937 fresli(:,i,10) = S11 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
938 fresli(:,i,11) = S22 * shp(i,intp) * WdetJ ! G*S_{2,2}*WdetJ
939 fresli(:,i,12) = S33 * shp(i,intp) * WdetJ ! G*S_{3,3}*WdetJ
940 fresli(:,i,13) = S12 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
941 fresli(:,i,14) = S13 * shp(i,intp) * WdetJ ! G*S_{1,3}*WdetJ
942 fresli(:,i,15) = S23 * shp(i,intp) * WdetJ ! G*S_{2,3}*WdetJ
944 fresli(:,i,22) = WdetJ*shp(i,intp) !Integral of filter kernel, G,
949 fresli(:,i,16) = S11 * strnrmi * shp(i,intp) * WdetJ
950 fresli(:,i,17) = S22 * strnrmi * shp(i,intp) * WdetJ
951 fresli(:,i,18) = S33 * strnrmi * shp(i,intp) * WdetJ
952 fresli(:,i,19) = S12 * strnrmi * shp(i,intp) * WdetJ
953 fresli(:,i,20) = S13 * strnrmi * shp(i,intp) * WdetJ
954 fresli(:,i,21) = S23 * strnrmi * shp(i,intp) * WdetJ
1265 do intp = 1, ngaussf ! Loop over quadrature points
1277 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
1278 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
1279 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
1280 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
1281 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
1282 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
1283 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
1284 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
1285 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
1316 wght=Qwtf(intp)
1324 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
1325 & + shgl(2,n,intp) * dxidx(:,2,1)
1326 & + shgl(3,n,intp) * dxidx(:,3,1))
1327 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
1328 & + shgl(2,n,intp) * dxidx(:,2,2)
1329 & + shgl(3,n,intp) * dxidx(:,3,2))
1330 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
1331 & + shgl(2,n,intp) * dxidx(:,2,3)
1332 & + shgl(3,n,intp) * dxidx(:,3,3))
1349 u1 = u1 + shp(i,intp)*yl(:,i,2)
1350 u2 = u2 + shp(i,intp)*yl(:,i,3)
1351 u3 = u3 + shp(i,intp)*yl(:,i,4)
1380 fresli(:,i,1) = WdetJ * u1 * shp(i,intp) !G * u1 * WdetJ
1381 fresli(:,i,2) = WdetJ * u2 * shp(i,intp) !G * u2 * WdetJ
1382 fresli(:,i,3) = WdetJ * u3 * shp(i,intp) !G * u3 * WdetJ
1384 fresli(:,i,4) = WdetJ * u1 * u1 * shp(i,intp) ! G*u1*u1*WdetJ
1385 fresli(:,i,5) = WdetJ * u2 * u2 * shp(i,intp) ! G*u2*u2*WdetJ
1386 fresli(:,i,6) = WdetJ * u3 * u3 * shp(i,intp) ! G*u3*u3*WdetJ
1387 fresli(:,i,7) = WdetJ * u1 * u2 * shp(i,intp) ! G*u1*u2*WdetJ
1388 fresli(:,i,8) = WdetJ * u1 * u3 * shp(i,intp) ! G*u1*u3*WdetJ
1389 fresli(:,i,9) = WdetJ * u2 * u3 * shp(i,intp) ! G*u2*u3*WdetJ
1391 fresli(:,i,10) = S11 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
1392 fresli(:,i,11) = S22 * shp(i,intp) * WdetJ ! G*S_{2,2}*WdetJ
1393 fresli(:,i,12) = S33 * shp(i,intp) * WdetJ ! G*S_{3,3}*WdetJ
1394 fresli(:,i,13) = S12 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
1395 fresli(:,i,14) = S13 * shp(i,intp) * WdetJ ! G*S_{1,3}*WdetJ
1396 fresli(:,i,15) = S23 * shp(i,intp) * WdetJ ! G*S_{2,3}*WdetJ
1399 fresli(:,i,16) = strnrmi*strnrmi*strnrmi*shp(i,intp)*WdetJ
1401 fresli(:,i,22) = WdetJ*shp(i,intp) !Integral of filter kernel, G,
1407 fresli(:,i,17) = S22 * strnrmi * shp(i,intp) * WdetJ
1408 fresli(:,i,18) = S33 * strnrmi * shp(i,intp) * WdetJ
1409 fresli(:,i,19) = S12 * strnrmi * shp(i,intp) * WdetJ
1410 fresli(:,i,20) = S13 * strnrmi * shp(i,intp) * WdetJ
1411 fresli(:,i,21) = S23 * strnrmi * shp(i,intp) * WdetJ