| /phasta/phSolver/incompressible/ |
| H A D | e3sts.f | 1 subroutine e3StsLhs( xl, lStsVec ) argument
13 & xl(npro,nenl,3) local
15 call e3StsDir( xl, lDir )
37 subroutine e3StsRes( xl, rl, lStsVec ) argument
46 real*8 xl(npro,nenl,3), rl(npro,nshl,ndof) local
49 call e3StsDir( xl, lDir )
74 subroutine e3StsDir( xl, lDir ) argument
82 real*8 xl(npro,nenl,3), lDir(npro,nshl,3) local
92 x12 = xl(e,2,1) - xl(e,1,1)
93 x13 = xl(e,3,1) - xl(e,1,1)
[all …]
|
| H A D | e3qvar.f | 2 & xl, g1yi, argument
38 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd),
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)
[all …]
|
| H A D | hessian.f | 98 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
121 call localx (x, xl, ien, nsd, 'gather ')
134 call local_gradient( yl(:,:,2:4), 3, shdrv, xl,
173 dimension gradul(npro,nshl,9), xl(npro,nenl,nsd),
196 call localx (x, xl, ien, nsd, 'gather ')
209 call local_gradient( gradul, 9, shdrv, xl,
288 subroutine local_gradient ( vector, vsize, shgl, xl, argument
299 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd),
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)
[all …]
|
| H A D | asigmr.f | 31 & xl(npro,nenl,nsd), dwl(npro,nenl),
60 call localx(x, xl, ien, nsd, 'gather ')
87 & shgl, xl, rl,
94 call e3StsRes ( xl, rl, lStsVec )
111 call timeseries(yl,xl,ien,sgn)
148 & xl(npro,nenl,nsd), local
166 call localx(x, xl, ien, nsd, 'gather ')
181 & shgl, xl, dwl,
|
| H A D | asiq.f | 32 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
54 call localx (x, xl, ien, nsd, 'gather ')
65 & xl, ql, rmassl,
100 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
112 call localx (x, xl, ien, nsd, 'gather ')
123 & xl, ql, rmassl,
|
| H A D | filters.f | 10 dimension x(numnp,3), xl(npro,nenl,3)
26 call localx (x, xl, ien, 3, 'gather ')
44 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
45 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
46 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
47 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
48 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
49 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
50 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
51 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
[all …]
|
| H A D | getdiff.f | 1 subroutine getDiff( dwl,yl, shape, xmudmi, xl, rmu, rho) argument
12 & xl(npro,nenl,nsd)
97 call EviscDESIC (xl,rmu,xmudmi)
104 subroutine EviscDESIC(xl,xmut,xmudmi) argument
107 real*8 xmut(npro),xl(npro,nenl,nsd),xmudmi(npro,ngauss) local
111 dx=maxval(xl(i,:,1))-minval(xl(i,:,1))
112 dy=maxval(xl(i,:,2))-minval(xl(i,:,2))
113 dz=maxval(xl(i,:,3))-minval(xl(i,:,3))
|
| H A D | e3.f | 2 & shgl, xl, rl, ql, argument
43 & xl(npro,nenl,nsd), dwl(npro,nenl),
76 & xl, ql, xmudmi,
94 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho)
99 & shdrv, xl,
164 & shgl, xl, dwl, argument
181 & xl(npro,nenl,nsd), rl(npro,nshl), local
211 call e3qlSclr (yl, dwl, shp, shgl, xl, ql, sgn)
232 & shdrv, xl, xmudmi,
|
| H A D | elmStats.f | 14 real*8, allocatable :: xl(:,:,:) local
30 allocate ( xl(npro,nenl,3) )
35 call localx ( x, xl, mien(iblk)%p, nsd, 'gather ' )
39 call e3StsLhs( xl, lStsVec )
46 deallocate ( xl )
|
| H A D | asiqGradV.f | 30 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
51 call localx (x, xl, ien, nsd, 'gather ')
58 & xl, ql,
|
| H A D | e3ivar.f | 2 & shgl, xl, argument
57 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd),
102 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
103 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
104 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3)
143 call e3metric( xl, shgl, dxidx,
292 dkei = dkei + shpfun(:,n) * (1.0-xl(:,n,2)**2) !u_ex^~ (in FEM space)
315 & shgl, xl, xmudmi, argument
330 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd),
378 call e3metric( xl, shgl, dxidx,
[all …]
|
| H A D | e3q.f | 2 & xl, ql, rmassl, argument
26 & xl(npro,nenl,nsd),
74 & xl, g1yi,
87 call getdiff(dwl, yl, shape, xmudmi, xl, rmu, tmp)
207 & xl, ql, rmassl, argument
221 & xl(npro,nenl,nsd),
259 & xl, gradT,
|
| H A D | e3ql.f | 2 & xl, ql, xmudmi, argument
27 & xl(npro,nenl,nsd), sgn(npro,nshl),
62 & xl, g1yi,
66 call getdiff(dwl, yl, shape, xmudmi, xl,rmu, tmp)
178 & xl, ql, sgn ) argument
193 & xl(npro,nenl,nsd), sgn(npro,nshl),
224 call e3qvarSclr (yl, shdrv, xl,
|
| H A D | e3qGradV.f | 2 & xl, ql, argument
25 & xl(npro,nenl,nsd),
55 & xl, g1yi,
|
| /phasta/phSolver/common/ |
| H A D | e3metric.f | 8 subroutine e3metric( xl, shgl, dxidx, argument
13 real*8 xl(npro,nenl,nsd), shgl(npro,nsd,nshl), local
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)
[all …]
|
| H A D | cmass.f | 20 subroutine cmass (shp, shgl, xl, em) argument
36 & xl(npro,nenl,nsd)
56 call e3metric( xl, shdrv,dxidx,shg,WdetJ)
73 subroutine cmassl (shp, shgl, shpf, shglf, xl, em, Qwtf) argument
90 & xl(npro,nenl,nsd)
120 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf)
155 call e3metric( xl, shdrv,dxidx,shg,WdetJ)
196 subroutine cmasstl (shp, shgl, shpf, shglf, xl, em, Qwtf) argument
213 & xl(npro,nenl,nsd)
237 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf)
[all …]
|
| H A D | getstrl.f | 6 dimension x(numnp,3), xl(npro,nenl,3)
16 call localx (x, xl, ien, 3, 'gather ')
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)
[all …]
|
| H A D | asithf.f | 6 dimension x(numnp,nsd), xl(npro,nenl,nsd)
19 call localx (x, xl, ien, 3, 'gather ')
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)
[all …]
|
| H A D | bardmc.f | 677 dimension x(numnp,3), xl(npro,nenl,3)
690 call localx (x, xl, ien, 3, 'gather ')
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)
[all …]
|
| H A D | hierarchic.f | 248 real*8, allocatable :: xl(:,:,:) local
274 allocate ( xl(npro,nenl,nsd ) )
280 call localx( x, xl, mien(iblk)%p, nsd, 'gather ')
282 call eval ( xl, ycl, yvl,
296 deallocate ( xl )
317 subroutine eval( xl, ycl, yvl, argument
328 & shgl(nsd,nshl,npts), xl(npro,nenl,nsd), local
351 call e3metric( xl, shdrv, dxidx,
|
| /phasta/phSolver/compressible/ |
| H A D | getdiff.f | 3 & xmudmi, xl) argument
45 & xmudmi(npro,ngauss), xl(npro,nenl,nsd),
56 xx(:)=xx(:) + shp(:,n) * xl(:,n,1)
257 call EviscDES (xl,xmut,xmudmi)
364 subroutine EviscDES(xl,xmut,xmudmi) argument
367 real*8 xmut(npro),xl(npro,nenl,nsd),xmudmi(npro,ngauss) local
371 dx=maxval(xl(i,:,1))-minval(xl(i,:,1))
372 dy=maxval(xl(i,:,2))-minval(xl(i,:,2))
373 dz=maxval(xl(i,:,3))-minval(xl(i,:,3))
|
| H A D | e3source.f | 5 & xl, shpfun, bcool)
34 & xl(npro,nenl,3), xx(npro,3)
66 xfind(i) = (xsor-minval(xl(i,:,1)))
67 & *(maxval(xl(i,:,1))-xsor)
82 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
83 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
84 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3)
230 & u2, u3, xl, elDwl)
292 & xl(npro,nenl,nsd)
300 dx=maxval(xl(i,:,1))-minval(xl(i,:,1))
[all …]
|
| H A D | asigmr.f | 30 & xl(npro,nenl,nsd), ytargetl(npro,nshl,nflow),
54 call localx(x, xl, ien, nsd, 'gather ')
76 & shgl, xl, rl, rml, xmudmi,
111 call timeseries(ycl,xl,ien,sgn)
152 & xl(npro,nenl,nsd), Diagl(npro,nshl),
168 call localx (x, xl, ien, nsd, 'gather ')
183 & sgn, shgl, xl,
|
| H A D | asiq.f | 34 dimension ycl(npro,nshl,ndof), xl(npro,nenl,nsd),
51 call localx(x, xl, ien, nsd, 'gather ')
59 & xl, ql, rmassl,
|
| H A D | asiflx.f | 23 & xl(npro,nenl,nsd),
45 call localx(x, xl, ien, nsd, 'gather ')
64 & shgl(lcsyst,:,1:nshl,:), xl, rl, rml,
|