| /phasta/phSolver/compressible/ |
| H A D | bc3lhs.f | 36 do iel = 1, npro 44 in = ien(iel,inod) 58 EGmass(iel,i,i1) = zero 59 EGmass(iel,i1,i) = zero 61 EGmass(iel,i1,i1) = one 75 EGmass(iel,i3,i) = EGmass(iel,i3,i) 76 & - BC(in,4) * EGmass(iel,i2,i) 77 EGmass(iel,i4,i) = EGmass(iel,i4,i) 78 & - BC(in,5) * EGmass(iel,i2,i) 81 EGmass(iel,i,i3) = EGmass(iel,i,i3) [all …]
|
| H A D | i3pcond.f | 37 iel = lcblk(1,iblk) 39 npro = lcblk(1,iblk+1) - iel 40 inum = iel + npro - 1 47 call itrPr1 (mien(iblk)%p, Binv(iel:inum,:,:), uBrg, 62 call itrPr1 (mien(iblk)%p, Binv(iel:inum,:,:), uBrg,
|
| H A D | au1gmr.f | 52 iel = lcblk(1,iblk) 54 npro = lcblk(1,iblk+1) - iel 55 inum = iel + npro - 1 60 call asAuGMR (mien(iblk)%p, EGmass(iel:inum,:,:), uBrg, 152 iel = lcblk(1,iblk) 154 npro = lcblk(1,iblk+1) - iel 155 inum = iel + npro - 1 160 call asAuGMRSclr (mien(iblk)%p, EGmassT(iel:inum,:,:), uBrg,
|
| H A D | elmgmr.f | 72 iel = lcblk(1,iblk) 81 npro = lcblk(1,iblk+1) - iel 132 iel = lcblk(1,iblk) 141 npro = lcblk(1,iblk+1) - iel 142 inum = iel + npro - 1 160 & qres, EGmass(iel:inum,:,:), 166 & EGmass(iel:inum,:,:) ) 183 iel = lcblkb(1,iblk) 193 npro = lcblkb(1,iblk+1) - iel 355 iel = lcblk(1,iblk) [all …]
|
| H A D | e3b.f | 326 do iel = 1, npro 331 iface = abs(iBCB(iel,2)) 333 flxID(1,iface) = flxID(1,iface) + WdetJb(iel)! measure area too 335 flxID(2,iface) = flxID(2,iface) - WdetJb(iel) * rou(iel) 337 & - ( tau1n(iel) - bnorm(iel,1)*pres(iel)) 338 & * WdetJb(iel) 340 & - ( tau2n(iel) - bnorm(iel,2)*pres(iel)) 341 & * WdetJb(iel) 343 & - ( tau3n(iel) - bnorm(iel,3)*pres(iel)) 344 & * WdetJb(iel)
|
| H A D | e3ql.f | 182 do iel=1,npro 186 rminv(iel,i,j) = 0.0 188 rminv(iel,i,i)=1.0 193 call ludcmp(rmass(iel,:,:),nshl,nshl,indx,d) 199 call lubksb(rmass(iel,:,:),nshl,nshl, 200 & indx,rminv(iel,:,j)) 208 do iel=1,npro 210 ql(iel,:,j) = matmul( rminv(iel,:,:),qrl(iel,:,j) )
|
| H A D | i3pre.f | 44 iel = lcblk(1,iblk) 46 npro = lcblk(1,iblk+1) - iel 47 n = iel - 1 48 inum = iel + npro - 1 53 call local (BDiag, BDiagl(iel:inum,:,:,:), abs(mien(iblk)%p), 186 iel = lcblk(1,iblk) 188 npro = lcblk(1,iblk+1) - iel 189 n = iel - 1 190 inum = iel + npro - 1 195 call local (Dinv, Dinvl(iel:inum,:), mien(iblk)%p,
|
| H A D | elmgmrpetsc.f | 86 iel = lcblk(1,iblk) 95 npro = lcblk(1,iblk+1) - iel 147 iel = lcblk(1,iblk) 156 npro = lcblk(1,iblk+1) - iel 157 inum = iel + npro - 1 219 iel = lcblkb(1,iblk) 229 npro = lcblkb(1,iblk+1) - iel 375 iel = lcblk(1,iblk) 383 npro = lcblk(1,iblk+1) - iel 384 inum = iel + npro - 1 [all …]
|
| H A D | elmmfg.f | 71 iel = lcblk(1,iblk) 80 npro = lcblk(1,iblk+1) - iel 130 iel = lcblk(1,iblk) 139 npro = lcblk(1,iblk+1) - iel 140 inum = iel + npro - 1 176 iel = lcblkb(1,iblk) 186 npro = lcblkb(1,iblk+1) - iel
|
| /phasta/phSolver/incompressible/ |
| H A D | bc3lhs.f | 38 do iel = 1, npro 46 in = abs(ien(iel,inod)) 84 xKebe(iel,iadj1,i,inod) = xKebe(iel,iadj1,i,inod) 85 & - BC(in,4) * xKebe(iel,irem1,i,inod) 86 xKebe(iel,iadj2,i,inod) = xKebe(iel,iadj2,i,inod) 87 & - BC(in,4) * xKebe(iel,irem2,i,inod) 88 xKebe(iel,iadj3,i,inod) = xKebe(iel,iadj3,i,inod) 89 & - BC(in,4) * xKebe(iel,irem3,i,inod) 104 xKebe(iel,iadj1,i,inod) = xKebe(iel,iadj1,i,inod) 105 & - BC(in,5) * xKebe(iel,irem1,i,inod) [all …]
|
| H A D | e3b.f | 139 do iel = 1, npro 144 iface = abs(iBCB(iel,2)) 154 flxID(1,iface) = flxID(1,iface) + WdetJb(iel)! measure area too 155 flxID(2,iface) = flxID(2,iface) - WdetJb(iel) * unm(iel) 157 & - ( tau1n(iel) - bnorm(iel,1)*pres(iel)) 158 & * WdetJb(iel) 160 & - ( tau2n(iel) - bnorm(iel,2)*pres(iel)) 161 & * WdetJb(iel) 163 & - ( tau3n(iel) - bnorm(iel,3)*pres(iel)) 164 & * WdetJb(iel) [all …]
|
| H A D | e3ql.f | 134 do iel=1,npro 137 lmassinv(iblock)%p(iel,i,j) = 0.0 139 lmassinv(iblock)%p(iel,i,i)=1.0 144 call ludcmp(rmass(iel,:,:),nshl,nshl,indx,d) 150 call lubksb(rmass(iel,:,:),nshl,nshl,indx, 151 & lmassinv(iblock)%p(iel,:,j)) 162 do iel=1,npro 165 ql(iel,:,j) = matmul( rminv(iel,:,:),qrl(iel,:,j) ) 281 do iel=1,npro 284 lmassinv(iblock)%p(iel,i,j) = 0.0 [all …]
|
| H A D | bflux.f | 69 iel = lcblkb(1,iblk) 76 npro = lcblkb(1,iblk+1) - iel 108 iel = lcblk(1,iblk) 113 npro = lcblk(1,iblk+1) - iel 165 iel = lcblkb(1,iblk) 171 npro = lcblkb(1,iblk+1) - iel 452 integer nshl, nshg, npro, npro2, i, iel local 459 do iel = 1, npro 460 map(iel) = map(iel) + mask( abs(ien(iel,i)) ) 465 do iel = 1, npro [all …]
|
| H A D | e3source.f | 24 do iel = 1, npro 26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2) 31 src(iel,1) = -xx(iel,2)/r * Stheta 32 src(iel,2) = xx(iel,1)/r * Stheta 33 src(iel,3) = dpdz 169 do iel = 1, npro 170 x = xx(iel,1) 171 y = xx(iel,2) 172 z = xx(iel,3) 194 src(iel,1) = 0.0 [all …]
|
| H A D | advLES.f | 568 iel = lcblk(1,iblk) !Element number where this block begins 569 npro = lcblk(1,iblk+1) - iel 573 inum = iel + npro - 1 676 iel = lcblk(1,iblk) !Element number where this block begins 677 npro = lcblk(1,iblk+1) - iel 681 inum = iel + npro - 1 686 & strl(iel:inum,:), shgl(lcsyst,:,1:nshl,:), 818 iel = lcblk(1,iblk) 819 npro = lcblk(1,iblk+1) - iel 821 inum = iel + npro - 1 [all …]
|
| H A D | elmgmr.f | 85 iel = lcblk(1,iblk) 94 npro = lcblk(1,iblk+1) - iel 143 iel = lcblk(1,iblk) 152 npro = lcblk(1,iblk+1) - iel 153 inum = iel + npro - 1 254 iel = lcblkb(1,iblk) 264 npro = lcblkb(1,iblk+1) - iel 424 iel = lcblk(1,iblk) 430 npro = lcblk(1,iblk+1) - iel 469 iel = lcblk(1,iblk) [all …]
|
| /phasta/phSolver/common/ |
| H A D | cmass.f | 32 integer ne, na, nb, nodlcla, nodlclb, iel local 58 do iel = 1, npro 61 shp2 = shape(iel,na) * shape(iel,nb) 62 em(iel,na,nb) = em(iel,na,nb) + shp2*WdetJ(iel) 85 integer ne, na, nb, nodlcla, nodlclb, iel local 158 do iel = 1, npro 161 shp2 = shape(iel,na) * shape(iel,nb) 162 em(iel,na,nb) = em(iel,na,nb) + shp2*WdetJ(iel) 182 do iel = 1, npro 184 em(iel,na,na) = em(iel,na,na)+flump*eml(iel,na) [all …]
|
| H A D | aveprep.f | 43 iel = lcblk(1,iblk) !Element number where this block begins 44 npro = lcblk(1,iblk+1) - iel 48 inum = iel + npro - 1 52 call getylist( ylist, ifathe(iel:inum,:), shp(lcsyst,1:nshl,:), 84 do iel = 1, npro 86 xmudmi(iel,intp) = cdelsq( lfathe(iel,intp) )* 87 & strl(iel,intp) 319 iel = lcblk(1,iblk) !Element number where this block begins 320 npro = lcblk(1,iblk+1) - iel 324 inum = iel + npro - 1 [all …]
|
| H A D | asbnabi.f | 105 do iel=1,npro 106 if (.not.btest(iBCB(iel,1),1)) then 107 Wfactor(iel) = zero ! we want zeros where we are not integrating 146 do iel=1,npro 147 if (iBCB(iel,2) .ne. 0) then 148 iface = iBCB(iel,2) 149 ndsurf(ienb(iel,1:nshlb))=iface 261 do iel=1,npro 262 if (.not.btest(iBCB(iel,1),1)) then 263 Wfactor(iel) = zero ! we want zeros where we are not integrating [all …]
|
| H A D | getdmc.f | 73 iel = lcblk(1,iblk) !Element number where this block begins 74 npro = lcblk(1,iblk+1) - iel 78 inum = iel + npro - 1 83 call asithf (y, x, strl(iel:inum,:), mien(iblk)%p, fres, 92 iel = lcblk(1,iblk) !Element number where this block begins 93 npro = lcblk(1,iblk+1) - iel 97 inum = iel + npro - 1 104 & strl(iel:inum,:), shgl(lcsyst,:,1:nshl,:), 386 iel = lcblk(1,iblk) 387 npro = lcblk(1,iblk+1) - iel [all …]
|
| H A D | ibcupdate.f | 14 do iel = 1, npro 15 if(btest(iBCB(iel,1),4)) then ! turbulence wall 17 iglobal = ienb(iel,inode)
|
| H A D | genblkPosix.f | 33 iel=1 61 lcblk(1,nelblk) = iel 95 iel=iel+npro 99 lcblk(1,nelblk+1) = iel
|
| H A D | genblkSyncIO.f | 33 iel=1 105 lcblk(1,nelblk) = iel 138 iel=iel+npro 144 lcblk(1,nelblk+1) = iel
|
| H A D | genblk.f | 33 iel=1 114 lcblk(1,nelblk) = iel 147 iel=iel+npro 153 lcblk(1,nelblk+1) = iel
|
| H A D | genbkbPosix.f | 40 iel=1 113 lcblkb(1,nelblb) = iel 146 iel=iel+npro 152 lcblkb(1,nelblb+1) = iel
|