/dports/cad/calculix-ccx/CalculiX/ccx_2.18/src/ |
H A D | normalsoninterface.f | 246 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 248 xsj2(1)=xsj2(1)/dd 249 xsj2(2)=xsj2(2)/dd 250 xsj2(3)=xsj2(3)/dd 274 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 276 xsj2(1)=xsj2(1)/dd 277 xsj2(2)=xsj2(2)/dd 278 xsj2(3)=xsj2(3)/dd 303 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 305 xsj2(1)=xsj2(1)/dd [all …]
|
H A D | normalsonsurface_robust.f | 307 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 309 xsj2(1)=xsj2(1)/dd 310 xsj2(2)=xsj2(2)/dd 311 xsj2(3)=xsj2(3)/dd 372 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 374 xsj2(1)=xsj2(1)/dd 375 xsj2(2)=xsj2(2)/dd 376 xsj2(3)=xsj2(3)/dd 437 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 439 xsj2(1)=xsj2(1)/dd [all …]
|
H A D | updatecontpen.f | 239 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 241 xsj2(1)=xsj2(1)/dd 242 xsj2(2)=xsj2(2)/dd 243 xsj2(3)=xsj2(3)/dd 267 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 269 xsj2(1)=xsj2(1)/dd 270 xsj2(2)=xsj2(2)/dd 271 xsj2(3)=xsj2(3)/dd 296 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 298 xsj2(1)=xsj2(1)/dd [all …]
|
H A D | normalsonsurface_se.f | 307 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 309 xsj2(1)=xsj2(1)/dd 310 xsj2(2)=xsj2(2)/dd 311 xsj2(3)=xsj2(3)/dd 372 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 374 xsj2(1)=xsj2(1)/dd 375 xsj2(2)=xsj2(2)/dd 376 xsj2(3)=xsj2(3)/dd 437 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 439 xsj2(1)=xsj2(1)/dd [all …]
|
H A D | nortanslav.f | 53 real*8 co(3,*),vold(0:mi(2),*),et,xi,xl2s(3,8),xsj2(3), local 159 call shape8q(xi,et,xl2s,xsj2,xs2,shp2,iflag) 161 call shape4q(xi,et,xl2s,xsj2,xs2,shp2,iflag) 167 dd=dsqrt(xsj2(1)*xsj2(1) + xsj2(2)*xsj2(2) 168 & + xsj2(3)*xsj2(3)) 169 xsj2(1)=xsj2(1)/dd 170 xsj2(2)=xsj2(2)/dd 171 xsj2(3)=xsj2(3)/dd 179 & +xsj2(1) 181 & +xsj2(2) [all …]
|
H A D | printoutface.f | 468 dd=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 469 & xsj2(3)*xsj2(3)) 471 tf(i1)=dvi*(t(i1,1)*xsj2(1)+t(i1,2)*xsj2(2)+ 472 & t(i1,3)*xsj2(3))-pres*xsj2(i1) 476 tn=(tf(1)*xsj2(1)+tf(2)*xsj2(2)+tf(3)*xsj2(3))/dd 502 dd=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 503 & xsj2(3)*xsj2(3)) 587 df(i1)=(t(i1,1)*xsj2(1)+ 588 & t(i1,2)*xsj2(2)+ 603 dd=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ [all …]
|
H A D | printoutfacefem.f | 43 & voldl(0:mi(2),8),xl2(3,8),xsj2(3), 215 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 217 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 294 dd=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 295 & xsj2(3)*xsj2(3)) 297 tf(i1)=dvi*(t(i1,1)*xsj2(1)+t(i1,2)*xsj2(2)+ 298 & t(i1,3)*xsj2(3))-pres*xsj2(i1) 302 tn=(tf(1)*xsj2(1)+tf(2)*xsj2(2)+tf(3)*xsj2(3))/dd 303 tt=dsqrt((tf(1)-tn*xsj2(1)/dd)**2+ 304 & (tf(2)-tn*xsj2(2)/dd)**2+ [all …]
|
H A D | springforc_n2f_th.f | 37 & elcon(0:ncmat_,ntmat_,*),pproj(3),xsj2(3),xs2(3,7),val, local 80 call shape8q(xi,et,pl,xsj2,xs2,shp2,iflag) 82 call shape4q(xi,et,pl,xsj2,xs2,shp2,iflag) 84 call shape6tri(xi,et,pl,xsj2,xs2,shp2,iflag) 88 call shape3tri(xi,et,pl,xsj2,xs2,shp2,iflag) 93 dm=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+xsj2(3)*xsj2(3)) 95 xn(i)=xsj2(i)/dm
|
H A D | e_c3d_em.f | 510 & shp2(4,k)*(h0l(2)*xsj2(3)-h0l(3)*xsj2(2))*weight 512 & shp2(4,k)*(h0l(3)*xsj2(1)-h0l(1)*xsj2(3))*weight 514 & shp2(4,k)*(h0l(1)*xsj2(2)-h0l(2)*xsj2(1))*weight 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)) 558 & -shp2(1,ii)*xsj2(2)) 567 & -shp2(2,jj)*xsj2(3)) 574 & -shp2(3,jj)*xsj2(1)) [all …]
|
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) 169 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 170 & xsj2(3)*xsj2(3))
|
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) 163 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 164 & xsj2(3)*xsj2(3))
|
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)
|
H A D | geocfd.f | 192 call shape4q(xi,et,xl2,xsj2,xs2,shp2,kflag) 196 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 197 & xsj2(3)*xsj2(3)) 222 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,kflag) 224 call shape4q(xi,et,xl2,xsj2,xs2,shp2,kflag) 229 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 230 & xsj2(3)*xsj2(3)) 258 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,kflag) 262 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 263 & xsj2(3)*xsj2(3)) [all …]
|
H A D | radmatrix.f | 45 & xl2(3,8),coords(3),dxsj2,temp,xi,et,weight,xsj2(3), 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) 245 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 246 & xsj2(3)*xsj2(3))
|
H A D | e_c3d_v1rhs.f | 756 xsj2(i1)=varf(index) 912 tt(1)=(t(1,1)*xsj2(1)+t(1,2)*xsj2(2)+t(1,3)*xsj2(3))* 914 tt(2)=(t(1,2)*xsj2(1)+t(2,2)*xsj2(2)+t(2,3)*xsj2(3))* 916 tt(3)=(t(1,3)*xsj2(1)+t(2,3)*xsj2(2)+t(3,3)*xsj2(3))* 934 tvn=tv(1)*xsj2(1)+tv(2)*xsj2(2)+tv(3)*xsj2(3) 940 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 941 & xsj2(3)*xsj2(3)) 1097 tvnk=tvk(1)*xsj2(1)+tvk(2)*xsj2(2)+tvk(3)*xsj2(3) 1098 tvnt=tvt(1)*xsj2(1)+tvt(2)*xsj2(2)+tvt(3)*xsj2(3) 1122 & (rhovel(1)*xsj2(1)+rhovel(2)*xsj2(2)+ [all …]
|
H A D | calcexternalwork.f | 37 real*8 co(3,*),xs2(3,7),pres,voldl2(3,8),xl2(3,8),xsj2(3), 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) 269 & (disp(1)*xsj2(1)+disp(2)*xsj2(2)+disp(3)*xsj2(3))
|
H A D | gendualcoeffs.f | 199 call shape8qtilde_lin(xis,ets,xl2s,xsj2,xs2, 203 call shape4q(xis,ets,xl2s,xsj2,xs2,shp2,iflag) 206 call shape6tritilde(xis,ets,xl2s,xsj2,xs2,shp2, 209 call shape6tritilde_lin(xis,ets,xl2s,xsj2,xs2, 213 call shape3tri(xis,ets,xl2s,xsj2,xs2,shp2,iflag) 215 dx=dsqrt(xsj2(1)**2+xsj2(2)**2+ 216 & xsj2(3)**2) 698 & (xis,ets,xl2s,xsj2,xs2,shp2s,ns, 702 & (xis,ets,xl2s,xsj2,xs2,shp2s,ns, 809 call dualshape8qtilde(xis,ets,xl2s,xsj2,xs2, [all …]
|
H A D | calcstabletimeincvol.f | 73 & xl2(3,9),xsj2(3),xs2(3,7),shp2(7,9),hmin,area, 392 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 394 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 396 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 400 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 403 area=weight*dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 404 & xsj2(3)*xsj2(3))
|
H A D | e_c3d_rhs_th.f | 46 & xl2(3,8),xsj2(3),shp2(7,8),vold(0:mi(2),*),xload(2,*), local 400 call shape8q(xi,et,xl2,xsj2,xs2,shp2,iflag) 402 call shape4q(xi,et,xl2,xsj2,xs2,shp2,iflag) 404 call shape6tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 406 call shape3tri(xi,et,xl2,xsj2,xs2,shp2,iflag) 409 dxsj2=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+ 410 & xsj2(3)*xsj2(3))
|
H A D | slavintpoints.f | 54 & xtri(2,6),xi,et,xsj2(3),xs2(3,2),shp2(7,8),anglesm 163 call shape8q(xi,et,xl2s,xsj2,xs2,shp2,iflag) 165 call shape4q(xi,et,xl2s,xsj2,xs2,shp2,iflag) 167 call shape6tri(xi,et,xl2s,xsj2,xs2,shp2,iflag) 169 call shape3tri(xi,et,xl2s,xsj2,xs2,shp2,iflag) 171 dd=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2) 172 & +xsj2(3)*xsj2(3)) 173 xsj2(1)=xsj2(1)/dd 174 xsj2(2)=xsj2(2)/dd 175 xsj2(3)=xsj2(3)/dd [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)
|
H A D | e_c3d.f | 1593 t(1)=stress(1)*xsj2(1)+stress(4)*xsj2(2)+ 1594 & stress(6)*xsj2(3) 1595 t(2)=stress(4)*xsj2(1)+stress(2)*xsj2(2)+ 1596 & stress(5)*xsj2(3) 1597 t(3)=stress(6)*xsj2(1)+stress(5)*xsj2(2)+ 1598 & stress(3)*xsj2(3) 1602 xsj2(k)=t(k) 1621 & *xsj2(1)*weight 1624 & *xsj2(2)*weight 1627 & *xsj2(3)*weight [all …]
|
H A D | genfirstactif.f | 58 & xs2(3,2),xsj2(3),tietol(3,*),adjust 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) 168 area=area+weight*dsqrt(xsj2(1)**2+xsj2(2)**2+ 169 & xsj2(3)**2)
|
H A D | springforc_n2f.f | 40 & elcon(0:ncmat_,ntmat_,*),pproj(3),xsj2(3),xs2(3,7),clear, local 233 call shape8q(xi,et,pl,xsj2,xs2,shp2,iflag) 235 call shape4q(xi,et,pl,xsj2,xs2,shp2,iflag) 237 call shape6tri(xi,et,pl,xsj2,xs2,shp2,iflag) 239 call shape3tri(xi,et,pl,xsj2,xs2,shp2,iflag) 244 dm=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2)+xsj2(3)*xsj2(3)) 246 xn(i)=xsj2(i)/dm
|
H A D | slavintmortar.f | 79 & xtri(2,6),xi,et,xsj2(3),xs2(3,7),shp2(7,8),anglesm 213 call shape8q(xi,et,xl2s,xsj2,xs2,shp2,iflag) 215 call shape4q(xi,et,xl2s,xsj2,xs2,shp2,iflag) 217 call shape6tri(xi,et,xl2s,xsj2,xs2,shp2,iflag) 219 call shape3tri(xi,et,xl2s,xsj2,xs2,shp2,iflag) 221 dd=dsqrt(xsj2(1)*xsj2(1)+xsj2(2)*xsj2(2) 222 & +xsj2(3)*xsj2(3)) 223 xsj2(1)=xsj2(1)/dd 224 xsj2(2)=xsj2(2)/dd 225 xsj2(3)=xsj2(3)/dd [all …]
|