Home
last modified time | relevance | path

Searched refs:u (Results 1 – 25 of 27) sorted by relevance

12

/phasta/svLS/
H A DGMRES.f58 REAL*8, ALLOCATABLE :: u(:,:,:), h(:,:), unCondU(:,:), y(:), c(:), local
61 ALLOCATE(h(ls%sD+1,ls%sD), u(dof,nNo,ls%sD+1), unCondU(dof,nNo),
70 u(:,:,1) = R
74 2 Val, X, u(:,:,1))
77 2 face, X, u(:,:,1))
79 u(:,:,1) = R - u(:,:,1)
82 unCondU = u(:,:,1)
84 2 face, unCondU, u(:,:,1))
87 err(1) = NORMV(dof, mynNo, commu, u(:,:,1))
100 u(:,:,1) = u(:,:,1)/err(1)
[all …]
/phasta/phSolver/compressible/
H A Dtnanq.f1 subroutine tnanq (u, n, arrname) argument
5 dimension u(nshg,n),rnan(2)
12 if (abs(u(i,j)).gt.1.0e10) nlarge=nlarge+1
13 if (u(i,j) .ne. u(i,j)) then
15 u(i,j)=9.876543e21
26 call write_restart(myrank,9876543,nshg,n,u,u)
H A Dtnanqe.f1 subroutine tnanqe (u, n, arrname) argument
5 dimension u(npro,n),rnan(2)
12 if(abs(u(i,j)).gt.1.0e10) nlarge=nlarge+1
13 IF (u(i,j) .ne. u(i,j)) then
16 u(i,j)=9.876543e21
25 write(8+myrank,245) (u(i,j), j=1,n)
H A De3eig1.f4 & dxidx, u, Q)
47 & u(npro), Q(npro,nflow,nflow)
56 where (u .ne. zero)
57 fact = one / u
175 Q(:,1,5) = fact * ((c + u) * c - eb1 * gamb)
176 Q(:,2,5) = -fact * (c + u * gamb) * Rcs(:,1,1)
177 Q(:,3,5) = -fact * (c + u * gamb) * Rcs(:,1,2)
178 Q(:,4,5) = -fact * (c + u * gamb) * Rcs(:,1,3)
181 Q(:,1,1) = fact * ((c - u) * c - eb1 * gamb)
182 Q(:,2,1) = fact * (c - u * gamb) * Rcs(:,1,1)
[all …]
H A Dget_h.c8 int intersect(double x[4][3],double u[],double points[4][3] );
12 get_h_(double xtmp[][4], double u[], double *h) in get_h_() argument
15 get_h_(double xtmp[][4], double u[], double *h) in get_h_()
18 get_h(double xtmp[][4], double u[], double *h) in get_h_()
21 void get_h(double xtmp[][4], double u[], double *h) in get_h_()
24 void get_h(double xtmp[][4], double u[], double *h) in get_h_()
27 void GET_H(double xtmp[][4], double u[], double *h) in get_h_()
37 *h = length(x,u); in get_h_()
40 double length(double x[4][3],double u[] ){ in length() argument
48 nint = intersect(x,u,points); in length()
[all …]
H A De3tau.f340 & u(npro), eb1(npro), Q(npro,5,5),
352 u = sqrt(u1**2+u2**2+u3**2)
374 & dxidx, u, Q)
381 call e3eig2 (u, c, a1,
802 & u(npro), eb1(npro), Q(npro,5,5),
814 u = sqrt(u1**2+u2**2+u3**2)
836 & dxidx, u, Q)
843 call e3eig2 (u, c, a1,
/phasta/phSolver/AMG/
H A Dramg_smoother.f7 & u,v)
21 real(kind=8),intent(inout),dimension(nshg) :: u local
42 call fLesSparseApNGtC(colm,rowp,lhsP,utmp4,u,nshg,nnz_tot)
43 call commIn(u,ilwork,1,iper,iBC,BC)
59 & r,u,v)
67 real(kind=8),intent(inout),dimension(Anshg) :: u local
76 u = 0
82 u(i) = u(i) - tmp
84 u(i) = damp_jacobi*u(i) + r(i)
85 u(i) = u(i)/Alhs(Acolm(i))
[all …]
H A Dramg_mls.f496 subroutine ramg_mls_calcPAv(level,u,v,colm,rowp,lhsK,lhsP, argument
513 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u local
522 call ramg_mls_sandw_post(u,aux2,level,colm,rowp,lhsK,lhsP,
563 subroutine ramg_mls_apply(u,v,r,level,colm,rowp,lhsK,lhsP, argument
581 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u local
598 call ramg_mls_apply_fwd(u,v,aux,level,colm,rowp,lhsK,lhsP,
601 call ramg_mls_apply_post(u,v,aux,level,colm,rowp,lhsK,lhsP,
605 u = v-1.1*u
610 subroutine ramg_mls_apply_fwd(u,v,r,level,colm,rowp,lhsK,lhsP, argument
628 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u local
[all …]
H A Dramg_tools.f107 & u,v,gcomm)
114 real(kind=8),intent(inout),dimension(gnshg) :: u local
119 u = 0
123 u(i) = u(i)+glhs(k)*v(j)
137 subroutine ramg_calcAv_g(level,u,v,colm,rowp,lhsK,lhsP, argument
154 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u local
164 call ramg_PPEAp(u,v,colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper)
169 u = 0
176 u(i) = u(i)+amg_A_lhs(level)%p(k,1)*v(j)
181 call ramg_commIn(u,level,ilwork,1,iper,iBC,BC)
[all …]
H A Dramg_cheby.f9 subroutine ramg_cheby_apply(u,v,r,level,colm,rowp,lhsK,lhsP, argument
27 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u local
63 u = chebyomega*(cy1+chebygamma*cz1)
65 u = chebyomega*(cy1-v+chebygamma*cz1)+v
99 u = epx
/phasta/phSolver/incompressible/
H A DitrPC.f65 subroutine itrPredict (yold, y, acold, ac, uold, u, iBC) argument
70 & u(nshg,nsd), yold(nshg,ndof), local
82 u(:,isd) = uold(:,isd) + Delt(itseq)*yold(:,isd) +
117 subroutine itrCorrect ( y, ac, u, solinc ,iBC) argument
122 & u(nshg,nsd), solinc(nshg,4) local
134 u(:,isd) = u(:,isd) +
210 & u, y, ac,
219 & u(nshg,nsd), uold(nshg,nsd), local
233 & + alfi * (u(:,1:3) - uold(:,1:3))
256 & y, ac, u )
[all …]
H A Dasbflx.f1 subroutine AsBFlx (u, y, ac, argument
19 & ac(nshg,ndofl), u(nshg,nsd),
52 call localx(u, ul, ienb, nsd, 'gather ')
H A Dwrite_hessian.c24 int u[6] = { 1, 10, 19, 13, 22, 25 }; in write_hessian() local
49 fprintf( uhess, "%f ", hessian[k+(u[j]-1)*(*nshg)] ); in write_hessian()
H A Dasbmfg.f1 subroutine AsBMFG (u, y, ac, x, argument
20 & ac(nshg,ndofl), u(nshg,nsd),
45 call localx(u, ul, ienb, nsd, 'gather ')
H A Dsolfar.f1 subroutine SolFlow(y, ac, u, argument
84 & u(nshg,nsd), uold(nshg,nsd), local
135 & u, y, ac,
341 subroutine SolSclr(y, ac, u, argument
393 & u(nshg,nsd), uold(nshg,nsd), local
430 & u, y, ac,
H A Dbflux.f1 subroutine Bflux ( y, ac, u, x, argument
30 & u(nshg,nsd), x(numnp,nsd) local
186 call AsBFlx (u, y,
H A Delmgmr.f1 subroutine ElmGMR (u, y, ac, x, argument
29 & u(nshg,nsd),
290 call AsBMFG (u, y,
/phasta/phSolver/common/
H A Dmpitools.f104 subroutine sumgat (u, n, summed) argument
110 dimension u(nshg,n), ilwork(nlwork)
113 summed = sum(u)
145 subroutine sumgatN (u, n, summed, nnp) argument
151 dimension u(nnp,n), ilwork(nlwork)
156 summed = sum(u) + 1.e-20
189 subroutine sumgatInt (u, n, summed ) argument
195 integer u(n), summed, sumvec local
200 summed = sum(u)
H A Dsettauw.f179 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
202 up=u/utau
214 write(*,*) u, y, rm
222 utau=sqrt(u/y*rm)
H A DslpwBC.f14 real*8 det,det3x3,tmp,u,v,w local
203 u=BC(3)
214 BC(4)=-wn1(1)/wn1(2)*u-wn1(3)/wn1(2)*w ! v is determined by u and w
216 BC(5)=-wn1(1)/wn1(3)*u-wn1(2)/wn1(3)*v ! w is determined by u and v
222 u=BC(3)
249 BC(4)=(c10*c7-c6*c11)/det-(c10*c4-c6*c8)/det*u ! v is determined by u
250 BC(5)=(c11*c5-c9*c7)/det-(c5*c8-c9*c4)/det*u ! w is determined by u
H A Dnewshape.cc614 for(int u=0;u<4;u++){ in HexShapeAndDrv() local
615 vrt[u]=face[f][u]; in HexShapeAndDrv()
/phasta/phSolver/testIC-DSIM/
H A Dgeom.sms166 3 24+f+u-4+1
167 3 24+g-9-5-u
173 0 34-q-x-l-u
H A Dgeom.sdm165 3 24+f+u-4+1
166 3 24+g-9-5-u
/phasta/phSolver/incompressible/e3source/
H A De3source.fRHSonly268 c phi, + |u - -----phi, | phi, - -----|(nu+phi)phi, |
281 c.... compute the global gradient of u
288 c i j indices match array where V is the velocity (u in our notes)
300 c a j u a i
413 c.... compute the global gradient of u
420 c i j indices match array where V is the velocity (u in our notes)
432 c a j u a i
523 cad d_{,t} + u * d_{,i} = sign(phi)
H A De3source.fimplicit268 c phi, + |u - -----phi, | phi, - -----|(nu+phi)phi, |
281 c.... compute the global gradient of u
288 c i j indices match array where V is the velocity (u in our notes)
300 c a j u a i
427 c.... compute the global gradient of u
434 c i j indices match array where V is the velocity (u in our notes)
446 c a j u a i
537 cad d_{,t} + u * d_{,i} = sign(phi)

12