| /phasta/phSolver/compressible/ |
| H A D | e3mtrx.f | 4 & u1, u2, u3, 55 & u1(npro), u2(npro), 104 A0(:,2,1) = drdp * u1 108 A0(:,2,5) = drdT * u1 123 A0(:,5,2) = rho * u1 132 A1(:,1,1) = drdp * u1 136 A1(:,1,5) = drdT * u1 138 A1(:,2,1) = drdp * u1 * u1 +1 139 A1(:,2,2) = two *rho * u1 142 A1(:,2,5) = drdT * u1 * u1 [all …]
|
| H A D | e3visc.f | 4 & u1, u2, u3, 45 & u1(npro), u2(npro), 72 stiff(:, 5, 2) = rlm2mu * u1 82 stiff(:, 5, 8) = rlm * u1 89 stiff(:, 5,14) = rlm * u1 97 stiff(:,10, 3) = rmu * u1 105 stiff(:,10, 7) = rmu * u1 122 stiff(:,15, 4) = rmu * u1 136 stiff(:,15,12) = rmu * u1 242 compK(:, 7) = f1 * T * rlm2mu * u1 + f2 * T * (rlm + rmu) * u2 [all …]
|
| H A D | e3conv.f | 4 & ei, rk, u1, argument 55 & rk(npro), u1(npro), 76 ri(:, 1) = (- u1) * rho 77 ri(:, 2) = (- u1) * rho * u1 - pres 78 ri(:, 3) = (- u1) * rho * u2 79 ri(:, 4) = (- u1) * rho * u3 80 ri(:, 5) = (- u1) * rho * (ei + rk) - u1 * pres 83 ri(:, 7) = (- u2) * rho * u1 89 ri(:,12) = (- u3) * rho * u1 358 & rho, u1, Sclr, argument [all …]
|
| H A D | e3tau.f | 2 & u1, u2, u3, 40 & cp(npro), u1(npro), 67 uh1= u1*u1*(gijd(:,3)*gijd(:,6)-gijd(:,5)*gijd(:,5)) 70 & + two *(u1*u2*(gijd(:,4)*gijd(:,5)-gijd(:,2)*gijd(:,6)) 71 & + u1*u3*(gijd(:,2)*gijd(:,5)-gijd(:,4)*gijd(:,3)) 82 h2o2u = u1*u1*gijd(:,1) 85 & +(u1*u2*gijd(:,2) 86 & + u1*u3*gijd(:,4) 163 & + (u1*(u1*gijd(:,1) + two*(u2*gijd(:,2)+u3*gijd(:,4))) 285 & u1, u2, u3, [all …]
|
| H A D | e3ql.f | 36 & u1(npro), u2(npro), 78 & cp, u1, u2, 102 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3) 104 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 105 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 118 qdi(:,8) = rlm * u2 * g1yi(:,2) + rmu * u1 * g1yi(:,3) 119 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 133 qdi(:,12) = rlm * u3 * g1yi(:,2) + rmu * u1 * g1yi(:,4) 135 & + rmu * u1 * g3yi(:,2) + rmu * u2 * g3yi(:,3)
|
| H A D | e3.f | 74 & u1(npro), u2(npro), 137 & u1, u2, u3, 152 & u1, u2, u3, 166 & ei, rk, u1, 180 & u1, u2, u3, 190 & rho, u1, u2, 204 & u1, u2, u3, 269 & u1, u2, u3, 278 & u1, u2, u3, 381 & u1(npro), u2(npro), [all …]
|
| H A D | e3q.f | 37 & u1(npro), u2(npro), 84 & cp, u1, u2, 111 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3) 113 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 114 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 127 qdi(:,8) = rlm * u2 * g1yi(:,2) + rmu * u1 * g1yi(:,3) 128 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 142 qdi(:,12) = rlm * u3 * g1yi(:,2) + rmu * u1 * g1yi(:,4) 144 & + rmu * u1 * g3yi(:,2) + rmu * u2 * g3yi(:,3)
|
| H A D | e3qvar.f | 4 & cp, u1, u2, u3) argument 47 & u1(npro), u2(npro), 61 u1 = zero 73 u1 = u1 + shp(:,n) * ycl(:,n,2)
|
| H A D | e3b.f | 71 & u1(npro), u2(npro), 132 & u1, u2, u3, 153 un = bnorm(:,1) * u1 + bnorm(:,2) * u2 + bnorm(:,3) * u3 168 F2 = rou * u1 + bnorm(:,1) * p 217 Fv5 = u1 * Fv2 + u2 * Fv3 + u3 * Fv4 437 dimension u1(npro), u2(npro), 485 & xlb, lnode, u1, 545 u1 = mytmp * g1yti 558 un = bnorm(:,1)*u1 + bnorm(:,2)*u2 + bnorm(:,3)*u3
|
| H A D | e3bvar.f | 4 & u1, u2, u3, rho, ei, argument 62 & u1(npro), u2(npro), 90 u1 = zero 99 u1 = u1 + shpb(:,nodlcl) * yl(:,nodlcl,2) 107 rk = pt5 * ( u1**2 + u2**2 + u3**2 ) 378 & u1, u2, u3, argument 422 & u1(npro), u2(npro), 446 u1 = zero 459 u1 = u1 + shpb(:,nodlcl) * ycl(:,nodlcl,2) 468 rk = pt5 * ( u1**2 + u2**2 + u3**2 )
|
| H A D | e3ivar.f | 9 & u1, u2, u3, argument 77 & u1(npro), u2(npro), 143 u1 = zero 152 u1 = u1 + shape(:,n) * ycl(:,n,2) 206 rk = pt5 * ( u1**2 + u2**2 + u3**2 ) 497 & rk, u1, u2, argument 567 & u1(npro), u2(npro), 589 u1 = zero 602 u1 = u1 + shape(:,n) * ycl( :,n,2) 637 rk = pt5 * ( u1**2 + u2**2 + u3**2 )
|
| H A D | e3source.f | 2 & rho, u1, u2, 29 & u1(npro), u2(npro), 54 src(:,5) = u1*src(:,2) + u2*src(:,3) + u3*src(:,4) 75 src(:,5) = u1*src(:,2) 188 src(:,5) = src(:,5) + (u1*sforce(:,1)+u2*sforce(:,2) 229 & ycl, shape, u1, 264 & g3yti (npro), u1 (npro), 484 u1 = mytmp(:) * g1yti(:)
|
| H A D | e3ls.f | 4 & u1, u2, u3, argument 65 & u1(npro), u2(npro), 191 & u1, u2, u3, 209 & u1, u2, u3, 777 & u1, u2, u3, argument 834 & u1(npro), u2(npro), 892 & u1, u2, u3,
|
| H A D | e3eig1.f | 2 & u1, u2, u3, a1, a2, 42 & c(npro), u1(npro), 58 Rcs(:,1,1) = u1 * fact
|
| /phasta/phSolver/incompressible/ |
| H A D | e3stab.f | 1 subroutine e3stab (rho, u1, u2, argument 30 dimension rho(npro), u1(npro), 78 3 + ( u1 * ( gijd(:,1) * u1 81 6 + u2 * ( gijd(:,4) * u1 84 9 + u3 * ( gijd(:,6) * u1 120 1 u1 * ( (gijd(:,2)*gijd(:,3)-gijd(:,5)*gijd(:,5)) * u1 142 3 u1 * ( gijd(:,1) * u1 145 6 + u2 * ( gijd(:,4) * u1 148 9 + u3 * ( gijd(:,6) * u1 156 unorm = (u1*u1+u2*u2+u3*u3)/velsq [all …]
|
| H A D | e3res.f | 1 subroutine e3Res ( u1, u2, u3, argument 41 dimension u1(npro), u2(npro), u3(npro), 139 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp - u1(:)*u1(:)*rho(:) 140 rGNa(:,1,2) = tmp1 - u1(:)*u2(:)*rho(:) 141 rGNa(:,1,3) = tmp3 - u1(:)*u3(:)*rho(:) 142 rGNa(:,2,1) = tmp1 - u1(:)*u2(:)*rho(:) 145 rGNa(:,3,1) = tmp3 - u1(:)*u3(:)*rho(:) 165 rGNa(:,1,1) = rGNa(:,1,1) + tmp1 * u1 168 rGNa(:,2,1) = rGNa(:,2,1) + tmp2 * u1 171 rGNa(:,3,1) = rGNa(:,3,1) + tmp3 * u1 [all …]
|
| H A D | e3ivar.f | 6 & u1, u2, u3, argument 63 & u1(npro), u2(npro), 82 u1 = zero 88 u1 = u1 + shpfun(:,n) * yl(:,n,2) 246 & aci, u1, u2, u3, Temp, rho, xx, 294 dkei = (u1(:)-dkei)**2 +u2(:)**2 ! u'^2+v'^2 318 & u1, u2, u3, argument 333 & u1(npro), u2(npro), 354 u1 = zero 361 u1 = u1 + shpfun(:,n) * yl(:,n,2) [all …]
|
| H A D | stats.f | 217 real*8 u1, u2, u3, r0, r1, r2, r3, r4, r5, t3, t4, t5 local 250 u1 = yAlpha(i,1) 257 r0 = res(1) + reg * u1 268 stsVelReg(i,1) = stsVelReg(i,1) + u1 272 r0 = res(1) * u1 + reg * u1 * u1 275 r3 = res(1) * u2 + res(2) * u1 + reg * u2 * u1 277 r5 = res(3) * u1 + res(1) * u3 + reg * u1 * u3 335 u1 = stsVec(i,1) 342 stsVel(i,1) = stsVel(i,1) + u1 346 stsVelSqr(i,1) = stsVelSqr(i,1) + u1*u1 [all …]
|
| H A D | e3.f | 55 & pres(npro), u1(npro), 103 & u1, u2, u3, 110 call e3stab (rho, u1, u2, 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 | 13 & u1(npro), u2(npro), 106 u1 = zero 116 u1 = u1 + shp(i,intp)*yl(:,i,2) 136 fresli(:,1) = WdetJ * u1 !G * u1 * WdetJ 140 fresli(:,4) = WdetJ * u1 * u1 ! G*u1*u1*WdetJ 143 fresli(:,7) = WdetJ * u1 * u2 ! G*u1*u2*WdetJ 144 fresli(:,8) = WdetJ * u1 * u3 ! G*u1*u3*WdetJ 186 & u1(npro), u2(npro), 412 & u1(npro), u2(npro), 506 u1 = zero [all …]
|
| H A D | getdiff.f | 215 u1=yl(e,n,2) 218 if((u1.eq.zero).and.(u2.eq.zero).and.(u3.eq.zero)) 301 u1=yl(e,n,2) 304 if((u1.eq.zero).and.(u2.eq.zero).and.(u3.eq.zero))
|
| H A D | e3source.f | 213 & u1, u2, u3, xl, argument 243 & rmu(npro), u1(npro), local 446 uMod(:,1) = u1(:) - tmp * gradS(:,1) 506 uMod(:,1) = u1(:) - zero 577 u1 = umod(:,1) ! These are for the RHS
|
| H A D | e3lhs.f | 1 subroutine e3LHS ( u1, u2, u3, argument 35 dimension u1(npro), u2(npro), u3(npro), 71 uBar(:,1) = tmp1 * u1(:) 115 tmp1 = tauM * ( u1 * shg(:,b,1) 124 t3(:,1) = t2(:,1) + tmp1 * u1 + tmp2 * rLui(:,1)
|
| /phasta/phSolver/common/ |
| H A D | asithf.f | 9 & u1(npro), u2(npro), 151 u1=zero 155 u1 = u1 + shp(i,intp)*ycl(:,i,2) 160 fresli(:,1) = fresli(:,22) * u1 !rho u1 * WdetJ 164 fresli(:,4) = fresli(:,1) * u1 !rho u1 u1 *WdetJ
|
| H A D | bardmc.f | 680 & u1(npro), u2(npro), 761 u1=zero 765 u1 = u1 + shp(i,intp)*yl(:,i,2) 773 fresli(:,i,1) = rho * u1 * shp(i,intp) !rho * u1 * WdetJ 779 fresli(:,i,6) = rho * u1 * u1 * shp(i,intp) !rho * u1 * u1 * WdetJ 782 fresli(:,i,9) = rho * u1 * u2 * shp(i,intp) !rho * u1 * u2 * WdetJ 783 fresli(:,i,10)= rho * u1 * u3 * shp(i,intp) !rho * u1 * u3 * WdetJ 834 & u1(npro), u2(npro), 999 u1=zero 1003 u1 = u1 + shp(i,intp)*yl(:,i,2) [all …]
|