43     lu_int i, k, ipivot, jpivot, pos;  in lu_residual_test()  local
60         ipivot = p[k];  in lu_residual_test()
61         rhs[ipivot] = d <= 0.0 ? 1.0 : -1.0;  in lu_residual_test()
62         lhs[ipivot] = rhs[ipivot] - d;  in lu_residual_test()
68         ipivot = pivotrow[k];  in lu_residual_test()
69         d = lhs[ipivot] /= row_pivot[ipivot];  in lu_residual_test()
70         for (pos = Ubegin[ipivot]; (i = Uindex[pos]) >= 0; pos++)  in lu_residual_test()
79         ipivot = pivotrow[k];  in lu_residual_test()
81         d = lhs[ipivot];  in lu_residual_test()
89         ipivot = pivotrow[k];  in lu_residual_test()
90         rhs[ipivot] -= lhs[ipivot];  in lu_residual_test()
102         ipivot = pivotrow[k];  in lu_residual_test()
104         for (pos = Ubegin[ipivot]; (i = Uindex[pos]) >= 0; pos++)  in lu_residual_test()
108         rhs[ipivot] = d <= 0.0 ? 1.0 : -1.0;  in lu_residual_test()
109         lhs[ipivot] = (rhs[ipivot] - d) / row_pivot[ipivot];  in lu_residual_test()
126         ipivot = pivotrow[k];  in lu_residual_test()
133         rhs[ipivot] -= d;  in lu_residual_test()
137         ipivot = pivotrow[k];  in lu_residual_test()
138         rhs[ipivot] -= lhs[ipivot];  in lu_residual_test()