42      2   K, U, KU)  argument
51       REAL*8, INTENT(OUT) :: KU(dof,nNo)  local
55       KU = 0D0
60                KU(1,i) = KU(1,i) + K(1,j)*U(1,colPtr(j))
67                KU(1,i) = KU(1,i) + K(1,j)*U(1,col) + K(2,j)*U(2,col)
68                KU(2,i) = KU(2,i) + K(3,j)*U(1,col) + K(4,j)*U(2,col)
75                KU(1,i) = KU(1,i) + K(1,j)*U(1,col) + K(2,j)*U(2,col)
77                KU(2,i) = KU(2,i) + K(4,j)*U(1,col) + K(5,j)*U(2,col)
79                KU(3,i) = KU(3,i) + K(7,j)*U(1,col) + K(8,j)*U(2,col)
87                KU(1,i) = KU(1,i) + K(1 ,j)*U(1,col) + K(2 ,j)*U(2,col)
89                KU(2,i) = KU(2,i) + K(5 ,j)*U(1,col) + K(6 ,j)*U(2,col)
91                KU(3,i) = KU(3,i) + K(9 ,j)*U(1,col) + K(10,j)*U(2,col)
93                KU(4,i) = KU(4,i) + K(13,j)*U(1,col) + K(14,j)*U(2,col)
104                   KU(l,i) = KU(l,i) + SUM(K(s:e,j)*U(:,col))
110       CALL COMMUV(dof, nNo, commu, cS, KU)
118      2   K, U, KU)  argument
127       REAL*8, INTENT(OUT) :: KU(nNo)  local
131       KU = 0D0
136                KU(i) = KU(i) + K(1,j)*U(1,colPtr(j))
143                KU(i) = KU(i) + K(1,j)*U(1,col) + K(2,j)*U(2,col)
150                KU(i) = KU(i) + K(1,j)*U(1,col) + K(2,j)*U(2,col)
158                KU(i) = KU(i) + K(1,j)*U(1,col) + K(2,j)*U(2,col)
165                KU(i) = KU(i) + SUM(K(:,j)*U(:,colPtr(j)))
170       CALL COMMUS(nNo, commu, cS, KU)
178      2   K, U, KU)  argument
187       REAL*8, INTENT(OUT) :: KU(dof,nNo)  local
191       KU = 0D0
196                KU(1,i) = KU(1,i) + K(1,j)*U(colPtr(j))
203                KU(1,i) = KU(1,i) + K(1,j)*U(col)
204                KU(2,i) = KU(2,i) + K(2,j)*U(col)
211                KU(1,i) = KU(1,i) + K(1,j)*U(col)
212                KU(2,i) = KU(2,i) + K(2,j)*U(col)
213                KU(3,i) = KU(3,i) + K(3,j)*U(col)
220                KU(1,i) = KU(1,i) + K(1,j)*U(col)
221                KU(2,i) = KU(2,i) + K(2,j)*U(col)
222                KU(3,i) = KU(3,i) + K(3,j)*U(col)
223                KU(4,i) = KU(4,i) + K(4,j)*U(col)
229                KU(:,i) = KU(:,i) + K(:,j)*U(colPtr(j))
234       CALL COMMUV(dof, nNo, commu, cS, KU)
242      2   K, U, KU)  argument
251       REAL*8, INTENT(OUT) :: KU(nNo)  local
255       KU = 0D0
258             KU(i) = KU(i) + K(j)*U(colPtr(j))
262       CALL COMMUS(nNo, commu, cS, KU)