/dports/cad/feappv/feappv-4.1i/elements/ |
H A D | fld2d1.f | 155 vl(1) = shp2(3,1,l)*ul(1,1,4) + shp2(3,2,l)*ul(1,2,4) 156 & + shp2(3,3,l)*ul(1,3,4) + shp2(3,4,l)*ul(1,4,4) 158 vl(2) = shp2(3,1,l)*ul(2,1,4) + shp2(3,2,l)*ul(2,2,4) 159 & + shp2(3,3,l)*ul(2,3,4) + shp2(3,4,l)*ul(2,4,4) 192 al(1) = cfac*(shp2(3,1,l)*ul(1,1,5) + shp2(3,2,l)*ul(1,2,5) 193 & + shp2(3,3,l)*ul(1,3,5) + shp2(3,4,l)*ul(1,4,5)) 195 al(2) = cfac*(shp2(3,1,l)*ul(2,1,5) + shp2(3,2,l)*ul(2,2,5) 196 & + shp2(3,3,l)*ul(2,3,5) + shp2(3,4,l)*ul(2,4,5)) 214 r1(1,i) = shp2(1,i,l)*sigv(1) + shp2(2,i,l)*sigv(4) 215 r1(2,i) = shp2(1,i,l)*sigv(4) + shp2(2,i,l)*sigv(2) [all …]
|
H A D | fld2d2.f | 104 call kine2m(shp2,ul,fi,df,ndf,nel,nen,xji,lint) 213 acc(i) = acc(i) + shp2(3,j,l)*ul(i,j,5) 223 ru(1,j) = shp2(1,j,l)*sigm(1) + shp2(2,j,l)*sigm(4) 224 ru(2,j) = shp2(1,j,l)*sigm(4) + shp2(2,j,l)*sigm(2) 242 bdb = shp2(1,i,l)*ru(1,j) + shp2(2,i,l)*ru(2,j) 279 bbd(1,jj) = shp2(1,i,l)*dd(1,jj) 280 & + shp2(2,i,l)*dd(4,jj) 283 bbd(2,jj) = shp2(2,i,l)*dd(2,jj) 284 & + shp2(1,i,l)*dd(4,jj) 288 dmshp = shp2(3,i,l)*dmass [all …]
|
H A D | sld2d1.f | 88 call strn2d(d,xl,ul,tl,shp2(1,1,l),ndf,ndm,nel, 101 shpr(j) = shp2(3,j,l)/xx 138 aj1 = shp2(1,j,l)*jac(l) 139 aj2 = shp2(2,j,l)*jac(l) 140 aj3 = shp2(3,j,l)*jac0 156 aj0 = shp2(3,j,l)*dv 163 s(j1 ,k1 ) = s(j1 ,k1 ) + bd11*shp2(1,k,l) 164 & + bd14*shp2(2,k,l) 167 s(j1 ,k1+1) = s(j1 ,k1+1) + bd12*shp2(2,k,l) 168 & + bd14*shp2(1,k,l) [all …]
|
H A D | sld2d3.f | 99 call shp2d(sg2(1,l),xl,shp2(1,1,l),jac(l),ndm,nel,ix,.false.) 104 r0(l) = shp2(3,1,l)*xl(1,1) + shp2(3,2,l)*xl(1,2) 105 & + shp2(3,3,l)*xl(1,3) + shp2(3,4,l)*xl(1,4) 108 shpr(j,l) = shp2(3,j,l)/r0(l) 146 call strn2d(d,xl,ul,tl,shp2(1,1,l),ndf,ndm,nel,rr,zz,ta,eps) 292 call rays2d(d,shp2(1,1,l),sig(1,l),dd(1,1,5,l),ul(1,1,4), 301 call resid2d(cfac,lfac,dvol(l),jac0,shp2(1,1,l),eps, 323 bd(i,1) = aa(i,1,l)*shp2(1,j,l) + aa(i,4,l)*shp2(2,j,l) 325 bd(i,2) = aa(i,2,l)*shp2(2,j,l) + aa(i,4,l)*shp2(1,j,l) 330 lms = shp2(3,j,l)*dmas [all …]
|
H A D | sld2d2.f | 73 call strn2m(d,xr1,shp2(1,1,l),xl,ul,tl,type,xr0,xz0, 105 call bmat2d(type,xr0,shp2(1,i,l),g(1,i),bbar(1,1,i)) 127 vl(k) = vl(k) + shp2(3,j,l)*ul(k,j,4) 136 w11 = shp2(3,i,l)*jac(l)*d(77)*d(4) 173 al(k) = al(k) + shp2(3,j,l)*ul(k,j,5) 192 p(1,i) = p(1,i) + (b1 - ac(1))*shp2(3,i,l)*jac(l) 193 p(2,i) = p(2,i) + (b2 - ac(2))*shp2(3,i,l)*jac(l) 212 w11 = shp2(3,i,l)*dl 217 w11 = shp2(3,i,l)*dv 225 s(i1+ii,j1+ii) = s(i1+ii,j1+ii) + w11*shp2(3,j,l) [all …]
|
H A D | convec3d.f | 61 real*8 shp2(3,9), xi(3,9) local 131 call shp2d (xi(1,l),xl,shp2,xsj,ndm,nel,ix,.true.) 134 uu = uu + shp2(3,i)*ul(1,i) 151 r(1,i) = r(1,i) - qq*shp2(3,i) 155 hh = tt*shp2(3,i) 158 s(i1,j1) = s(i1,j1) + hh*shp2(3,j)
|
H A D | interp2d.f | 33 call shp2d(sg2(1,l),xl,shp2(1,1,l),jac(l),ndm,nel,ix,flag) 37 call shptri(el2(1,l), nel, shp2(1,1,l)) 40 call trishp(el2(1,l),xl,ndm,nel-4,jac(l),shp2(1,1,l))
|
H A D | shell3d.f | 212 call pzero(shp2,12*lint) 283 shp2(i,j,l) = shp2(i,j,l) - dv2 576 shp2(i,j,1) = shp2(i,j,1) - gshp2(i,j)/dv 644 shp2(i,j,l) = shp2(i,j,l) - gshp2(i,j)/dv 1076 shp2(i,j,l) = shp2(i,j,l) - dv2 1184 shp21 = shp2(1,i,l) 1185 shp22 = shp2(2,i,l) 1186 shp23 = shp2(3,i,l) 1257 y2n = - shp2(2,j,l) 1373 shp2(i,j,1) = shp2(i,j,1) - gshp2(i,j)/dv [all …]
|
/dports/cad/calculix-ccx/CalculiX/ccx_2.18/src/ |
H A D | springdamp_n2f.f | 72 call shape8q(xi,et,pl,xm,xs2,shp2,iflag) 96 b1(i,j)=shp2(4,j)*xs2(i,1)-shp2(1,j)*al(i) 97 b2(i,j)=shp2(4,j)*xs2(i,2)-shp2(2,j)*al(i) 126 dal(j,j,i)=dal(j,j,i)-shp2(4,i) 205 xmu(1,2,k)=shp2(1,k)*xs2(3,2)-shp2(2,k)*xs2(3,1) 206 xmu(2,3,k)=shp2(1,k)*xs2(1,2)-shp2(2,k)*xs2(1,1) 207 xmu(3,1,k)=shp2(1,k)*xs2(2,2)-shp2(2,k)*xs2(2,1) 309 dal(j,j,i)=-shp2(4,i) 355 shp2(1,nope)=0.d0 356 shp2(2,nope)=0.d0 [all …]
|
H A D | e_c3d_em.f | 543 & +shp2(4,jj)*(shp2(3,ii)*xsj2(2) 544 & -shp2(2,ii)*xsj2(3)) 550 & +shp2(4,jj)*(shp2(1,ii)*xsj2(3) 551 & -shp2(3,ii)*xsj2(1)) 557 & +shp2(4,jj)*(shp2(2,ii)*xsj2(1) 558 & -shp2(1,ii)*xsj2(2)) 566 & +shp2(4,ii)*(shp2(3,jj)*xsj2(2) 567 & -shp2(2,jj)*xsj2(3)) 573 & +shp2(4,ii)*(shp2(1,jj)*xsj2(3) 574 & -shp2(3,jj)*xsj2(1)) [all …]
|
H A D | evalshapefunc.f | 25 real*8 xl2(3,8),xs2(3,2),shp2(7,8),p(3),xsj2(3), local 33 call shape8q(xil,etl,xl2,xsj2,xs2,shp2,iflag) 35 call shape4q(xil,etl,xl2,xsj2,xs2,shp2,iflag) 37 call shape6tri(xil,etl,xl2,xsj2,xs2,shp2,iflag) 39 call shape3tri(xil,etl,xl2,xsj2,xs2,shp2,iflag) 43 p(j)=p(j)+xl2(j,i)*shp2(4,i)
|
H A D | printoutface.f | 44 & shp2(7,8),vold(0:mi(2),*),xi,et,xsj,temp,xi3d,et3d,ze3d,weight, 355 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 357 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 359 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 361 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 369 coords(j1)=coords(j1)+shp2(4,i1)*xl2(j1,i1) 537 & shp2(4,i1)*stn(1,konl(ifaceq(i1,ig))) 539 & shp2(4,i1)*stn(2,konl(ifaceq(i1,ig))) 541 & shp2(4,i1)*stn(3,konl(ifaceq(i1,ig))) 543 & shp2(4,i1)*stn(4,konl(ifaceq(i1,ig))) [all …]
|
H A D | advecstiff.f | 42 & xl2(3,8),xsj2(3),shp2(7,9),coords(3),xs2(3,7),dxsj2,areaj, local 154 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 156 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 158 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 160 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 169 temp=temp+tl2(j)*shp2(4,j) 180 coords(k)=coords(k)+xl2(k,j)*shp2(4,j) 199 shp2(4,nope)=-1.d0 202 sref2=sref*shp2(4,j) 205 s(k,j)=s(k,j)+sref2*shp2(4,k)
|
H A D | advecforc.f | 42 & xl2(3,8),xsj2(3),shp2(7,9),coords(3),xs2(3,7),dxsj2,areaj, local 160 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 162 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 164 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 166 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 175 temp=temp+tl2(j)*shp2(4,j) 186 coords(k)=coords(k)+xl2(k,j)*shp2(4,j) 203 tnl(j)=tnl(j)-tnlref*shp2(4,j)
|
H A D | springstiff_n2f_th.f | 34 real*8 xl(3,10),ratio(9),q(3),val,shp2(7,9),ak(5), 82 call shape8q(xi,et,pl,xm,xs2,shp2,iflag) 84 call shape4q(xi,et,pl,xm,xs2,shp2,iflag) 86 call shape6tri(xi,et,pl,xm,xs2,shp2,iflag) 90 call shape3tri(xi,et,pl,xm,xs2,shp2,iflag) 203 s(k,nope)=-shp2(4,k)*constant 207 s(i,j)=shp2(4,i)*shp2(4,j)*constant
|
H A D | springstiff_n2f.f | 310 call shape8q(xi,et,pl,xm,xs2,shp2,iflag) 334 b1(i,j)=shp2(4,j)*xs2(i,1)-shp2(1,j)*al(i) 335 b2(i,j)=shp2(4,j)*xs2(i,2)-shp2(2,j)*al(i) 364 dal(j,j,i)=dal(j,j,i)-shp2(4,i) 499 xmu(1,2,k)=shp2(1,k)*xs2(3,2)-shp2(2,k)*xs2(3,1) 500 xmu(2,3,k)=shp2(1,k)*xs2(1,2)-shp2(2,k)*xs2(1,1) 501 xmu(3,1,k)=shp2(1,k)*xs2(2,2)-shp2(2,k)*xs2(2,1) 619 dal(j,j,i)=-shp2(4,i) 745 shp2(1,nope)=0.d0 746 shp2(2,nope)=0.d0 [all …]
|
H A D | e_c3d.f | 205 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1514 shptil2(1,i1)=shp2(1,i1) 1515 shptil2(2,i1)=shp2(2,i1) 1516 shptil2(3,i1)=shp2(3,i1) 1517 shptil2(4,i1)=shp2(4,i1) 1658 shptil2(1,i1)=shp2(1,i1) 1659 shptil2(2,i1)=shp2(2,i1) 1660 shptil2(3,i1)=shp2(3,i1) 1661 shptil2(4,i1)=shp2(4,i1) 1817 & shp2(1,jj))) [all …]
|
H A D | dload.f | 118 real*8 xl2(3,8),pres(8),xi,et,xsj2(3),xs2(3,7),shp2(7,8), local 252 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 254 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 256 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 258 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 265 f=f+pres(j)*shp2(4,j)
|
H A D | extrapolatecontact.f | 44 & co(3,*),xi,et,vold(0:mi(2),*),xs2(3,7),xsj2(3),shp2(7,8), 207 call shape4q(xi,et,pl,xsj2,xs2,shp2,iflag) 209 aa(j+1-i,k)=shp2(4,k) 215 call shape4q(xi,et,pl,xsj2,xs2,shp2,iflag) 218 aa(k,l)=aa(k,l)+shp2(4,k)*shp2(4,l) 223 bb(k,l)=bb(k,l)+shp2(4,k)*yi(l,1,j) 331 call shape3tri(xi,et,pl,xsj2,xs2,shp2,iflag) 333 cc(j+1-i,k)=shp2(4,k) 339 call shape3tri(xi,et,pl,xsj2,xs2,shp2,iflag) 342 cc(k,l)=cc(k,l)+shp2(4,k)*shp2(4,l) [all …]
|
H A D | calcexternalwork.f | 38 & shp2(7,8),vold(0:mi(2),*),xi,et,weight,voldprev(0:mi(2),*), local 248 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 250 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 252 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 254 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 262 disp(j1)=disp(j1)+shp2(4,i1)*voldl2(j1,i1)
|
H A D | e_c3d_cs_se.f | 77 & xl2(3,9),xsj2(3),shp2(7,9),vold(0:mi(2),*),xload(2,*), local 207 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 246 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 1434 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1436 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1659 & shp2(2,jj))+xsj2(2)* 1661 & shp2(3,jj))+xsj2(3)* 1663 & shp2(1,jj))) 1689 & shp2(2,jj))+xsj2(2)* 1691 & shp2(3,jj))+xsj2(3)* [all …]
|
H A D | radmatrix.f | 46 & vold(0:mi(2),*),co(3,*),shp2(7,8),xs2(3,7), 236 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 238 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 240 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 242 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 250 temp=temp+tl2(j)*shp2(4,j) 263 coords(k)=coords(k)+xl2(k,j)*shp2(4,j)
|
H A D | resultsem.f | 40 & elconloc(21),xi,et,ze,xsj,t1l,dtime,weight,xsj2(3),shp2(7,9), 421 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 423 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 425 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 427 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 438 phi=phi+shp2(4,m1)*vl2(5,m1) 448 & (shp2(l,m1)*xsj2(mm)-shp2(mm,m1)*xsj2(l)) 462 vkl(k,m3)=vkl(k,m3)+shp2(m3,m1)*vl2(k,m1) 474 & shp2(4,m1)*(vkl(mm,l)-vkl(l,mm))*xsj2(k)
|
H A D | e_c3d_se.f | 77 & xl2(3,9),xsj2(3),shp2(7,9),vold(0:mi(2),*),xload(2,*), local 246 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 284 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 1574 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1576 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1789 & shp2(2,jj))+xsj2(2)* 1791 & shp2(3,jj))+xsj2(3)* 1793 & shp2(1,jj))) 1819 & shp2(2,jj))+xsj2(2)* 1821 & shp2(3,jj))+xsj2(3)* [all …]
|
H A D | e_c3d_duds.f | 80 & xl2(3,9),xsj2(3),shp2(7,9),vold(0:mi(2),*),xload(2,*), 251 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 289 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 1598 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1600 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 1813 & shp2(2,jj))+xsj2(2)* 1815 & shp2(3,jj))+xsj2(3)* 1817 & shp2(1,jj))) 1843 & shp2(2,jj))+xsj2(2)* 1845 & shp2(3,jj))+xsj2(3)* [all …]
|