| /phasta/phSolver/common/ |
| H A D | cmass.f | 34 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT), 36 & xl(npro,nenl,nsd) 38 dimension shape(npro,nshl), shdrv(npro,nsd,nshl), 39 & sgn(npro,nshl), dxidx(npro,nsd,nsd), 40 & shg(npro,nshl,nsd), 87 & shp(nshl,MAXQPT), shgl(nsd,nshl,MAXQPT), 88 & shpf(nshl,MAXQPT), shglf(nsd,nshl,MAXQPT), 90 & xl(npro,nenl,nsd) 92 dimension shape(npro,nshl), shdrv(npro,nsd,nshl), 93 & sgn(npro,nshl), dxidx(npro,nsd,nsd), [all …]
|
| H A D | asbnabi.f | 15 dimension xlb(npro,nenl,nsd), bnorm(npro,nsd), 16 & rl(npro,nshl,nsd), WdetJb(npro), 19 dimension x(numnp,nsd), 20 & shpb(nshl,ngaussb), shglb(nsd,nshl,ngaussb), 25 & shpfun(npro,nshl), shdrv(npro,nsd,nshl) 28 dimension dxdxib(npro,nsd,nsd), temp(npro), 31 & v1(npro,nsd), v2(npro,nsd) 42 call localx(x, xlb, ienb, nsd, 'gather ') 173 dimension xlb(npro,nenl,nsd), 177 dimension x(numnp,nsd), [all …]
|
| H A D | mod_wallData.f90 | 29 real*8 x(numnp,nsd) 34 real*8 eNorm(nsd),en(nsd),wn(nsd),aedg(nsd),bedg(nsd) 48 allocate(wNorm(nshg,nsd)) 86 allocate(wNormTmp(nshg,nsd)) 163 real*8 :: wNorm(nshg, nsd) 164 real*8 :: eNorm(nsd), en(nsd), wn(nsd)
|
| H A D | errorsmooth.f | 10 & shgl(MAXTOP,nsd,maxsh,MAXQPT), 12 & shglb(MAXTOP,nsd,maxsh,MAXQPT) 50 allocate (tmpshgl(nsd,nshl,MAXQPT)) 133 dimension rerr(nshg,10), x(numnp,nsd), 135 & shgl(nsd,nshl,maxsh), 141 dimension rerrl(npro,nshl,10), xl(npro,nenl,nsd), 145 & shdrv(npro,nsd,nshl), WdetJ(npro), 146 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd) 161 call localx(x, xl, ien, nsd, 'gather ')
|
| H A D | hierarchic.f | 37 dimension shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 39 & shdrv(npro,nsd,nshl) 67 dimension shp(nshl,ngaussb), shgl(nsd,nshl,ngaussb), 69 & shdrv(npro,nsd,nshl) 244 & shp(nshl,npts), shgl(nsd,nshl,npts), 245 & intpnt(3,npts), x(numnp,nsd) 274 allocate ( xl(npro,nenl,nsd ) ) 280 call localx( x, xl, mien(iblk)%p, nsd, 'gather ') 327 & shdrv(npro,nsd,nshl), shp(nshl,npts), 328 & shgl(nsd,nshl,npts), xl(npro,nenl,nsd), [all …]
|
| H A D | qpbc.f | 20 do i=1,idflx/nsd 45 idflow = (nflow-1)*nsd 72 idflow=(nflow-1)*nsd 84 do i=1, nsd 97 do i=1,idflx/nsd 131 dimension rmass(nshg), qres(nshg,nsd), 141 call commu (qres , ilwork, nsd , 'in ') 168 do i=1, nsd 173 call commu (qres, ilwork, nsd, 'out')
|
| H A D | findslpw.f | 6 real*8 x(numnp,nsd) 11 real*8 elnrm(nsd),wn1(nsd),wn2(nsd),wn3(nsd),aedg(nsd),bedg(nsd) 15 real*8 AA(nsd),BB(nsd),tmp 20 allocate(wlnorm(nshg,nsd,9))
|
| H A D | e3metric.f | 13 real*8 xl(npro,nenl,nsd), shgl(npro,nsd,nshl), 14 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd), 17 real*8 dxdxi(npro,nsd,nsd), tmp(npro)
|
| /phasta/phSolver/incompressible/ |
| H A D | e3bvar.f | 52 & shpb(npro,nshl), shglb(npro,nsd,nshl), 53 & xlb(npro,nenl,nsd), 56 & WdetJb(npro), bnorm(npro,nsd), 63 & acl(npro,nshl,ndof), ul(npro,nshl,nsd), 64 & vdot(npro,nsd), rlKwall(npro,nshlb,nsd) 67 & gl3yi(npro,ndof), dxdxib(npro,nsd,nsd), 68 & dxidxb(npro,nsd,nsd), temp(npro), 71 & v1(npro,nsd), v2(npro,nsd), 72 & v3(npro,nsd), 73 & rotnodallocal(npro,nsd,nsd), [all …]
|
| H A D | asbflx.f | 18 dimension y(nshg,ndofl), x(numnp,nsd), 19 & ac(nshg,ndofl), u(nshg,nsd), 21 & shglb(nsd,nshl,ngaussb), 25 & flxLHS(nshg,1), flxnrm(nshg,nsd) 27 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 29 & flhsl(npro,nshl,1), fnrml(npro,nshl,nsd), 31 & ul(npro,nshl,nsd), acl(npro,nshl,ndofl) 51 call localx(x, xlb, ienb, nsd, 'gather ') 52 call localx(u, ul, ienb, nsd, 'gather ') 87 do i = 1, nsd [all …]
|
| H A D | asiq.f | 27 dimension y(nshg,ndof), x(numnp,nsd), 28 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 32 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 54 call localx (x, xl, ien, nsd, 'gather ') 95 dimension y(nshg,ndof), x(numnp,nsd), 96 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 98 & qres(nshg,nsd), rmass(nshg) 100 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 101 & ql(npro,nshl,nsd), rmassl(npro,nshl) 112 call localx (x, xl, ien, nsd, 'gather ') [all …]
|
| H A D | e3ql.f | 26 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 27 & xl(npro,nenl,nsd), sgn(npro,nshl), 33 & g3yi(npro,ndof), shg(npro,nshl,nsd), 34 & dxidx(npro,nsd,nsd), WdetJ(npro), 37 & qrl(npro,nshl,nsd*nsd) 39 dimension qdi(npro,nsd*nsd), shape(npro,nshl), 40 & shdrv(npro,nsd,nshl), indx(nshl), 192 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 193 & xl(npro,nenl,nsd), sgn(npro,nshl), 194 & ql(npro,nshl,nsd) [all …]
|
| H A D | e3.f | 42 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 43 & xl(npro,nenl,nsd), dwl(npro,nenl), 52 & g3yi(npro,ndof), shg(npro,nshl,nsd), 53 & aci(npro,3), dxidx(npro,nsd,nsd), 57 & rLui(npro,nsd), uBar(npro,nsd), 59 & shpfun(npro,nshl), shdrv(npro,nsd,nshl), 180 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 181 & xl(npro,nenl,nsd), rl(npro,nshl), 182 & ql(npro,nshl,nsd), xSebe(npro,nshl,nshl), 187 real*8 gradS(npro,nsd), shg(npro,nshl,nsd), [all …]
|
| H A D | asbmfg.f | 19 dimension y(nshg,ndofl), x(numnp,nsd), 20 & ac(nshg,ndofl), u(nshg,nsd), 22 & shglb(nsd,nshl,ngaussb), 27 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 29 & ul(npro,nshl,nsd), acl(npro,nshl,ndofl) 44 call localx(x, xlb, ienb, nsd, 'gather ') 45 call localx(u, ul, ienb, nsd, 'gather ') 93 dimension y(nshg,ndofl), x(numnp,nsd), 95 & shglb(nsd,nshl,*), 100 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd), [all …]
|
| H A D | e3qvar.f | 38 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd), 40 & g3yi(npro,nflow), shg(npro,nshl,nsd), 41 & dxidx(npro,nsd,nsd), WdetJ(npro) 45 dimension tmp(npro), dxdxi(npro,nsd,nsd) 156 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd), 157 & dxidx(npro,nsd,nsd), WdetJ(npro), 158 & gradT(npro,nsd) 162 real*8 shg(npro,nshl,nsd)
|
| H A D | errsmooth.f | 10 & shgl(MAXTOP,nsd,maxsh,MAXQPT), 12 & shglb(MAXTOP,nsd,maxsh,MAXQPT), 13 & x(numnp,nsd) 51 allocate (tmpshgl(nsd,nshl,MAXQPT)) 134 dimension rerr(nshg,10), x(numnp,nsd), 136 & shgl(nsd,nshl,maxsh), 142 dimension rerrl(npro,nshl,10), xl(npro,nenl,nsd), 146 & shdrv(npro,nsd,nshl), WdetJ(npro), 147 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd) 162 call localx(x, xl, ien, nsd, 'gather ')
|
| H A D | hessian.f | 10 & x(numnp,nsd), iBC(nshg), 14 & shgl(MAXTOP,nsd,maxsh,MAXQPT), 16 & shglb(MAXTOP,nsd,maxsh,MAXQPT) 92 dimension y(nshg,ndof), x(numnp,nsd), 93 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 95 & shdrv(npro,nsd,nshl), shape( npro, nshl ), 98 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 99 & ql(npro,nshl,9), dxidx(npro,nsd,nsd), 121 call localx (x, xl, ien, nsd, 'gather ') 167 dimension gradu(nshg,9), x(numnp,nsd), [all …]
|
| H A D | asigmr.f | 23 & x(numnp,nsd), 24 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 31 & xl(npro,nenl,nsd), dwl(npro,nenl), 60 call localx(x, xl, ien, nsd, 'gather ') 141 & x(numnp,nsd), 142 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 144 & res(nshg), qres(nshg,nsd) 148 & xl(npro,nenl,nsd), 149 & rl(npro,nshl), ql(npro,nshl,nsd), 166 call localx(x, xl, ien, nsd, 'gather ') [all …]
|
| H A D | e3q.f | 25 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 26 & xl(npro,nenl,nsd), 33 & g3yi(npro,nflow), shg(npro,nshl,nsd), 34 & dxidx(npro,nsd,nsd), WdetJ(npro), 40 & shdrv(npro,nsd,nshl), shpsum(npro) 220 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 221 & xl(npro,nenl,nsd), 222 & ql(npro,nshl,nsd), rmassl(npro,nshl) 226 dimension gradT(npro,nsd), 227 & dxidx(npro,nsd,nsd), WdetJ(npro) [all …]
|
| H A D | genlmass.f | 8 real*8 x(numnp,nsd) 11 & shgl(MAXTOP,nsd,maxsh,MAXQPT) 41 allocate (tmpshgl(nsd,nshl,MAXQPT)) 74 real*8 x(numnp,nsd), 75 & shp(nshl,maxsh), shgl(nsd,nshl,ngauss), 81 real*8 xl(npro,nenl,nsd), WdetJ(npro), 83 & locmass(npro,nshl), shg(npro,nshl,nsd), 85 & dxidx(npro,nsd,nsd), shdrv(npro,nsd,nshl) 100 call localx(x, xl, ien, nsd, 'gather ')
|
| H A D | e3qGradV.f | 24 & shp(nshl,ngauss), shgl(nsd,nshl,ngauss), 25 & xl(npro,nenl,nsd), 31 & g3yi(npro,nflow), shg(npro,nshl,nsd), 32 & dxidx(npro,nsd,nsd), WdetJ(npro) 36 & shdrv(npro,nsd,nshl)
|
| /phasta/svLS/ |
| H A D | NSSOLVER.f | 56 INTEGER i, j, k, iB, iBB, nB, nsd local 62 nsd = dof - 1 65 ALLOCATE(Rm(nsd,nNo), Rc(nNo), Rmi(nsd,nNo), Rci(nNo), 66 2 U(nsd,nNo,iB), P(nNo,iB), MU(nsd,nNo,nB), MP(nNo,nB), 69 Rmi = Ri(1:nsd,:) 76 eps = SQRT(NORMV(nsd, mynNo, commu, Rm)**2D0 95 CALL GMRES(nFaces, nsd, nNo, nnz, mynNo, commu, cS, face, 98 CALL SPARMULVS(nsd, nNo, nnz, commu, cS, 102 CALL CGRAD(nFaces, nsd, nNo, nnz, mynNo, commu, cS, face, 105 CALL SPARMULSV(nsd, nNo, nnz, commu, cS, [all …]
|
| /phasta/phSolver/compressible/ |
| H A D | asbflx.f | 15 dimension y(nshg,ndof), x(numnp,nsd), 17 & shglb(nsd,nshl,ngaussb), 21 & flxLHS(numnp,1), flxnrm(numnp,nsd) 23 dimension ycl(npro,nshl,ndof), xlb(npro,nenl,nsd), 25 & flhsl(npro,nshl,1), fnrml(npro,nshl,nsd), 47 call localx(x, xlb,ienb, nsd, 'gather ') 86 & min(1, invflx(abs(ienb(:,mnodeb(n,lelCat,nsd))))) 92 do i = 1, nsd 101 call local (flxnrm, fnrml, ienb, nsd, 'scatter ')
|
| H A D | asbmfg.f | 15 dimension y(nshg,ndofl), x(numnp,nsd), 17 & shglb(nsd,nshl,ngaussb), 22 dimension ycl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 40 call localx(x, xlb, ienb, nsd, 'gather ') 82 dimension y(nshg,ndofl), x(numnp,nsd), 84 & shglb(nsd,nshl,maxsh), 89 dimension ycl(npro,nshl,ndofl), xlb(npro,nenl,nsd), 103 call localx (x, xlb, ienb, nsd, 'gather ')
|
| H A D | asiflx.f | 16 & x(numnp,nsd), 18 & shgl(MAXTOP,nsd,maxsh,MAXQPT), 23 & xl(npro,nenl,nsd), 26 dimension ql(npro,nshl,(nflow-1)*nsd) 45 call localx(x, xl, ien, nsd, 'gather ') 62 if (nsd .eq. 3) then
|