Home
last modified time | relevance | path

Searched refs:tmp1 (Results 1 – 25 of 27) sorted by relevance

12

/phasta/phSolver/incompressible/
H A De3lhs.f46 & tmp1(npro), tmp2(npro),
59 tmp1 = tlW * rho
66 uBar(:,1) = tmp1 * uBar(:,1)
67 uBar(:,2) = tmp1 * uBar(:,2)
68 uBar(:,3) = tmp1 * uBar(:,3)
71 uBar(:,1) = tmp1 * u1(:)
72 uBar(:,2) = tmp1 * u2(:)
73 uBar(:,3) = tmp1 * u3(:)
87 tmp1 = tsFct * shpfun(:,aa) * shpfun(:,b)
88 tmp2 = tmp1 + t1(:,1) * shpfun(:,aa)
[all …]
H A Dftools.f336 real*8 tmp1, tmp2, tmp3, tmp4 local
350 tmp1 = 0
352 tmp1 = tmp1 + x(i,1) * y(i)
354 c(1) = tmp1
358 tmp1 = 0
361 tmp1 = tmp1 + x(i,1) * y(i)
364 c(1) = tmp1
369 tmp1 = 0
373 tmp1 = tmp1 + x(i,1) * y(i)
377 c(1) = tmp1
[all …]
H A Dgenlmass.f178 real*8 ac(nshg,ndof), res(nshg,4), tmp,tmp1 local
187 tmp1=flmpl*almi
195 tmp=gmass(n)*tmp1*rho(n)
202 tmp1=flmpr
205 rho=rho*gmass*tmp1 ! reuse rho
220 real*8 ac(nshg), res(nshg), tmp, tmp1 local
229 tmp1=flmpl*almi
237 tmp=gmass(n)*tmp1*rho(n)
242 tmp1=flmpr
244 rho=rho*gmass*tmp1 ! reuse rho
H A De3res.f55 real*8 tmp1(npro), tmp2(npro), tmp3(npro), local
103 tmp1 = rmu * ( g2yi(:,2) + g1yi(:,3) )
126 rGNa(:,1,2) = tmp1
128 rGNa(:,2,1) = tmp1
140 rGNa(:,1,2) = tmp1 - u1(:)*u2(:)*rho(:)
142 rGNa(:,2,1) = tmp1 - u1(:)*u2(:)*rho(:)
161 tmp1 = tauM * rLui(:,1)
165 rGNa(:,1,1) = rGNa(:,1,1) + tmp1 * u1
166 rGNa(:,1,2) = rGNa(:,1,2) + tmp1 * u2
167 rGNa(:,1,3) = rGNa(:,1,3) + tmp1 * u3
[all …]
H A De3stab.f401 & tmp1(npro), tmp2(npro), local
443 tmp1(:) = c1 * dxidx(:,1,1)+c2 *(dxidx(:,2,1)+dxidx(:,3,1))
446 gijd(:,1) = dxidx(:,1,1) * tmp1
450 tmp1(:) = c1 * dxidx(:,1,2)+c2 *(dxidx(:,2,2)+dxidx(:,3,2))
453 gijd(:,2) = dxidx(:,1,2) * tmp1
457 gijd(:,4) = dxidx(:,1,1) * tmp1
461 tmp1(:) = c1 * dxidx(:,1,3)+c2 *(dxidx(:,2,3)+dxidx(:,3,3))
464 gijd(:,3) = dxidx(:,1,3) * tmp1
468 gijd(:,5) = dxidx(:,1,2) * tmp1
472 gijd(:,6) = dxidx(:,1,1) * tmp1
[all …]
H A DlesSparse.f266 real*8 tmp1, tmp2, tmp3, pisave
281 tmp1 = 0
288 tmp1 = tmp1
305 q(i,1) = q(i,1) + tmp1
425 real*8 tmp1, tmp2, tmp3, tmp4, pisave
440 tmp1 = 0
449 tmp1 = tmp1
472 q(i,1) = q(i,1) + tmp1
H A De3source.f255 & dP2Inv, sixth, tmp, tmp1, p, dp, local
417 tmp1 = ( (1.0 + saCw3P6) * tmp ) ** sixth
418 fw = g * tmp1
419 fwp = gp * tmp1 * saCw3P6 * tmp
422 tmp1 = saCw1 * dP2Inv*SclrNN
424 srcRat(e)= prat - ep*fw * tmp1
429 if(fw .gt. zero) tmp=tmp-two*tmp1*fw
430 if(fwp .gt. zero) tmp=tmp-tmp1*SclrNN*fwp
H A Dgetdiff.f352 double precision epsInv, Rey, Ret, RetInv, tmp1, fmuKE
375 tmp1 = exp(-0.0165*Rey)
376 fmuKE = (1. -tmp1) ** 2 * (1.+20.5*RetInv) ! fmu of Lam-Bremhorst
H A De3ivar.f68 dimension tmp(npro), tmp1(npro), dkei(npro), dist2w(npro)
263 tmp1=shpfun(:,ia)*Qwt(lcsyst,intp)
265 & tmp1(:)*rLui(:,1)*rLui(:,1)
267 & tmp1(:)*rLui(:,2)*rLui(:,2)
269 & tmp1(:)*rLui(:,3)*rLui(:,3)
H A De3bvar.f104 dimension tmp1(npro)
662 tmp1 = tsFctvw * shpb(:,aa) * shpb(:,b)
666 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp1
667 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp1
668 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp1
H A Dfilters.f2204 tmp1 = MINVAL(cdelsq)
2207 call MPI_REDUCE (tmp1, tmp3, 1,MPI_DOUBLE_PRECISION,
2211 tmp1=tmp3
2215 write(34,*)lstep,tmp1,tmp2
2956 tmp1 = MINVAL(cdelsq)
2959 call MPI_REDUCE (tmp1, tmp3, 1,MPI_DOUBLE_PRECISION,
2963 tmp1=tmp3
2967 write(34,*)lstep,tmp1,tmp2
H A DadvLES.f3019 tmp1 = MINVAL(cdelsq)
3022 call MPI_REDUCE (tmp1, tmp3, 1,MPI_DOUBLE_PRECISION,
3026 tmp1=tmp3
3030 write(34,*)lstep,tmp1,tmp2
3776 tmp1 = MINVAL(cdelsq)
3779 call MPI_REDUCE (tmp1, tmp3, 1,MPI_DOUBLE_PRECISION,
3783 tmp1=tmp3
3787 write(34,*)lstep,tmp1,tmp2
/phasta/phSolver/common/
H A Dnewshape.cc187 double tmp1 =2*xi[k]*tmp0; in Pyr_eB() local
188 psi =0.125*(tmp1*tmp1-1)*(1+sign[m]*(2*xi[m]*tmp0))*(1-xi[2]); in Pyr_eB()
211 double tmp1 =2*xi[k]*tmp0; in dPeBdxi() local
213 dPsidxi =0.5*tmp1*(1+sign[m]*(2*xi[m]*tmp0)); in dPeBdxi()
215 dPsidxi =0.25*(tmp1*tmp1-1)*sign[m]; in dPeBdxi()
217 dPsidxi =0.25*(tmp1*tmp1)*(1+sign[m]*(2*xi[m]*tmp0))-0.125*(tmp1*tmp1-1); in dPeBdxi()
231 double tmp1 =2*xi[0]*tmp0; in Pyr_fB() local
234 psi =0.125*(1-tmp1*tmp1)*(1-tmp2*tmp2)*(1-xi[1]); in Pyr_fB()
238 double tmp1 =2*xi[m]*tmp0; in Pyr_fB() local
240 psi =0.125*(1+sign[k]*2*xi[k]*tmp0)*(1-tmp1*tmp1)*(1-xi[1]*xi[1]); in Pyr_fB()
[all …]
H A Delm3keps.f130 real*8 tmp1, tmp2 , tmp0 local
239 tmp1 = exp(-0.0165*Rey)
241 fmuKE = (1. -tmp1) ** 2 * (1.+20.5*RetInv) ! fmu of Lam-Bremhorst
252 fmuKE_k = 0.033*(1.-tmp1)*(1.+20.5*RetInv)*Rey_k*tmp1
253 & -20.5 *(1.-tmp1)**2 * Ret_k*RetInv**2
255 fmuKE_eps= -20.5*(1.-tmp1)**2* Ret_eps*RetInv**2
395 tmp1 = min( srcJac(e,1) * srcJac(e,4),
401 if ( tmp2 .gt. tmp1 ) then
402 tmp2 = sqrt(tmp1/tmp2)
H A Dsolvecon.f185 & h_prime(npro), tmp1(npro),
233 tmp1 =zero
261 tmp1(i)=-h_prime(i)*(sclr(i)-sclrtmp(i))*dtgl
268 v_lambdal1(:,i)=v_lambdal1(:,i)+ shape(:,i)*WdetJ*tmp1
H A Dgetdmc.f363 tmp1 = MINVAL(cdelsq)
366 call MPI_REDUCE (tmp1, tmp3, 1,MPI_DOUBLE_PRECISION,
370 tmp1=tmp3
374 write(34,*)lstep,tmp1,tmp2
H A Dbardmc.f529 tmp1 = MINVAL(cdelsq)
532 call MPI_REDUCE (tmp1, tmp3, 1,MPI_DOUBLE_PRECISION,
536 tmp1=tmp3
540 write(34,*)lstep,tmp1,tmp2
/phasta/phSolver/compressible/
H A Dshuffle.f11 dimension restmp (nshg,nflow), tmp(nshg), tmp1(nshg,nflow)
22 tmp1(:,:) = restmp(:,:)
24 restmp(:,i+1) = tmp1(:,i)
26 restmp(:,1) = tmp1(:,4)
H A Dsparseap.f23 real*8 tmp1, tmp2, tmp3, tmp4, tmp5 local
50 tmp1 = 0.
59 tmp1 = tmp1 + lhsK(1 ,k)*p(j,1)
87 q(i,1) = q(i,1) + tmp1
H A De3tau.f1402 & tmp1(npro), tmp2(npro),
1443 tmp1(:) = c1 * dxidx(:,1,1) + c2 *(dxidx(:,2,1)+dxidx(:,3,1))
1446 gijd(:,1) = dxidx(:,1,1) * tmp1
1450 tmp1(:) = c1 * dxidx(:,1,2) + c2 *(dxidx(:,2,2)+dxidx(:,3,2))
1453 gijd(:,2) = dxidx(:,1,1) * tmp1
1457 gijd(:,3) = dxidx(:,1,2) * tmp1
1461 tmp1(:) = c1 * dxidx(:,1,3) + c2 *(dxidx(:,2,3)+dxidx(:,3,3))
1464 gijd(:,4) = dxidx(:,1,1) * tmp1
1468 gijd(:,5) = dxidx(:,1,2) * tmp1
1472 gijd(:,6) = dxidx(:,1,3) * tmp1
H A De3ls.f86 real*8 rerrl(npro,nshl,6), tmp(npro), tmp1(npro) local
162 tmp1=shape(:,ia)*Qwt(lcsyst,intp)
164 & tmp1(:)*rLyi(:,1)*rLyi(:,1)
166 & tmp1(:)*rLyi(:,2)*rLyi(:,2)
168 & tmp1(:)*rLyi(:,3)*rLyi(:,3)
H A De3bdg.f27 dimension tmp1(npro), tmp2(npro)
331 tmp1 = shp(:,i) * WdetJ * almi/gami/alfi*dtgl
335 & + tmp1*(
2623 dimension tmp1(npro), tmp2(npro)
2718 tmp1 = shp(:,i) * WdetJ * almi/gami/alfi*dtgl
2720 BDiagtl(:,i) = BDiagtl(:,i) + tmp1*(
/phasta/phSolver/AMG/
H A Dramg_coarse.f35 real(kind=8) :: tmp1,tmp2,tmp3,tmp4 local
323 tmp1 = 0
332 tmp1 = tmp1 + min(amg_la(k),zero)
339 if ( abs(tmp1).gt.1e-5) tmp1 = tmp3/tmp1
343 amg_la(k) = amg_la(k)*tmp1
/phasta/phSolver/incompressible/e3source/
H A De3source.fRHSonly255 & dP2Inv, sixth, tmp, tmp1, p, dp,
379 tmp1 = ( (1.0 + saCw3P6) * tmp ) ** sixth
380 fw = g * tmp1
381 fwp = gp * tmp1 * saCw3P6 * tmp
384 tmp1 = saCw1 * dP2Inv*SclrNN
386 srcRat(e)= prat - ep*fw * tmp1
391 if(fw .gt. zero) tmp=tmp-two*tmp1*fw
392 if(fwp .gt. zero) tmp=tmp-tmp1*SclrNN*fwp
H A De3source.fimplicit255 & dP2Inv, sixth, tmp, tmp1, p, dp,
393 tmp1 = ( (1.0 + saCw3P6) * tmp ) ** sixth
394 fw = g * tmp1
395 fwp = gp * tmp1 * saCw3P6 * tmp
398 tmp1 = saCw1 * dP2Inv*SclrNN
400 srcRat(e)= prat - ep*fw * tmp1
405 if(fw .gt. zero) tmp=tmp-two*tmp1*fw
406 if(fwp .gt. zero) tmp=tmp-tmp1*SclrNN*fwp

12