| /phasta/phSolver/compressible/ |
| H A D | e3eig1.f | 49 dimension Rcs(npro,nsd,nsd), fact(npro),
57 fact = one / u
58 Rcs(:,1,1) = u1 * fact
59 Rcs(:,1,2) = u2 * fact
60 Rcs(:,1,3) = u3 * fact
71 fact = a1 * Rcs(:,1,1) + a2 * Rcs(:,1,2) + a3 * Rcs(:,1,3)
72 a1 = a1 - fact * Rcs(:,1,1)
73 a2 = a2 - fact * Rcs(:,1,2)
74 a3 = a3 - fact * Rcs(:,1,3)
75 fact = a1**2 + a2**2 + a3**2
[all …]
|
| H A D | e3massl.f | 24 dimension fact(npro,5), temp(npro),
44 fact(:,1) = shpij * (temp + bcool*spongeContinuity)
45 fact(:,2) = shpij * (temp + bcool*spongeMomentum1)
46 fact(:,3) = shpij * (temp + bcool*spongeMomentum2)
47 fact(:,4) = shpij * (temp + bcool*spongeMomentum3)
48 fact(:,5) = shpij * (temp + bcool*spongeEnergy)
54 & + fact(:,1) * A0(:,1,jdof)
56 & + fact(:,2) * A0(:,2,jdof)
58 & + fact(:,3) * A0(:,3,jdof)
60 & + fact(:,4) * A0(:,4,jdof)
[all …]
|
| H A D | e3juel.f | 31 & ub(npro,nflow), fact(npro),
50 fact=WdetJ/(Qwt(lcsyst,intp)*15.0d0)
53 fddt=fact*fct1
104 yl(:,i,1) = fact*(acl(:,i,1)+ub(:,1))
105 yl(:,i,2) = fact*(acl(:,i,2)+ub(:,2))
106 yl(:,i,3) = fact*(acl(:,i,3)+ub(:,3))
107 yl(:,i,4) = fact*(acl(:,i,4)+ub(:,4))
108 yl(:,i,5) = fact*(acl(:,i,5)+ub(:,5))
|
| H A D | e3dc.f | 261 fact = one
262 if (ipord .eq. 2) fact = 0.9
263 if (ipord .eq. 3) fact = 0.75
276 DC(:,intp)=max(zero,(fact*sqrt(raLS*gnorm))-(rtLS*gnorm))
310 fact = one
311 if (ipord .eq. 2) fact = pt5
324 DC(:,intp) = min( max(zero,fact * sqrt(raLS * gnorm)-
576 fact = one
577 if (ipord .eq. 2) fact = 0.9
578 if (ipord .eq. 3) fact = 0.75
[all …]
|
| H A D | e3ls.f | 81 & A3tauA0(npro,nflow,nflow), fact(npro),
735 fact = shape(:,j) * WdetJ * almi/gami/alfi*dtgl
743 & fact * Atau(:,idof,1)
745 & fact * Atau(:,idof,2)
747 & fact * Atau(:,idof,3)
749 & fact * Atau(:,idof,4)
751 & fact * Atau(:,idof,5)
849 & A3tautA0(npro), fact(npro)
1054 fact = shape(:,jb) * WdetJ
1056 EGmasst(:,ia,jb) = EGmasst(:,ia,jb) + fact * Ataut(:)
|
| H A D | e3tau.f | 50 & fact(npro), h2o2u(npro), giju(npro,6),
61 fact = gijd(:,1) * gijd(:,3) * gijd(:,6)
78 uh1= two * sqrt(uh1/fact)
101 fact=pt5*rho*uh1/rmu
106 tau(:,1)=pt5*uh1*min(one,fact)*taucfct
135 fact=h2o2u*h2o2u*rk*celt/rmu
137 tau(:,2)=min(tau(:,2),fact)
138 tau(:,3)=min(tau(:,3),fact*rmu/con)*temper
168 fact=sqrt(tau(:,2))
169 tau(:,1)=pt125*fact/(rho*(gijd(:,1)+gijd(:,3)+gijd(:,6)))*taucfct
[all …]
|
| H A D | getthm.f | 131 fact = datmat(1,1,2)/datmat(1,1,1)
285 fact = datmat(1,1,2)/datmat(1,1,1)
|
| /phasta/phSolver/incompressible/ |
| H A D | e3stab.f | 38 dimension gijd(npro,6), fact(npro), rnu(npro),
97 fact = sqrt(tauM)
100 tauC =rho* pt125*fact/(gijd(:,1)+gijd(:,2)+gijd(:,3))*ff
101 tauM = one/fact
107 fact = gijd(:,1) * gijd(:,2) * gijd(:,3)
126 tauM=fact/taum ! here we have (u_i g^{ij} u^j)^{-1} approx 4/u^2h^2
141 fact =
159 taubar= dtsfsq/( dt*dt + .01*unorm/fact) ! never gets above (C_1 20*u_inf/h)^2
168 1 * ( four*taubar + fact
178 fact=sqrt(tauM)
[all …]
|
| H A D | e3dc.f | 539 fact = one
540 if (ipord .eq. 2) fact = 0.9
541 if (ipord .eq. 3) fact = 0.75
545 dcFct(:)=max(zero,(fact*sqrt(dcFct(:)))-(tauS(:)*dcFct(:)))
562 fact = one
563 if (ipord .eq. 2) fact = 0.9
565 dcFct(:) = min( max(zero, (fact * sqrt(dcFct(:)) -
|
| /phasta/phSolver/common/ |
| H A D | genscale.f | 315 fact = 1.0
373 ymapped(nen1(1:npin),5)= correction * fact*(ubar2(:,5)+fluc2(:,5))
|