| /phasta/phSolver/compressible/ |
| H A D | e3mtrx.f | 4 & u1, u2, u3, 56 & u3(npro), fact1(npro), 116 A0(:,4,1) = drdp * u3 120 A0(:,4,5) = drdT * u3 125 A0(:,5,4) = rho * u3 150 A1(:,4,1) = drdp * u1 * u3 151 A1(:,4,2) = rho * u3 154 A1(:,4,5) = drdT * u1 * u3 159 A1(:,5,4) = rho * u1 * u3 182 A2(:,4,1) = drdp * u2 * u3 [all …]
|
| H A D | e3visc.f | 4 & u1, u2, u3, 46 & u3(npro), rho(npro), 74 stiff(:, 5, 4) = rmu * u3 88 stiff(:, 5,12) = rmu * u3 107 stiff(:,10, 9) = rmu * u3 114 stiff(:,10,13) = rmu * u3 121 stiff(:,15, 2) = rlm * u3 128 stiff(:,15, 8) = rlm * u3 138 stiff(:,15,14) = rlm2mu * u3 243 & + f3 * T * rmu * u1 + f4 * T * (rlm + rmu) * u3 [all …]
|
| H A D | e3tau.f | 2 & u1, u2, u3, 41 & u2(npro), u3(npro), 69 & + u3*u3*(gijd(:,1)*gijd(:,3)-gijd(:,2)*gijd(:,2)) 71 & + u1*u3*(gijd(:,2)*gijd(:,5)-gijd(:,4)*gijd(:,3)) 72 & + u2*u3*(gijd(:,4)*gijd(:,2)-gijd(:,1)*gijd(:,5))) 84 & + u3*u3*gijd(:,6) 86 & + u1*u3*gijd(:,4) 87 & + u2*u3*gijd(:,5))*two + 1.0e-15 !FIX FOR INVALID MESHES 163 & + (u1*(u1*gijd(:,1) + two*(u2*gijd(:,2)+u3*gijd(:,4))) 164 & + u2*(u2*gijd(:,3) + two*u3*gijd(:,5)) [all …]
|
| H A D | e3conv.f | 5 & u2, u3, rLyi, argument 56 & u2(npro), u3(npro), 79 ri(:, 4) = (- u1) * rho * u3 85 ri(:, 9) = (- u2) * rho * u3 88 ri(:,11) = (- u3) * rho 89 ri(:,12) = (- u3) * rho * u1 90 ri(:,13) = (- u3) * rho * u2 91 ri(:,14) = (- u3) * rho * u3 - pres 92 ri(:,15) = (- u3) * rho * (ei + rk) - u3 * pres 359 & u2, u3, rLyti, argument [all …]
|
| H A D | e3ql.f | 37 & u3(npro), rmu(npro), 79 & u3 ) 103 & + rmu * u3 * g1yi(:,4) 105 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 120 & + rmu * u3 * g2yi(:,4) 121 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 133 qdi(:,12) = rlm * u3 * g1yi(:,2) + rmu * u1 * g1yi(:,4) 134 & + rlm * u3 * g2yi(:,3) + rmu * u2 * g2yi(:,4) 136 & + rlm2mu * u3 * g3yi(:,4)
|
| H A D | e3.f | 75 & u3(npro), A0DC(npro,4), 137 & u1, u2, u3, 152 & u1, u2, u3, 167 & u2, u3, rLyi, 180 & u1, u2, u3, 191 & u3, pres, sforce, 204 & u1, u2, u3, 269 & u1, u2, u3, 278 & u1, u2, u3, 382 & u3(npro) [all …]
|
| H A D | e3q.f | 38 & u3(npro), rmu(npro), 85 & u3 ) 112 & + rmu * u3 * g1yi(:,4) 114 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 129 & + rmu * u3 * g2yi(:,4) 130 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 142 qdi(:,12) = rlm * u3 * g1yi(:,2) + rmu * u1 * g1yi(:,4) 143 & + rlm * u3 * g2yi(:,3) + rmu * u2 * g2yi(:,4) 145 & + rlm2mu * u3 * g3yi(:,4)
|
| H A D | e3qvar.f | 4 & cp, u1, u2, u3) argument 48 & u3(npro), pres(npro) 63 u3 = zero 75 u3 = u3 + shp(:,n) * ycl(:,n,4)
|
| H A D | e3b.f | 72 & u3(npro), 132 & u1, u2, u3, 153 un = bnorm(:,1) * u1 + bnorm(:,2) * u2 + bnorm(:,3) * u3 170 F4 = rou * u3 + bnorm(:,3) * p 217 Fv5 = u1 * Fv2 + u2 * Fv3 + u3 * Fv4 438 & u3(npro), 486 & u2, u3, g1yti, 547 u3 = mytmp * g3yti 558 un = bnorm(:,1)*u1 + bnorm(:,2)*u2 + bnorm(:,3)*u3
|
| H A D | e3bvar.f | 4 & u1, u2, u3, rho, ei, argument 63 & u3(npro), rho(npro), 92 u3 = zero 101 u3 = u3 + shpb(:,nodlcl) * yl(:,nodlcl,4) 107 rk = pt5 * ( u1**2 + u2**2 + u3**2 ) 378 & u1, u2, u3, argument 423 & u3(npro), rho(npro), 448 u3 = zero 461 u3 = u3 + shpb(:,nodlcl) * ycl(:,nodlcl,4) 468 rk = pt5 * ( u1**2 + u2**2 + u3**2 )
|
| H A D | e3ivar.f | 9 & u1, u2, u3, argument 78 & u3(npro), divqi(npro,nflow), 145 u3 = zero 154 u3 = u3 + shape(:,n) * ycl(:,n,4) 206 rk = pt5 * ( u1**2 + u2**2 + u3**2 ) 498 & u3, shg, dwl, argument 568 & u3(npro), divqi(npro,nflow-1), 591 u3 = zero 604 u3 = u3 + shape(:,n) * ycl( :,n,4) 637 rk = pt5 * ( u1**2 + u2**2 + u3**2 )
|
| H A D | e3source.f | 3 & u3, pres, sforce, 30 & u3(npro), rho(npro), 54 src(:,5) = u1*src(:,2) + u2*src(:,3) + u3*src(:,4) 189 & + u3*sforce(:,3))*rho(:) 230 & u2, u3, xl, elDwl) 265 & u2 (npro), u3 (npro), 486 u3 = mytmp(:) * g3yti(:)
|
| H A D | e3ls.f | 4 & u1, u2, u3, argument 66 & u3(npro), rk(npro), 191 & u1, u2, u3, 209 & u1, u2, u3, 777 & u1, u2, u3, argument 835 & u3(npro), rk(npro), 892 & u1, u2, u3,
|
| H A D | e3eig1.f | 2 & u1, u2, u3, a1, a2, 43 & u2(npro), u3(npro), 60 Rcs(:,1,3) = u3 * fact
|
| /phasta/phSolver/incompressible/ |
| H A D | e3stab.f | 2 & u3, dxidx, rLui, argument 31 & u2(npro), u3(npro), 80 5 + gijd(:,6) * u3 ) 83 8 + gijd(:,5) * u3 ) 84 9 + u3 * ( gijd(:,6) * u1 86 1 + gijd(:,3) * u3 ) ) ) 122 3 + two * (gijd(:,4)*gijd(:,5)-gijd(:,6)*gijd(:,2)) * u3) 124 3 + two * (gijd(:,4)*gijd(:,6)-gijd(:,1)*gijd(:,5)) * u3) 125 1 + u3 * ( (gijd(:,1)*gijd(:,2)-gijd(:,4)*gijd(:,4)) * u3) 144 5 + gijd(:,6) * u3 ) [all …]
|
| H A D | e3res.f | 1 subroutine e3Res ( u1, u2, u3, argument 41 dimension u1(npro), u2(npro), u3(npro), 141 rGNa(:,1,3) = tmp3 - u1(:)*u3(:)*rho(:) 144 rGNa(:,2,3) = tmp2 - u3(:)*u2(:)*rho(:) 145 rGNa(:,3,1) = tmp3 - u1(:)*u3(:)*rho(:) 146 rGNa(:,3,2) = tmp2 - u3(:)*u2(:)*rho(:) 147 rGNa(:,3,3) = two * rmu * g3yi(:,4) + tmp - u3(:)*u3(:)*rho(:) 167 rGNa(:,1,3) = rGNa(:,1,3) + tmp1 * u3 170 rGNa(:,2,3) = rGNa(:,2,3) + tmp2 * u3 173 rGNa(:,3,3) = rGNa(:,3,3) + tmp3 * u3 [all …]
|
| H A D | stats.f | 217 real*8 u1, u2, u3, r0, r1, r2, r3, r4, r5, t3, t4, t5 local 252 u3 = yAlpha(i,3) 259 r2 = res(3) + reg * u3 270 stsVelReg(i,3) = stsVelReg(i,3) + u3 274 r2 = res(3) * u3 + reg * u3 * u3 276 r4 = res(2) * u3 + res(3) * u2 + reg * u3 * u2 277 r5 = res(3) * u1 + res(1) * u3 + reg * u1 * u3 337 u3 = stsVec(i,3) 344 stsVel(i,3) = stsVel(i,3) + u3 348 stsVelSqr(i,3) = stsVelSqr(i,3) + u3*u3 [all …]
|
| H A D | e3ivar.f | 6 & u1, u2, u3, argument 64 & u3(npro), divqi(npro,nflow-1+isurf), 84 u3 = zero 90 u3 = u3 + shpfun(:,n) * yl(:,n,4) 246 & aci, u1, u2, u3, Temp, rho, xx, 318 & u1, u2, u3, argument 334 & u3(npro), divS(npro), 356 u3 = zero 363 u3 = u3 + shpfun(:,n) * yl(:,n,4) 458 & aci, u1, u2, u3, Temp, rho, x, [all …]
|
| H A D | e3.f | 56 & u2(npro), u3(npro), 103 & u1, u2, u3, 111 & u3, dxidx, rLui, 117 call e3Res ( u1, u2, u3, 129 call e3LHS ( u1, u2, u3, 190 & u1(npro), u2(npro), u3(npro), local 235 & u1, u2, u3,
|
| H A D | filters.f | 14 & u3(npro), 108 u3 = zero 118 u3 = u3 + shp(i,intp)*yl(:,i,4) 138 fresli(:,3) = WdetJ * u3 !G * u3 * WdetJ 142 fresli(:,6) = WdetJ * u3 * u3 ! G*u3*u3*WdetJ 144 fresli(:,8) = WdetJ * u1 * u3 ! G*u1*u3*WdetJ 145 fresli(:,9) = WdetJ * u2 * u3 ! G*u2*u3*WdetJ 187 & u3(npro), dxdxi(npro,nsd,nsd), 413 & u3(npro), 508 u3 = zero [all …]
|
| H A D | getdiff.f | 217 u3=yl(e,n,4) 218 if((u1.eq.zero).and.(u2.eq.zero).and.(u3.eq.zero)) 303 u3=yl(e,n,4) 304 if((u1.eq.zero).and.(u2.eq.zero).and.(u3.eq.zero))
|
| H A D | e3source.f | 213 & u1, u2, u3, xl, argument 244 & u2(npro), u3(npro), local 448 uMod(:,3) = u3(:) - tmp * gradS(:,3) 508 uMod(:,3) = u3(:) - zero 579 u3 = umod(:,3)
|
| H A D | e3lhs.f | 1 subroutine e3LHS ( u1, u2, u3, argument 35 dimension u1(npro), u2(npro), u3(npro), 73 uBar(:,3) = tmp1 * u3(:) 117 & + u3 * shg(:,b,3) )*rho 126 t3(:,3) = t2(:,3) + tmp1 * u3 + tmp2 * rLui(:,3)
|
| /phasta/phSolver/common/ |
| H A D | asithf.f | 10 & u3(npro), dxdxi(npro,nsd,nsd), 153 u3=zero 157 u3 = u3 + shp(i,intp)*ycl(:,i,4) 162 fresli(:,3) = fresli(:,22) * u3 !rho u3 * WdetJ 166 fresli(:,6) = fresli(:,3) * u3 !rho u3 u3 *WdetJ 168 fresli(:,8) = fresli(:,1) * u3 !rho u1 u3 *WdetJ 169 fresli(:,9) = fresli(:,2) * u3 !rho u2 u3 *WdetJ
|
| H A D | bardmc.f | 681 & u3(npro), rho(npro), 763 u3=zero 767 u3 = u3 + shp(i,intp)*yl(:,i,4) 775 fresli(:,i,3) = rho * u3 * shp(i,intp) !rho * u3 * WdetJ 781 fresli(:,i,8) = rho * u3 * u3 * shp(i,intp) !rho * u3 * u3 * WdetJ 783 fresli(:,i,10)= rho * u1 * u3 * shp(i,intp) !rho * u1 * u3 * WdetJ 784 fresli(:,i,11)= rho * u2 * u3 * shp(i,intp) !rho * u2 * u3 * WdetJ 835 & u3(npro), dxdxi(npro,nsd,nsd), 1001 u3=zero 1005 u3 = u3 + shp(i,intp)*yl(:,i,4) [all …]
|