| /phasta/phSolver/compressible/ |
| H A D | au1gmr.f | 1 subroutine Au1GMR (EGmass, uBrg, ilwork,iBC,iper ) argument
22 dimension EGmass(numel,nedof,nedof), uBrg(nshg,nflow),
30 call commu (uBrg, ilwork, nflow , 'out')
36 uBrg(:,j)=uBrg(iper(:),j)
42 & call rotabc(uBrg(1,2), iBC, 'out')
60 call asAuGMR (mien(iblk)%p, EGmass(iel:inum,:,:), uBrg,
65 uBrg = uBtmp
71 & call rotabc(uBrg(1,2), iBC, 'in ')
77 call commu (uBrg , ilwork, nflow , 'in ')
94 uBrg(isgbeg:isgend,:) = zero
[all …]
|
| H A D | au2mfg.f | 3 & uBrg, BDiag, iBC, BC,
38 & uBrg(nshg,nflow), BDiag(nshg,nflow,nflow),
65 uBrg = ypre + eps * Dy
67 call i3LU (BDiag, uBrg, 'backward')
69 call yshuffle(uBrg, 'old2new ')
71 call itrBC (uBrg, uBrg, iBC, BC, iper, ilwork)
73 call itrRes (uBrg, y,
89 uBrg = ypre - eps * Dy
91 call i3LU (BDiag, uBrg, 'backward')
93 call yshuffle(uBrg, 'old2new ')
[all …]
|
| H A D | au1mfg.f | 2 & rmes, res, uBrg,
42 & res(nshg,nflow), uBrg(nshg,nflow),
64 uBrg = ypre + eGMRES * uBrg
66 call i3LU (BDiag, uBrg, 'backward')
68 call yshuffle(uBrg, 'old2new ')
70 call itrBC (uBrg, uBrg, iBC, BC, iper, ilwork)
72 call itrRes (uBrg, y,
83 uBrg = ( uBtmp - rmes ) / eGMRES
|
| H A D | solmfg.f | 54 & uBrg(nshg,nflow,Kspace+1), ytmp(nshg,nflow)
117 uBrg(:,:,1) = res
195 uBrg(:,:,1) = temp
217 uBrg(:,:,1) = uBrg(:,:,1) / unorm
226 temp = uBrg(:,:,iKs)
236 uBrg(:,:,iKs+1) = temp ! u_{i+1}= J u_i In Johan Thesis p 15c
251 temp = uBrg(:,:,iKs+1) * uBrg(:,:,1) ! {u_{i+1}*u_1} vector
258 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) - beta * uBrg(:,:,jK-1)
260 temp = uBrg(:,:,iKs+1) * uBrg(:,:,jK) !{u_{i+1}*u_j} vector
279 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) / unorm ! normalize the next Krylov
[all …]
|
| H A D | solgmr.f | 56 & uBrg(nshg,nflow,Kspace+1)
120 uBrg(:,:,1) = res
175 uBrg(:,:,1) = temp
196 uBrg(:,:,1) = uBrg(:,:,1) / unorm
203 uBrg(:,:,iKs+1) = uBrg(:,:,iKs)
211 call Au1GMR ( EGmass, uBrg(:,:,iKs+1), ilwork, iBC,iper)
215 call bc3per (iBC, uBrg(:,:,iKs+1), iper, ilwork, nflow)
228 temp = uBrg(:,:,iKs+1) * uBrg(:,:,1) ! {u_{i+1}*u_1} vector
235 uBrg(:,:,iKs+1) = uBrg(:,:,iKs+1) - beta * uBrg(:,:,jK-1)
237 temp = uBrg(:,:,iKs+1) * uBrg(:,:,jK) !{u_{i+1}*u_j} vector
[all …]
|
| H A D | itrpr1.f | 1 subroutine itrPr1 (ien, Binv, uBrg, uBtmp, code) argument
23 dimension Binv(npro,nedof,nedof), uBrg(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 | i3pcond.f | 1 subroutine i3Pcond ( Binv, uBrg, ilwork, code ) argument
24 dimension Binv(numel,nedof,nedof), uBrg(nshg,nflow)
47 call itrPr1 (mien(iblk)%p, Binv(iel:inum,:,:), uBrg,
62 call itrPr1 (mien(iblk)%p, Binv(iel:inum,:,:), uBrg,
|
| H A D | itrfdi.f | 2 & rmes, uBrg, BDiag,
42 & uBrg(nshg,nflow), BDiag(nshg,nflow,nflow),
97 ytmp = ypre + epsSD * uBrg
113 ytmp = ypre - epsSD * uBrg
152 & rmes, uBrg, BDiag,
192 & uBrg(nshg,nflow), BDiag(nshg,nflow,nflow),
232 ytmp = ypre + epsSD * uBrg
244 ytmp = ypre - epsSD * uBrg
|
| H A D | asaugmr.f | 1 subroutine AsAuGMR (ien, EGmass, uBrg, uBtmp ) argument
20 & uBrg(nshg,nflow), uBtmp(nshg,nflow)
26 call localt (uBrg, uBrgl, abs(ien), nflow, 'gather ')
83 subroutine AsAuGMRSclr (ien, EGmass, uBrg, uBtmp ) argument
104 & uBrg(nshg), uBtmp(nshg)
111 call localtSclr(uBrg, uBrgl, ien, 'gather ')
|