| /phasta/phSolver/incompressible/ |
| H A D | e3lhs.f | 5 & shpfun, shg, xKebe, argument 39 & shpfun(npro,nshl),shg(npro,nshl,3) 80 t1(:,1) = uBar(:,1) * shg(:,b,1) 81 & + uBar(:,2) * shg(:,b,2) 82 & + uBar(:,3) * shg(:,b,3) 104 t1(:,1) = tauC * shg(:,b,1) 105 t1(:,2) = tauC * shg(:,b,2) 106 t1(:,3) = tauC * shg(:,b,3) 110 t2(:,1) = rmu * shg(:,b,1) 111 t2(:,2) = rmu * shg(:,b,2) [all …]
|
| H A D | e3ivar.f | 4 & g3yi, shg, dxidx, argument 60 & shg(npro,nshl,nsd), dxidx(npro,nsd,nsd), 144 & shg, WdetJ) 153 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1) 154 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2) 155 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3) 156 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4) 158 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * yl(:,n,1) 159 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) 160 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) [all …]
|
| H A D | e3qvar.f | 3 & g2yi, g3yi, shg, argument 40 & g3yi(npro,nflow), shg(npro,nshl,nsd), 107 shg(:,n,1) = shgl(:,1,n) * dxidx(:,1,1) + 110 shg(:,n,2) = shgl(:,1,n) * dxidx(:,1,2) + 113 shg(:,n,3) = shgl(:,1,n) * dxidx(:,1,3) + 122 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2) 123 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3) 124 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4) 126 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) 127 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) [all …]
|
| H A D | e3source.f | 211 & dwl, shape_funct, shg, argument 241 & shape_funct(npro,nshl), shg(npro,nshl,nsd), local 293 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2) 294 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3) 295 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4) 297 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2) 298 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3) 299 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4) 301 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2) 302 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3) [all …]
|
| H A D | e3res.f | 6 & shpfun, shg, src, argument 46 & shpfun(npro,nshl),shg(npro,nshl,nsd), src(npro,nsd), 232 & * ( shg(:,aa,1) * uBar(:,1) 233 & + shg(:,aa,2) * uBar(:,2) 234 & + shg(:,aa,3) * uBar(:,3) ) 240 & + shg(:,aa,1) * rGNa(:,1,1) 241 & + shg(:,aa,2) * rGNa(:,1,2) 242 & + shg(:,aa,3) * rGNa(:,1,3) ) 245 & + shg(:,aa,1) * rGNa(:,2,1) 246 & + shg(:,aa,2) * rGNa(:,2,2) [all …]
|
| H A D | e3.f | 52 & g3yi(npro,ndof), shg(npro,nshl,nsd), 101 & g3yi, shg, dxidx, 122 & shpfun, shg, src, 133 & shpfun, shg, xKebe, 187 real*8 gradS(npro,nsd), shg(npro,nshl,nsd), local 234 & shg, dxidx, WdetJ, 265 & shpfun, shg, srcR, 275 & shpfun, shg, srcL,
|
| H A D | filters.f | 19 & shgl(nsd,nshl,ngauss), shg(npro,nshl,nsd), 91 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1) 94 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2) 97 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3) 120 S11 = S11 + shg(:,i,1)*yl(:,i,2) 121 S22 = S22 + shg(:,i,2)*yl(:,i,3) 122 S33 = S33 + shg(:,i,3)*yl(:,i,4) 124 S12 = S12 + shg(:,i,2)*yl(:,i,2) 125 & +shg(:,i,1)*yl(:,i,3) 126 S13 = S13 + shg(:,i,3)*yl(:,i,2) [all …]
|
| H A D | e3q.f | 33 & g3yi(npro,nflow), shg(npro,nshl,nsd), 75 & g2yi, g3yi, shg, 168 g1yti(:) = g1yti(:) + shg(:,n,1) * yl(:,n,6) 169 g2yti(:) = g2yti(:) + shg(:,n,2) * yl(:,n,6) 170 g3yti(:) = g3yti(:) + shg(:,n,3) * yl(:,n,6)
|
| /phasta/phSolver/compressible/ |
| H A D | e3ivar.f | 5 & shg, dxidx, WdetJ, argument 70 & shg(npro,nshl,nsd), dxidx(npro,nsd,nsd), 238 & shg, WdetJ) 260 g1yi(:,1) = g1yi(:,1) + shg(:,1,1) * yl(:,1,1) 261 & + shg(:,2,1) * yl(:,2,1) 262 & + shg(:,3,1) * yl(:,3,1) 263 & + shg(:,4,1) * yl(:,4,1) 265 g1yi(:,2) = g1yi(:,2) + shg(:,1,1) * yl(:,1,2) 266 & + shg(:,2,1) * yl(:,2,2) 267 & + shg(:,3,1) * yl(:,3,2) [all …]
|
| H A D | e3wmlt.f | 1 subroutine e3wmlt (shp, shg, WdetJ, argument 54 & shg(npro,nshl,nsd), 76 & shg(:,i,1) * ri(:, 1) + shg(:,i,2) * ri(:, 6) 77 & + shg(:,i,3) * ri(:,11) ) 79 & shg(:,i,1) * ri(:, 2) + shg(:,i,2) * ri(:, 7) 80 & + shg(:,i,3) * ri(:,12) ) 82 & shg(:,i,1) * ri(:, 3) + shg(:,i,2) * ri(:, 8) 83 & + shg(:,i,3) * ri(:,13) ) 85 & shg(:,i,1) * ri(:, 4) + shg(:,i,2) * ri(:, 9) 86 & + shg(:,i,3) * ri(:,14) ) [all …]
|
| H A D | e3qvar.f | 3 & shg, dxidx, WdetJ, T, argument 44 & g3yi(npro,nflow), shg(npro,nshl,nsd), 104 & shg, WdetJ) 120 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * ycl(:,n,1) 121 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * ycl(:,n,2) 122 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * ycl(:,n,3) 123 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * ycl(:,n,4) 124 g1yi(:,5) = g1yi(:,5) + shg(:,n,1) * ycl(:,n,5) 126 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * ycl(:,n,1) 127 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * ycl(:,n,2) [all …]
|
| H A D | e3bdg.f | 1 subroutine e3bdg (shp, shg, WdetJ, argument 14 & shg(npro,nshl,nsd), WdetJ(npro), 41 & shg(:,j,1) * A1(:,1,1) + 42 & shg(:,j,2) * A2(:,1,1) + 43 & shg(:,j,3) * A3(:,1,1) 47 & shg(:,j,1) * A1(:,1,2) 54 & +shg(:,j,2) * A2(:,1,3) 61 & shg(:,j,3) * A3(:,1,4) 65 & shg(:,j,1) * A1(:,1,5) + 66 & shg(:,j,2) * A2(:,1,5) + [all …]
|
| H A D | e3bdg_nd.f | 1 subroutine e3bdg_nd (shp, shg, WdetJ, argument 14 & shg(npro,nshl,nsd), WdetJ(npro), 42 & shg(:,j,1) * A1(:,1,1) + 43 & shg(:,j,2) * A2(:,1,1) + 44 & shg(:,j,3) * A3(:,1,1) 48 & shg(:,j,1) * A1(:,1,2) 55 & +shg(:,j,2) * A2(:,1,3) 62 & shg(:,j,3) * A3(:,1,4) 66 & shg(:,j,1) * A1(:,1,5) + 67 & shg(:,j,2) * A2(:,1,5) + [all …]
|
| H A D | vorticity.f | 99 real*8 shg(npro,nshl,nsd) local 123 & shg, WdetJ) 131 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2) ! du/dx 132 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) ! du/dy 133 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) ! du/dz 135 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3) ! dv/dx 136 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) ! dv/dy 137 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) ! dv/dz 139 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4) ! dw/dx 140 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4) ! dw/dy [all …]
|
| H A D | e3conv.f | 7 & shg, shape, WdetJ ) argument 59 & shg(npro,nshl,nsd), shape(npro,nshl), 197 fact1 = WdetJ * shg(:,j,1) 198 fact2 = WdetJ * shg(:,j,2) 199 fact3 = WdetJ * shg(:,j,3) 361 & shg, shape, WdetJ) argument 408 & shg(npro,nshl,nsd), shape(npro,nshl), 466 & WdetJ * shg(:,j,1) * A1t(:) 467 & + WdetJ * shg(:,j,2) * A2t(:) 468 & + WdetJ * shg(:,j,3) * A3t(:)
|
| H A D | e3.f | 64 & g3yi(npro,nflow), shg(npro,nshl,nsd), 133 & shg, dxidx, WdetJ, 169 & shg, shape, WdetJ) 208 & divqi, shape, shg, 266 call e3bdg(shape, shg, WdetJ, 275 call e3bdg_nd(shape, shg, WdetJ, 292 call e3wmlt (shape, shg, WdetJ, 376 dimension shg(npro,nshl,nsd) 437 & u3, shg, dwl, 473 & shg, shape, WdetJ) [all …]
|
| H A D | e3q.f | 34 & g3yi(npro,nflow), shg(npro,nshl,nsd), 82 & g2yi, g3yi, shg, 210 g1yti(:) = g1yti(:) + shg(:,n,1) * ycl(:,n,6) 211 g2yti(:) = g2yti(:) + shg(:,n,2) * ycl(:,n,6) 212 g3yti(:) = g3yti(:) + shg(:,n,3) * ycl(:,n,6)
|
| /phasta/phSolver/common/ |
| H A D | getstrl.f | 11 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd), 88 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1) 91 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2) 94 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3) 103 fresli(:,j) = fresli(:,j)+shg(:,i,ig)*yl(:,i,iv) 111 fresli(:,13) = fresli(:,13)+shg(:,i,2)*yl(:,i,2) 112 & +shg(:,i,1)*yl(:,i,3) 113 fresli(:,14) = fresli(:,14)+shg(:,i,3)*yl(:,i,2) 114 & +shg(:,i,1)*yl(:,i,4) 115 fresli(:,15) = fresli(:,15)+shg(:,i,3)*yl(:,i,3) [all …]
|
| H A D | cmass.f | 40 & shg(npro,nshl,nsd), 56 call e3metric( xl, shdrv,dxidx,shg,WdetJ) 94 & shg(npro,nshl,nsd), Qwtf(ngaussf), 120 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf) 136 shg = zero 155 call e3metric( xl, shdrv,dxidx,shg,WdetJ) 217 & shg(npro,nshl,nsd), Qwtf(ngaussf), 237 call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf) 253 shg = zero 272 call e3metric( xl, shdrv,dxidx,shg,WdetJ) [all …]
|
| H A D | asithf.f | 12 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd), 97 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1) 100 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2) 103 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3) 112 fresli(:,j) = fresli(:,j)+shg(:,i,ig)*ycl(:,i,iv) 120 fresli(:,13) = fresli(:,13)+shg(:,i,2)*ycl(:,i,2) 121 & +shg(:,i,1)*ycl(:,i,3) 122 fresli(:,14) = fresli(:,14)+shg(:,i,3)*ycl(:,i,2) 123 & +shg(:,i,1)*ycl(:,i,4) 124 fresli(:,15) = fresli(:,15)+shg(:,i,3)*ycl(:,i,3) [all …]
|
| H A D | hierarchic.f | 329 & shg(npro,nshl,nsd), gradV(npro,nsd,nsd), local 352 & shg, tmp) 356 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * ycl(:,n,2) 357 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * ycl(:,n,3) 358 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * ycl(:,n,4) 360 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * ycl(:,n,2) 361 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * ycl(:,n,3) 362 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * ycl(:,n,4) 364 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * ycl(:,n,2) 365 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * ycl(:,n,3) [all …]
|
| H A D | e3metric.f | 9 & shg, WdetJ) argument 14 & dxidx(npro,nsd,nsd), shg(npro,nshl,nsd), local 68 shg(:,n,1) = shgl(:,1,n) * dxidx(:,1,1) + 71 shg(:,n,2) = shgl(:,1,n) * dxidx(:,1,2) + 74 shg(:,n,3) = shgl(:,1,n) * dxidx(:,1,3) +
|
| H A D | bardmc.f | 683 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd), 837 & shgl(nsd,nshl,maxsh), shg(npro,nshl,nsd), 945 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1) 948 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2) 951 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3) 960 fresli(:,j) = fresli(:,j)+shg(:,i,ig)*yl(:,i,iv) 968 fresli(:,13) = fresli(:,13)+shg(:,i,2)*yl(:,i,2) 969 & +shg(:,i,1)*yl(:,i,3) 970 fresli(:,14) = fresli(:,14)+shg(:,i,3)*yl(:,i,2) 971 & +shg(:,i,1)*yl(:,i,4) [all …]
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 211 & dwl, shape_funct, shg, 241 & shape_funct(npro,nshl), shg(npro,nshl,nsd), 289 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2) 290 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3) 291 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4) 293 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2) 294 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3) 295 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4) 297 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2) 298 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3) [all …]
|
| H A D | e3source.fRHSonly | 211 & dwl, shape_funct, shg, 241 & shape_funct(npro,nshl), shg(npro,nshl,nsd), 289 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2) 290 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3) 291 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4) 293 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2) 294 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3) 295 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4) 297 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2) 298 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3) [all …]
|