| /phasta/phSolver/incompressible/ |
| H A D | e3qGradV.f | 2 & xl, ql, argument 26 & ql(npro,nshl,nsdsq) 63 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*g1yi(:,2 ) ! du/dx 64 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*g2yi(:,2 ) ! du/dy 65 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*g3yi(:,2 ) ! du/dz 67 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*g1yi(:,3 ) ! dv/dx 68 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*g2yi(:,3 ) ! dv/dy 69 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*g3yi(:,3 ) ! dv/dz 71 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*g1yi(:,4 ) ! dw/dx 72 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*g2yi(:,4 ) ! dw/dy [all …]
|
| H A D | e3q.f | 2 & xl, ql, rmassl, argument 27 & ql(npro,nshl,idflx), rmassl(npro,nshl), 112 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*qdi(:,1 ) 113 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*qdi(:,2 ) 114 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*qdi(:,3 ) 116 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*qdi(:,4 ) 117 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*qdi(:,5 ) 118 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*qdi(:,6 ) 120 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*qdi(:,7 ) 121 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*qdi(:,8 ) [all …]
|
| H A D | asiq.f | 33 & ql(npro,nshl,idflx), rmassl(npro,nshl), 61 ql = zero 65 & xl, ql, rmassl, 71 call local (qres, ql, ien, idflx, 'scatter ') 101 & ql(npro,nshl,nsd), rmassl(npro,nshl) 119 ql = zero 123 & xl, ql, rmassl, 129 call local (qres, ql, ien, nsd, 'scatter ')
|
| H A D | asiqGradV.f | 31 & ql(npro,nshl,nsdsq) 55 ql = zero 58 & xl, ql, 64 call local (qres, ql, ien, nsdsq, 'scatter ')
|
| H A D | e3ivar.f | 7 & ql, rLui, src, argument 65 & ql(npro,nshl,idflx), rLui(npro,nsd), 177 divqi(:,1) = divqi(:,1) + shg(:,n,1)*ql(:,n,1 ) 178 & + shg(:,n,2)*ql(:,n,4 ) 179 & + shg(:,n,3)*ql(:,n,7 ) 181 divqi(:,2) = divqi(:,2) + shg(:,n,1)*ql(:,n,2 ) 182 & + shg(:,n,2)*ql(:,n,5 ) 183 & + shg(:,n,3)*ql(:,n,8) 185 divqi(:,3) = divqi(:,3) + shg(:,n,1)*ql(:,n,3 ) 186 & + shg(:,n,2)*ql(:,n,6 ) [all …]
|
| H A D | asigmr.f | 33 & ql(npro,nshl,idflx) 61 call local (qres, ql, ien, idflx, 'gather ') 88 & ql, xKebe, xGoC, xmudmi, 149 & rl(npro,nshl), ql(npro,nshl,nsd), local 169 call local (qres, ql, ien, nsd, 'gather ') 182 & rl, ql, xSebe,
|
| H A D | e3.f | 2 & shgl, xl, rl, ql, argument 44 & rl(npro,nshl,nflow), ql(npro,nshl,idflx) 76 & xl, ql, xmudmi, 104 & ql, rLui, src, 165 & rl, ql, xSebe, argument 182 & ql(npro,nshl,nsd), xSebe(npro,nshl,nshl), local 211 call e3qlSclr (yl, dwl, shp, shgl, xl, ql, sgn) 236 & ql, rLS, SrcR,
|
| H A D | hessian.f | 99 & ql(npro,nshl,9), dxidx(npro,nsd,nsd), 125 ql = zero 141 ql(:,i,j) = ql(:,i,j)+shape(:,i)*WdetJ*gradul(:,j) 151 call local (gradu, ql, ien, 9, 'scatter ') 174 & ql(npro,nshl,27), dxidx(npro,nsd,nsd), 200 ql = zero 216 ql(:,i,j)=ql(:,i,j)+shape(:,i)*WdetJ*uhessl(:,j ) 225 call local (uhess, ql, ien, 27, 'scatter ')
|
| H A D | e3ql.f | 2 & xl, ql, xmudmi, 28 & ql(npro,nshl,idflx), xmudmi(npro,ngauss) 165 ql(iel,:,j) = matmul( rminv(iel,:,:),qrl(iel,:,j) ) 178 & xl, ql, sgn ) 194 & ql(npro,nshl,nsd) 311 ql(iel,:,j) = matmul( rminv(iel,:,:),qrl(iel,:,j) )
|
| H A D | advLES.f | 26 dimension qres(nshg,nsd*nsd), ql(npro,nshl,nsd*nsd), 62 call local (qres, ql, ien, nsd*nsd, 'gather ') 67 & xl, ql, xmudmi, 250 divqi(:,1) = divqi(:,1) + shg(:,n,1)*ql(:,n,1 ) 251 & + shg(:,n,2)*ql(:,n,4 ) 252 & + shg(:,n,3)*ql(:,n,7 ) 254 divqi(:,2) = divqi(:,2) + shg(:,n,1)*ql(:,n,2 ) 255 & + shg(:,n,2)*ql(:,n,5 ) 256 & + shg(:,n,3)*ql(:,n,8) 258 divqi(:,3) = divqi(:,3) + shg(:,n,1)*ql(:,n,3 ) [all …]
|
| /phasta/phSolver/compressible/ |
| H A D | e3q.f | 2 & xl, ql, rmassl, argument 28 & ql(npro,nshl,idflx), rmassl(npro,nshl), 154 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*qdi(:,1 ) 155 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*qdi(:,2 ) 156 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*qdi(:,3 ) 157 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*qdi(:,4 ) 158 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*qdi(:,5 ) 159 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*qdi(:,6 ) 160 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*qdi(:,7 ) 161 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*qdi(:,8 ) [all …]
|
| H A D | asiq.f | 35 & ql(npro,nshl,idflx), rmassl(npro,nshl), 55 ql = zero 59 & xl, ql, rmassl, 65 call local (qres, ql, ien, idflx,'scatter ')
|
| H A D | e3ivar.f | 10 & ql, divqi, sgn, tmp, argument 79 & ql(npro,nshl,idflx), 378 divqi(:,1) = divqi(:,1) + shg(:,n,1)*ql(:,n,1 ) 379 & + shg(:,n,2)*ql(:,n,5 ) 380 & + shg(:,n,3)*ql(:,n,9 ) 382 divqi(:,2) = divqi(:,2) + shg(:,n,1)*ql(:,n,2 ) 383 & + shg(:,n,2)*ql(:,n,6 ) 384 & + shg(:,n,3)*ql(:,n,10) 386 divqi(:,3) = divqi(:,3) + shg(:,n,1)*ql(:,n,3 ) 387 & + shg(:,n,2)*ql(:,n,7 ) [all …]
|
| H A D | asiflx.f | 26 dimension ql(npro,nshl,(nflow-1)*nsd) 58 ql = zero !only used now in elmmfg 65 & xmudmi, BDiagl, ql, sgn, rlsl, EGmassd)
|
| H A D | asimfg.f | 31 & ql(npro,nshl,idflx) 54 & call local (qres, ql, ien, idflx, 'gather ') 79 & BDiagl, ql, sgn, rlsl, EGmassd,
|
| H A D | asires.f | 25 & rml(npro,nshl,nflow), ql(npro,nshl,(nflow-1)*nsd) 80 & xmudmi, BDiagl, ql, sgn, rlsl, EGmassd,
|
| H A D | asigmr.f | 33 & ql(npro,nshl,idflx) 55 call local (qres, ql, ien, idflx, 'gather ') 77 & BDiagl, ql, sgn, rlsl, EGmass,
|
| H A D | e3ql.f | 2 & xl, ql, xmudmi, 28 & ql(npro,nshl,idflx), xmudmi(npro,ngauss) 210 ql(iel,:,j) = matmul( rminv(iel,:,:),qrl(iel,:,j) )
|
| H A D | e3.f | 3 & BDiagl, ql, sgn, rlsl, EGmass, argument 55 & rl(npro,nshl,nflow), ql(npro,nshl,idflx), 101 call e3ql (ycl, shp, shgl, xl, ql, xmudmi, sgn) 138 & ql, divqi, sgn,
|