phSolver/incompressible/filters.f

11       dimension ien(npro,nshl),        yl(npro,nshl,5),
25       call local (y,      yl,     ien,    5,  'gather  ')
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)
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)
127      &                       +shg(:,i,1)*yl(:,i,4)
128             S23 = S23 + shg(:,i,3)*yl(:,i,3)
129      &                       +shg(:,i,2)*yl(:,i,4)
184       dimension ien(npro,nshl),        yl(npro,nshl,ndof),
203       call local (y,      yl,     ien,    5,  'gather  ')
410       dimension ien(npro,nshl),        yl(npro,nshl,5),
425       call local (y,      yl,     ien,    5,  'gather  ')
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)
520             S11 = S11 + shg(:,i,1)*yl(:,i,2)
521             S22 = S22 + shg(:,i,2)*yl(:,i,3)
522             S33 = S33 + shg(:,i,3)*yl(:,i,4)
524             S12 = S12 + shg(:,i,2)*yl(:,i,2)
525      &                       +shg(:,i,1)*yl(:,i,3)
526             S13 = S13 + shg(:,i,3)*yl(:,i,2)
527      &                       +shg(:,i,1)*yl(:,i,4)
528             S23 = S23 + shg(:,i,3)*yl(:,i,3)
529      &                       +shg(:,i,2)*yl(:,i,4)
600       dimension ien(npro,nshl),        yl(npro,nshl,ndof),
616       call local (y,      yl,     ien,    5,  'gather  ')
788       dimension ien(npro,nshl),        yl(npro,nshl,5),
804       call local (y,      yl,     ien,    5,  'gather  ')
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)
899             S11 = S11 + shg(:,i,1)*yl(:,i,2)
900             S22 = S22 + shg(:,i,2)*yl(:,i,3)
901             S33 = S33 + shg(:,i,3)*yl(:,i,4)
903             S12 = S12 + shg(:,i,2)*yl(:,i,2)
904      &                       +shg(:,i,1)*yl(:,i,3)
905             S13 = S13 + shg(:,i,3)*yl(:,i,2)
906      &                       +shg(:,i,1)*yl(:,i,4)
907             S23 = S23 + shg(:,i,3)*yl(:,i,3)
908      &                       +shg(:,i,2)*yl(:,i,4)
1242       dimension ien(npro,nshl),        yl(npro,nshl,5),
1258       call local (y,      yl,     ien,    5,  'gather  ')
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)
1353             S11 = S11 + shg(:,i,1)*yl(:,i,2)
1354             S22 = S22 + shg(:,i,2)*yl(:,i,3)
1355             S33 = S33 + shg(:,i,3)*yl(:,i,4)
1357             S12 = S12 + shg(:,i,2)*yl(:,i,2)
1358      &                       +shg(:,i,1)*yl(:,i,3)
1359             S13 = S13 + shg(:,i,3)*yl(:,i,2)
1360      &                       +shg(:,i,1)*yl(:,i,4)
1361             S23 = S23 + shg(:,i,3)*yl(:,i,3)
1362      &                       +shg(:,i,2)*yl(:,i,4)