9   #define PetscSparseDensePlusDot_no_function(sum, r, xv, xi, nnz) \  argument
15           sum += *xv++ * r[*xi++]; \
17           sum += *xv++ * r[*xi++]; \
19           sum += *xv++ * r[*xi++]; \
23           sum += xv[0] * r[xi[0]] + xv[1] * r[xi[1]] + xv[2] * r[xi[2]] + xv[3] * r[xi[3]]; \
34   #define PetscSparseDensePlusDot_no_function(sum, r, xv, xi, nnz) \  argument
40         sum += (xv[__i] * r[__i1] + xv[__i + 1] * r[__i2]); \
42       if (nnz & 0x1) sum += xv[__i] * r[xi[__i]]; \
46   #define PetscSparseDensePlusDot_no_function(sum, r, xv, xi, nnz) \  argument
49       for (__i = 0; __i < nnz; __i++) sum += xv[__i] * r[xi[__i]]; \
61   PetscScalar       *z, x1, sum;  in MatMult_SeqSBAIJ_1_ushort()  local
86     sum  = 0.0;  in MatMult_SeqSBAIJ_1_ushort()
90       sum = v[0] * x1; /* diagonal term */  in MatMult_SeqSBAIJ_1_ushort()
100         sum += vj * x[ibt];           /* (strict upper triangular part of A)*x  */  in MatMult_SeqSBAIJ_1_ushort()
107         sum += vj * x[ibt]; /* (strict upper triangular part of A)*x  */  in MatMult_SeqSBAIJ_1_ushort()
110     z[i] += sum;  in MatMult_SeqSBAIJ_1_ushort()
129   PetscScalar       *x, *t, sum;  in MatSOR_SeqSBAIJ_ushort()  local
172       sum = b[i] * d / omega;  in MatSOR_SeqSBAIJ_ushort()
174       PetscSparseDensePlusDot_no_function(sum, b, v, vj, nz);  in MatSOR_SeqSBAIJ_ushort()
176       PetscSparseDensePlusDot(sum, b, v, vj, nz);  in MatSOR_SeqSBAIJ_ushort()
178       x[i] = sum;  in MatSOR_SeqSBAIJ_ushort()
205           sum = b[i];  in MatSOR_SeqSBAIJ_ushort()
209             for (__i = 0; __i < nz; __i++) sum -= v[-__i] * x[vj[-__i]];  in MatSOR_SeqSBAIJ_ushort()
211           x[i] = omega * sum * aidiag[i];  in MatSOR_SeqSBAIJ_ushort()
221           sum = t[i];  in MatSOR_SeqSBAIJ_ushort()
225           PetscSparseDenseMinusDot(sum, x, v, vj, nz);  in MatSOR_SeqSBAIJ_ushort()
226           x[i] = (1 - omega) * x[i] + omega * sum * aidiag[i];  in MatSOR_SeqSBAIJ_ushort()
256         sum = t[i];  in MatSOR_SeqSBAIJ_ushort()
257         while (nz1--) sum -= (*v1++) * x[*vj1++];  in MatSOR_SeqSBAIJ_ushort()
258         x[i] = (1 - omega) * x[i] + omega * sum * aidiag[i];  in MatSOR_SeqSBAIJ_ushort()
286         sum = t[i];  in MatSOR_SeqSBAIJ_ushort()
287         while (nz--) sum -= x[*vj++] * (*v++);  in MatSOR_SeqSBAIJ_ushort()
288         x[i] = (1 - omega) * x[i] + omega * sum * aidiag[i];  in MatSOR_SeqSBAIJ_ushort()