| /phasta/phSolver/common/ |
| H A D | genbc1.f | 29 dimension tmp(nshg) 46 tmp = BCtmp(:,4)**2 + BCtmp(:,5)**2 + BCtmp(:,6)**2 47 BC(:,3) = tmp * BCtmp(:,7) / BCtmp(:,4) 55 tmp = BCtmp(:,4)**2 + BCtmp(:,5)**2 + BCtmp(:,6)**2 56 BC(:,3) = tmp * BCtmp(:,7) / BCtmp(:,5) 76 tmp = sqrt (BCtmp(:, 4)**2 + BCtmp(:, 5)**2 78 BCtmp(:, 4) = BCtmp(:, 4) / tmp 79 BCtmp(:, 5) = BCtmp(:, 5) / tmp 80 BCtmp(:, 6) = BCtmp(:, 6) / tmp 81 BCtmp(:, 7) = BCtmp(:, 7) * tmp [all …]
|
| H A D | e3metric.f | 17 real*8 dxdxi(npro,nsd,nsd), tmp(npro) local 44 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 47 dxidx(:,1,1) = dxidx(:,1,1) * tmp 48 dxidx(:,1,2) = dxidx(:,1,2) * tmp 49 dxidx(:,1,3) = dxidx(:,1,3) * tmp 51 & - dxdxi(:,2,1) * dxdxi(:,3,3)) * tmp 53 & - dxdxi(:,3,1) * dxdxi(:,1,3)) * tmp 55 & - dxdxi(:,1,1) * dxdxi(:,2,3)) * tmp 57 & - dxdxi(:,2,2) * dxdxi(:,3,1)) * tmp 59 & - dxdxi(:,1,1) * dxdxi(:,3,2)) * tmp [all …]
|
| H A D | dtn.f | 23 real*8 BC(nshg,ndofBC),y(nshg,ndof),tmp(nsclr) local 29 tmp(j)=y(i,5+j) 31 call Dirichlet2Neumann(nsclr, itype, tmp) 39 BC(i,6+j)=-tmp(j) 46 subroutine Dirichlet2Neumann_faux(nscalar, itype, tmp) argument 64 real*8 tmp(nscalar) local 79 tmp(i) = 0.0d0 91 tmp(i) = D(i) * ( tmp(i) - 0.0 ) / distance * units 92 if(itype.eq.2) tmp(i) = tmp(i) / 2.0d+00 141 subroutine Dirichlet2Neumann(nscalar, itype, tmp) argument [all …]
|
| H A D | rotabc.f | 21 & tmp(nshg) 27 tmp = global(:,1)*acs(:,1) - global(:,2)*acs(:,2) 29 global(:,1) = tmp 33 tmp = global(:,1)*acs(:,1) + global(:,2)*acs(:,2) 35 global(:,1) = tmp
|
| H A D | getstrl.f | 13 dimension tmp(npro), fresli(npro,24) 56 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 59 dxidx(:,1,1) = dxidx(:,1,1) * tmp 60 dxidx(:,1,2) = dxidx(:,1,2) * tmp 61 dxidx(:,1,3) = dxidx(:,1,3) * tmp 63 & - dxdxi(:,2,1) * dxdxi(:,3,3)) * tmp 65 & - dxdxi(:,3,1) * dxdxi(:,1,3)) * tmp 67 & - dxdxi(:,1,1) * dxdxi(:,2,3)) * tmp 69 & - dxdxi(:,2,2) * dxdxi(:,3,1)) * tmp 71 & - dxdxi(:,1,1) * dxdxi(:,3,2)) * tmp [all …]
|
| H A D | get_time.c | 17 uint64_t tmp = ((*last)-(*first)); in get_max_time_diff() local 23 assert(tmp > 0); in get_max_time_diff() 24 MPI_Allreduce(&tmp, &iresult, 1, MPI_UINT64_T, MPI_MAX, MPI_COMM_WORLD); in get_max_time_diff() 28 tmp = *c_last + (UINT64_MAX - *c_first); in get_max_time_diff() 32 tmp = ((*c_last)-(*c_first)); in get_max_time_diff() 34 MPI_Allreduce(&tmp, &iresult, 1, MPI_UINT64_T, MPI_MAX, MPI_COMM_WORLD); in get_max_time_diff()
|
| H A D | sfID2np.f | 66 integer, allocatable :: nodeMap(:), tmp(:) local 68 allocate(tmp(nshg)) 69 call sfID2np(isurfID, nSurfNP, tmp) 72 nodeMap = tmp(1:nSurfNP) 81 deallocate(tmp)
|
| H A D | asithf.f | 16 dimension tmp(npro) 62 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 65 dxidx(:,1,1) = dxidx(:,1,1) * tmp 66 dxidx(:,1,2) = dxidx(:,1,2) * tmp 67 dxidx(:,1,3) = dxidx(:,1,3) * tmp 69 & - dxdxi(:,2,1) * dxdxi(:,3,3)) * tmp 71 & - dxdxi(:,3,1) * dxdxi(:,1,3)) * tmp 73 & - dxdxi(:,1,1) * dxdxi(:,2,3)) * tmp 75 & - dxdxi(:,2,2) * dxdxi(:,3,1)) * tmp 77 & - dxdxi(:,1,1) * dxdxi(:,3,2)) * tmp [all …]
|
| H A D | qpbc.f | 17 & tmpvec(nshg,4), tmp(nshg) 76 tmp = sqrt(qres(:,idflow+1)**2 80 if (tmp(i) .lt. 0.0001) tmp(i) = 0.0001 82 tmp = one/tmp 85 qres(:,idflow+i) = tmp*qres(:,idflow+i)
|
| /phasta/phSolver/compressible/ |
| H A D | shuffle.f | 11 dimension restmp (nshg,nflow), tmp(nshg), tmp1(nshg,nflow) 15 tmp(:)=restmp(:,1) ! copying the res of continuity 17 restmp(:, 4) = tmp(:) 42 dimension bdiagtmp (nshg,nflow,nflow), tmp(nshg,nflow) 44 tmp(:,:) = bdiagtmp(:,:,1) ! reshuffling 1st column with 46 bdiagtmp(:,:,4) = tmp(:,:) 49 tmp(:,:) = bdiagtmp(:,1,:) ! reshuffling 1st row with 51 bdiagtmp(:,4,:) = tmp (:,:)
|
| H A D | e3ql.f | 46 real*8 tmp(npro) local 144 tmp = Qwt(lcsyst,intp) 146 tmp = tmp*(three/four) 150 qrl(:,i,1 ) = qrl(:,i,1 )+ shape(:,i)*tmp*qdi(:,1 ) 151 qrl(:,i,2 ) = qrl(:,i,2 )+ shape(:,i)*tmp*qdi(:,2 ) 152 qrl(:,i,3 ) = qrl(:,i,3 )+ shape(:,i)*tmp*qdi(:,3 ) 153 qrl(:,i,4 ) = qrl(:,i,4 )+ shape(:,i)*tmp*qdi(:,4 ) 154 qrl(:,i,5 ) = qrl(:,i,5 )+ shape(:,i)*tmp*qdi(:,5 ) 155 qrl(:,i,6 ) = qrl(:,i,6 )+ shape(:,i)*tmp*qdi(:,6 ) 156 qrl(:,i,7 ) = qrl(:,i,7 )+ shape(:,i)*tmp*qdi(:,7 ) [all …]
|
| H A D | itrbc.f | 30 real*8 tmp(nshg), y1(nshg),q1(nshg) local 148 call getthm (y1, y(:,5), tmp, 149 & rk, q1, tmp, 150 & tmp, tmp, tmp, 151 & tmp, tmp, tmp, 152 & tmp, tmp)
|
| H A D | e3qvar.f | 41 & shp(npro,nshl), tmp(npro), 96 & tmp, rho, tmp, 97 & tmp, tmp, tmp, 98 & cp, tmp, tmp, 99 & tmp, tmp)
|
| H A D | e3bdg.f | 28 dimension tmp(npro) 38 tmp=shp(:,j)*WdetJ 40 & + tmp * ( 46 & + tmp * ( 52 & + tmp * ( 58 & + tmp * ( 64 & + tmp * ( 70 & + tmp * ( 76 & + tmp * ( 82 & + tmp * ( [all …]
|
| H A D | e3juel.f | 30 & ei(npro), tmp(npro), 60 & tmp, tmp, tmp, 61 & tmp, tmp, tmp, 62 & tmp, tmp)
|
| H A D | e3bdg_nd.f | 29 dimension tmp(npro) 39 tmp=shp(:,j)*WdetJ 41 & + tmp * ( 47 & + tmp * ( 53 & + tmp * ( 59 & + tmp * ( 65 & + tmp * ( 71 & + tmp * ( 77 & + tmp * ( 83 & + tmp * ( [all …]
|
| H A D | solgmr.f | 271 tmp = Rcos(jK) * HBrg(jK, iKs) + 275 HBrg(jK, iKs) = tmp 278 tmp = sqrt(HBrg(iKs,iKs)**2 + HBrg(iKs+1,iKs)**2) 279 Rcos(iKs) = HBrg(iKs, iKs) / tmp 280 Rsin(iKs) = HBrg(iKs+1,iKs) / tmp 281 HBrg(iKs, iKs) = tmp 286 tmp = Rcos(iKs) * eBrg(iKs) + Rsin(iKs) * eBrg(iKs+1) 288 eBrg(iKs) = tmp 646 tmp = Rcos(jK) * HBrg(jK, iKs) + 650 HBrg(jK, iKs) = tmp [all …]
|
| H A D | e3ivar.f | 10 & ql, divqi, sgn, tmp, argument 84 dimension tmp(npro), dxdxi(npro,nsd,nsd), sgn(npro,nshl) 122 & tmp, tmp, tmp, 123 & tmp, tmp, tmp, 124 & tmp, tmp) 221 & h, tmp, cv, 223 & tmp, tmp) 578 dimension tmp(npro), dxdxi(npro,nsd,nsd) 643 & rk, rho, tmp, 644 & tmp, tmp, tmp, [all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | e3ql.f | 44 real*8 tmp(npro) local 66 call getdiff(dwl, yl, shape, xmudmi, xl,rmu, tmp) 90 tmp = Qwt(lcsyst,intp) 92 tmp = tmp*(three/four) 99 qrl(:,i,1 ) = qrl(:,i,1 )+ shape(:,i)*tmp*qdi(:,1 ) 100 qrl(:,i,2 ) = qrl(:,i,2 )+ shape(:,i)*tmp*qdi(:,2 ) 101 qrl(:,i,3 ) = qrl(:,i,3 )+ shape(:,i)*tmp*qdi(:,3 ) 103 qrl(:,i,4 ) = qrl(:,i,4 )+ shape(:,i)*tmp*qdi(:,4 ) 104 qrl(:,i,5 ) = qrl(:,i,5 )+ shape(:,i)*tmp*qdi(:,5 ) 105 qrl(:,i,6 ) = qrl(:,i,6 )+ shape(:,i)*tmp*qdi(:,6 ) [all …]
|
| H A D | e3source.f | 255 & dP2Inv, sixth, tmp, tmp1, p, dp, local 390 tmp = 1.0 / ( chiP3 + saCv1P3 ) 391 fv1 = chiP3 * tmp 392 fv1p = 3.0 * viscInv * saCv1P3 * chiP2 * tmp**2 394 tmp = 1.0 / (1.0 + chi * fv1) 395 fv2 = 1.0 - chi * tmp 396 fv2p = (chiP2 * fv1p - viscInv) * tmp**2 399 tmp = SclrNN * k2d2Inv !eOkdP2 400 st = s + fv2 * tmp 401 stp = ep*k2d2Inv * fv2 + tmp*fv2p [all …]
|
| H A D | e3qvar.f | 45 dimension tmp(npro), dxdxi(npro,nsd,nsd) 72 tmp = one / ( dxidx(:,1,1) * dxdxi(:,1,1) 75 dxidx(:,1,1) = dxidx(:,1,1) * tmp 76 dxidx(:,1,2) = dxidx(:,1,2) * tmp 77 dxidx(:,1,3) = dxidx(:,1,3) * tmp 79 & - dxdxi(:,2,1) * dxdxi(:,3,3)) * tmp 81 & - dxdxi(:,3,1) * dxdxi(:,1,3)) * tmp 83 & - dxdxi(:,1,1) * dxdxi(:,2,3)) * tmp 85 & - dxdxi(:,2,2) * dxdxi(:,3,1)) * tmp 87 & - dxdxi(:,1,1) * dxdxi(:,3,2)) * tmp [all …]
|
| H A D | e3lhs.f | 47 & tmp(npro), tlW(npro) 137 tmp = t3(:,1) * shg(:,aa,1) 142 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp 145 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp 148 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp 170 tmp = t3(:,1) * shg(:,aa,1) 174 tmp1 = tmp 180 tmp1 = tmp 186 tmp1 = tmp 293 & tmp(npro), dcFct(npro), giju(npro,6) [all …]
|
| /phasta/svLS/ |
| H A D | NSSOLVER.f | 197 REAL*8, ALLOCATABLE :: tmp(:) local 200 2 Gt(nsd,nnz), tmp((nsd+1)*(nsd+1))) 204 tmp = Val(:,i) 206 mK(1,i) = tmp(1) 207 mK(2,i) = tmp(2) 208 mK(3,i) = tmp(4) 209 mK(4,i) = tmp(5) 211 mG(1,i) = tmp(3) 212 mG(2,i) = tmp(6) 214 mD(1,i) = tmp(7) [all …]
|
| H A D | NORM.f | 50 REAL*8 tmp, NORMV local 78 CALL MPI_ALLREDUCE(NORMV, tmp, 1, mpreal, MPI_SUM, 80 NORMV = tmp 99 REAL*8 tmp, NORMS local 107 CALL MPI_ALLREDUCE(NORMS, tmp, 1, mpreal, MPI_SUM, 109 NORMS = tmp
|
| H A D | DOT.f | 50 REAL*8 tmp, DOTV local 79 CALL MPI_ALLREDUCE(DOTV, tmp, 1, mpreal, MPI_SUM, commu%comm, 82 DOTV = tmp 98 REAL*8 tmp, DOTS local 107 CALL MPI_ALLREDUCE(DOTS, tmp, 1, mpreal, MPI_SUM, commu%comm, 110 DOTS = tmp
|