40 A.rowComputeNonZeroElements();
45 int i, j, numRows = A.rows();
49 for (
int k = 0; k < numIterations; k++)
52 for (i = 0; i < numRows; i++)
57 for (
int h = 0; h < A.m_rowNonZeroElements1[i].size(); h++)
59 j = A.m_rowNonZeroElements1[i][h];
62 delta += A(i, j) * x[j];
68 for (j = 0; j < i; j++)
69 delta += A(i, j) * x[j];
70 for (j = i + 1; j < numRows; j++)
71 delta += A(i, j) * x[j];
76 x[i] = (b[i] - delta) / aDiag;
79 if (limitDependency[i] >= 0)
81 s = x[limitDependency[i]];
97#ifdef VERBOSE_PRINTF_RESIDUAL