| /phasta/phSolver/incompressible/ |
| H A D | e3sts.f | 1 subroutine e3StsLhs( xl, lStsVec ) argument 13 & xl(npro,nenl,3) local 15 call e3StsDir( xl, lDir ) 37 subroutine e3StsRes( xl, rl, lStsVec ) argument 46 real*8 xl(npro,nenl,3), rl(npro,nshl,ndof) local 49 call e3StsDir( xl, lDir ) 74 subroutine e3StsDir( xl, lDir ) argument 82 real*8 xl(npro,nenl,3), lDir(npro,nshl,3) local 92 x12 = xl(e,2,1) - xl(e,1,1) 93 x13 = xl(e,3,1) - xl(e,1,1) [all …]
|
| H A D | e3qvar.f | 2 & xl, g1yi, argument 38 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd), 53 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n) 54 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n) 55 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n) 56 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n) 57 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n) 58 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n) 59 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n) 60 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n) [all …]
|
| H A D | hessian.f | 98 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 121 call localx (x, xl, ien, nsd, 'gather ') 134 call local_gradient( yl(:,:,2:4), 3, shdrv, xl, 173 dimension gradul(npro,nshl,9), xl(npro,nenl,nsd), 196 call localx (x, xl, ien, nsd, 'gather ') 209 call local_gradient( gradul, 9, shdrv, xl, 288 subroutine local_gradient ( vector, vsize, shgl, xl, argument 299 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd), 313 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n) 314 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n) [all …]
|
| H A D | asigmr.f | 31 & xl(npro,nenl,nsd), dwl(npro,nenl), 60 call localx(x, xl, ien, nsd, 'gather ') 87 & shgl, xl, rl, 94 call e3StsRes ( xl, rl, lStsVec ) 111 call timeseries(yl,xl,ien,sgn) 148 & xl(npro,nenl,nsd), local 166 call localx(x, xl, ien, nsd, 'gather ') 181 & shgl, xl, dwl,
|
| H A D | asiq.f | 32 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 54 call localx (x, xl, ien, nsd, 'gather ') 65 & xl, ql, rmassl, 100 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 112 call localx (x, xl, ien, nsd, 'gather ') 123 & xl, ql, rmassl,
|
| H A D | filters.f | 10 dimension x(numnp,3), xl(npro,nenl,3) 26 call localx (x, xl, ien, 3, 'gather ') 44 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp) 45 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp) 46 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp) 47 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp) 48 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp) 49 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp) 50 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp) 51 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp) [all …]
|
| H A D | getdiff.f | 1 subroutine getDiff( dwl,yl, shape, xmudmi, xl, rmu, rho) argument 12 & xl(npro,nenl,nsd) 97 call EviscDESIC (xl,rmu,xmudmi) 104 subroutine EviscDESIC(xl,xmut,xmudmi) argument 107 real*8 xmut(npro),xl(npro,nenl,nsd),xmudmi(npro,ngauss) local 111 dx=maxval(xl(i,:,1))-minval(xl(i,:,1)) 112 dy=maxval(xl(i,:,2))-minval(xl(i,:,2)) 113 dz=maxval(xl(i,:,3))-minval(xl(i,:,3))
|
| H A D | e3.f | 2 & shgl, xl, rl, ql, argument 43 & xl(npro,nenl,nsd), dwl(npro,nenl), 76 & xl, ql, xmudmi, 94 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho) 99 & shdrv, xl, 164 & shgl, xl, dwl, argument 181 & xl(npro,nenl,nsd), rl(npro,nshl), local 211 call e3qlSclr (yl, dwl, shp, shgl, xl, ql, sgn) 232 & shdrv, xl, xmudmi,
|
| H A D | elmStats.f | 14 real*8, allocatable :: xl(:,:,:) local 30 allocate ( xl(npro,nenl,3) ) 35 call localx ( x, xl, mien(iblk)%p, nsd, 'gather ' ) 39 call e3StsLhs( xl, lStsVec ) 46 deallocate ( xl )
|
| H A D | asiqGradV.f | 30 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd), 51 call localx (x, xl, ien, nsd, 'gather ') 58 & xl, ql,
|
| H A D | e3ivar.f | 2 & shgl, xl, argument 57 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd), 102 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1) 103 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2) 104 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3) 143 call e3metric( xl, shgl, dxidx, 292 dkei = dkei + shpfun(:,n) * (1.0-xl(:,n,2)**2) !u_ex^~ (in FEM space) 315 & shgl, xl, xmudmi, argument 330 & shgl(npro,nsd,nshl), xl(npro,nenl,nsd), 378 call e3metric( xl, shgl, dxidx, [all …]
|
| H A D | e3q.f | 2 & xl, ql, rmassl, argument 26 & xl(npro,nenl,nsd), 74 & xl, g1yi, 87 call getdiff(dwl, yl, shape, xmudmi, xl, rmu, tmp) 207 & xl, ql, rmassl, argument 221 & xl(npro,nenl,nsd), 259 & xl, gradT,
|
| H A D | e3ql.f | 2 & xl, ql, xmudmi, argument 27 & xl(npro,nenl,nsd), sgn(npro,nshl), 62 & xl, g1yi, 66 call getdiff(dwl, yl, shape, xmudmi, xl,rmu, tmp) 178 & xl, ql, sgn ) argument 193 & xl(npro,nenl,nsd), sgn(npro,nshl), 224 call e3qvarSclr (yl, shdrv, xl,
|
| H A D | e3qGradV.f | 2 & xl, ql, argument 25 & xl(npro,nenl,nsd), 55 & xl, g1yi,
|
| /phasta/phSolver/common/ |
| H A D | e3metric.f | 8 subroutine e3metric( xl, shgl, dxidx, argument 13 real*8 xl(npro,nenl,nsd), shgl(npro,nsd,nshl), local 25 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n) 26 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n) 27 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n) 28 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n) 29 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n) 30 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n) 31 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n) 32 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n) [all …]
|
| H A D | cmass.f | 20 subroutine cmass (shp, shgl, xl, em) argument 36 & xl(npro,nenl,nsd) 56 call e3metric( xl, shdrv,dxidx,shg,WdetJ) 73 subroutine cmassl (shp, shgl, shpf, shglf, xl, em, Qwtf) argument 90 & xl(npro,nenl,nsd) 120 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf) 155 call e3metric( xl, shdrv,dxidx,shg,WdetJ) 196 subroutine cmasstl (shp, shgl, shpf, shglf, xl, em, Qwtf) argument 213 & xl(npro,nenl,nsd) 237 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf) [all …]
|
| H A D | getstrl.f | 6 dimension x(numnp,3), xl(npro,nenl,3) 16 call localx (x, xl, ien, 3, 'gather ') 37 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp) 38 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp) 39 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp) 40 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp) 41 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp) 42 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp) 43 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp) 44 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp) [all …]
|
| H A D | asithf.f | 6 dimension x(numnp,nsd), xl(npro,nenl,nsd) 19 call localx (x, xl, ien, 3, 'gather ') 43 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp) 44 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp) 45 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp) 46 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp) 47 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp) 48 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp) 49 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp) 50 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp) [all …]
|
| H A D | bardmc.f | 677 dimension x(numnp,3), xl(npro,nenl,3) 690 call localx (x, xl, ien, 3, 'gather ') 713 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp) 714 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp) 715 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp) 716 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp) 717 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp) 718 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp) 719 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp) 720 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp) [all …]
|
| H A D | hierarchic.f | 248 real*8, allocatable :: xl(:,:,:) local 274 allocate ( xl(npro,nenl,nsd ) ) 280 call localx( x, xl, mien(iblk)%p, nsd, 'gather ') 282 call eval ( xl, ycl, yvl, 296 deallocate ( xl ) 317 subroutine eval( xl, ycl, yvl, argument 328 & shgl(nsd,nshl,npts), xl(npro,nenl,nsd), local 351 call e3metric( xl, shdrv, dxidx,
|
| /phasta/phSolver/compressible/ |
| H A D | getdiff.f | 3 & xmudmi, xl) argument 45 & xmudmi(npro,ngauss), xl(npro,nenl,nsd), 56 xx(:)=xx(:) + shp(:,n) * xl(:,n,1) 257 call EviscDES (xl,xmut,xmudmi) 364 subroutine EviscDES(xl,xmut,xmudmi) argument 367 real*8 xmut(npro),xl(npro,nenl,nsd),xmudmi(npro,ngauss) local 371 dx=maxval(xl(i,:,1))-minval(xl(i,:,1)) 372 dy=maxval(xl(i,:,2))-minval(xl(i,:,2)) 373 dz=maxval(xl(i,:,3))-minval(xl(i,:,3))
|
| H A D | e3source.f | 5 & xl, shpfun, bcool) 34 & xl(npro,nenl,3), xx(npro,3) 66 xfind(i) = (xsor-minval(xl(i,:,1))) 67 & *(maxval(xl(i,:,1))-xsor) 82 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1) 83 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2) 84 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3) 230 & u2, u3, xl, elDwl) 292 & xl(npro,nenl,nsd) 300 dx=maxval(xl(i,:,1))-minval(xl(i,:,1)) [all …]
|
| H A D | asigmr.f | 30 & xl(npro,nenl,nsd), ytargetl(npro,nshl,nflow), 54 call localx(x, xl, ien, nsd, 'gather ') 76 & shgl, xl, rl, rml, xmudmi, 111 call timeseries(ycl,xl,ien,sgn) 152 & xl(npro,nenl,nsd), Diagl(npro,nshl), 168 call localx (x, xl, ien, nsd, 'gather ') 183 & sgn, shgl, xl,
|
| H A D | asiq.f | 34 dimension ycl(npro,nshl,ndof), xl(npro,nenl,nsd), 51 call localx(x, xl, ien, nsd, 'gather ') 59 & xl, ql, rmassl,
|
| H A D | asiflx.f | 23 & xl(npro,nenl,nsd), 45 call localx(x, xl, ien, nsd, 'gather ') 64 & shgl(lcsyst,:,1:nshl,:), xl, rl, rml,
|