| /phasta/phSolver/compressible/ |
| H A D | e3ls.f | 60 dimension A1(npro,nflow,nflow), A2(npro,nflow,nflow), 61 & A3(npro,nflow,nflow), cv(npro), 62 & A0(npro,nflow,nflow), rho(npro), 64 & dui(npro,nflow), aci(npro,nflow), 67 & rLyi(npro,nflow), dxidx(npro,nsd,nsd), 70 & ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1)), 71 & divqi(npro,nflow-1), stiff(npro,3*nflow,3*nflow), 79 dimension rLymi(npro,nflow), Atau(npro,nflow,nflow), 80 & A1tauA0(npro,nflow,nflow), A2tauA0(npro,nflow,nflow), 81 & A3tauA0(npro,nflow,nflow), fact(npro), [all …]
|
| H A D | au2mfg.f | 34 dimension y(nshg,ndof), ypre(nshg,nflow), 36 & rmes(nshg,nflow), ytmp(nshg,nflow), 37 & res(nshg,nflow), Dy(nshg,nflow), 38 & uBrg(nshg,nflow), BDiag(nshg,nflow,nflow), 40 & iper(nshg), Dy2(nshg,nflow) 42 dimension uBtmp1(nshg,nflow), uBtmp2(nshg,nflow), 56 call sumgat (Dy2, nflow, summed, ilwork)
|
| H A D | asimfg.f | 22 & mater(npro), res(nshg,nflow), 23 & rmes(nshg,nflow), BDiag(nshg,nflow,nflow), 29 & rl(npro,nshl,nflow), rml(npro,nshl,nflow), 30 & BDiagl(npro,nshl,nflow,nflow), 87 call local (res, rl, ien, nflow, 'scatter ') 88 call local (rmes, rml, ien, nflow, 'scatter ') 99 & call local (BDiag, BDiagl, ien, nflow*nflow, 'scatter ')
|
| H A D | itrfdi.f | 39 dimension y(nshg,ndof), ypre(nshg,nflow), 41 & rmes(nshg,nflow), 42 & uBrg(nshg,nflow), BDiag(nshg,nflow,nflow), 47 dimension ytmp(nshg,nflow), rtmp(nshg,nflow) 82 call sumgat (rtmp, nflow, summed, ilwork) 134 call sumgat (rtmp, nflow, summed, ilwork) 191 & rmes(nshg,nflow), 192 & uBrg(nshg,nflow), BDiag(nshg,nflow,nflow), 197 dimension ytmp(nshg,ndof), rtmp(nshg,nflow) 219 call sumgat (rtmp, nflow, summed, ilwork) [all …]
|
| H A D | asigmr.f | 24 & mater(npro), res(nshg,nflow), 25 & rmes(nshg,nflow), BDiag(nshg,nflow,nflow), 30 & xl(npro,nenl,nsd), ytargetl(npro,nshl,nflow), 31 & rl(npro,nshl,nflow), rml(npro,nshl,nflow), 32 & BDiagl(npro,nshl,nflow,nflow), 58 & call localy (ytarget, ytargetl, ien, nflow, 'gather ') 84 call local (res, rl, ien, nflow, 'scatter ') 94 do j = 1, nflow 95 i0 = (i - 1) * nflow + j 96 do k = 1, nflow [all …]
|
| H A D | asires.f | 16 dimension y(nshg,nflow), yc(nshg,ndofl), 17 & x(numnp,nsd), ac(nshg,nflow), 21 & rmes(nshg,nflow) 23 dimension yl(npro,nshl,nflow), ycl(npro,nshl,ndofl), 24 & xl(npro,nenl,nsd), acl(npro,nshl,nflow), 25 & rml(npro,nshl,nflow), ql(npro,nshl,(nflow-1)*nsd) 53 call localy(y, yl, ien, nflow, 'gather ') 57 call localy(ac, acl, ien, nflow, 'gather ') 89 call local (rmes, rml, ien, nflow, 'scatter ')
|
| H A D | itrpr1.f | 23 dimension Binv(npro,nedof,nedof), uBrg(nshg,nflow), 24 & uBrgl(npro,nshl*nflow), ien(npro,nshl), 25 & uBtmp(nshg,nflow) 35 call localt (uBrg, uBrgl, abs(ien), nflow, 'gather ' ) 43 call localt (uBrg, uBrgl, abs(ien), nflow, 'globaliz') 55 call localt (uBrg, uBrgl, abs(ien), nflow, 'gather ') 63 call localt (uBrg, uBrgl, abs(ien), nflow, 'globaliz')
|
| H A D | au1mfg.f | 39 dimension y(nshg,ndof), ypre(nshg,nflow), 41 & rmes(nshg,nflow), ytmp(nshg,nflow), 42 & res(nshg,nflow), uBrg(nshg,nflow), 43 & BDiag(nshg,nflow,nflow), iBC(nshg), 46 dimension uBtmp(nshg,nflow), tmpBC(nshg),
|
| H A D | e3wmlt.f | 55 & WdetJ(npro), ri(npro,nflow*(nsd+1)), 56 & rmi(npro,nflow*(nsd+1)), stiff(npro,3*nflow,3*nflow), 57 & rl(npro,nshl,nflow), rml(npro,nshl,nflow), 61 & shg3(npro), stif1(npro,nflow,nflow), 62 & stif2(npro,nflow,nflow), stif3(npro,nflow,nflow) 155 j0 = nflow * (j - 1) 165 do jdof = 1, nflow 166 do idof = 1, nflow 167 idof2 = idof + nflow 168 jdof2 = jdof + nflow [all …]
|
| H A D | solmfg.f | 45 & res(nshg,nflow), 46 & BDiag(nshg,nflow,nflow), 52 dimension Dy(nshg,nflow), rmes(nshg,nflow), 53 & ypre(nshg,nflow), temp(nshg,nflow), 54 & uBrg(nshg,nflow,Kspace+1), ytmp(nshg,nflow) 68 if(idiff >= 1) idflx= idflx + (nflow-1) * nsd 125 ypre(:,:) = y(:,1:nflow) ! ypre is the pre-conditioned, 143 call sumgat (temp, nflow, summed, ilwork) 252 call sumgat (temp, nflow, beta, ilwork) ! sum vector=(u_{i+1},u_1) 261 call sumgat (temp, nflow, beta, ilwork) ! sum vector=(u_{i+1},u_j)
|
| H A D | solgmr.f | 46 & res(nshg,nflow), 47 & BDiag(nshg,nflow,nflow), 54 dimension Dy(nshg,nflow), rmes(nshg,nflow), 55 & temp(nshg,nflow), 56 & uBrg(nshg,nflow,Kspace+1) 71 if(idiff >= 1) idflx= idflx + (nflow-1) * nsd 126 call sumgat (temp, nflow, summed, ilwork) 166 call bc3per (iBC, temp, iper, ilwork, nflow) 180 call sumgat (temp, nflow, summed, ilwork) 215 call bc3per (iBC, uBrg(:,:,iKs+1), iper, ilwork, nflow) [all …]
|
| H A D | asiflx.f | 20 & mater(npro), flxres(numnp,nflow) 24 & rl(npro,nshl,nflow), BDiagl(npro,nshl,nflow,nflow) 26 dimension ql(npro,nshl,(nflow-1)*nsd) 30 dimension EGmass(npro,nedof,nedof), rml(npro,nshl,nflow) 70 call local (flxres, rl, ien, nflow, 'scatter ')
|
| H A D | shuffle.f | 11 dimension restmp (nshg,nflow), tmp(nshg), tmp1(nshg,nflow) 42 dimension bdiagtmp (nshg,nflow,nflow), tmp(nshg,nflow)
|
| H A D | asbres.f | 15 dimension y(nshg,nflow), x(numnp,nsd), 20 & rmes(nshg,nflow) 22 dimension yl(npro,nshl,nflow), xlb(npro,nenl,nsd), 23 & ycl(npro,nshl,ndof), rml(npro,nshl,nflow) 36 call localy(y, yl, ienb, nflow, 'gather ') 51 call local(rmes, rml, ienb, nflow, 'scatter ')
|
| H A D | e3.f | 50 dimension yl(npro,nshl,nflow), ycl(npro,nshl,ndof), 55 & rl(npro,nshl,nflow), ql(npro,nshl,idflx), 56 & rml(npro,nshl,nflow), xmudmi(npro,ngauss), 57 & BDiagl(npro,nshl,nflow,nflow), 59 & ytargetl(npro,nshl,nflow) 63 dimension g1yi(npro,nflow), g2yi(npro,nflow), 64 & g3yi(npro,nflow), shg(npro,nshl,nsd), 65 & divqi(npro,nflow), tau(npro,5) 80 dimension A0(npro,nflow,nflow), A1(npro,nflow,nflow), 81 & A2(npro,nflow,nflow), A3(npro,nflow,nflow) [all …]
|
| H A D | e3massr.f | 27 dimension aci(npro,nflow), dui(npro,nflow), 28 & ri(npro,nflow*(nsd+1)), 29 & rmi(npro,nflow*(nsd+1)), A0(npro,nflow,nflow)
|
| H A D | e3conv.f | 49 dimension g1yi(npro,nflow), g2yi(npro,nflow), 50 & g3yi(npro,nflow), 51 & A1(npro,nflow,nflow), 52 & A2(npro,nflow,nflow), A3(npro,nflow,nflow), 57 & rLyi(npro,nflow), ri(npro,nflow*(nsd+1)), 58 & rmi(npro,nflow*(nsd+1)), EGmass(npro,nedof,nedof), 64 dimension AiNbi(npro,nflow,nflow), fact1(npro), 193 j0 = nflow * (j - 1) 312 i0 = nflow * (i - 1) 316 do jdof = 1, nflow
|
| H A D | asaugmr.f | 20 & uBrg(nshg,nflow), uBtmp(nshg,nflow) 26 call localt (uBrg, uBrgl, abs(ien), nflow, 'gather ') 32 do i = 1, nflow*nshl, nflow 33 do j = 1, nflow*nshl, nflow 74 call localt (uBtmp, ubBgl, abs(ien), nflow, 'scatter ')
|
| H A D | i3pre.f | 24 dimension BDtmp(nshg,nflow,nflow), 29 dimension BDiagl(numel,nshape,nflow,nflow), 30 & BDiag(nshg,nflow,nflow) 35 call commu (BDiag , ilwork, nflow*nflow , 'out') 54 & nflow*nflow, 'gather ' ) 59 i = (inode - 1) * nflow ! EGmass dof offset 91 i = (inode - 1) * nflow ! EGmass dof offset
|
| H A D | e3dc.f | 34 dimension g1yi(npro,nflow), g2yi(npro,nflow), 35 & g3yi(npro,nflow), A0(npro,5,5), 38 & ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1)), 39 & stiff(npro,3*nflow,3*nflow),dtmp(npro) 42 dimension ggyi(npro,nflow), gAgyi(npro,15), 411 nflow2=two*nflow 412 do j = 1, nflow 413 do i = 1, nflow 424 stiff(:,i,j+nflow) = stiff(:,i,j+nflow) 436 stiff(:,i+nflow,j) = stiff(:,i+nflow,j) [all …]
|
| H A D | itrPC.f | 132 real*8 y(nshg,ndof), ac(nshg,ndof), Dy(nshg,nflow) 197 acold(:,1:nflow) = acold(:,1:nflow) 198 & + (ac(:,1:nflow)-acold(:,1:nflow))*fct2 199 yold(:,1:nflow) = yold(:,1:nflow) 200 & + (y(:,1:nflow)-yold(:,1:nflow))*fct3
|
| H A D | sparseap.f | 20 real*8 lhsK(nflow*nflow,nnz_tot) 21 real*8 p(nshg,nflow), q(nshg,nflow) 28 call commu (p, ilwork, nflow , 'out') 33 do j=1,nflow 110 call commu (p , ilwork, nflow , 'in ')
|
| H A D | asbmfg.f | 20 & res(nshg,nflow), rmes(nshg,nflow) 23 & rl(npro,nshl,nflow), 24 & rml(npro,nshl,nflow), 58 call local(res, rl, ienb, nflow, 'scatter ') 60 & call local(rmes, rml, ienb, nflow, 'scatter ')
|
| /phasta/phSolver/common/ |
| H A D | getvel.f | 11 & velo(numnp,nflow), velf(nfath,nflow), 12 & velft(nfath,nflow), 13 & velbar(nfath,nflow), 25 if(nflow.eq.5) velo(:,5)=y(:,5) 39 if(nflow.eq.5) then 78 velf(ifathi,1:nflow) = velf(ifathi,1:nflow) 79 & + velo(i,1:nflow) 87 call drvAllreduce(velf, velft,nfath*nflow) 104 if(nflow.eq.5) velft(:,5) = velft(:,5) * rinvsons(:)
|
| H A D | genscale.f | 25 & x(numnp,nsd), velbarR(nfint,nflow) 26 dimension ifath(numnp), velbarl(nelint,nshl,nflow), 29 & xsi(nsd), yintl(nelint,nshl,nflow), 30 & flucl(nelint,nshl,nflow), 31 & ubarintl(nelint,nshl,nflow), 32 & fluc1(npin,nflow), fluc2(npin,nflow), 33 & ubar1(npin,nflow), ubar2(npin,nflow) 37 real*8 freestream(nflow) 60 if(nflow.gt.4) then 61 do j = 5, nflow [all …]
|