Home
last modified time | relevance | path

Searched refs:y (Results 1 – 25 of 73) sorted by relevance

123

/phasta/phSolver/incompressible/
H A Ditrbc.f1 subroutine itrBC (y,ac, iBC, BC, iper, ilwork) argument
23 dimension y(nshg,nflow), iBC(nshg),
32 & y(:,i)=min(ylimit(3,i),max(ylimit(2,i),y(:,i)))
41 y(:,1) = BC(:,3) - BC(:,4) * y(:,2)
42 & - BC(:,5) * y(:,3)
48 y(:,2) = BC(:,3) - BC(:,4) * y(:,1)
49 & - BC(:,5) * y(:,3)
55 y(:,1) = BC(:,3) - BC(:,4) * y(:,3)
56 y(:,2) = BC(:,5) - BC(:,6) * y(:,3)
62 y(:,3) = BC(:,3) - BC(:,4) * y(:,1)
[all …]
H A DitrPC.f7 subroutine itrSetup ( y, acold ) argument
11 real*8 y(nshg,ndof), acold(nshg,ndof) local
44 y=y-(one-alfi)*Delt(1)*acold
65 subroutine itrPredict (yold, y, acold, ac, uold, u, iBC) argument
69 real*8 y(nshg,ndof), ac(nshg,ndof), local
77 y = yold
91 y = yold + (one - gami)/Dtgl * acold
96 y = yold+alfi*Delt(itseq)*acold
104 y = yold+fct1*(yold-y)
105 ac = acold*fct2+(y-yold)*fct3
[all …]
H A Dftools.f267 subroutine flesDaxpy ( x, y, a, m, n ) argument
271 real*8 x(n,m), y(n,m) local
276 y(i,j) = y(i,j) + a * x(i,j)
287 subroutine flesDxpay ( x, y, a, m, n ) argument
291 real*8 x(n,m), y(n,m) local
296 y(i,j) = a * y(i,j) + x(i,j)
328 subroutine fMtxBlkDot2( x, y, c, m, n ) argument
334 real*8 x(n,m), y(n), c(m) local
352 tmp1 = tmp1 + x(i,1) * y(i)
361 tmp1 = tmp1 + x(i,1) * y(i)
[all …]
H A Delmgmr.f1 subroutine ElmGMR (u, y, ac, x, argument
28 dimension y(nshg,ndof), ac(nshg,ndof),
102 call AsIqGradV (y, x,
108 call AsIq (y, x,
169 call AsIGMR (y, ac,
290 call AsBMFG (u, y,
330 call ElmpvsQ (res,y,-1.0d0)
368 subroutine ElmGMRSclr (y, ac, x, argument
388 dimension y(nshg,ndof), ac(nshg,ndof),
437 call AsIqSclr (y, x,
[all …]
H A Dasiq.f1 subroutine AsIq (y, x, shp, argument
27 dimension y(nshg,ndof), x(numnp,nsd),
53 call localy(y, yl, ien, ndof, 'gather ')
88 subroutine AsIqSclr (y, x, shp, argument
95 dimension y(nshg,ndof), x(numnp,nsd),
111 call localy(y, yl, ien, ndof, 'gather ')
H A Ditrdrv.f1 subroutine itrdrv (y, ac, argument
55 real*8 y(nshg,ndof), ac(nshg,ndof), local
154 yold = y
244 call itrSetup ( y, acold )
256 call calcRCRic(y,nsrflistRCR,numRCRSrfs)
362 call itrPredict(yold, y, acold, ac , uold, u, iBC)
363 call itrBC (y, ac, iBC, BC, iper,ilwork)
367 call itrBCSclr (y, ac, iBC, BC, iper, ilwork)
370 call itrBCSclr (y, ac, iBC, BC, iper, ilwork)
389 call SolFlow(y, ac, u,
[all …]
H A Dsoldir.f1 subroutine SolDir (y, ac, yold, acold, argument
18 dimension y(nshg,ndof), ac(nshg,ndof),
45 call itrYAlpha( yold, acold, y, ac, yAlpha, acAlpha)
123 call itrCorrect ( y, ac, solinc, iBC)
124 call itrBC (y, ac, iBC, BC, iper, ilwork)
128 call rstatic (res, y, solinc)
H A Dasbmfg.f1 subroutine AsBMFG (u, y, ac, x, argument
19 dimension y(nshg,ndofl), x(numnp,nsd),
42 call localy(y, yl, ienb, ndofl, 'gather ')
87 subroutine AsBSclr (y, x, shpb, shglb, argument
93 dimension y(nshg,ndofl), x(numnp,nsd),
112 call localy(y, yl, ienb, ndofl, 'gather ')
/phasta/phSolver/compressible/
H A Ditrbc.f1 subroutine itrBC (y,ac, iBC, BC, iper, ilwork) argument
25 dimension y(nshg,nflow), iBC(nshg),
40 locmax=maxloc(y(:,i))
41 locmin=minloc(y(:,i))
42 ymax=maxval(y(:,i))
43 ymin=minval(y(:,i))
47 ypc=y(j,i)
48 y(j,i)=min(ylimit(3,i),max(ylimit(2,i),y(j,i)))
49 if(ypc.ne.y(j,i) )limitcount(i)=limitcount(i)+one
63 y(:,5) = BC(:,2)
[all …]
H A DitrPC.f6 subroutine itrSetup ( y, acold ) argument
10 real*8 y(nshg,ndof), acold(nshg,ndof) local
32 & y=y-(one-alfi)*Delt(1)*CFLfl(1)*acold
57 & y, ac ) argument
62 & y(nshg,ndof), ac(nshg,ndof) local
71 y = yold
79 y = yold+alfi/Dtgl*acold*(one-gami)
80 call itrBC (y, ac, iBC, BC, iper,ilwork)
83 call genscale(y, x, iBC)
92 y = yold+alfi/Dtgl*acold
[all …]
H A Ditrdrv.f1 subroutine itrdrv (y, ac, uold, x, argument
45 dimension y(nshg,ndof), ac(nshg,ndof),
148 yold = y
197 ybar(:,1:ndof) = y(:,1:ndof)
199 ybar(:,ndof+idof) = y(:,idof)*y(:,idof)
201 ybar(:,ndof+6) = y(:,1)*y(:,2)
202 ybar(:,ndof+7) = y(:,1)*y(:,3)
203 ybar(:,ndof+8) = y(:,2)*y(:,3)
252 call itrSetup ( y, acold)
342 yold(:,6:ndof)=y(:,6:ndof)*tdecay
[all …]
H A Dbc3res.f1 subroutine bc3Res (y, iBC, BC, res, iper, ilwork) argument
23 dimension y(nshg,ndof), iBC(nshg),
46 res(:,2) = res(:,2) -y(:,1)*res(:,1)
47 res(:,3) = res(:,3) -y(:,2)*res(:,1)
48 res(:,4) = res(:,4) -y(:,3)*res(:,1)
50 & (gamma*Rgas/gamma1*y(:,5)
51 & + pt5 * ( y(:,1)**2 + y(:,2)**2 + y(:,3)**2))*res(:,1)
207 subroutine bc3ResSclr (y, iBC, BC, rest, argument
230 dimension y(nshg,ndof), iBC(nshg),
H A Dvorticity.f1 subroutine vortGLB (y, x, shp, shgl, ilwork, vortG) argument
7 real*8 y(nshg,ndof), x(numnp,nsd) local
48 call vortELE(y,x,tmpshp,tmpshgl,mien(iblk)%p,vortIntG,lpmassG)
78 subroutine vortELE(y,x,shp1,shgl1,ien,vortIntG,lpmassG) argument
83 real*8 y(nshg,ndof), x(numnp,nsd) local
108 call localy (y, yl, ien, ndof, 'gather ')
192 subroutine PRcalcDuct (y) argument
197 real*8 y(nshg,ndof) local
207 Mach(:)=sqrt(y(:,1)**2+y(:,2)**2+y(:,3)**2)/sqrt(1.4*287*y(:,5))
208 PresRec(:)=y(:,4)*(1+0.2*Mach(:)**2)**3.5/Ptotal_contra
[all …]
H A Ditrfdi.f1 subroutine itrFDI (ypre, y, ac, x, argument
39 dimension y(nshg,ndof), ypre(nshg,nflow),
69 call itrRes (ytmp, y,
103 call itrRes (ytmp, y,
120 call itrRes (ytmp, y,
151 subroutine itrFDISclr (y, ypre, x, argument
189 dimension y(nshg,ndof), ypre(nshg,ndof),
208 call itrRes (ytmp, y,
234 call itrRes (ytmp, y,
246 call itrRes (ytmp, y,
H A DrstatCheck.f1 subroutine rstatCheck (res, ilwork,y,ac) argument
27 real*8 y(nshg,ndof),ac(nshg,ndof) local
74 call restar('out ',y,res) ! 'res' is used instead of 'ac'
89 subroutine rstatCheckSclr (rest, ilwork,y,ac) argument
114 real*8 y(nshg,ndof),ac(nshg,ndof) local
166 call restar('out ',y,ac)
H A Dasbmfg.f1 subroutine AsBMFG (y, x, shpb, shglb, argument
15 dimension y(nshg,ndofl), x(numnp,nsd),
39 call localy(y, ycl, ienb, ndofl, 'gather ')
68 subroutine AsBMFGSclr (y, x, shpb, shglb, argument
82 dimension y(nshg,ndofl), x(numnp,nsd),
102 call localy (y, ycl, ienb, ndofl, 'gather ')
/phasta/phSolver/common/
H A DINIprofile.f1 subroutine INIprofile(BC,y,x) argument
5 real*8 y(nshg,ndof) local
9 call TakeBCfromIC_ScaleInlVelX(BC,y,x)
13 call setInitial(x,y)
19 subroutine TakeBCfromIC_ScaleInlVelX(BC,y,x) argument
28 real*8 y(nshg,ndof) local
38 if(y(nn,1).gt.u_max)u_max=y(nn,1)
50 BC(inlf(:),2) = y(inlf(:),5) ! set temp
51 BC(inlf(:),3) = y(inlf(:),1)/u_max_all*inletVelX ! set and scale x velocity
55 BC(inlf(:),7) = y(inlf(:),6) ! set scalar_1
[all …]
H A Dproces.f24 dimension y(nshg,ndof),
58 call gendat (y, ac, point2x,
94 call restar ('out ', y ,ac)
101 call initSponge( y,point2x)
136 y(i,1:3)=(xi*bcinterp(iupper,4:6)
138 y(i,4)=(xi*bcinterp(iupper,2)
140 y(i,5)=(xi*bcinterp(iupper,3)
183 … if(y(nn,6) .eq. 0) cycle !don't change wall boundary conditions, should be iTurbWall == 1
186 y(nn,6) = EVrampMax
191 y(nn,6) = c0 + (point2x(nn,1) - x1)*(
[all …]
H A Dgenini.f1 subroutine genini (iBC, BC, y, ac, iper, ilwork, argument
33 & BC(nshg,ndofBC), y(nshg,ndof),
56 call restar ('in ', y, ac)
60 call INIprofile(BC,y,x)
72 call getvel (y, ilwork, iBC,
142 call itrBC (y, ac, iBC, BC, iper, ilwork)
147 call itrBCSclr (y, ac, iBC, BC, iper, ilwork)
150 call itrBCSclr (y, ac, iBC, BC, iper, ilwork)
154 call genscale(y, x, iBC)
164 write (iecho,1100) n, (y(n,i),i=1,ndof)
H A Dgenscale.f1 subroutine genscale(y, x, iBC) argument
24 dimension y(numnp,ndof), iBC(numnp),
42 ymapped(:,2:4)=y(:,1:3)
43 ymapped(:,1)=y(:,4)
44 ymapped(:,5)=y(:,5)
57 yintl(i,:,j+1) = y(ien2D(i,:),j)
59 yintl(i,:,1) = y(ien2D(i,:),4)
62 yintl(i,:,j) = y(ien2D(i,:),j)
383 y(:,1) = ymapped(:,2)
384 y(:,2) = ymapped(:,3)
[all …]
H A Dgetvel.f1 subroutine getvel (y,ilwork, iBC, argument
14 & y(numnp,ndof),
21 velo(:,1)=y(:,4)
22 velo(:,2)=y(:,1)
23 velo(:,3)=y(:,2)
24 velo(:,4)=y(:,3)
25 if(nflow.eq.5) velo(:,5)=y(:,5)
H A Dsettauw.f1 subroutine settauw (y, x, BC, argument
16 dimension y(nshg,ndofl), x(numnp,nsd),
35 ull(:)=y(otwn(nodw),1:3)
63 ull(:)=y(otwn(nodw),1:3)
104 & y(nodw,1)**2 !
105 & +y(nodw,2)**2 ! flow is auto-|| at boundary
106 & +y(nodw,3)**2 !
179 function utaul(u,y,rm)
185 function utau(u,y,rm,nodw,x)
187 real*8 u,err,utau,yrmi,y,rm,yp,up,kup,f,dfds,rat,efac local
[all …]
H A Dasbwmod.f1 subroutine asbwmod (y, ac, x, BC, iBC, argument
13 dimension y(nshg,ndof), x(numnp, nsd),
22 call settauw (y, x,
30 call itrBCSclr (y, ac, iBC, BC, iper, ilwork)
H A Ddtn.f20 subroutine DtN(iBC,BC,y) argument
23 real*8 BC(nshg,ndofBC),y(nshg,ndof),tmp(nsclr) local
29 tmp(j)=y(i,5+j)
197 real*8 x,y, z(3,2) local
215 y = tmp(2)
229 if( y .gt. table(2,0,0,itype) .or.
230 & y .lt. table(2,0,numdata(2,itype)-1,itype) ) then
231 write(*,*) 'Sorry, concentration 2 asked for: ', y
242 j = int ( (y - table(2,0,0,itype) ) / dx(2,itype))
251 y = table(2,0,0,itype)
[all …]
/phasta/svLS/
H A DGMRES.f58 REAL*8, ALLOCATABLE :: u(:,:,:), h(:,:), unCondU(:,:), y(:), c(:), local
62 2 y(ls%sD), c(ls%sD), s(ls%sD), err(ls%sD+1))
140 y = err(1:i)
143 y(j) = y(j) - h(j,k)*y(k)
145 y(j) = y(j)/h(j,j)
149 X = X + u(:,:,j)*y(j)
180 REAL*8, ALLOCATABLE :: u(:,:,:), h(:,:), X(:,:), y(:), c(:), s(:), local
184 2 y(ls%sD), c(ls%sD), s(ls%sD), err(ls%sD+1))
243 y = err(1:i)
246 y(j) = y(j) - h(j,k)*y(k)
[all …]

123