Home
last modified time | relevance | path

Searched refs:x (Results 1 – 25 of 103) sorted by relevance

12345

/phasta/phSolver/incompressible/
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)
307 subroutine flesInv ( x, m, n ) argument
311 real*8 x(n,m) local
315 if ( x(i,j) .ne. 0 ) x(i,j) = 1. / x(i,j)
328 subroutine fMtxBlkDot2( x, y, c, m, n ) argument
[all …]
H A DBCprofile.f77 subroutine BCprofileInit(vbc_prof,x) argument
81 real*8 vbc_prof(nshg,3), x(numnp,nsd) local
93 vbc_prof(kk,2)=-1.6*1E7*(x(kk,1)-x1)*(x(kk,1)-x2)*
94 & (x(kk,3)-z1)*(x(kk,3)-z2)
106 987 format(i6,3(2x,e14.7))
H A Dwrite_hessian.c10 int x; in write_hessian() local
28 for( x=0; x < *nshg; x++ ) fscanf(idmap,"%d\n", map+x ); in write_hessian()
H A Dasiq.f1 subroutine AsIq (y, x, shp, argument
27 dimension y(nshg,ndof), x(numnp,nsd),
54 call localx (x, xl, ien, nsd, 'gather ')
88 subroutine AsIqSclr (y, x, shp, argument
95 dimension y(nshg,ndof), x(numnp,nsd),
112 call localx (x, xl, ien, nsd, 'gather ')
H A DelmStats.f1 subroutine elmStatsLhs( x, iBC, iper, ilwork ) argument
11 real*8 x(numnp,3) local
35 call localx ( x, xl, mien(iblk)%p, nsd, 'gather ' )
79 subroutine elmStatsRes( y, ac, x, shp, shgl, argument
88 real*8 y(nshg,ndof), ac(nshg,ndof), x(numnp,nsd), local
107 call ElmGMR (y, ac, x,
H A Dasbmfg.f1 subroutine AsBMFG (u, y, ac, x, argument
19 dimension y(nshg,ndofl), x(numnp,nsd),
44 call localx(x, xlb, ienb, nsd, 'gather ')
87 subroutine AsBSclr (y, x, shpb, shglb, argument
93 dimension y(nshg,ndofl), x(numnp,nsd),
113 call localx(x, xlb, ienb, nsd, 'gather ')
/phasta/phSolver/common/
H A Drenum.f1 subroutine renum_cyl(x) argument
30 dimension x(numnp,nsd), nrin(numnp), nula(numnp),
46 xcyl(:,1) = sqrt(x(:,1)*x(:,1) + x(:,2)*x(:,2))
50 if ((x(i,1).eq.0).and.(x(i,2).eq.0)) then
54 xcyl(i,2) = atan2(x(i,2),x(i,1))
57 xcyl(:,3) = x(:,3)
67 if((x(nen1(i),1).eq.0).and.(x(nen1(i),2).eq.0)) then
152 xts1 = x(nrint(jj),1) !*cos(xcyl(nrint(jj),2)+tolerence)
153 xts2 = x(nrint(jj),2) !*sin(xcyl(nrint(jj),2)+tolerence)
154 xts3 = x(nrint(jj),3) + plandist
[all …]
H A Dgendat.f1 subroutine gendat (y, ac, x, iBC, BC, argument
27 & x(numnp,nsd), iBC(nshg),
66 call setInitial_Duct(x) ! in INIprofile.f
74 call xyzbound(x)
81 write (iecho,1100) n, (x(n,i),i=1,nsd)
136 … call setBlowing_Duct3(x,BC,iTurbWall) ! in setBlowing_Duct3.f, fixed jet inlet velocity
145 call findWallNorm(x,iBC,ilwork,iper)
147 call setSuction_Duct3(x, BC, y, ilwork)
209 call aveprep(shp,x)
229 & nsons, x,
[all …]
H A Dmod_wallData.f9021 subroutine findWallNorm(x,iBC,ilwork,iper) argument
29 real*8 x(numnp,nsd) local
62 aedg(:)=x(ienb(j+1),:)-x(nn,:)
63 bedg(:)=x(ienb(j-1),:)-x(nn,:)
65 aedg(:)=x(ienb(j+1),:)-x(nn,:)
66 bedg(:)=x(ienb(nenbl),:)-x(nn,:)
68 aedg(:)=x(ienb(1),:)-x(nn,:)
69 bedg(:)=x(ienb(j-1),:)-x(nn,:)
119 if(((dabs(dabs(x(nn,3))-zwall)).lt.0.000000001d00)
126 & ", x = (", x(nn,1), x(nn,2), x(nn,3), ")",
[all …]
H A DINIprofile.f1 subroutine INIprofile(BC,y,x) argument
6 real*8 x(nshg,nsd) local
9 call TakeBCfromIC_ScaleInlVelX(BC,y,x)
13 call setInitial(x,y)
19 subroutine TakeBCfromIC_ScaleInlVelX(BC,y,x) argument
29 real*8 x(nshg,nsd) local
68 subroutine setInitial(x,y) argument
76 real*8 x(nshg,nsd) local
92 xinlet=x(inlf(1),1)
98 xoutlet=x(outf(1),1)
[all …]
H A Deqn_plane.f1 subroutine eqn_plane(x, iBC) argument
30 dimension x(numnp,nsd), nrin(numnp),
78 xA(:) = x(ienb(inum,1),:)
79 xB(:) = x(ienb(inum,2),:)
80 xC(:) = x(ienb(inum,3),:)
81 xD(:) = x(ienb(inum,4),:)
154 erreur(k) = aR - xnrml*x(ien(i,k),1)
155 & - ynrml*x(ien(i,k),2) - znrml*x(ien(i,k),3)
180 xintl(i,:,1) = x(ien2D(i,:),1)
181 xintl(i,:,2) = x(ien2D(i,:),2)
[all …]
H A Dsettauw.f1 subroutine settauw (y, x, BC, argument
16 dimension y(nshg,ndofl), x(numnp,nsd),
46 & (x(otwn(nodw),1)-x(nodw,1))*wnrm(nodw,1)
47 & +(x(otwn(nodw),2)-x(nodw,2))*wnrm(nodw,2)
48 & +(x(otwn(nodw),3)-x(nodw,3))*wnrm(nodw,3)
54 ut=utau(ub,dw,rm,nodw,x(nodw,:)) ! find utau
108 BC(nodw,7)=savarw(up,rm,saCv1,nodw,x(nodw,:))
185 function utau(u,y,rm,nodw,x)
191 real*8 x(3) local
212 write(*,*)'utau failed to converge at ',nodw,x(1),x(2),x(3)
[all …]
H A Dfindslpw.f1 subroutine findslpw(x,ilwork,iper,iBC) argument
6 real*8 x(numnp,nsd) local
50 aedg(:)=x(ienb(j+1),:)-x(nn,:)
51 bedg(:)=x(ienb(j-1),:)-x(nn,:)
53 aedg(:)=x(ienb(j+1),:)-x(nn,:)
54 bedg(:)=x(ienb(nenbl),:)-x(nn,:)
56 aedg(:)=x(ienb(1),:)-x(nn,:)
57 bedg(:)=x(ienb(j-1),:)-x(nn,:)
H A Dgenini.f2 & ifath, velbar, nsons,x, argument
34 & ac(nshg,ndof), x(numnp,nsd),
60 call INIprofile(BC,y,x)
82 call initBCt( x, iBC, BC)
83 call BCint(lstep*Delt(1),shp,shgl,shpb,shglb,x, BC, iBC)
154 call genscale(y, x, iBC)
175 & ' Number ',6x,6('dof',i1,:,10x))
176 1100 format(1p,2x,i5,5x,5(e12.5,2x))
H A Dgenbc.f1 subroutine genBC (iBC, BC, x, ilwork, iper) argument
31 & x(numnp,nsd), ilwork(nlwork),
75 call genwnm (iBC, BCtmp, x, ilwork, iper, nsurf)
79 call genotwn (x,BCtmp, iBC, nsurf)
102 call findslpw(x,ilwork,iper,iBC)
147 & ' Number ',5x,6('BC',i1,:,10x))
148 1100 format(1p,2x,i5,3x,6(e12.5,1x))
161 subroutine genwnm (iBC, BCtmp, x, ilwork, iper, nsurf) argument
183 & x(numnp,nsd), ilwork(nlwork)
245 e1 = x(ienb(2),:)-x(ienb(1),:)
[all …]
H A DBCprofile2.f28 subroutine setInlet_Duct(x,BC,iTurbWall) argument
35 real*8 x(nshg,nsd) local
54 xcoor=x(nn,1)
55 ycoor=x(nn,2)
56 zcoor=x(nn,3)
H A DValType.cc59 double x = atof(str.c_str()); in get_double() local
60 return x; in get_double()
73 double *x = new double[n]; in get_double_array() local
75 x[i] = vec[i]; in get_double_array()
77 return x; in get_double_array()
H A Dbctint.f66 subroutine initSponge( y,x) argument
71 real*8 y(nshg,nflow), x(numnp,3) local
102 radsqr=x(id,2)**2+x(id,1)**2
106 zval=max(x(id,3),zstart)
128 subroutine initBCt( x, iBC, BC ) argument
133 real*8 x(numnp,nsd), BC(nshg,ndofBC), rj1,rj2,rj3,rj4,distd,epsd local
196 dd= distds(x1,x2,x3,x(i,1),x(i,2),x(i,3))
297 subroutine BCint(timel,shp,shgl,shpb,shglb,x,BC,iBC) argument
304 real*8 x(numnp,nsd), local
339 real*8 x1,y1,z1,x2,y2,z2,x,y,z local
[all …]
H A DsetBlowing_Duct3.f3 subroutine setBlowing_Duct3(x,BC,iTurbWall) argument
12 real*8 x(nshg,nsd) local
42 xcoor=x(nn,1)
43 ycoor=x(nn,2)
44 zcoor=x(nn,3)
/phasta/phSolver/compressible/
H A Dtimer.f160 1100 format(1x,
162 1200 format(1x,a8,4x,i7,5x,f11.2,8x,f7.2)
163 1300 format(1x,a8,16x,f11.2,8x,f7.2)
167 2100 format(1x,'No. of floating-point operations : ',i20,/)
168 2200 format(1x,'Flop rate with gather/scatter : ',f5.1,
170 2300 format(1x,'Flop rate without gather/scatter :',f5.1,
172 2400 format(1x,'No. of bytes gathered : ',i20,/)
173 2500 format(1x,'No. of bytes scattered : ',i20,/)
174 2600 format(1x,'Gather transfer rate : ',f6.2, ' Mbytes/s',/)
175 2700 format(1x,'Scatter transfer rate :',f6.2, ' Mbytes/s')
H A Dget_h.c8 int intersect(double x[4][3],double u[],double points[4][3] );
30 double x[4][3]; in get_h_() local
35 x[i][j] = xtmp[j][i]; in get_h_()
37 *h = length(x,u); in get_h_()
40 double length(double x[4][3],double u[] ){ in length()
48 nint = intersect(x,u,points); in length()
97 int intersect(double x[4][3],double u[],double points[4][3] ){ in intersect()
112 xc[i] = 0.25*(x[0][i]+x[1][i]+x[2][i]+x[3][i]); in intersect()
119 xbar[j][k] = x[f_index[face][j]][k]; in intersect()
H A Ditrfdi.f1 subroutine itrFDI (ypre, y, ac, x, argument
40 & x(numnp,nsd), ac(nshg,ndof),
70 & x, shp,
104 & x, shp,
121 & x, shp,
151 subroutine itrFDISclr (y, ypre, x, argument
190 & x(numnp,nsd),
209 & x, shp,
235 & x, shp,
247 & x, shp,
H A Ditrdrv.f1 subroutine itrdrv (y, ac, uold, x, argument
47 & x(numnp,nsd), iBC(nshg),
215 call vortGLB(yold, x, shp, shgl, ilwork, vortG)
301 call initBCprofileScale(vbc_prof,BC,yold,x)
357 & shpb, shglb, x, BC, iBC)
371 & ifath, x)
377 & nsons, ifath, x)
382 & iper, ilwork, x)
392 call asbwmod(yold, acold, x, BC, iBC,
463 & x,
[all …]
H A Dd2wall.f2 subroutine D2wall (dwall, x) argument
17 dimension x(numnp,nsd), dwall(numnp),
54 xwi(nwalli,node,1:3)=x(mienb(iblk)%p(j,node),:)
112 distance = ( x(i,1) - xw(j,k,1) )**2
113 & +( x(i,2) - xw(j,k,2) )**2
114 & +( x(i,3) - xw(j,k,3) )**2
H A Dasbmfg.f1 subroutine AsBMFG (y, x, shpb, shglb, argument
15 dimension y(nshg,ndofl), x(numnp,nsd),
40 call localx(x, xlb, ienb, nsd, 'gather ')
68 subroutine AsBMFGSclr (y, x, shpb, shglb, argument
82 dimension y(nshg,ndofl), x(numnp,nsd),
103 call localx (x, xlb, ienb, nsd, 'gather ')

12345