Lines Matching refs:dxidx
39 & sgn(npro,nshl), dxidx(npro,nsd,nsd),
56 call e3metric( xl, shdrv,dxidx,shg,WdetJ)
93 & sgn(npro,nshl), dxidx(npro,nsd,nsd),
120 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf)
134 dxidx = zero
155 call e3metric( xl, shdrv,dxidx,shg,WdetJ)
216 & sgn(npro,nshl), dxidx(npro,nsd,nsd),
237 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf)
251 dxidx = zero
272 call e3metric( xl, shdrv,dxidx,shg,WdetJ)
316 subroutine e3metricf( xl, shgl, dxidx, argument
322 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd), local
346 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3)
348 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3)
350 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3)
352 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1)
353 & + dxidx(:,1,2) * dxdxi(:,2,1)
354 & + dxidx(:,1,3) * dxdxi(:,3,1) )
355 dxidx(:,1,1) = dxidx(:,1,1) * tmp
356 dxidx(:,1,2) = dxidx(:,1,2) * tmp
357 dxidx(:,1,3) = dxidx(:,1,3) * tmp
358 dxidx(:,2,1) = (dxdxi(:,2,3) * dxdxi(:,3,1)
360 dxidx(:,2,2) = (dxdxi(:,1,1) * dxdxi(:,3,3)
362 dxidx(:,2,3) = (dxdxi(:,2,1) * dxdxi(:,1,3)
364 dxidx(:,3,1) = (dxdxi(:,2,1) * dxdxi(:,3,2)
366 dxidx(:,3,2) = (dxdxi(:,3,1) * dxdxi(:,1,2)
368 dxidx(:,3,3) = (dxdxi(:,1,1) * dxdxi(:,2,2)
378 shg(:,n,1) = shgl(:,1,n) * dxidx(:,1,1) +
379 & shgl(:,2,n) * dxidx(:,2,1) +
380 & shgl(:,3,n) * dxidx(:,3,1)
381 shg(:,n,2) = shgl(:,1,n) * dxidx(:,1,2) +
382 & shgl(:,2,n) * dxidx(:,2,2) +
383 & shgl(:,3,n) * dxidx(:,3,2)
384 shg(:,n,3) = shgl(:,1,n) * dxidx(:,1,3) +
385 & shgl(:,2,n) * dxidx(:,2,3) +
386 & shgl(:,3,n) * dxidx(:,3,3)