| /phasta/phSolver/common/ |
| H A D | shp10t.f | 1 subroutine shp10T (QuadPt, nQpt, shp, shgl, wght) argument 39 dimension QuadPt(4,*), shp(nen,*), 50 shp(1,m) = one - QuadPt(1,m) - QuadPt(2,m) - QuadPt(3,m) 51 shp(2,m) = QuadPt(1,m) 52 shp(3,m) = QuadPt(2,m) 53 shp(4,m) = QuadPt(3,m)
|
| H A D | shp4t.f | 1 subroutine shp4T (QuadPt, nQpt, shp, shgl, wght) argument 39 dimension QuadPt(4,*), shp(nen,*), 48 shp(1,m) = one - QuadPt(1,m) - QuadPt(2,m) - QuadPt(3,m) 49 shp(2,m) = QuadPt(1,m) 50 shp(3,m) = QuadPt(2,m) 51 shp(4,m) = QuadPt(3,m)
|
| H A D | genshp.f | 1 subroutine genshp (shp, shgl, nshp, nblk) argument 13 dimension shp(MAXTOP,maxsh,MAXQPT), 32 call shpTet(ipord,Qpt(1,1:3,i),shp(1,:,i),shgl(1,:,:,i)) 41 call shphex (ipord, Qpt(2,1:3,i),shp(2,:,i), 49 call shp6w (ipord,Qpt(3,1:3,i),shp(3,:,i), 57 call shppyr (ipord,Qpt(5,1:3,i),shp(5,:,i),shgl(5,:,:,i))
|
| H A D | hierarchic.f | 30 subroutine getshp(shp, shgl, sgn, shape, shdrv) argument 37 dimension shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 43 shape(:,i) = shp(i,intp) 50 shape(:,i) = sgn(:,i) * shp(i,intp) 60 subroutine getshpb(shp, shgl, sgn, shape, shdrv) argument 67 dimension shp(nshl,ngaussb), shgl(nsd,nshl,ngaussb), 73 shape(:,i) = shp(i,intp) 80 shape(:,i) = sgn(:,i) * shp(i,intp) 244 & shp(nshl,npts), shgl(nsd,nshl,npts), local 258 call shpTet(ipord,intpnt(:,i),shp(:,i),shgl(:,:,i)) [all …]
|
| H A D | aveprep.f | 27 subroutine aveprep(shp,x) argument 34 dimension shp(MAXTOP,maxsh,MAXQPT) 52 call getylist( ylist, ifathe(iel:inum,:), shp(lcsyst,1:nshl,:), 108 subroutine getxnut (fres, xden, xnum, ien, shp) argument 113 & shp(nshl,maxsh), fresli(npro,22), 132 fresli(:,j) = fresli(:,j) + shp(i,intp)*fresl(:,i,j) 208 subroutine getylist (ylist, lfathe, shp, x, ien) argument 212 dimension lfathe(npro,maxsh), shp(nshl,maxsh), 230 yint = yint + xl(nel,n,2)*shp(n,intp) 231 xint = xint + xl(nel,n,1)*shp(n,intp) [all …]
|
| H A D | bardmc.f | 40 subroutine bardmc (y, shgl, shp, argument 71 & shgl(MAXTOP,nsd,maxsh,MAXQPT), shp(MAXTOP,maxsh,MAXQPT), 251 & shp(lcsyst,1:nshl,:),Qwtf(lcsyst,1:ngaussf)) 270 & shp(lcsyst,1:nshl,:)) 561 & mxmudmi(iblk)%p,cdelsq,shp(lcsyst,1:nshl,:)) 672 subroutine hfilter (y, x, ien, hfres, shgl, shp, Qwtf) argument 684 & shp(nshl,maxsh), Qwtf(ngaussf), 756 rho = rho+shp(i,intp)*yl(:,i,1) !density at qpt 765 u1 = u1 + shp(i,intp)*yl(:,i,2) 766 u2 = u2 + shp(i,intp)*yl(:,i,3) [all …]
|
| H A D | scatnu.f | 1 subroutine scatnu (ien, strl, xmudmi, xnut, shp) argument 6 & xmudmi(npro,ngauss), shp(nshl,ngauss) 16 & *shp(in,int)
|
| H A D | errorsmooth.f | 2 & shp, shgl, iBC) argument 9 dimension shp(MAXTOP,maxsh,MAXQPT), 52 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 108 subroutine smooth (rerr, x, shp, argument 134 & shp(nshl,maxsh), 180 call getshp(shp, shgl, sgn,
|
| H A D | cmass.f | 20 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,
|
| /phasta/phSolver/compressible/ |
| H A D | e3wmlt.f | 1 subroutine e3wmlt (shp, shg, WdetJ, argument 53 dimension shp(npro,nshl), 102 & + shp(:,i) * rmi(:,16) ) 106 & + shp(:,i) * rmi(:,17) ) 110 & + shp(:,i) * rmi(:,18) ) 114 & + shp(:,i) * rmi(:,19) ) 118 & + shp(:,i) * rmi(:,20) ) 129 rl(:,i,2) = rl(:,i,2) + shp(:,i) * WdetJ * ri(:,17) 130 rl(:,i,3) = rl(:,i,3) + shp(:,i) * WdetJ * ri(:,18) 131 rl(:,i,4) = rl(:,i,4) + shp(:,i) * WdetJ * ri(:,19) [all …]
|
| H A D | vorticity.f | 1 subroutine vortGLB (y, x, shp, shgl, ilwork, vortG) argument 9 real*8 shp(MAXTOP,maxsh,MAXQPT), local 40 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 86 real*8 shp(npro,nshl), shgl(npro,nsd,nshl) local 120 call getshp (shp2, shgl2, sgn, shp, shgl ) 143 alph2(:) = alph2(:) + shp(:,n)*shp(:,n)*WdetJ(:) 151 lpmass(:,n) = lpmass(:,n) + WdetJ(:)*shp(:,n)*shp(:,n) 155 & *WdetJ(:)*shp(:,n) 159 & *WdetJ(:)*shp(:,n) 163 & *WdetJ(:)*shp(:,n) [all …]
|
| H A D | e3qvar.f | 1 subroutine e3qvar (ycl, shp, shgl, rho, argument 41 & shp(npro,nshl), tmp(npro), 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) 89 Sclr = Sclr + shp(:,n) * ycl(:,n,isc)
|
| H A D | itrfdi.f | 4 & ilwork, shp, shgl, 50 dimension shp(MAXTOP,maxsh,MAXQPT), 70 & x, shp, 104 & x, shp, 121 & x, shp, 154 & ilwork, shp, shgl, wght, 209 & x, shp, 235 & x, shp, 247 & x, shp,
|
| H A D | itrdrv.f | 3 & iper, ilwork, shp, argument 54 dimension shp(MAXTOP,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 D | asiq.f | 1 subroutine AsIq (y, x, shp, argument 27 & shp(nshl,ngauss), 58 call e3q (ycl, shp, shgl,
|
| H A D | au2mfg.f | 5 & shp, shgl, 45 dimension shp(MAXTOP,maxsh,MAXQPT), 74 & x, shp, 98 & x, shp,
|
| H A D | asigmr.f | 2 & shp, shgl, ien, argument 21 & shp(nshl,MAXQPT), 75 call e3 (ycl, ycl, acl, shp, 125 & shp, shgl, ien, argument 143 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT), 182 & dwl, elDwl, shp,
|
| /phasta/phSolver/incompressible/ |
| H A D | filters.f | 1 subroutine hfilterC (y, x, ien, hfres, shgl, shp, Qwtf) argument 20 & shp(nshl,ngauss), 116 u1 = u1 + shp(i,intp)*yl(:,i,2) 117 u2 = u2 + shp(i,intp)*yl(:,i,3) 118 u3 = u3 + shp(i,intp)*yl(:,i,4) 178 & hfres, shgl, shp, Qwtf) argument 190 & shp(nshl,ngauss), 304 fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1) ! hat{u1} 305 fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2) ! hat{u2} 306 fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3) ! hat{u3} [all …]
|
| H A D | asiq.f | 1 subroutine AsIq (y, x, shp, argument 28 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 64 call e3q (yl, dwl, shp, shgl, 88 subroutine AsIqSclr (y, x, shp, argument 96 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 122 call e3qSclr (yl, dwl, shp, shgl,
|
| H A D | errsmooth.f | 2 & shp, shgl, iBC) argument 9 dimension shp(MAXTOP,maxsh,MAXQPT), 53 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 109 subroutine smooth (rerr, x, shp, argument 135 & shp(nshl,maxsh), 181 call getshp(shp, shgl, sgn,
|
| H A D | asigmr.f | 2 & shp, shgl, ien, argument 24 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 86 call e3 (yl, acl, dwl, shp, 126 & shp, shgl, ien, argument 142 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 180 call e3Sclr (yl, acl, shp,
|
| H A D | e3.f | 1 subroutine e3 (yl, acl, dwl, shp, argument 42 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 75 call e3ql (yl, dwl, shp, shgl, 89 call getshp(shp, shgl, sgn, 163 subroutine e3Sclr (yl, acl, shp, 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 D | asiqGradV.f | 1 subroutine AsIqGradV (y, x, shp, argument 26 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 57 call e3qGradV (yl, shp, shgl,
|
| H A D | genlmass.f | 1 subroutine genlmass (x, shp,shgl) argument 10 real*8 shp(MAXTOP,maxsh,MAXQPT), local 42 tmpshp(1:nshl,:) = shp(lcsyst,1:nshl,:) 60 subroutine AsImass (x, shp, argument 75 & shp(nshl,maxsh), shgl(nsd,nshl,ngauss), local 114 call getshp(shp, shgl, sgn,
|
| H A D | hessian.f | 2 & shp, shgl, iBC, argument 13 dimension shp(MAXTOP,maxsh,MAXQPT), 42 & shp(lcsyst,1:nshl,:), 72 & shp(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 )
|