phSolver/incompressible/filters.f

1       subroutine hfilterC (y, x, ien, hfres, shgl, shp, Qwtf)  argument
20      &          shp(nshl,ngauss),
116             u1 = u1 + shp(i,intp)*yl(:,i,2)
117             u2 = u2 + shp(i,intp)*yl(:,i,3)
118             u3 = u3 + shp(i,intp)*yl(:,i,4)
178      &     hfres, shgl, shp, Qwtf)  argument
190      &          shp(nshl,ngauss),
304          fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1)  ! hat{u1}
305          fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2)  ! hat{u2}
306          fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3)  ! hat{u3}
308          fresli(:,4) = fresli(:,4) + shp(i,intp)*hfresl(:,i,1)*
310          fresli(:,5) = fresli(:,5) + shp(i,intp)*hfresl(:,i,2)*
312          fresli(:,6) = fresli(:,6) + shp(i,intp)*hfresl(:,i,3)*
314          fresli(:,7) = fresli(:,7) + shp(i,intp)*hfresl(:,i,1)*
316          fresli(:,8) = fresli(:,8) + shp(i,intp)*hfresl(:,i,1)*
318          fresli(:,9) = fresli(:,9) + shp(i,intp)*hfresl(:,i,2)*
321          fresli(:,10) = fresli(:,10) + shp(i,intp)*hfresl(:,i,10)  ! hat{S11}
322          fresli(:,11) = fresli(:,11) + shp(i,intp)*hfresl(:,i,11)  ! hat{S22}
323          fresli(:,12) = fresli(:,12) + shp(i,intp)*hfresl(:,i,12)  ! hat{S33}
324          fresli(:,13) = fresli(:,13) + shp(i,intp)*hfresl(:,i,13)  ! hat{S12}
325          fresli(:,14) = fresli(:,14) + shp(i,intp)*hfresl(:,i,14)  ! hat{S13}
326          fresli(:,15) = fresli(:,15) + shp(i,intp)*hfresl(:,i,15)  ! hat{S23}
360          fresli(:,16) =fresli(:,16)+shp(i,intp)*SS11(:,i)! |hat{S}|*hat{S11}
361          fresli(:,17) =fresli(:,17)+shp(i,intp)*SS22(:,i)! |hat{S}|*hat{S22}
362          fresli(:,18) =fresli(:,18)+shp(i,intp)*SS33(:,i)! |hat{S}|*hat{S33}
363          fresli(:,19) =fresli(:,19)+shp(i,intp)*SS12(:,i)! |hat{S}|*hat{S12}
364          fresli(:,20) =fresli(:,20)+shp(i,intp)*SS13(:,i)! |hat{S}|*hat{S13}
365          fresli(:,21) =fresli(:,21)+shp(i,intp)*SS23(:,i)! |hat{S}|*hat{S23}
404       subroutine hfilterCC (y, x, ien, hfres, shgl, shp, Qwtf)  argument
419      &          shp(nshl,ngauss),
516             u1 = u1 + shp(i,intp)*yl(:,i,2)
517             u2 = u2 + shp(i,intp)*yl(:,i,3)
518             u3 = u3 + shp(i,intp)*yl(:,i,4)
594      &     hfres, shgl, shp, Qwtf)  argument
606      &          shp(nshl,ngauss),
717          fresli(:,1) = fresli(:,1) + shp(i,intp)*hfresl(:,i,1)  ! hat{u1}
718          fresli(:,2) = fresli(:,2) + shp(i,intp)*hfresl(:,i,2)  ! hat{u2}
719          fresli(:,3) = fresli(:,3) + shp(i,intp)*hfresl(:,i,3)  ! hat{u3}
721          fresli(:,4) = fresli(:,4) + shp(i,intp)*hfresl(:,i,4) ! hat{u1*u1}
722          fresli(:,5) = fresli(:,5) + shp(i,intp)*hfresl(:,i,5) ! hat{u2*u2}
723          fresli(:,6) = fresli(:,6) + shp(i,intp)*hfresl(:,i,6) ! hat{u3*u3}
724          fresli(:,7) = fresli(:,7) + shp(i,intp)*hfresl(:,i,7) ! hat{u1*u2}
725          fresli(:,8) = fresli(:,8) + shp(i,intp)*hfresl(:,i,8) ! hat{u1*u3}
726          fresli(:,9) = fresli(:,9) + shp(i,intp)*hfresl(:,i,9) ! hat{u2*u3}
728          fresli(:,10) = fresli(:,10) + shp(i,intp)*hfresl(:,i,10)  ! hat{S11}
729          fresli(:,11) = fresli(:,11) + shp(i,intp)*hfresl(:,i,11)  ! hat{S22}
730          fresli(:,12) = fresli(:,12) + shp(i,intp)*hfresl(:,i,12)  ! hat{S33}
731          fresli(:,13) = fresli(:,13) + shp(i,intp)*hfresl(:,i,13)  ! hat{S12}
732          fresli(:,14) = fresli(:,14) + shp(i,intp)*hfresl(:,i,14)  ! hat{S13}
733          fresli(:,15) = fresli(:,15) + shp(i,intp)*hfresl(:,i,15)  ! hat{S23}
735          fresli(:,16) = fresli(:,16) +shp(i,intp)*hfresl(:,i,17)! hat{S11*|S|}
736          fresli(:,17) = fresli(:,17) +shp(i,intp)*hfresl(:,i,18)! hat{S22*|S|}
737          fresli(:,18) = fresli(:,18) +shp(i,intp)*hfresl(:,i,19)! hat{S33*|S|}
738          fresli(:,19) = fresli(:,19) +shp(i,intp)*hfresl(:,i,20)! hat{S12*|S|}
739          fresli(:,20) = fresli(:,20) +shp(i,intp)*hfresl(:,i,21)! hat{S13*|S|}
740          fresli(:,21) = fresli(:,21) +shp(i,intp)*hfresl(:,i,22)! hat{S23*|S|}
782       subroutine hfilterBB (y, x, ien, hfres, shgl, shp, Qwtf)  argument
797      &          shp(nshl,ngauss),
895             u1 = u1 + shp(i,intp)*yl(:,i,2)
896             u2 = u2 + shp(i,intp)*yl(:,i,3)
897             u3 = u3 + shp(i,intp)*yl(:,i,4)
926             fresli(:,i,1) = WdetJ * u1 * shp(i,intp) !G * u1 * WdetJ
927             fresli(:,i,2) = WdetJ * u2 * shp(i,intp) !G * u2 * WdetJ
928             fresli(:,i,3) = WdetJ * u3 * shp(i,intp) !G * u3 * WdetJ
930             fresli(:,i,4) = WdetJ * u1 * u1 * shp(i,intp) ! G*u1*u1*WdetJ
931             fresli(:,i,5) = WdetJ * u2 * u2 * shp(i,intp) ! G*u2*u2*WdetJ
932             fresli(:,i,6) = WdetJ * u3 * u3 * shp(i,intp) ! G*u3*u3*WdetJ
933             fresli(:,i,7) = WdetJ * u1 * u2 * shp(i,intp) ! G*u1*u2*WdetJ
934             fresli(:,i,8) = WdetJ * u1 * u3 * shp(i,intp) ! G*u1*u3*WdetJ
935             fresli(:,i,9) = WdetJ * u2 * u3 * shp(i,intp) ! G*u2*u3*WdetJ
937             fresli(:,i,10) = S11 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
938             fresli(:,i,11) = S22 * shp(i,intp) * WdetJ ! G*S_{2,2}*WdetJ
939             fresli(:,i,12) = S33 * shp(i,intp) * WdetJ ! G*S_{3,3}*WdetJ
940             fresli(:,i,13) = S12 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
941             fresli(:,i,14) = S13 * shp(i,intp) * WdetJ ! G*S_{1,3}*WdetJ
942             fresli(:,i,15) = S23 * shp(i,intp) * WdetJ ! G*S_{2,3}*WdetJ
944             fresli(:,i,22) = WdetJ*shp(i,intp) !Integral of filter kernel, G,
949             fresli(:,i,16) = S11 * strnrmi * shp(i,intp) * WdetJ
950             fresli(:,i,17) = S22 * strnrmi * shp(i,intp) * WdetJ
951             fresli(:,i,18) = S33 * strnrmi * shp(i,intp) * WdetJ
952             fresli(:,i,19) = S12 * strnrmi * shp(i,intp) * WdetJ
953             fresli(:,i,20) = S13 * strnrmi * shp(i,intp) * WdetJ
954             fresli(:,i,21) = S23 * strnrmi * shp(i,intp) * WdetJ
972      &                  shp,         iper,       ilwork,      argument
999      &          shgl(MAXTOP,nsd,maxsh,MAXQPT), shp(MAXTOP,maxsh,MAXQPT),
1053      &                   shp(lcsyst,1:nshl,:),
1090      &                   shgl(lcsyst,:,1:nshl,:), shp(lcsyst,1:nshl,:),
1236       subroutine resSij (y, x, ien, hfres, shgl, shp, Qwtf)  argument
1251      &          shp(nshl,ngauss),
1349             u1 = u1 + shp(i,intp)*yl(:,i,2)
1350             u2 = u2 + shp(i,intp)*yl(:,i,3)
1351             u3 = u3 + shp(i,intp)*yl(:,i,4)
1380             fresli(:,i,1) = WdetJ * u1 * shp(i,intp) !G * u1 * WdetJ
1381             fresli(:,i,2) = WdetJ * u2 * shp(i,intp) !G * u2 * WdetJ
1382             fresli(:,i,3) = WdetJ * u3 * shp(i,intp) !G * u3 * WdetJ
1384             fresli(:,i,4) = WdetJ * u1 * u1 * shp(i,intp) ! G*u1*u1*WdetJ
1385             fresli(:,i,5) = WdetJ * u2 * u2 * shp(i,intp) ! G*u2*u2*WdetJ
1386             fresli(:,i,6) = WdetJ * u3 * u3 * shp(i,intp) ! G*u3*u3*WdetJ
1387             fresli(:,i,7) = WdetJ * u1 * u2 * shp(i,intp) ! G*u1*u2*WdetJ
1388             fresli(:,i,8) = WdetJ * u1 * u3 * shp(i,intp) ! G*u1*u3*WdetJ
1389             fresli(:,i,9) = WdetJ * u2 * u3 * shp(i,intp) ! G*u2*u3*WdetJ
1391             fresli(:,i,10) = S11 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
1392             fresli(:,i,11) = S22 * shp(i,intp) * WdetJ ! G*S_{2,2}*WdetJ
1393             fresli(:,i,12) = S33 * shp(i,intp) * WdetJ ! G*S_{3,3}*WdetJ
1394             fresli(:,i,13) = S12 * shp(i,intp) * WdetJ ! G*S_{1,1}*WdetJ
1395             fresli(:,i,14) = S13 * shp(i,intp) * WdetJ ! G*S_{1,3}*WdetJ
1396             fresli(:,i,15) = S23 * shp(i,intp) * WdetJ ! G*S_{2,3}*WdetJ
1399             fresli(:,i,16) = strnrmi*strnrmi*strnrmi*shp(i,intp)*WdetJ
1401             fresli(:,i,22) = WdetJ*shp(i,intp) !Integral of filter kernel, G,
1407             fresli(:,i,17) = S22 * strnrmi * shp(i,intp) * WdetJ
1408             fresli(:,i,18) = S33 * strnrmi * shp(i,intp) * WdetJ
1409             fresli(:,i,19) = S12 * strnrmi * shp(i,intp) * WdetJ
1410             fresli(:,i,20) = S13 * strnrmi * shp(i,intp) * WdetJ
1411             fresli(:,i,21) = S23 * strnrmi * shp(i,intp) * WdetJ
1639       subroutine stdfdmc (y,      shgl,      shp,   argument
1678      &          shgl(MAXTOP,nsd,maxsh,MAXQPT), shp(MAXTOP,maxsh,MAXQPT),
1767      &               shp(lcsyst,1:nshl,:))
2237      &        mxmudmi(iblk)%p,cdelsq,shp(lcsyst,1:nshl,:))
2326       subroutine widefdmc (y,      shgl,      shp,   argument
2363       dimension shgl(MAXTOP,nsd,maxsh,MAXQPT), shp(MAXTOP,maxsh,MAXQPT),
2542      &               shp(lcsyst,1:nshl,:))
2989      &        mxmudmi(iblk)%p,cdelsq,shp(lcsyst,1:nshl,:))