Lines Matching refs:i
235 do i = 1, npro
237 shpfun(i,n) = shp(n,intp)
611 i = iper(j)
612 if (i .ne. j) then
613 fres(i,:) = fres(i,:) + fres(j,:)
618 i = iper(j)
619 if (i .ne. j) then
659 do i = 1, 21
660 call sparseCG (fres(:,i), pfres(:,i), lhsG,
668 do i = 1, 21
669 fres(:,i) = pfres(:,i)
732 i = iper(j)
733 if (i .ne. j) then
781 do i = 1, numnp
782 xnuml = xnuml + xnum(i)
783 xdenl = xdenl + xden(i)
797 do i = 1, numnp
798 xnumt = xnumt + xnum(i)
799 xdent = xdent + xden(i)
1018 i = iper(j)
1019 if (i .ne. j) then
1020 hfres(i,:) = hfres(i,:) + hfres(j,:)
1024 i = iper(j)
1025 if (i .ne. j) then
1026 hfres(j,:) = hfres(i,:)
1074 i = iper(j)
1075 if (i .ne. j) then
1076 fres(i,:) = fres(i,:) + fres(j,:)
1081 i = iper(j)
1082 if (i .ne. j) then
1083 fres(j,:) = fres(i,:)
1158 i = iper(j)
1159 if (i .ne. j) then
1206 do i = 1,nshg
1207 xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
1208 xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
1235 do i = 1,nfath
1236 cdelsq1(i) = xnuder(i,1)
1244 do i = 1,nfath
1245 cdelsq1(i) = xnude(i,1)
1251 do i = 1, nfath
1252 write(22,*)i, xnuder(i,1)
1255 do i = 1, nfath
1256 write(22,*)i, xnude(i,1)
1262 do i = 1, nfath
1263 if (cdelsq1(i) .lt. zero) then
1264 cdelsq1(i) = zero
1470 i = iper(j)
1471 if (i .ne. j) then
1472 stress(i,:) = stress(i,:) + stress(j,:)
1473 diss(i,:) = diss(i,:) + diss(j,:)
1474 vol(i) = vol(i) + vol(j)
1479 i = iper(j)
1480 if (i .ne. j) then
1481 stress(j,:) = stress(i,:)
1482 diss(j,:) = diss(i,:)
1483 vol(j) = vol(i)
1492 do i = 1, 9
1493 stress(:,i) = stress(:,i)*vol(:)
1495 do i = 1, 3
1496 diss(:,i) = diss(:,i)*vol(:)
1501 do i = 1, 9
1502 xave(:,i) = stress(:,i)
1510 i = iper(j)
1511 if (i .ne. j) then
1544 do i = 1,nshg
1545 xaveg(ifath(i),:) = xaveg(ifath(i),:) + xave(i,:)
1586 do i = 1, nfath
1590 write(380,*)xavegt(i,10),xavegt(i,11),xavegt(i,12)
1722 i = iper(j)
1723 if (i .ne. j) then
1724 fres(i,:) = fres(i,:) + fres(j,:)
1728 i = iper(j)
1729 if (i .ne. j) then
1730 fres(j,:) = fres(i,:)
1906 i = iper(j)
1907 if (i .ne. j) then
1908 stress(i,:) = stress(i,:) + stress(j,:)
1909 diss(i,:) = diss(i,:) + diss(j,:)
1910 vol(i) = vol(i) + vol(j)
1915 i = iper(j)
1916 if (i .ne. j) then
1917 stress(j,:) = stress(i,:)
1918 diss(j,:) = diss(i,:)
1919 vol(j) = vol(i)
1928 do i = 1, 9
1929 stress(:,i) = stress(:,i)*vol(:)
1931 do i = 1, 3
1932 diss(:,i) = diss(:,i)*vol(:)
1937 do i = 1, 9
1938 xave(:,i) = stress(:,i)
1946 i = iper(j)
1947 if (i .ne. j) then
1980 do i = 1,nshg
1981 xaveg(ifath(i),:) = xaveg(ifath(i),:) + xave(i,:)
2018 do i = 1, nfath
2022 write(381,*)xavegt(i,10),xavegt(i,11),xavegt(i,12)
2037 do i = 1, nfath
2038 if (stabdis(i) .gt. zero) then
2039 cdel(i) = (two*xavegr(i,11)-stabdis(i))/xavegr(i,12)
2040 xnume(i) = two*xavegr(i,11)-stabdis(i)
2041 xdeno(i) = xavegr(i,12)
2043 xnume(i) = rnum(i)
2044 xdeno(i) = rden(i)
2048 do i = 1, nfath
2049 if (stabdis(i) .gt. zero) then
2050 cdel(i) = (two*xaveg(i,11)-stabdis(i))/xaveg(i,12)
2051 xnume(i) = two*xaveg(i,11)-stabdis(i)
2052 xdeno(i) = xaveg(i,12)
2054 xnume(i) = rnum(i)
2055 xdeno(i) = rden(i)
2060 do i = 1, nfath
2061 if (xnume(i) .lt. zero) then
2062 xnume(i) = rnum(i)
2063 xdeno(i) = rden(i)
2067 do i = 1, nshg
2068 numNden(i,1) = xnume(ifath(i))
2069 numNden(i,2) = xdeno(ifath(i))
2187 i = iper(j)
2188 if (i .ne. j) then
2189 fres(i,:) = fres(i,:) + fres(j,:)
2193 i = iper(j)
2194 if (i .ne. j) then
2195 fres(j,:) = fres(i,:)
2255 i = iper(j)
2256 if (i .ne. j) then
2315 do i = 1,nshg
2316 xnude(ifath(i),1) = xnude(ifath(i),1) + xfac(i,1)
2317 xnude(ifath(i),2) = xnude(ifath(i),2) + xfac(i,2)
2318 xnude(ifath(i),3) = xnude(ifath(i),3) + xfac(i,3)
2319 xnude(ifath(i),4) = xnude(ifath(i),4) + xfac(i,4)
2320 xnude(ifath(i),5) = xnude(ifath(i),5) + xfac(i,5)
2350 do i = 1,nfath
2351 xa(i,1) = two*cdelsq1(i)*xnuder(i,1) +
2352 & xnuder(i,2)
2353 xa(i,2) = four*cdelsq1(i)*xnuder(i,3) +
2354 & xnuder(i,4)
2355 xa(i,3) = two*cdelsq1(i)*xnuder(i,5)
2373 do i = 1,nfath
2374 xa(i,1) = two*cdelsq1(i)*xnude(i,1) +
2375 & xnude(i,2)
2376 xa(i,2) = four*cdelsq1(i)*xnude(i,3) +
2377 & xnude(i,4)
2378 xa(i,3) = two*cdelsq1(i)*xnude(i,5)
2392 do i = 1, nfath
2394 xdisc = xa(i,2)**2 - four*xa(i,1)*xa(i,3)
2398 xlamb1(i) = fwr1
2399 xlamb2(i) = fwr1
2406 xlamb1(i) = xa(i,2) / (two*xa(i,1))
2407 xlamb2(i) = xa(i,2) / (two*xa(i,1))
2411 xlamb1(i)= ( xa(i,2) + sqrt( xa(i,2)**2 - four*xa(i,1)*xa(i,3) ) )
2412 & / (two*xa(i,1))
2413 xlamb2(i)= ( xa(i,2) - sqrt( xa(i,2)**2 - four*xa(i,1)*xa(i,3) ) )
2414 & / (two*xa(i,1))
2419 do i = 1, nshg
2420 fwr2(i) = xlamb1(ifath(i))
2421 fwr3(i) = xlamb2(ifath(i))
2425 do i = 1, nfath
2426 write(23,*)i,xlamb1(i), xlamb2(i)
2433 do i = 1, nfath
2434 xkap(i) = cdelsq1(i) / xlamb2(i)
2435 xa(i,1) = two*xkap(i)*xnuder(i,1)
2436 xa(i,2) = four*xkap(i)*xnuder(i,3) - xnuder(i,2)
2437 xa(i,3) = two*xkap(i)*xnuder(i,5) - xnuder(i,4)
2439 xlamb1(i)= ( xa(i,2) + sqrt( xa(i,2)**2 - four*xa(i,1)*xa(i,3)
2440 & ) )/ (two*xa(i,1))
2441 xlamb2(i)= ( xa(i,2) - sqrt( xa(i,2)**2 - four*xa(i,1)*xa(i,3)
2442 & ) )/ (two*xa(i,1))
2448 do i = 1, nfath
2449 write(255,*)i, xlamb1(i), xlamb2(i)
2598 i = iper(j)
2599 if (i .ne. j) then
2600 fres(i,:) = fres(i,:) + fres(j,:)
2604 i = iper(j)
2605 if (i .ne. j) then
2606 fres(j,:) = fres(i,:)
2687 i = iper(j)
2688 if (i .ne. j) then
2753 do i = 1,nshg
2754 xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
2755 xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
2757 ynude(ifath(i),1) = ynude(ifath(i),1) + xfac(i,1)
2758 ynude(ifath(i),2) = ynude(ifath(i),2) + xfac(i,2)
2759 ynude(ifath(i),3) = ynude(ifath(i),3) + xfac(i,3)
2760 ynude(ifath(i),4) = ynude(ifath(i),4) + xfac(i,4)
2761 ynude(ifath(i),5) = ynude(ifath(i),5) + xfac(i,5)
2763 xm(ifath(i),1) = xm(ifath(i),1) + xmij(i,1)
2764 xm(ifath(i),2) = xm(ifath(i),2) + xmij(i,2)
2765 xm(ifath(i),3) = xm(ifath(i),3) + xmij(i,3)
2766 xm(ifath(i),4) = xm(ifath(i),4) + xmij(i,4)
2767 xm(ifath(i),5) = xm(ifath(i),5) + xmij(i,5)
2768 xm(ifath(i),6) = xm(ifath(i),6) + xmij(i,6)
2770 xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
2771 xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
2772 xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
2773 xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
2774 xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
2775 xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)
2777 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
2778 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
2779 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
2780 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
2781 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
2782 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
2784 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
2785 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
2786 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
2787 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
2788 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
2789 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
2791 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
2792 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
2793 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
2795 snorm(ifath(i)) = snorm(ifath(i)) + strnrm(i)
2813 do i = 1, nfath
2814 ynuder(i,6) = ( ynuder(i,4) - fwr1*ynuder(i,2) ) /
2815 & ( two*ynuder(i,5) - four*fwr1*ynuder(i,3)
2816 & + two*fwr1*fwr1*ynuder(i,1) )
2915 do i = 1, nfath
2916 ynude(i,6) = ( ynude(i,4) - fwr1*ynude(i,2) ) /
2917 & ( two*ynude(i,5) - four*fwr1*ynude(i,3)
2918 & + fwr1*fwr1*ynude(i,1) )
2969 do i = 1, nfath
2970 write(365,*)xmodcomp(i,1),xmodcomp(i,2),xmodcomp(i,3),
2971 & xmodcomp(i,4),xmodcomp(i,5)
2973 write(366,*)xmcomp(i,1),xmcomp(i,2),xmcomp(i,3)
2974 write(367,*)xmcomp(i,4),xmcomp(i,5),xmcomp(i,6)
2976 write(368,*)xlcomp(i,1),xlcomp(i,2),xlcomp(i,3)
2977 write(369,*)xlcomp(i,4),xlcomp(i,5),xlcomp(i,6)
2979 write(370,*)xl1comp(i,1),xl1comp(i,2),xl1comp(i,3)
2980 write(371,*)xl1comp(i,4),xl1comp(i,5),xl1comp(i,6)
2982 write(372,*)xl2comp(i,1),xl2comp(i,2),xl2comp(i,3)
2983 write(373,*)xl2comp(i,4),xl2comp(i,5),xl2comp(i,6)
2985 write(374,*)xnd(i,1),xnd(i,2),scomp(i)
2986 write(375,*)ucomp(i,1),ucomp(i,2),ucomp(i,3)
3102 do i = 1,numnp
3103 xnude(ifath(i),2) = xnude(ifath(i),2) + strnrm(i)
3119 do i = 1,numnp
3120 xnude(ifath(i),1) = xnude(ifath(i),1) + x(i,1)
3121 xnude(ifath(i),2) = xnude(ifath(i),2) + x(i,2)
3132 do i=1,nfath ! cdelsq * |S|
3133 write(444,*) xnuder(i,1),xnuder(i,2),xnutf(i)
3262 i = iper(j)
3263 if (i .ne. j) then
3264 hfres(i,:) = hfres(i,:) + hfres(j,:)
3268 i = iper(j)
3269 if (i .ne. j) then
3270 hfres(j,:) = hfres(i,:)
3321 i = iper(j)
3322 if (i .ne. j) then
3323 fres(i,:) = fres(i,:) + fres(j,:)
3328 i = iper(j)
3329 if (i .ne. j) then
3330 fres(j,:) = fres(i,:)
3433 i = iper(j)
3434 if (i .ne. j) then
3499 do i = 1,nshg
3500 xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
3501 xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
3503 ynude(ifath(i),1) = ynude(ifath(i),1) + xfac(i,1)
3504 ynude(ifath(i),2) = ynude(ifath(i),2) + xfac(i,2)
3505 ynude(ifath(i),3) = ynude(ifath(i),3) + xfac(i,3)
3506 ynude(ifath(i),4) = ynude(ifath(i),4) + xfac(i,4)
3507 ynude(ifath(i),5) = ynude(ifath(i),5) + xfac(i,5)
3509 xm(ifath(i),1) = xm(ifath(i),1) + xmij(i,1)
3510 xm(ifath(i),2) = xm(ifath(i),2) + xmij(i,2)
3511 xm(ifath(i),3) = xm(ifath(i),3) + xmij(i,3)
3512 xm(ifath(i),4) = xm(ifath(i),4) + xmij(i,4)
3513 xm(ifath(i),5) = xm(ifath(i),5) + xmij(i,5)
3514 xm(ifath(i),6) = xm(ifath(i),6) + xmij(i,6)
3516 xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
3517 xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
3518 xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
3519 xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
3520 xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
3521 xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)
3523 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
3524 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
3525 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
3526 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
3527 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
3528 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
3530 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
3531 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
3532 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
3533 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
3534 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
3535 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
3537 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
3538 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
3539 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
3541 snorm(ifath(i)) = snorm(ifath(i)) + strnrm(i)
3559 do i = 1, nfath
3560 ynuder(i,6) = ( ynuder(i,4) - fwr1*ynuder(i,2) ) /
3561 & ( two*ynuder(i,5) - four*fwr1*ynuder(i,3)
3562 & + two*fwr1*fwr1*ynuder(i,1) )
3661 do i = 1, nfath
3662 ynude(i,6) = ( ynude(i,4) - fwr1*ynude(i,2) ) /
3663 & ( two*ynude(i,5) - four*fwr1*ynude(i,3)
3664 & + fwr1*fwr1*ynude(i,1) )
3713 do i = 1, nfath
3714 write(365,*)xmodcomp(i,1),xmodcomp(i,2),xmodcomp(i,3),
3715 & xmodcomp(i,4),xmodcomp(i,5)
3717 write(366,*)xmcomp(i,1),xmcomp(i,2),xmcomp(i,3)
3718 write(367,*)xmcomp(i,4),xmcomp(i,5),xmcomp(i,6)
3720 write(368,*)xlcomp(i,1),xlcomp(i,2),xlcomp(i,3)
3721 write(369,*)xlcomp(i,4),xlcomp(i,5),xlcomp(i,6)
3723 write(370,*)xl1comp(i,1),xl1comp(i,2),xl1comp(i,3)
3724 write(371,*)xl1comp(i,4),xl1comp(i,5),xl1comp(i,6)
3726 write(372,*)xl2comp(i,1),xl2comp(i,2),xl2comp(i,3)
3727 write(373,*)xl2comp(i,4),xl2comp(i,5),xl2comp(i,6)
3729 write(374,*)xnd(i,1),xnd(i,2),scomp(i)
3730 write(375,*)ucomp(i,1),ucomp(i,2),ucomp(i,3)
3859 do i = 1,numnp
3860 xnude(ifath(i),2) = xnude(ifath(i),2) + strnrm(i)
3876 do i = 1,numnp
3877 xnude(ifath(i),1) = xnude(ifath(i),1) + x(i,1)
3878 xnude(ifath(i),2) = xnude(ifath(i),2) + x(i,2)
3889 do i=1,nfath ! cdelsq * |S|
3890 write(444,*) xnuder(i,1),xnuder(i,2),xnutf(i)