| /phasta/phSolver/incompressible/ ! |
| H A D | e3q.f | 34 & dxidx(npro,nsd,nsd), WdetJ(npro), 76 & dxidx, WdetJ ) 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 | e3qGradV.f | 32 & dxidx(npro,nsd,nsd), WdetJ(npro) 57 & dxidx, WdetJ ) 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 | filters.f | 12 & fresl(npro,16), WdetJ(npro), 84 WdetJ(:) = wght / tmp(:) 136 fresli(:,1) = WdetJ * u1 !G * u1 * WdetJ 137 fresli(:,2) = WdetJ * u2 !G * u2 * WdetJ 138 fresli(:,3) = WdetJ * u3 !G * u3 * WdetJ 140 fresli(:,4) = WdetJ * u1 * u1 ! G*u1*u1*WdetJ 141 fresli(:,5) = WdetJ * u2 * u2 ! G*u2*u2*WdetJ 142 fresli(:,6) = WdetJ * u3 * u3 ! G*u3*u3*WdetJ 143 fresli(:,7) = WdetJ * u1 * u2 ! G*u1*u2*WdetJ 144 fresli(:,8) = WdetJ * u1 * u3 ! G*u1*u3*WdetJ [all …]
|
| H A D | e3.f | 54 & WdetJ(npro), rho(npro), 102 & WdetJ, rho, pres, 118 & uBar, aci, WdetJ, 130 & uBar, WdetJ, rho, 189 & dxidx(npro,nsd,nsd), WdetJ(npro), local 234 & shg, dxidx, WdetJ, 264 & WdetJ, rLS, tauS, 274 & WdetJ, rLS, tauS,
|
| H A D | e3lhs.f | 2 & uBar, WdetJ, rho, argument 36 & uBar(npro,3), WdetJ(npro), rho(npro), 58 tlW = lhsFct * WdetJ 63 tsFct = lhmFct * WdetJ * rho 65 tauBar = lhsFct * WdetJ * tauBar 272 & WdetJ, rLS, tauS, argument 282 & WdetJ(npro), rLS(npro), rho(npro), local 304 tauS = lhsFct * WdetJ * tauS 305 kptmp = lhsFct * WdetJ * diffus 306 tsFct = almi * WdetJ * (one - flmpl) [all …]
|
| H A D | e3qvar.f | 4 & dxidx, WdetJ ) argument 41 & dxidx(npro,nsd,nsd), WdetJ(npro) 91 WdetJ = Qwt(lcsyst,intp)/ tmp 149 & gradT, dxidx, WdetJ ) argument 157 & dxidx(npro,nsd,nsd), WdetJ(npro), local 166 & shg, WdetJ)
|
| H A D | hessian.f | 100 & WdetJ(npro), rmassl(npro,nshl) 135 & gradul , dxidx, WdetJ ) 141 ql(:,i,j) = ql(:,i,j)+shape(:,i)*WdetJ*gradul(:,j) 144 rmassl(:,i) = rmassl(:,i) + shape(:,i)*WdetJ 175 & WdetJ(npro), rmassl(npro, nshl) 210 & uhessl , dxidx, WdetJ ) 216 ql(:,i,j)=ql(:,i,j)+shape(:,i)*WdetJ*uhessl(:,j ) 219 rmassl(:,i) = rmassl(:,i) + shape(:,i)*WdetJ 289 & gradient, dxidx, WdetJ ) argument 301 & dxidx(npro,nsd,nsd), WdetJ(npro) [all …]
|
| H A D | e3res.f | 2 & uBar, aci, WdetJ, argument 42 & uBar(npro,nsd), aci(npro,nsd), WdetJ(npro), 231 rl(:,aa,4) = rl(:,aa,4) + WdetJ 238 rl(:,aa,1) = rl(:,aa,1) - WdetJ 243 rl(:,aa,2) = rl(:,aa,2) - WdetJ 248 rl(:,aa,3) = rl(:,aa,3) - WdetJ 270 & WdetJ, rLS, tauS, argument 279 & WdetJ(npro), rLS(npro), rho(npro), local 334 rl(:,aa) = rl(:,aa) - WdetJ
|
| H A D | errsmooth.f | 146 & shdrv(npro,nsd,nshl), WdetJ(npro), 185 & shg, WdetJ) 195 & + shape(:,i)*WdetJ*error(:,j) 198 rmassl(:,i) = rmassl(:,i) + shape(:,i)*WdetJ
|
| H A D | genlmass.f | 81 real*8 xl(npro,nenl,nsd), WdetJ(npro), local 121 & shg, WdetJ) 129 & + shape(:,aa)*shape(:,aa)*WdetJ 135 fmstot= fmstot + WdetJ ! intregral d Omega
|
| /phasta/phSolver/compressible/ ! |
| H A D | e3wmlt.f | 1 subroutine e3wmlt (shp, shg, WdetJ, argument 55 & WdetJ(npro), ri(npro,nflow*(nsd+1)), 75 rl(:,i,1) = rl(:,i,1) + WdetJ * ( 78 rl(:,i,2) = rl(:,i,2) + WdetJ * ( 81 rl(:,i,3) = rl(:,i,3) + WdetJ * ( 84 rl(:,i,4) = rl(:,i,4) + WdetJ * ( 87 rl(:,i,5) = rl(:,i,5) + WdetJ * ( 99 rml(:,i,1) = rml(:,i,1) + WdetJ * ( 103 rml(:,i,2) = rml(:,i,2) + WdetJ * ( 107 rml(:,i,3) = rml(:,i,3) + WdetJ * ( [all …]
|
| H A D | e3q.f | 35 & dxidx(npro,nsd,nsd), WdetJ(npro), 83 & dxidx, WdetJ, T, 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 | e3massl.f | 1 subroutine e3massl (bcool, shape, WdetJ, A0, argument 22 & WdetJ(npro), EGmass(npro,nedof,nedof) 29 temp = WdetJ * almi/gami/alfi*Dtgl 30 bcool = WdetJ * bcool ! note that bcool is not used after this 80 subroutine e3masslSclr (shape, WdetJ, A0t, argument 101 & WdetJ(npro), EGmasst(npro,nshape,nshape) 108 temp = WdetJ * almi/gami/alfi*Dtgl 121 & - shpij*WdetJ(:) * srcp(:)
|
| H A D | vorticity.f | 92 real*8 WdetJ(npro) local 123 & shg, WdetJ) 143 alph2(:) = alph2(:) + shp(:,n)*shp(:,n)*WdetJ(:) 146 volum(:) = volum(:) + WdetJ(:) 151 lpmass(:,n) = lpmass(:,n) + WdetJ(:)*shp(:,n)*shp(:,n) 155 & *WdetJ(:)*shp(:,n) 159 & *WdetJ(:)*shp(:,n) 163 & *WdetJ(:)*shp(:,n) 174 & *WdetJ(:)*shp(:,n)
|
| H A D | e3.f | 67 dimension dxidx(npro,nsd,nsd), WdetJ(npro) 133 & shg, dxidx, WdetJ, 169 & shg, shape, WdetJ) 209 & EGmass, stiff, WdetJ, 245 call e3juel (yl, acl, Sclr, A0, WdetJ, rl, rml) 255 call e3massl (bcool,shape, WdetJ, A0, EGmass) 266 call e3bdg(shape, shg, WdetJ, 275 call e3bdg_nd(shape, shg, WdetJ, 292 call e3wmlt (shape, shg, WdetJ, 378 dimension dxidx(npro,nsd,nsd), WdetJ(npro) [all …]
|
| H A D | e3conv.f | 7 & shg, shape, WdetJ ) argument 60 & WdetJ(npro) 197 fact1 = WdetJ * shg(:,j,1) 198 fact2 = WdetJ * shg(:,j,2) 199 fact3 = WdetJ * shg(:,j,3) 361 & shg, shape, WdetJ) argument 409 & WdetJ(npro) 466 & WdetJ * shg(:,j,1) * A1t(:) 467 & + WdetJ * shg(:,j,2) * A2t(:) 468 & + WdetJ * shg(:,j,3) * A3t(:)
|
| H A D | e3bdg.f | 1 subroutine e3bdg (shp, shg, WdetJ, argument 14 & shg(npro,nshl,nsd), WdetJ(npro), 38 tmp=shp(:,j)*WdetJ 197 tmp = 0.4d0 * WdetJ * Dtgl/Qwt(lcsyst,intp)/three*almi/gami/alfi 223 tmp = WdetJ * (Dtgl * ff + bcool) 331 tmp1 = shp(:,i) * WdetJ * almi/gami/alfi*dtgl 612 tmp = WdetJ * shg(:,j,1) * shg(:,j,1) 873 tmp = WdetJ * shg(:,j,2) * shg(:,j,2) 1134 tmp = WdetJ * shg(:,j,3) * shg(:,j,3) 1617 tmp = WdetJ * shg(:,j,1) * shg(:,j,2) [all …]
|
| H A D | e3qvar.f | 3 & shg, dxidx, WdetJ, T, argument 45 & dxidx(npro,nsd,nsd), WdetJ(npro), 104 & shg, WdetJ)
|
| H A D | e3juel.f | 2 & WdetJ, rl, rml) argument 26 & WdetJ(npro), A0(npro,nflow,nflow), 50 fact=WdetJ/(Qwt(lcsyst,intp)*15.0d0)
|
| /phasta/phSolver/common/ ! |
| H A D | cmass.f | 41 & WdetJ(npro) 56 call e3metric( xl, shdrv,dxidx,shg,WdetJ) 62 em(iel,na,nb) = em(iel,na,nb) + shp2*WdetJ(iel) 95 & WdetJ(npro) 120 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf) 123 eml(:,i) = eml(:,i) + shape(:,i)*WdetJ(:) 135 WdetJ = zero 155 call e3metric( xl, shdrv,dxidx,shg,WdetJ) 162 em(iel,na,nb) = em(iel,na,nb) + shp2*WdetJ(iel) 218 & WdetJ(npro) [all …]
|
| H A D | asithf.f | 8 & fresl(npro,24), WdetJ(npro), 83 WdetJ = wght / tmp 141 fresli(:,10) = fresli(:,10) * WdetJ ! u_{1,1}*WdetJ 142 fresli(:,11) = fresli(:,11) * WdetJ ! u_{2,2}*WdetJ 143 fresli(:,12) = fresli(:,12) * WdetJ ! u_{3,3}*WdetJ 144 fresli(:,13) = fresli(:,13) * WdetJ ! (1/2)*(u_{1,2}+u_{2,1})*WdetJ 145 fresli(:,14) = fresli(:,14) * WdetJ ! (1/2)*(u_{1,3}+u_{3,1})*WdetJ 146 fresli(:,15) = fresli(:,15) * WdetJ ! (1/2)*(u_{2,3}+u_{3,2})*WdetJ 148 fresli(:,22) = fresli(:,22) * WdetJ !rho * WdetJ 178 fresli(:,23) = WdetJ ! Integral of 1 over the element
|
| H A D | bardmc.f | 679 & fresl(npro,nshl,11), WdetJ(npro), 752 WdetJ(:) = wght / tmp(:) 759 rho = rho * WdetJ !rho * WdetJ 777 fresli(:,i,5) = WdetJ*shp(i,intp) !Integral of filter kernel 833 & fresl(npro,33), WdetJ(npro), 909 WdetJ = wght / tmp 927 fresli(:,28) = WdetJ(:) * fresli(:,25) * 929 fresli(:,29) = WdetJ(:) * fresli(:,26) * 931 fresli(:,30) = WdetJ(:) * fresli(:,27) * 933 fresli(:,31) = WdetJ(:) * fresli(:,25) * [all …]
|
| H A D | solvecon.f | 176 & dxidx(npro,nsd,nsd), WdetJ(npro) 240 & shg, WdetJ) 268 v_lambdal1(:,i)=v_lambdal1(:,i)+ shape(:,i)*WdetJ*tmp1 269 v_lambdal2(:,i)=v_lambdal2(:,i)+ shape(:,i)*WdetJ*tmp2 271 hprimel(:,i) = hprimel(:,i) + shape(:,i)*WdetJ*h_prime 272 rmassl(:,i) =rmassl(:,i) + shape(:,i)*WdetJ
|
| H A D | errorsmooth.f | 145 & shdrv(npro,nsd,nshl), WdetJ(npro), 184 & shg, WdetJ) 194 & + shape(:,i)*WdetJ*error(:,j) 197 rmassl(:,i) = rmassl(:,i) + shape(:,i)*WdetJ
|
| H A D | e3metric.f | 9 & shg, WdetJ) argument 15 & WdetJ(npro) local 63 WdetJ = Qwt(lcsyst,intp) / tmp
|