Home
last modified time | relevance | path

Searched refs:e (Results 1 – 25 of 48) sorted by relevance

12

/phasta/phSolver/incompressible/
H A De3sts.f83 integer e local
90 do e = 1, npro
92 x12 = xl(e,2,1) - xl(e,1,1)
93 x13 = xl(e,3,1) - xl(e,1,1)
94 x14 = xl(e,4,1) - xl(e,1,1)
95 x23 = xl(e,3,1) - xl(e,2,1)
96 x24 = xl(e,4,1) - xl(e,2,1)
97 x34 = xl(e,4,1) - xl(e,3,1)
99 y12 = xl(e,2,2) - xl(e,1,2)
100 y13 = xl(e,3,2) - xl(e,1,2)
[all …]
H A Dtimeseries.f18 integer e, founde local
144 do e = 1, npro
145 if ((abs(zi0(e,1)).lt.(one+tolpt)).and.
146 & (abs(zi0(e,2)).lt.(one+tolpt)).and.
147 & (abs(zi0(e,3)).lt.(one+tolpt))) then ! got the element
149 call shphex (ipord, zi0(e,:),shape(:),
152 founde=e
164 do e = 1, npro
165 detaij(e) = al(e,2,1)*(-al(e,3,2)*al(e,4,3) +
166 & al(e,4,2)*al(e,3,3)) + al(e,2,2)*
[all …]
H A Dgetdiff.f13 integer n, e
137 integer n, e local
205 integer e, n
209 do e = 1, npro
215 u1=yl(e,n,2)
216 u2=yl(e,n,3)
217 u3=yl(e,n,4)
235 evisc = evisc + yl(e,n,6)
240 rmu(e) = rmu(e) + abs(evisc)
265 Turb = Turb + shape(e,n) * yl(e,n,6)
[all …]
H A De3source.f267 integer e,n local
323 do e = 1, npro
324 SclrNN= max(Sclr(e),zero)
348 dx=maxval(xl(e,:,1))-minval(xl(e,:,1))
349 dy=maxval(xl(e,:,2))-minval(xl(e,:,2))
350 dz=maxval(xl(e,:,3))-minval(xl(e,:,3))
355 dwall(e) = min(dmax,dwall(e))
360 & gradV(e,1,1)**2+gradV(e,1,2)**2+gradV(e,1,3)**2
361 & + gradV(e,2,1)**2+gradV(e,2,2)**2+gradV(e,2,3)**2
362 & + gradV(e,3,1)**2+gradV(e,3,2)**2+gradV(e,3,3)**2 )
[all …]
/phasta/phSolver/compressible/
H A Dtimeseries.f18 integer e, founde local
144 do e = 1, npro
145 if ((abs(zi0(e,1)).lt.(one+tolpt)).and.
146 & (abs(zi0(e,2)).lt.(one+tolpt)).and.
147 & (abs(zi0(e,3)).lt.(one+tolpt))) then ! got the element
149 call shphex (ipord, zi0(e,:),shape(:),
152 founde=e
164 do e = 1, npro
165 detaij(e) = al(e,2,1)*(-al(e,3,2)*al(e,4,3) +
166 & al(e,4,2)*al(e,3,3)) + al(e,2,2)*
[all …]
H A Dgetdiff.f51 integer n, e local
185 do e = 1, npro
196 u1=ycl(e,n,2)
197 u2=ycl(e,n,3)
198 u3=ycl(e,n,4)
215 evisc = evisc + ycl(e,n,6)
220 xmut(e)= abs(evisc)
245 savar = savar + shp(e,n) * ycl(e,n,6)
247 xki = rho(e)*abs(savar)/rmu(e)
250 xmut(e) = fv1*abs(savar)*rho(e)
/phasta/phSolver/common/
H A Delm3keps.f128 integer e local
172 do e = 1, npro
174 nuInv = rho(e)/rmu(e)
175 k = abs(kay(e))
176 if (k.lt.1.e-32) k=0
177 eps = abs(epsilon(e))
178 y = dwi(e)
193 if ( k .gt.1.e-32 ) then
208 if ( eps .gt.1.e-32) then
214 ss = gradV(e,1,1) ** 2
[all …]
H A Dfillsparse.f12 integer aa, b, c, e, i, k, n local
26 do e = 1, npro ! loop over the elements
28 i = ien(e,aa) ! finds the global equation number or
37 k = sparseloc( row(c), n, ien(e,b) ) + c-1
44 lhsK(1,k) = lhsK(1,k) + xKebe(e,1,aa,b)
45 lhsK(2,k) = lhsK(2,k) + xKebe(e,2,aa,b)
46 lhsK(3,k) = lhsK(3,k) + xKebe(e,3,aa,b)
47 lhsK(4,k) = lhsK(4,k) + xKebe(e,4,aa,b)
48 lhsK(5,k) = lhsK(5,k) + xKebe(e,5,aa,b)
49 lhsK(6,k) = lhsK(6,k) + xKebe(e,6,aa,b)
[all …]
H A Delem-search.f32 integer e, elmt, result, count, counting, find local
75 do e = 1, nelint
79 if (zi0(e,1).lt.(one+tolpt).and.
80 & zi0(e,1).gt.(zero-tolpt)) then
84 if (zi0(e,2).lt.(one+tolpt).and.
85 & zi0(e,2).gt.(zero-tolpt)) then
89 if (zi0(e,3).lt.(one+tolpt).and.
90 & zi0(e,3).gt.(zero-tolpt)) then
94 if ((1-zi0(e,1)-zi0(e,2)-zi0(e,3)).lt.(one+tolpt).and.
95 & (1-zi0(e,1)-zi0(e,2)-zi0(e,3)).gt.(zero-tolpt)) then
[all …]
H A Dfillsparse.c26 int e,i,j,aa; //following along with fillsparse.f in fillsparsecpetscs() local
28 for(e=0;e<npro;e++) in fillsparsecpetscs()
30 for(aa=0;aa<nshl;aa++) locat[aa]=ieng[e+npro*aa]-1; in fillsparsecpetscs()
34 mb[nshl*i + j] = EGmass[e + npro*(i + nshl*j)]; in fillsparsecpetscs()
51 int e,i,j,aa; //following along with fillsparse.f in fillsparsecpetscc() local
54 for(e=0;e<npro;e++) in fillsparsecpetscc()
56 for(aa=0;aa<nshl;aa++) locat[aa]=ieng[e+npro*aa]-1; in fillsparsecpetscc()
60 mb[nedof*i + j] = EGmass[e + npro*(i + nedof*j)]; in fillsparsecpetscc()
H A DBCprofile2.f63 if(ycoor<(-0.587375+1.0e-4))then
64 ry=(ycoor+0.587375)/1.0e-4
65 elseif(ycoor>(0.587375-1.0e-4))then
66 ry=(0.587375-ycoor)/1.0e-4
71 if(zcoor<(-0.587375+1.0e-4))then
72 rz=(zcoor+0.587375)/1.0e-4
73 elseif(zcoor>(0.587375-1.0e-4))then
74 rz=(0.587375-zcoor)/1.0e-4
H A DsetBlowing_Duct3.f52 if(zcoor<(-0.05715+1.0e-3))then
53 rz=(zcoor+0.05715)/1.0e-3
54 elseif(zcoor>(0.05715-1.0e-3))then
55 rz=(0.05715-zcoor)/1.0e-3
86 BC(nn,7) = 1.825e-5 ! set varSA = 0
/phasta/phSolver/Test_IC/
H A Dgeom.sms20 0 6.12323399573677e-17 0.5 1.039777937572983e-312
22 0.5 -1 0.5 1.039777937572983e-312
24 0 6.12323399573677e-17 0 1.039777937572983e-312
26 0.5 -1 0 1.039777937572983e-312
28 0.5 1 0.5 1.039777937572983e-312
30 0.5 1 0 1.039777937572983e-312
32 1 6.12323399573677e-17 0.5 1.039777937572983e-312
34 1 6.12323399573677e-17 0 1.039777937572983e-312
36 0.5 6.12323399573677e-17 0 1.039777937572983e-312 0.5 0
38 0.5 6.12323399573677e-17 0.5 1.039777937572983e-312 0.5 0
H A Dgeom.bdr20 0 6.12323399573677e-17 0.5 0
24 0 6.12323399573677e-17 0 0
32 1 6.12323399573677e-17 0.5 0
34 1 6.12323399573677e-17 0 0
36 0.5 6.12323399573677e-17 0 0 0 0
38 0.5 6.12323399573677e-17 0.5 0 0 0
/phasta/phSolver/Test_C/
H A Dgeom.sms20 0 6.12323399573677e-17 0.5 1.039777937572983e-312
22 0.5 -1 0.5 1.039777937572983e-312
24 0 6.12323399573677e-17 0 1.039777937572983e-312
26 0.5 -1 0 1.039777937572983e-312
28 0.5 1 0.5 1.039777937572983e-312
30 0.5 1 0 1.039777937572983e-312
32 1 6.12323399573677e-17 0.5 1.039777937572983e-312
34 1 6.12323399573677e-17 0 1.039777937572983e-312
36 0.5 6.12323399573677e-17 0.5 1.039777937572983e-312 0.5 0
38 0.5 6.12323399573677e-17 0 1.039777937572983e-312 0.5 0
H A Dgeom.bdr20 0 6.12323399573677e-17 0.5 0
24 0 6.12323399573677e-17 0 0
32 1 6.12323399573677e-17 0.5 0
34 1 6.12323399573677e-17 0 0
36 0.5 6.12323399573677e-17 0.5 0 0 0
38 0.5 6.12323399573677e-17 0 0 0 0
H A Dsolver.inp36 Viscosity: 1.0e-2 # fills datmat (2 values REQUIRED if iLset=1)
54 Tolerance on Momentum Equations: 1.0e-4 # epstol(1)
55 Tolerance on ACUSIM Pressure Projection: 1.0e-2 # prestol
/phasta/phSolver/testIC-DSIM/
H A Dgeom.sdm104 0 6.12323399573677e-17 0 0
110 1 6.12323399573677e-17 0 0
112 0 6.12323399573677e-17 0.5 0
118 1 6.12323399573677e-17 0.5 0
120 0.5 6.12323399573677e-17 0 0 0 0
122 0.5 6.12323399573677e-17 0.5 0 0 0
139 a 12 5 e
140 a 12 e 7
149 1 22 e h
153 4 22 a e
[all …]
H A Dgeom.sms104 0 6.12323399573677e-17 0 0
110 1 6.12323399573677e-17 0 0
112 0 6.12323399573677e-17 0.5 0
118 1 6.12323399573677e-17 0.5 0
120 0.5 6.12323399573677e-17 0 0 0 0
122 0.5 6.12323399573677e-17 0.5 0 0 0
139 a 12 5 e
140 a 12 e 7
149 1 23 e h
153 4 23 a e
[all …]
/phasta/svLS/
H A DINCOMMU.f50 INTEGER i, j, k, s, e, ierr, nTasks, tF, stat(MPI_STATUS_SIZE) local
69 e = s + cS(i)%n - 1
71 CALL MPI_IRECV(rTmp(:,s:e), cS(i)%n*dof, mpreal, i-1,
74 CALL MPI_ISEND(R(:,s:e), cS(i)%n*dof, mpreal, i-1,
86 e = s + cS(i)%blN(j) - 1
87 R(:,s:e) = R(:,s:e) + rTmp(:,k:k+e-s)
97 e = s + cS(i)%blN(j) - 1
98 rTmp(:,k:k+e-s) = R(:,s:e)
106 e = s + cS(i)%n - 1
109 CALL MPI_IRECV(R(:,s:e), cS(i)%n*dof, mpreal, i-1,
[all …]
H A DLHS.f52 INTEGER i, j, k, a, Ac, ai, s, e, nTasks, tF, maxnNo, ierr, local
84 e = rowPtr(Ac+1) - 1
85 DO i=s, e
94 lhs%rowPtr(2,Ac) = e
219 e = s + lhs%cS(i)%n - 1
221 CALL MPI_RECV(part(s:e), lhs%cS(i)%n, mpint, i-1,
225 DO j=s, e
239 DO j=s, e
255 CALL MPI_SEND(ltg(s:e), lhs%cS(i)%n, mpint, i-1,
H A DSPARMUL.f53 INTEGER i, j, l, col, s, e local
102 e = l*dof
103 s = e - dof + 1
104 KU(l,i) = KU(l,i) + SUM(K(s:e,j)*U(:,col))
/phasta/phSolver/incompressible/e3source/
H A De3source.fimplicit263 integer e
319 do e = 1, npro
320 SclrNN= max(Sclr(e),zero)
325 dx=maxval(xl(e,:,1))-minval(xl(e,:,1))
326 dy=maxval(xl(e,:,2))-minval(xl(e,:,2))
327 dz=maxval(xl(e,:,3))-minval(xl(e,:,3))
333 qfac = sqrt(gradV(e,1,1)**2+gradV(e,1,2)**2+gradV(e,1,3)**2
334 & +gradV(e,2,1)**2+gradV(e,2,2)**2+gradV(e,2,3)**2
335 & +gradV(e,3,1)**2+gradV(e,3,2)**2+gradV(e,3,3)**2)
336 qfac = max(qfac,1.0e-10)
[all …]
H A De3source.fRHSonly263 integer e
319 do e = 1, npro
320 SclrNN= max(Sclr(e),zero)
323 dx=maxval(xl(e,:,1))-minval(xl(e,:,1))
324 dy=maxval(xl(e,:,2))-minval(xl(e,:,2))
325 dz=maxval(xl(e,:,3))-minval(xl(e,:,3))
331 qfac = sqrt(gradV(e,1,1)**2+gradV(e,1,2)**2+gradV(e,1,3)**2
332 & +gradV(e,2,1)**2+gradV(e,2,2)**2+gradV(e,2,3)**2
333 & +gradV(e,3,1)**2+gradV(e,3,2)**2+gradV(e,3,3)**2)
334 rd=SclrNN*saKappaP2Inv/(dwall(e)**2*qfact)
[all …]
/phasta/phSolver/compressible/test/
H A Drunphasta.cmake26 ${COMPARE_EXE} ${CASEDIR} ${CASEDIR}_ref ${IS_SYNCIO} 1e-6)

12