Home
last modified time | relevance | path

Searched refs:shgl (Results 1 – 25 of 61) sorted by relevance

123

/phasta/phSolver/common/
H A Dshp10t.f1 subroutine shp10T (QuadPt, nQpt, shp, shgl, wght) argument
40 & shgl(nsd,nen,*), wght(*)
57 shgl(1,1,m) = -pt5
58 shgl(2,1,m) = -pt5
59 shgl(3,1,m) = -pt5
60 shgl(1,2,m) = pt5
61 shgl(2,2,m) = zero
62 shgl(3,2,m) = zero
63 shgl(1,3,m) = zero
64 shgl(2,3,m) = pt5
[all …]
H A Dshp4t.f1 subroutine shp4T (QuadPt, nQpt, shp, shgl, wght) argument
40 & shgl(nsd,nen,*), wght(*)
55 shgl(1,1,m) = -pt5
56 shgl(2,1,m) = -pt5
57 shgl(3,1,m) = -pt5
58 shgl(1,2,m) = pt5
59 shgl(2,2,m) = zero
60 shgl(3,2,m) = zero
61 shgl(1,3,m) = zero
62 shgl(2,3,m) = pt5
[all …]
H A De3metric.f8 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 Dgetstrl.f1 subroutine getstrl( y, x, ien, strnrm, shgl, shp ) argument
11 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd),
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 Dcmass.f20 subroutine cmass (shp, shgl, xl, em) argument
34 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT),
51 call getshp(shp, shgl, sgn,
73 subroutine cmassl (shp, shgl, shpf, shglf, xl, em, Qwtf) argument
87 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT),
149 call getshp(shp, shgl, sgn,
196 subroutine cmasstl (shp, shgl, shpf, shglf, xl, em, Qwtf) argument
210 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT),
266 call getshp(shp, shgl, sgn,
316 subroutine e3metricf( xl, shgl, dxidx, argument
[all …]
H A Dgenshp.f1 subroutine genshp (shp, shgl, nshp, nblk) argument
14 & shgl(MAXTOP,nsd,maxsh,MAXQPT)
32 call shpTet(ipord,Qpt(1,1:3,i),shp(1,:,i),shgl(1,:,:,i))
34 shgl(1,:,1:nshl,1:nint(lcsyst)) =
35 & shgl(1,:,1:nshl,1:nint(lcsyst))/two
42 & shgl(2,:,:,i))
50 & shgl(3,:,:,i))
57 call shppyr (ipord,Qpt(5,1:3,i),shp(5,:,i),shgl(5,:,:,i))
H A Dasithf.f1 subroutine asithf (y, x, strnrm, ien, fres, shgl, shp, Qwtf) argument
12 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd),
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 Dbardmc.f40 subroutine bardmc (y, shgl, shp, argument
71 & shgl(MAXTOP,nsd,maxsh,MAXQPT), shp(MAXTOP,maxsh,MAXQPT),
250 & fres, hfres, shgl(lcsyst,:,1:nshl,:),
269 & strl(iel:inum,:), shgl(lcsyst,:,1:nshl,:),
672 subroutine hfilter (y, x, ien, hfres, shgl, shp, Qwtf) argument
683 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd),
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)
[all …]
H A Dhierarchic.f30 subroutine getshp(shp, shgl, sgn, shape, shdrv) argument
37 dimension shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
45 shdrv(:,j,i) = shgl(j,i,intp)
52 shdrv(:,j,i) = shgl(j,i,intp)*sgn(:,i)
60 subroutine getshpb(shp, shgl, sgn, shape, shdrv) argument
67 dimension shp(nshl,ngaussb), shgl(nsd,nshl,ngaussb),
75 shdrv(:,j,i) = shgl(j,i,intp)
82 shdrv(:,j,i) = shgl(j,i,intp)*sgn(:,i)
244 & shp(nshl,npts), shgl(nsd,nshl,npts), local
258 call shpTet(ipord,intpnt(:,i),shp(:,i),shgl(:,:,i))
[all …]
H A Derrorsmooth.f2 & shp, shgl, iBC) argument
10 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
53 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:)
109 & shgl, ien, argument
135 & shgl(nsd,nshl,maxsh),
180 call getshp(shp, shgl, sgn,
/phasta/phSolver/incompressible/
H A De3qvar.f1 subroutine e3qvar (yl, shgl, 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 Dhessian.f2 & shp, shgl, iBC, argument
14 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
43 & shgl(lcsyst,:,1:nshl,:),
73 & shgl(lcsyst,:,1:nshl,:),
87 subroutine velocity_gradient ( y, x, shp, shgl, argument
93 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
133 call getshp( shp, shgl, sgn, shape, shdrv )
162 subroutine velocity_hessian ( gradu, x, shp, shgl, argument
168 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
208 call getshp( shp, shgl, sgn, shape, shdrv )
[all …]
H A Dfilters.f1 subroutine hfilterC (y, x, ien, hfres, shgl, shp, Qwtf) argument
19 & shgl(nsd,nshl,ngauss), shg(npro,nshl,nsd),
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 Dasiq.f2 & shgl, ien, xmudmi, argument
28 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
64 call e3q (yl, dwl, shp, shgl,
89 & shgl, ien, qres, argument
96 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
122 call e3qSclr (yl, dwl, shp, shgl,
H A Derrsmooth.f2 & shp, shgl, iBC) argument
10 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
54 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:)
110 & shgl, ien, argument
136 & shgl(nsd,nshl,maxsh),
181 call getshp(shp, shgl, sgn,
H A Dasigmr.f2 & shp, shgl, ien, argument
24 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
87 & shgl, xl, rl,
126 & shp, shgl, ien, argument
142 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
181 & shgl, xl, dwl,
H A De3.f2 & shgl, xl, rl, ql, argument
42 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
75 call e3ql (yl, dwl, shp, shgl,
89 call getshp(shp, shgl, sgn,
164 & shgl, xl, dwl, argument
180 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), local
211 call e3qlSclr (yl, dwl, shp, shgl, xl, ql, sgn)
222 call getshp(shp, shgl, sgn,
H A DasiqGradV.f2 & shgl, ien, argument
26 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss),
57 call e3qGradV (yl, shp, shgl,
H A Dgenlmass.f1 subroutine genlmass (x, shp,shgl) argument
11 & shgl(MAXTOP,nsd,maxsh,MAXQPT) local
43 tmpshgl(:,1:nshl,:) = shgl(lcsyst,:,1:nshl,:)
61 & shgl, ien, argument
75 & shp(nshl,maxsh), shgl(nsd,nshl,ngauss), local
114 call getshp(shp, shgl, sgn,
/phasta/phSolver/compressible/
H A Ditrfdi.f4 & ilwork, shp, shgl,
51 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
71 & shgl, iBC,
105 & shgl, iBC,
122 & shgl, iBC,
154 & ilwork, shp, shgl, wght,
210 & shgl, wght, iBC,
236 & shgl, wght, iBC,
248 & shgl, wght, iBC,
H A Ditrdrv.f4 & shgl, shpb, shglb, argument
55 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
215 call vortGLB(yold, x, shp, shgl, ilwork, vortG)
356 if(itvn.gt.0) call BCint((lstep+1)*Delt(1), shp, shgl,
369 call getdmc (yold, shgl, shp,
375 call bardmc (yold, shgl, shp,
381 call projdmc (yold, shgl, shp,
469 & shp, shgl,
485 & shp, shgl,
503 & shp, shgl,
[all …]
H A Dasiq.f2 & shgl, ien, xmudmi, argument
28 & shgl(nsd,nshl,ngauss),
58 call e3q (ycl, shp, shgl,
H A Dau2mfg.f5 & shp, shgl,
46 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
75 & shgl, iBC,
99 & shgl, iBC,
H A Dasigmr.f2 & shp, shgl, ien, argument
22 & shgl(nsd,nshl,MAXQPT),
76 & shgl, xl, rl, rml, xmudmi,
125 & shp, shgl, ien, argument
143 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT),
183 & sgn, shgl, xl,
H A Dau1mfg.f5 & shgl, shpb,
51 & shgl(MAXTOP,nsd,maxsh,MAXQPT),
74 & shgl, iBC,

123