| /phasta/phSolver/incompressible/ |
| H A D | ftools.f | 336 real*8 tmp1, tmp2, tmp3, tmp4 local
359 tmp2 = 0
362 tmp2 = tmp2 + x(i,2) * y(i)
365 c(2) = tmp2
370 tmp2 = 0
374 tmp2 = tmp2 + x(i,2) * y(i)
378 c(2) = tmp2
384 tmp2 = 0
389 tmp2 = tmp2 + x(i,2) * y(i)
394 c(2) = tmp2
[all …]
|
| H A D | e3res.f | 55 real*8 tmp1(npro), tmp2(npro), tmp3(npro), local
104 tmp2 = rmu * ( g3yi(:,3) + g2yi(:,4) )
130 rGNa(:,2,3) = tmp2
132 rGNa(:,3,2) = tmp2
144 rGNa(:,2,3) = tmp2 - u3(:)*u2(:)*rho(:)
146 rGNa(:,3,2) = tmp2 - u3(:)*u2(:)*rho(:)
162 tmp2 = tauM * rLui(:,2)
168 rGNa(:,2,1) = rGNa(:,2,1) + tmp2 * u1
169 rGNa(:,2,2) = rGNa(:,2,2) + tmp2 * u2
170 rGNa(:,2,3) = rGNa(:,2,3) + tmp2 * u3
[all …]
|
| H A D | e3lhs.f | 46 & tmp1(npro), tmp2(npro),
88 tmp2 = tmp1 + t1(:,1) * shpfun(:,aa)
94 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp2
95 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp2
96 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp2
120 tmp2 = tauBar * ( rLui(:,1) * shg(:,b,1)
124 t3(:,1) = t2(:,1) + tmp1 * u1 + tmp2 * rLui(:,1)
125 t3(:,2) = t2(:,2) + tmp1 * u2 + tmp2 * rLui(:,2)
126 t3(:,3) = t2(:,3) + tmp1 * u3 + tmp2 * rLui(:,3)
292 dimension t1(npro,3), tmp1(npro), tmp2(npro),
[all …]
|
| H A D | e3stab.f | 401 & tmp1(npro), tmp2(npro), local
444 tmp2(:) = c1 * dxidx(:,2,1)+c2 *(dxidx(:,1,1)+dxidx(:,3,1))
447 1 + dxidx(:,2,1) * tmp2
451 tmp2(:) = c1 * dxidx(:,2,2)+c2 *(dxidx(:,1,2)+dxidx(:,3,2))
454 1 + dxidx(:,2,2) * tmp2
458 1 + dxidx(:,2,1) * tmp2
462 tmp2(:) = c1 * dxidx(:,2,3)+c2 *(dxidx(:,1,3)+dxidx(:,3,3))
465 1 + dxidx(:,2,3) * tmp2
469 1 + dxidx(:,2,2) * tmp2
473 1 + dxidx(:,2,1) * tmp2
[all …]
|
| H A D | lesSparse.f | 266 real*8 tmp1, tmp2, tmp3, pisave
282 tmp2 = 0
292 tmp2 = tmp2
306 q(i,2) = q(i,2) + tmp2
425 real*8 tmp1, tmp2, tmp3, tmp4, pisave
441 tmp2 = 0
453 tmp2 = tmp2
473 q(i,2) = q(i,2) + tmp2
|
| H A D | filters.f | 2205 tmp2 = MAXVAL(cdelsq)
2209 call MPI_REDUCE (tmp2, tmp4, 1, MPI_DOUBLE_PRECISION,
2212 tmp2=tmp4
2215 write(34,*)lstep,tmp1,tmp2
2957 tmp2 = MAXVAL(cdelsq)
2961 call MPI_REDUCE (tmp2, tmp4, 1, MPI_DOUBLE_PRECISION,
2964 tmp2=tmp4
2967 write(34,*)lstep,tmp1,tmp2
|
| H A D | advLES.f | 3020 tmp2 = MAXVAL(cdelsq)
3024 call MPI_REDUCE (tmp2, tmp4, 1, MPI_DOUBLE_PRECISION,
3027 tmp2=tmp4
3030 write(34,*)lstep,tmp1,tmp2
3777 tmp2 = MAXVAL(cdelsq)
3781 call MPI_REDUCE (tmp2, tmp4, 1, MPI_DOUBLE_PRECISION,
3784 tmp2=tmp4
3787 write(34,*)lstep,tmp1,tmp2
|
| /phasta/phSolver/compressible/ |
| H A D | sparseap.f | 23 real*8 tmp1, tmp2, tmp3, tmp4, tmp5 local
51 tmp2 = 0.
64 tmp2 = tmp2 + lhsK(2 ,k)*p(j,1)
88 q(i,2) = q(i,2) + tmp2
|
| H A D | e3bdg.f | 27 dimension tmp1(npro), tmp2(npro)
227 tmp2 = (shp(:,j)*shp(:,j)) * tmp
229 BDiagl(:,j,2,2) = BDiagl(:,j,2,2) + tmp2 * A0(:,2,2)
230 BDiagl(:,j,3,3) = BDiagl(:,j,3,3) + tmp2 * A0(:,3,3)
231 BDiagl(:,j,4,4) = BDiagl(:,j,4,4) + tmp2 * A0(:,4,4)
232 BDiagl(:,j,1,5) = BDiagl(:,j,1,5) + tmp2 * A0(:,1,5)
233 BDiagl(:,j,2,5) = BDiagl(:,j,2,5) + tmp2 * A0(:,2,5)
234 BDiagl(:,j,3,5) = BDiagl(:,j,3,5) + tmp2 * A0(:,3,5)
235 BDiagl(:,j,4,5) = BDiagl(:,j,4,5) + tmp2 * A0(:,4,5)
236 BDiagl(:,j,1,1) = BDiagl(:,j,1,1) + tmp2 * A0(:,1,1)
[all …]
|
| H A D | e3bdg_nd.f | 28 dimension tmp1(npro), tmp2(npro)
228 tmp2 = (shp(:,j)*shp(:,j)) * tmp
230 BDiagl(:,j,2,2) = BDiagl(:,j,2,2) + tmp2 * A0(:,2,2)
231 BDiagl(:,j,3,3) = BDiagl(:,j,3,3) + tmp2 * A0(:,3,3)
232 BDiagl(:,j,4,4) = BDiagl(:,j,4,4) + tmp2 * A0(:,4,4)
233 BDiagl(:,j,1,5) = BDiagl(:,j,1,5) + tmp2 * A0(:,1,5)
234 BDiagl(:,j,2,5) = BDiagl(:,j,2,5) + tmp2 * A0(:,2,5)
235 BDiagl(:,j,3,5) = BDiagl(:,j,3,5) + tmp2 * A0(:,3,5)
236 BDiagl(:,j,4,5) = BDiagl(:,j,4,5) + tmp2 * A0(:,4,5)
237 BDiagl(:,j,1,1) = BDiagl(:,j,1,1) + tmp2 * A0(:,1,1)
[all …]
|
| H A D | e3tau.f | 1402 & tmp1(npro), tmp2(npro),
1444 tmp2(:) = c1 * dxidx(:,2,1) + c2 *(dxidx(:,1,1)+dxidx(:,3,1))
1447 1 + dxidx(:,2,1) * tmp2
1451 tmp2(:) = c1 * dxidx(:,2,2) + c2 *(dxidx(:,1,2)+dxidx(:,3,2))
1454 1 + dxidx(:,2,1) * tmp2
1458 1 + dxidx(:,2,2) * tmp2
1462 tmp2(:) = c1 * dxidx(:,2,3) + c2 *(dxidx(:,1,3)+dxidx(:,3,3))
1465 1 + dxidx(:,2,1) * tmp2
1469 1 + dxidx(:,2,2) * tmp2
1473 1 + dxidx(:,2,3) * tmp2
|
| /phasta/AcuStat/ |
| H A D | create_AcuStatInput.sh | 27 cat tmp | awk -F dat '{print $2}' | sed 's/^\.//g' | sed 's/\.1//g' | sort -n > tmp2
46 done < tmp2
47 rm tmp2
|
| /phasta/phSolver/common/ |
| H A D | elm3keps.f | 130 real*8 tmp1, tmp2 , tmp0 local
398 tmp2 = srcJac(e,2) * srcJac(e,3)
401 if ( tmp2 .gt. tmp1 ) then
402 tmp2 = sqrt(tmp1/tmp2)
403 srcJac(e,2) = tmp2 * srcJac(e,2)
404 srcJac(e,3) = tmp2 * srcJac(e,3)
|
| H A D | newshape.cc | 232 double tmp2 =2*xi[2]*tmp0; in Pyr_fB() local
234 psi =0.125*(1-tmp1*tmp1)*(1-tmp2*tmp2)*(1-xi[1]); in Pyr_fB()
253 double tmp2 =2*xi[2]*tmp0; in dPfBdxi() local
257 dPsidxi =-0.5*tmp1*(1-tmp2*tmp2); in dPfBdxi()
260 dPsidxi = -0.125* (1+tmp1*tmp1+tmp2*tmp2-3*tmp1*tmp1*tmp2*tmp2); in dPfBdxi()
263 dPsidxi =-0.5*(1-tmp1*tmp1)*tmp2; in dPfBdxi()
276 double tmp2 =xi[m]*tmp0; in dPfBdxi() local
277 dPsidxi =(1+xi[1])*(0.25*sign[k]*tmp1*(1-4*tmp2*tmp2)- in dPfBdxi()
278 (1+2*sign[k]*tmp1)*tmp2*tmp2)- in dPfBdxi()
279 0.25*(1+2*sign[k]*tmp1)*(1-4*tmp2*tmp2)*xi[1]; in dPfBdxi()
|
| H A D | solvecon.f | 186 & tmp2(npro),
234 tmp2 =zero
262 tmp2(i)=h_prime(i)**2
269 v_lambdal2(:,i)=v_lambdal2(:,i)+ shape(:,i)*WdetJ*tmp2
|
| H A D | getdmc.f | 364 tmp2 = MAXVAL(cdelsq)
368 call MPI_REDUCE (tmp2, tmp4, 1, MPI_DOUBLE_PRECISION,
371 tmp2=tmp4
374 write(34,*)lstep,tmp1,tmp2
|
| H A D | bardmc.f | 530 tmp2 = MAXVAL(cdelsq)
534 call MPI_REDUCE (tmp2, tmp4, 1, MPI_DOUBLE_PRECISION,
537 tmp2=tmp4
540 write(34,*)lstep,tmp1,tmp2
|
| /phasta/phSolver/AMG/ |
| H A D | ramg_coarse.f | 35 real(kind=8) :: tmp1,tmp2,tmp3,tmp4 local
324 tmp2 = 0
333 tmp2 = tmp2 + max(amg_la(k),zero)
340 if ( abs(tmp2).gt.1e-5) tmp2 = tmp4/tmp2
345 amg_la(k) = amg_la(k)*tmp2
|