Lines Matching refs:gijd
49 & gijd(npro,6), uh1(npro),
53 call e3gijd( dxidx, gijd ) ! both diagonal taus need this
61 fact = gijd(:,1) * gijd(:,3) * gijd(:,6)
62 & - gijd(:,1) * gijd(:,5) * gijd(:,5)
63 & - gijd(:,3) * gijd(:,4) * gijd(:,4)
64 & - gijd(:,6) * gijd(:,2) * gijd(:,2)
65 & + gijd(:,2) * gijd(:,4) * gijd(:,5) * two
67 uh1= u1*u1*(gijd(:,3)*gijd(:,6)-gijd(:,5)*gijd(:,5))
68 & + u2*u2*(gijd(:,1)*gijd(:,6)-gijd(:,4)*gijd(:,4))
69 & + u3*u3*(gijd(:,1)*gijd(:,3)-gijd(:,2)*gijd(:,2))
70 & + two *(u1*u2*(gijd(:,4)*gijd(:,5)-gijd(:,2)*gijd(:,6))
71 & + u1*u3*(gijd(:,2)*gijd(:,5)-gijd(:,4)*gijd(:,3))
72 & + u2*u3*(gijd(:,4)*gijd(:,2)-gijd(:,1)*gijd(:,5)))
82 h2o2u = u1*u1*gijd(:,1)
83 & + u2*u2*gijd(:,3)
84 & + u3*u3*gijd(:,6)
85 & +(u1*u2*gijd(:,2)
86 & + u1*u3*gijd(:,4)
87 & + u2*u3*gijd(:,5))*two + 1.0e-15 !FIX FOR INVALID MESHES
163 & + (u1*(u1*gijd(:,1) + two*(u2*gijd(:,2)+u3*gijd(:,4)))
164 & + u2*(u2*gijd(:,3) + two*u3*gijd(:,5))
165 & + u3*u3*gijd(:,6)))
166 & +fff*rmu**2*(gijd(:,1)**2 + gijd(:,3)**2 + gijd(:,6)**2 +
167 & two*(gijd(:,2)**2 + gijd(:,4)**2 + gijd(:,5)**2))
169 tau(:,1)=pt125*fact/(rho*(gijd(:,1)+gijd(:,3)+gijd(:,6)))*taucfct
233 gijdu(:,1)=gijd(:,1)
234 gijdu(:,2)=gijd(:,3)
235 gijdu(:,3)=gijd(:,6)
236 gijdu(:,4)=gijd(:,2)
237 gijdu(:,5)=gijd(:,4)
238 gijdu(:,6)=gijd(:,5)
334 & gijd(npro,6),
691 call e3gijd( dxidx, gijd )
693 gijdu(:,1)=gijd(:,1)
694 gijdu(:,2)=gijd(:,3)
695 gijdu(:,3)=gijd(:,6)
696 gijdu(:,4)=gijd(:,2)
697 gijdu(:,5)=gijd(:,4)
698 gijdu(:,6)=gijd(:,5)
796 & gijd(npro,6),
1156 gijdu(:,1)=gijd(:,1)
1157 gijdu(:,2)=gijd(:,3)
1158 gijdu(:,3)=gijd(:,6)
1159 gijdu(:,4)=gijd(:,2)
1160 gijdu(:,5)=gijd(:,4)
1161 gijdu(:,6)=gijd(:,5)
1258 & gijd(npro,6), uh1(npro),
1265 call e3gijd( dxidx, gijd )
1272 fact = gijd(:,1) * gijd(:,3) * gijd(:,6)
1273 & - gijd(:,1) * gijd(:,5) * gijd(:,5)
1274 & - gijd(:,3) * gijd(:,4) * gijd(:,4)
1275 & - gijd(:,6) * gijd(:,2) * gijd(:,2)
1276 & + gijd(:,2) * gijd(:,4) * gijd(:,5) * two
1278 uh1= u1*u1*(gijd(:,3)*gijd(:,6)-gijd(:,5)*gijd(:,5))
1279 & + u2*u2*(gijd(:,1)*gijd(:,6)-gijd(:,4)*gijd(:,4))
1280 & + u3*u3*(gijd(:,1)*gijd(:,3)-gijd(:,2)*gijd(:,2))
1281 & + two *(u1*u2*(gijd(:,4)*gijd(:,5)-gijd(:,2)*gijd(:,6))
1282 & + u1*u3*(gijd(:,2)*gijd(:,5)-gijd(:,4)*gijd(:,3))
1283 & + u2*u3*(gijd(:,4)*gijd(:,2)-gijd(:,1)*gijd(:,5)))
1295 h2o2u = u1*u1*gijd(:,1)
1296 & + u2*u2*gijd(:,3)
1297 & + u3*u3*gijd(:,6)
1298 & +(u1*u2*gijd(:,2)
1299 & + u1*u3*gijd(:,4)
1300 & + u2*u3*gijd(:,5))*two + 1.0e-15 !FIX FOR INVALID MESHES
1345 gijdu(:,1)=gijd(:,1)
1346 gijdu(:,2)=gijd(:,3)
1347 gijdu(:,3)=gijd(:,6)
1348 gijdu(:,4)=gijd(:,2)
1349 gijdu(:,5)=gijd(:,4)
1350 gijdu(:,6)=gijd(:,5)
1397 subroutine e3gijd( dxidx, gijd ) argument
1401 real*8 dxidx(npro,nsd,nsd), gijd(npro,6), local
1410 gijd(:,1) = dxidx(:,1,1) * dxidx(:,1,1)
1414 gijd(:,2) = dxidx(:,1,1) * dxidx(:,1,2)
1418 gijd(:,3) = dxidx(:,1,2) * dxidx(:,1,2)
1422 gijd(:,4) = dxidx(:,1,1) * dxidx(:,1,3)
1426 gijd(:,5) = dxidx(:,1,2) * dxidx(:,1,3)
1430 gijd(:,6) = dxidx(:,1,3) * dxidx(:,1,3)
1446 gijd(:,1) = dxidx(:,1,1) * tmp1
1453 gijd(:,2) = dxidx(:,1,1) * tmp1
1457 gijd(:,3) = dxidx(:,1,2) * tmp1
1464 gijd(:,4) = dxidx(:,1,1) * tmp1
1468 gijd(:,5) = dxidx(:,1,2) * tmp1
1472 gijd(:,6) = dxidx(:,1,3) * tmp1
1481 gijd(:,1) = dxidx(:,1,1) * dxidx(:,1,1)
1485 gijd(:,2) = dxidx(:,1,1) * dxidx(:,1,2)
1489 gijd(:,3) = dxidx(:,1,2) * dxidx(:,1,2)
1493 gijd(:,4) = dxidx(:,1,1) * dxidx(:,1,3)
1497 gijd(:,5) = dxidx(:,1,2) * dxidx(:,1,3)
1501 gijd(:,6) = dxidx(:,1,3) * dxidx(:,1,3)