Home
last modified time | relevance | path

Searched refs:n (Results 1 – 25 of 94) sorted by relevance

1234

/phasta/phSolver/incompressible/
H A De3ivar.f86 do n = 1, nshl
87 pres = pres + shpfun(:,n) * yl(:,n,1)
88 u1 = u1 + shpfun(:,n) * yl(:,n,2)
89 u2 = u2 + shpfun(:,n) * yl(:,n,3)
90 u3 = u3 + shpfun(:,n) * yl(:,n,4)
94 do n = 1, nshl
95 Temp = Temp + shpfun(:,n) * yl(:,n,5)
101 do n = 1,nenl
102 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
103 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
[all …]
H A De3qvar.f52 do n = 1, nenl
53 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n)
54 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n)
55 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n)
56 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n)
57 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n)
58 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n)
59 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n)
60 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n)
61 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(:,3,n)
[all …]
H A DlesSparse.f30 integer n, k
38 do n = 1, nshg
39 k = sparseloc( rowp(colm(n)), colm(n+1)-colm(n), n )
40 & + colm(n)-1
42 flowdiag(n,1) = lhsK(1,k)
43 flowdiag(n,2) = lhsK(5,k)
44 flowdiag(n,3) = lhsK(9,k)
46 flowdiag(n,4) = lhsP(4,k)
50 do n = 1, nshg ! rowsum put on the diagonal instead of diag entry
51 do k=colm(n),colm(n+1)-1
[all …]
H A Dgetdiff.f13 integer n, e
50 do n = 1, nshl
51 Sclr = Sclr + shape(:,n) * yl(:,n,isc)
137 integer n, e local
205 integer e, n
214 do n = 1, nshl
215 u1=yl(e,n,2)
216 u2=yl(e,n,3)
217 u3=yl(e,n,4)
220 wallnode(n)=.true.
[all …]
H A Dftools.f52 subroutine fMtxVdimVecMult( a, b, c, na, nb, nc, m, n ) argument
57 integer na, nb, nc, m, n local
58 real*8 a(n,na), b(n,nb), c(n,nc)
68 do i = 1, n
74 do i = 1, n
81 do i = 1, n
89 do i = 1, n
98 do i = 1, n
113 subroutine flesZero ( a, m, n ) argument
116 integer m, n, i, j local
[all …]
H A De3source.f267 integer e,n local
288 do n = 1, nshl
293 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2)
294 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3)
295 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4)
297 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2)
298 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3)
299 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4)
301 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2)
302 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3)
[all …]
H A Dhessian.f312 do n = 1, nenl
313 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n)
314 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n)
315 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n)
316 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n)
317 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n)
318 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n)
319 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n)
320 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n)
321 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(:,3,n)
[all …]
H A Dgenlmass.f181 integer n, k local
192 do n = 1, nshg
193 k = sparseloc( rowp(colm(n)), colm(n+1)-colm(n), n )
194 & + colm(n)-1
195 tmp=gmass(n)*tmp1*rho(n)
223 integer n, k local
234 do n = 1, nshg
235 k = sparseloc( rowp(colm(n)), colm(n+1)-colm(n), n )
236 & + colm(n)-1
237 tmp=gmass(n)*tmp1*rho(n)
/phasta/phSolver/compressible/
H A De3ivar.f96 do n = 1, nshl
97 dui(:,1) = dui(:,1) + shape(:,n) * yl(:,n,1) ! p
98 dui(:,2) = dui(:,2) + shape(:,n) * yl(:,n,2) ! u1
99 dui(:,3) = dui(:,3) + shape(:,n) * yl(:,n,3) ! u2
100 dui(:,4) = dui(:,4) + shape(:,n) * yl(:,n,4) ! u3
101 dui(:,5) = dui(:,5) + shape(:,n) * yl(:,n,5) ! T
114 do n = 1, nshl
115 Sclr = Sclr + shape(:,n) * ycl(:,n,isc)
147 do n = 1, nshl
151 pres = pres + shape(:,n) * ycl(:,n,1)
[all …]
H A Di3pre.f47 n = iel - 1
67 EGmass(n+iv,i+1,j) = EGmass(n+iv,i+1,j)
69 EGmass(n+iv,i+2,j) = (EGmass(n+iv,i+2,j)
70 & - BDiagl(n+iv,inode,2,1) * EGmass(n+iv,i+1,j))
72 EGmass(n+iv,i+3,j) = (EGmass(n+iv,i+3,j)
73 & - BDiagl(n+iv,inode,3,1) * EGmass(n+iv,i+1,j)
74 & - BDiagl(n+iv,inode,3,2) * EGmass(n+iv,i+2,j))
76 EGmass(n+iv,i+4,j) = (EGmass(n+iv,i+4,j)
77 & - BDiagl(n+iv,inode,4,1) * EGmass(n+iv,i+1,j)
78 & - BDiagl(n+iv,inode,4,2) * EGmass(n+iv,i+2,j)
[all …]
H A De3qvar.f68 do n = 1, nshl
72 pres = pres + shp(:,n) * ycl(:,n,1)
73 u1 = u1 + shp(:,n) * ycl(:,n,2)
74 u2 = u2 + shp(:,n) * ycl(:,n,3)
75 u3 = u3 + shp(:,n) * ycl(:,n,4)
76 T = T + shp(:,n) * ycl(:,n,5)
88 do n = 1, nshl
89 Sclr = Sclr + shp(:,n) * ycl(:,n,isc)
114 do n = 1, nshl
120 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * ycl(:,n,1)
[all …]
H A De3bvar.f95 do n = 1, nshlb
96 nodlcl = lnode(n)
114 do n = 1, nshlb
115 temp = temp + shpb(:,n) * ycl(:,n,isc)
133 do n = 1, nenl
134 dxdxib(:,1,1) = dxdxib(:,1,1) + xlb(:,n,1) * shglb(:,1,n)
135 dxdxib(:,1,2) = dxdxib(:,1,2) + xlb(:,n,1) * shglb(:,2,n)
136 dxdxib(:,1,3) = dxdxib(:,1,3) + xlb(:,n,1) * shglb(:,3,n)
137 dxdxib(:,2,1) = dxdxib(:,2,1) + xlb(:,n,2) * shglb(:,1,n)
138 dxdxib(:,2,2) = dxdxib(:,2,2) + xlb(:,n,2) * shglb(:,2,n)
[all …]
H A Dvorticity.f130 do n=1,nshl
131 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2) ! du/dx
132 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) ! du/dy
133 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) ! du/dz
135 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3) ! dv/dx
136 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) ! dv/dy
137 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) ! dv/dz
139 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4) ! dw/dx
140 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4) ! dw/dy
141 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4) ! dw/dz
[all …]
/phasta/phSolver/common/
H A De3metric.f24 do n = 1, nenl
25 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n)
26 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n)
27 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n)
28 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n)
29 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n)
30 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n)
31 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n)
32 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n)
33 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(:,3,n)
[all …]
H A Dlocal.f1 subroutine local (global, rlocal, ientmp, n, code) argument
23 dimension global(nshg,n), rlocal(npro,nshl,n),
44 do j = 1, n
54 gbytes = gbytes + n*nshl*npro
67 do j = 1, n
79 sbytes = sbytes + n*nshl*npro
80 flops = flops + n*nshl*npro
93 do j = 1, n
114 subroutine localx (global, rlocal, ien, n, code) argument
136 dimension global(numnp,n), rlocal(npro,nenl,n),
[all …]
H A Dgetstrl.f36 do n = 1, nenl
37 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
38 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
39 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
40 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
41 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
42 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
43 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
44 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
45 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
H A Dhierarchic.f103 do n = 1, nenbl
104 lnode(n) = n
108 do n=1,3*nem
109 lnode(nenbl+n) = nenbl+1+n
114 do n=1,nfm
115 lnode(3+3*nem+n) = 4+6*nem+n
122 do n = 1, nenbl
123 lnode(n) = n
127 do n=1,4*nem
128 lnode(nenbl+n) = 8+n
[all …]
H A Dasithf.f42 do n = 1, nenl
43 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
44 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
45 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
46 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
47 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
48 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
49 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
50 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
51 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
H A Df3lhs.f101 do n = 1, nenl ! when changed to nshl solution degraded ipord 10
102 flhsl(:,n,1) = flhsl(:,n,1) + WdetJb * shape(:,n)
105 fnrml(:,n,1) = fnrml(:,n,1) + WdetJb * bnorm(:,1)
106 & * shape(:,n)
107 fnrml(:,n,2) = fnrml(:,n,2) + WdetJb * bnorm(:,2)
108 & * shape(:,n)
109 fnrml(:,n,3) = fnrml(:,n,3) + WdetJb * bnorm(:,3)
110 & * shape(:,n)
128 do n = 1, nshl
129 temp = temp + flhsl(:,n,1)
[all …]
H A Dcmass.f332 do n = 1, nenl
333 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n)
334 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n)
335 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n)
336 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n)
337 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n)
338 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n)
339 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n)
340 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n)
341 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(:,3,n)
[all …]
H A Dcommu.f1 subroutine commu (global, ilwork, n, code) argument
37 dimension global(nshg,n),
38 & rtemp(maxfront*n,maxtask),
57 if (n .eq. 1) then ! like a scalar
59 elseif (n .eq. nsd) then ! like the normal vectors
61 elseif (n .eq. ndof) then ! res, y, ac, krylov vectors....
63 elseif (n .eq. nflow*nflow) then ! bdiag
65 elseif (n .eq. (nflow-1)*nsd) then ! qres
67 elseif (n .eq. nflow) then
69 elseif (n .eq. 24 ) then
[all …]
H A Dbardmc.f712 do n = 1, nenl
713 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
714 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
715 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
716 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
717 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
718 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
719 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
720 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
721 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
[all …]
/phasta/AcuStat/src/
H A DcommuMax.f1 subroutine commuMax (global, ilwork, n, code) argument
37 dimension global(nshg,n),
38 & rtemp(maxfront*n,maxtask),
53 if (n .eq. 1) then ! like a scalar
55 elseif (n .eq. nsd) then ! like the normal vectors
57 elseif (n .eq. ndof) then ! res, y, ac, krylov vectors....
59 elseif (n .eq. nflow*nflow) then ! bdiag
61 elseif (n .eq. (nflow-1)*nsd) then ! qres
63 elseif (n .eq. nflow) then
65 elseif (n .eq. 24 ) then
[all …]
/phasta/M2NFixBnd/src/
H A DcommuMax.f1 subroutine commuMax (global, ilwork, n, code) argument
37 dimension global(nshg,n),
38 & rtemp(maxfront*n,maxtask),
53 if (n .eq. 1) then ! like a scalar
55 elseif (n .eq. nsd) then ! like the normal vectors
57 elseif (n .eq. ndof) then ! res, y, ac, krylov vectors....
59 elseif (n .eq. nflow*nflow) then ! bdiag
61 elseif (n .eq. (nflow-1)*nsd) then ! qres
63 elseif (n .eq. nflow) then
65 elseif (n .eq. 24 ) then
[all …]
/phasta/M2N/src/
H A Dcommu.f1 subroutine commu (global, ilwork, n, code) argument
37 dimension global(nshg,n),
38 & rtemp(maxfront*n,maxtask),
51 if (n .eq. 1) then ! like a scalar
53 elseif (n .eq. nsd) then ! like the normal vectors
55 elseif (n .eq. ndof) then ! res, y, ac, krylov vectors....
57 elseif (n .eq. nflow*nflow) then ! bdiag
59 elseif (n .eq. (nflow-1)*nsd) then ! qres
61 elseif (n .eq. nflow) then
63 elseif (n .eq. 24 ) then
[all …]

1234