Lines Matching refs:gijd
38 dimension gijd(npro,6), fact(npro), rnu(npro),
44 call e3gijd( dxidx, gijd )
78 3 + ( u1 * ( gijd(:,1) * u1
79 4 + gijd(:,4) * u2
80 5 + gijd(:,6) * u3 )
81 6 + u2 * ( gijd(:,4) * u1
82 7 + gijd(:,2) * u2
83 8 + gijd(:,5) * u3 )
84 9 + u3 * ( gijd(:,6) * u1
85 a + gijd(:,5) * u2
86 1 + gijd(:,3) * u3 ) ) )
88 3 * ( gijd(:,1) ** 2
89 4 + gijd(:,2) ** 2
90 5 + gijd(:,3) ** 2
92 7 * ( gijd(:,4) ** 2
93 8 + gijd(:,5) ** 2
94 9 + gijd(:,6) ** 2 )
100 tauC =rho* pt125*fact/(gijd(:,1)+gijd(:,2)+gijd(:,3))*ff
107 fact = gijd(:,1) * gijd(:,2) * gijd(:,3)
108 & - gijd(:,2) * gijd(:,6) * gijd(:,6)
109 & - gijd(:,1) * gijd(:,5) * gijd(:,5)
110 & - gijd(:,3) * gijd(:,4) * gijd(:,4)
111 & + gijd(:,6) * gijd(:,4) * gijd(:,5) * two
120 1 u1 * ( (gijd(:,2)*gijd(:,3)-gijd(:,5)*gijd(:,5)) * u1
121 2 + two * (gijd(:,5)*gijd(:,6)-gijd(:,4)*gijd(:,3)) * u2
122 3 + two * (gijd(:,4)*gijd(:,5)-gijd(:,6)*gijd(:,2)) * u3)
123 1 + u2 * ( (gijd(:,1)*gijd(:,3)-gijd(:,6)*gijd(:,6)) * u2
124 3 + two * (gijd(:,4)*gijd(:,6)-gijd(:,1)*gijd(:,5)) * u3)
125 1 + u3 * ( (gijd(:,1)*gijd(:,2)-gijd(:,4)*gijd(:,4)) * u3)
142 3 u1 * ( gijd(:,1) * u1
143 4 + gijd(:,4) * u2
144 5 + gijd(:,6) * u3 )
145 6 + u2 * ( gijd(:,4) * u1
146 7 + gijd(:,2) * u2
147 8 + gijd(:,5) * u3 )
148 9 + u3 * ( gijd(:,6) * u1
149 a + gijd(:,5) * u2
150 1 + gijd(:,3) * u3 )
170 3 * ( gijd(:,1) ** 2
171 4 + gijd(:,2) ** 2
172 5 + gijd(:,3) ** 2
174 7 * ( gijd(:,4) ** 2
175 8 + gijd(:,5) ** 2
176 9 + gijd(:,6) ** 2 ) )
186 tauM=(gijd(:,1)+gijd(:,2)+gijd(:,3))/unorm ! here we have 4/u^2h^2
202 3 u1 * ( gijd(:,1) * u1
203 4 + gijd(:,4) * u2
204 5 + gijd(:,6) * u3 )
205 6 + u2 * ( gijd(:,4) * u1
206 7 + gijd(:,2) * u2
207 8 + gijd(:,5) * u3 )
208 9 + u3 * ( gijd(:,6) * u1
209 a + gijd(:,5) * u2
210 1 + gijd(:,3) * u3 )
230 3 * ( gijd(:,1) ** 2
231 4 + gijd(:,2) ** 2
232 5 + gijd(:,3) ** 2
234 7 * ( gijd(:,4) ** 2
235 8 + gijd(:,5) ** 2
236 9 + gijd(:,6) ** 2 ) )
247 fact = gijd(:,1) * gijd(:,2) * gijd(:,3)
248 & - gijd(:,2) * gijd(:,6) * gijd(:,6)
249 & - gijd(:,1) * gijd(:,5) * gijd(:,5)
250 & - gijd(:,3) * gijd(:,4) * gijd(:,4)
251 & + gijd(:,6) * gijd(:,4) * gijd(:,5) * two
260 1 u1 * ( (gijd(:,2)*gijd(:,3)-gijd(:,5)*gijd(:,5)) * u1
261 2 + two * (gijd(:,5)*gijd(:,6)-gijd(:,4)*gijd(:,3)) * u2
262 3 + two * (gijd(:,4)*gijd(:,5)-gijd(:,6)*gijd(:,2)) * u3)
263 1 + u2 * ( (gijd(:,1)*gijd(:,3)-gijd(:,6)*gijd(:,6)) * u2
264 3 + two * (gijd(:,4)*gijd(:,6)-gijd(:,1)*gijd(:,5)) * u3)
265 1 + u3 * ( (gijd(:,1)*gijd(:,2)-gijd(:,4)*gijd(:,4)) * u3)
280 3 u1 * ( gijd(:,1) * u1
281 4 + gijd(:,4) * u2
282 5 + gijd(:,6) * u3 )
283 6 + u2 * ( gijd(:,4) * u1
284 7 + gijd(:,2) * u2
285 8 + gijd(:,5) * u3 )
286 9 + u3 * ( gijd(:,6) * u1
287 a + gijd(:,5) * u2
288 1 + gijd(:,3) * u3 )
299 tauBar = rLui(:,1) * ( gijd(:,1) * rLui(:,1)
300 & + gijd(:,4) * rLui(:,2)
301 & + gijd(:,6) * rLui(:,3) )
302 & + rLui(:,2) * ( gijd(:,4) * rLui(:,1)
303 & + gijd(:,2) * rLui(:,2)
304 & + gijd(:,5) * rLui(:,3) )
305 & + rLui(:,3) * ( gijd(:,6) * rLui(:,1)
306 & + gijd(:,5) * rLui(:,2)
307 & + gijd(:,3) * rLui(:,3) )
338 real*8 gijd(npro,6), tauM(npro) local
343 call e3gijd( dxidx, gijd )
364 3 + ( ui(:,1) * ( gijd(:,1) * ui(:,1)
365 4 + gijd(:,4) * ui(:,2)
366 5 + gijd(:,6) * ui(:,3) )
367 6 + ui(:,2) * ( gijd(:,4) * ui(:,1)
368 7 + gijd(:,2) * ui(:,2)
369 8 + gijd(:,5) * ui(:,3) )
370 9 + ui(:,3) * ( gijd(:,6) * ui(:,1)
371 a + gijd(:,5) * ui(:,2)
372 1 + gijd(:,3) * ui(:,3) ) ) )
374 3 * ( gijd(:,1) ** 2
375 4 + gijd(:,2) ** 2
376 5 + gijd(:,3) ** 2
378 7 * ( gijd(:,4) ** 2
379 8 + gijd(:,5) ** 2
380 9 + gijd(:,6) ** 2 ) )
396 subroutine e3gijd( dxidx, gijd ) argument
400 real*8 dxidx(npro,nsd,nsd), gijd(npro,6), local
410 gijd(:,1) = dxidx(:,1,1) * dxidx(:,1,1)
414 gijd(:,4) = dxidx(:,1,1) * dxidx(:,1,2)
418 gijd(:,2) = dxidx(:,1,2) * dxidx(:,1,2)
422 gijd(:,5) = dxidx(:,1,2) * dxidx(:,1,3)
426 gijd(:,6) = dxidx(:,1,1) * dxidx(:,1,3)
430 gijd(:,3) = dxidx(:,1,3) * dxidx(:,1,3)
446 gijd(:,1) = dxidx(:,1,1) * tmp1
453 gijd(:,2) = dxidx(:,1,2) * tmp1
457 gijd(:,4) = dxidx(:,1,1) * tmp1
464 gijd(:,3) = dxidx(:,1,3) * tmp1
468 gijd(:,5) = dxidx(:,1,2) * tmp1
472 gijd(:,6) = dxidx(:,1,1) * tmp1
500 real*8 gijd(npro,6), giju(npro,6), local
513 call e3gijd( dxidx, gijd )
539 3 + uMod(:,1) * ( gijd(:,1) * uMod(:,1)
540 4 + gijd(:,4) * uMod(:,2)
541 5 + gijd(:,6) * uMod(:,3) )
542 6 + uMod(:,2) * ( gijd(:,4) * uMod(:,1)
543 7 + gijd(:,2) * uMod(:,2)
544 8 + gijd(:,5) * uMod(:,3) )
545 9 + uMod(:,3) * ( gijd(:,6) * uMod(:,1)
546 a + gijd(:,5) * uMod(:,2)
547 1 + gijd(:,3) * uMod(:,3) )
549 3 * ( gijd(:,1) ** 2
550 4 + gijd(:,2) ** 2
551 5 + gijd(:,3) ** 2
553 7 * ( gijd(:,4) ** 2
554 8 + gijd(:,5) ** 2
555 9 + gijd(:,6) ** 2 ) )
565 fact = one/(gijd(:,1) * gijd(:,2) * gijd(:,3)
566 & - gijd(:,2) * gijd(:,6) * gijd(:,6)
567 & - gijd(:,1) * gijd(:,5) * gijd(:,5)
568 & - gijd(:,3) * gijd(:,4) * gijd(:,4)
569 & + gijd(:,6) * gijd(:,4) * gijd(:,5) * two)
574 giju(:,1) = fact * (gijd(:,2)*gijd(:,3)
575 & - gijd(:,5)**2)
576 giju(:,2) = fact * (gijd(:,1)*gijd(:,3)
577 & - gijd(:,6)**2)
578 giju(:,3) = fact * (gijd(:,1)*gijd(:,2)
579 & - gijd(:,4)**2)
580 giju(:,4) = fact * (gijd(:,5)*gijd(:,6)
581 & - gijd(:,4)*gijd(:,3) )
582 giju(:,5) = fact * (gijd(:,4)*gijd(:,6)
583 & - gijd(:,1)*gijd(:,5) )
584 giju(:,6) = fact * (gijd(:,4)*gijd(:,5)
585 & - gijd(:,6)*gijd(:,2) )