Lines Matching refs:dxdxi
305 dimension tmp(npro), dxdxi(npro,nsd,nsd)
310 dxdxi = zero
313 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n)
314 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n)
315 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n)
316 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n)
317 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n)
318 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n)
319 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n)
320 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n)
321 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(:,3,n)
326 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3)
327 & - dxdxi(:,3,2) * dxdxi(:,2,3)
328 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3)
329 & - dxdxi(:,1,2) * dxdxi(:,3,3)
330 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3)
331 & - dxdxi(:,1,3) * dxdxi(:,2,2)
332 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1)
333 & + dxidx(:,1,2) * dxdxi(:,2,1)
334 & + dxidx(:,1,3) * dxdxi(:,3,1) )
338 dxidx(:,2,1) = (dxdxi(:,2,3) * dxdxi(:,3,1)
339 & - dxdxi(:,2,1) * dxdxi(:,3,3)) * tmp
340 dxidx(:,2,2) = (dxdxi(:,1,1) * dxdxi(:,3,3)
341 & - dxdxi(:,3,1) * dxdxi(:,1,3)) * tmp
342 dxidx(:,2,3) = (dxdxi(:,2,1) * dxdxi(:,1,3)
343 & - dxdxi(:,1,1) * dxdxi(:,2,3)) * tmp
344 dxidx(:,3,1) = (dxdxi(:,2,1) * dxdxi(:,3,2)
345 & - dxdxi(:,2,2) * dxdxi(:,3,1)) * tmp
346 dxidx(:,3,2) = (dxdxi(:,3,1) * dxdxi(:,1,2)
347 & - dxdxi(:,1,1) * dxdxi(:,3,2)) * tmp
348 dxidx(:,3,3) = (dxdxi(:,1,1) * dxdxi(:,2,2)
349 & - dxdxi(:,1,2) * dxdxi(:,2,1)) * tmp