| /phasta/phSolver/compressible/ |
| H A D | e3visc.f | 2 & dxidx, 44 & Sclr(npro), dxidx(npro,nsd,nsd), 148 f1 = dxidx(:,1,1) * dxidx(:,1,1) + 149 & dxidx(:,2,1) * dxidx(:,2,1) + 150 & dxidx(:,3,1) * dxidx(:,3,1) 151 f2 = dxidx(:,1,1) * dxidx(:,1,2) + 152 & dxidx(:,2,1) * dxidx(:,2,2) + 153 & dxidx(:,3,1) * dxidx(:,3,2) 154 f3 = dxidx(:,1,2) * dxidx(:,1,2) + 155 & dxidx(:,2,2) * dxidx(:,2,2) + [all …]
|
| H A D | e3eig1.f | 4 & dxidx, u, Q) 46 & eb1(npro), dxidx(npro,nsd,nsd), 111 a1 = Rcs(:,1,1) * dxidx(:,1,1) + Rcs(:,1,2) * dxidx(:,1,2) + 112 & Rcs(:,1,3) * dxidx(:,1,3) 113 a2 = Rcs(:,2,1) * dxidx(:,1,1) + Rcs(:,2,2) * dxidx(:,1,2) + 114 & Rcs(:,2,3) * dxidx(:,1,3) 115 a3 = Rcs(:,3,1) * dxidx(:,1,1) + Rcs(:,3,2) * dxidx(:,1,2) + 116 & Rcs(:,3,3) * dxidx(:,1,3) 117 dxidx(:,1,1) = a1 118 dxidx(:,1,2) = a2 [all …]
|
| H A D | e3tau.f | 3 & con, dxidx, rLyi, 42 & dxidx(npro,nsd,nsd), rk(npro), 53 call e3gijd( dxidx, gijd ) ! both diagonal taus need this 287 & con, dxidx, rLyi, 327 & dxidx(npro,nsd,nsd), rk(npro), 374 & dxidx, u, Q) 691 call e3gijd( dxidx, gijd ) 748 & con, dxidx, rLyi, 789 & dxidx(npro,nsd,nsd), rk(npro), 836 & dxidx, u, Q) [all …]
|
| H A D | e3.f | 67 dimension dxidx(npro,nsd,nsd), WdetJ(npro) 133 & shg, dxidx, WdetJ, 178 & dxidx, 192 & dui, dxidx, ytargetl, 205 & rLyi, dxidx, tau, 378 dimension dxidx(npro,nsd,nsd), WdetJ(npro) 432 & dxidx, Sclr, 495 & rLyti, dxidx, raLSt,
|
| H A D | e3source.f | 4 & dui, dxidx, ytargetl, 40 & dxidx( npro, nsd, nsd), xfind( npro ), delta(npro), rat 62 delta(:) = 0.5*sqrt(dxidx(:,1,1)*dxidx(:,1,1) ! 1/dx 63 . +dxidx(:,2,1)*dxidx(:,2,1) 64 . +dxidx(:,3,1)*dxidx(:,3,1))
|
| H A D | e3qvar.f | 3 & shg, dxidx, WdetJ, T, argument 45 & dxidx(npro,nsd,nsd), WdetJ(npro), 103 call e3metric( xl, shgl, dxidx,
|
| H A D | e3ivar.f | 5 & shg, dxidx, WdetJ, argument 70 & shg(npro,nshl,nsd), dxidx(npro,nsd,nsd), 237 call e3metric( xl, shgl, dxidx, 493 & dxidx, Sclr, argument 560 & shg(npro,nshl,nsd), dxidx(npro,nsd,nsd), 659 call e3metric( xl, shgl, dxidx,
|
| H A D | e3ls.f | 5 & rLyi, dxidx, tau, argument 67 & rLyi(npro,nflow), dxidx(npro,nsd,nsd), 192 & con, dxidx, rLyi, 211 & con, dxidx, rLyi, 778 & rLyti, dxidx, raLSt, argument 836 & rLyti(npro), dxidx(npro,nsd,nsd), 893 & dxidx, rLyti, rLymti,
|
| /phasta/phSolver/incompressible/ |
| H A D | e3stab.f | 2 & u3, dxidx, rLui, argument 32 & dxidx(npro,nsd,nsd), 44 call e3gijd( dxidx, gijd ) 329 subroutine e3uBar (rho, ui, dxidx, argument 335 & dxidx(npro,nsd,nsd), rLui(npro,nsd), local 343 call e3gijd( dxidx, gijd ) 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) 411 & + dxidx(:,2,1) * dxidx(:,2,1) [all …]
|
| H A D | e3qvar.f | 4 & dxidx, WdetJ ) argument 41 & dxidx(npro,nsd,nsd), WdetJ(npro) 66 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 68 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 70 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 72 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 73 & + dxidx(:,1,2) * dxdxi(:,2,1) 74 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 75 dxidx(:,1,1) = dxidx(:,1,1) * tmp 76 dxidx(:,1,2) = dxidx(:,1,2) * tmp [all …]
|
| H A D | hessian.f | 99 & ql(npro,nshl,9), dxidx(npro,nsd,nsd), 135 & gradul , dxidx, WdetJ ) 174 & ql(npro,nshl,27), dxidx(npro,nsd,nsd), 210 & uhessl , dxidx, WdetJ ) 289 & gradient, dxidx, WdetJ ) argument 301 & dxidx(npro,nsd,nsd), WdetJ(npro) 326 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 328 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 330 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 332 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) [all …]
|
| H A D | filters.f | 18 & dxdxi(npro,nsd,nsd), dxidx(npro,nsd,nsd), 57 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 59 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 61 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 63 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 64 & + dxidx(:,1,2) * dxdxi(:,2,1) 65 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 66 dxidx(:,1,1) = dxidx(:,1,1) * tmp 67 dxidx(:,1,2) = dxidx(:,1,2) * tmp 68 dxidx(:,1,3) = dxidx(:,1,3) * tmp [all …]
|
| H A D | e3.f | 53 & aci(npro,3), dxidx(npro,nsd,nsd), 101 & g3yi, shg, dxidx, 111 & u3, dxidx, rLui, 189 & dxidx(npro,nsd,nsd), WdetJ(npro), local 234 & shg, dxidx, WdetJ, 244 call e3StabSclr (uMod, dxidx, tauS,
|
| H A D | e3qGradV.f | 32 & dxidx(npro,nsd,nsd), WdetJ(npro) 57 & dxidx, WdetJ )
|
| H A D | e3ivar.f | 4 & g3yi, shg, dxidx, argument 60 & shg(npro,nshl,nsd), dxidx(npro,nsd,nsd), 143 call e3metric( xl, shgl, dxidx, 317 & shg, dxidx, WdetJ, argument 331 & shg(npro,nshl,nsd), dxidx(npro,nsd,nsd), 378 call e3metric( xl, shgl, dxidx, 466 call e3uBar(rho, src, dxidx, rLui, rmu, uBar) 493 & shpfun, shg, yl, dxidx,
|
| H A D | e3q.f | 34 & dxidx(npro,nsd,nsd), WdetJ(npro), 76 & dxidx, WdetJ ) 227 & dxidx(npro,nsd,nsd), WdetJ(npro) 260 & dxidx, WdetJ )
|
| H A D | e3ql.f | 34 & dxidx(npro,nsd,nsd), WdetJ(npro), 64 & dxidx, WdetJ ) 198 dimension dxidx(npro,nsd,nsd), WdetJ(npro), 225 & gradT, dxidx, WdetJ )
|
| H A D | errsmooth.f | 147 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd) 184 call e3metric( xl, shdrv, dxidx,
|
| H A D | advLES.f | 21 & dxdxi(npro,nsd,nsd), dxidx(npro,nsd,nsd), 133 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 135 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 137 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 139 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 140 & + dxidx(:,1,2) * dxdxi(:,2,1) 141 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 142 dxidx(:,1,1) = dxidx(:,1,1) * tmp 143 dxidx(:,1,2) = dxidx(:,1,2) * tmp 144 dxidx(:,1,3) = dxidx(:,1,3) * tmp [all …]
|
| /phasta/phSolver/common/ |
| H A D | e3metric.f | 8 subroutine e3metric( xl, shgl, dxidx, argument 14 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd), local 38 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 40 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 42 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 44 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 45 & + dxidx(:,1,2) * dxdxi(:,2,1) 46 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 47 dxidx(:,1,1) = dxidx(:,1,1) * tmp 48 dxidx(:,1,2) = dxidx(:,1,2) * tmp [all …]
|
| H A D | getstrl.f | 10 & strnrm(npro,maxsh), dxidx(npro,nsd,nsd), 50 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 52 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 54 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 56 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 57 & + dxidx(:,1,2) * dxdxi(:,2,1) 58 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 59 dxidx(:,1,1) = dxidx(:,1,1) * tmp 60 dxidx(:,1,2) = dxidx(:,1,2) * tmp 61 dxidx(:,1,3) = dxidx(:,1,3) * tmp [all …]
|
| H A D | cmass.f | 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) [all …]
|
| H A D | asithf.f | 11 & strnrm(npro,maxsh), dxidx(npro,nsd,nsd), 56 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 58 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 60 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 62 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 63 & + dxidx(:,1,2) * dxdxi(:,2,1) 64 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 65 dxidx(:,1,1) = dxidx(:,1,1) * tmp 66 dxidx(:,1,2) = dxidx(:,1,2) * tmp 67 dxidx(:,1,3) = dxidx(:,1,3) * tmp [all …]
|
| H A D | bardmc.f | 682 & dxdxi(npro,nsd,nsd), dxidx(npro,nsd,nsd), 726 dxidx(:,1,1) = dxdxi(:,2,2) * dxdxi(:,3,3) 728 dxidx(:,1,2) = dxdxi(:,3,2) * dxdxi(:,1,3) 730 dxidx(:,1,3) = dxdxi(:,1,2) * dxdxi(:,2,3) 732 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 733 & + dxidx(:,1,2) * dxdxi(:,2,1) 734 & + dxidx(:,1,3) * dxdxi(:,3,1) ) 735 dxidx(:,1,1) = dxidx(:,1,1) * tmp 736 dxidx(:,1,2) = dxidx(:,1,2) * tmp 737 dxidx(:,1,3) = dxidx(:,1,3) * tmp [all …]
|
| H A D | errorsmooth.f | 146 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd) 183 call e3metric( xl, shdrv, dxidx,
|