phSolver/incompressible/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)
52             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
183       dimension x(numnp,nsd),            xl(npro,nenl,nsd)
204       call localx (x,      xl,     ien,    3,  'gather  ')
235             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
236             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
237             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
238             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
239             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
240             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
241             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
242             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
243             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
409       dimension x(numnp,3),            xl(npro,nenl,3)
426       call localx (x,      xl,     ien,    3,  'gather  ')
444             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
445             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
446             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
447             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
448             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
449             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
450             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
451             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
452             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
599       dimension x(numnp,nsd),            xl(npro,nenl,nsd)
617       call localx (x,      xl,     ien,    3,  'gather  ')
648             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
649             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
650             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
651             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
652             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
653             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
654             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
655             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
656             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
787       dimension x(numnp,3),            xl(npro,nenl,3)
805       call localx (x,      xl,     ien,    3,  'gather  ')
823             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
824             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
825             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
826             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
827             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
828             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
829             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
830             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
831             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
1241       dimension x(numnp,3),            xl(npro,nenl,3)
1259       call localx (x,      xl,     ien,    3,  'gather  ')
1277             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
1278             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
1279             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
1280             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
1281             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
1282             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
1283             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
1284             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
1285             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
1668      &          xm(nfath,6),           xl(nfath,6),
1932       xl    = zero
1955          xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
1956          xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
1957          xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
1958          xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
1959          xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
1960          xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)
1992          call drvAllreduce(xl, xlr,6*nfath)
2078          xl(:,m) = xl(:,m)/nsons(:)
2126             xlcomp(:,:) = xlcomp(:,:)+xl(:,:)
2353       dimension xm(nfath,6),           xl(nfath,6),
2673       xl    = zero
2696          xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
2697          xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
2698          xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
2699          xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
2700          xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
2701          xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)
2733          call drvAllreduce(xl, xlr,6*nfath)
2820          xl(:,m) = xl(:,m)/nsons(:)
2867             xlcomp(:,:) = xlcomp(:,:)+xl(:,:)