Searched refs:gijd (Results 1 – 2 of 2) sorted by relevance
| /phasta/phSolver/incompressible/ |
| H A D | e3stab.f | 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
[all …]
|
| /phasta/phSolver/compressible/ |
| H A D | e3tau.f | 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))
[all …]
|