18   PetscScalar       *x, *xk, *xj, *xk_tmp, *t;  in MatSolve_SeqSBAIJ_N_inplace()  local
30   xk = t;  in MatSolve_SeqSBAIJ_N_inplace()
33     for (k1 = 0; k1 < bs; k1++) *xk++ = b[idx + k1];  in MatSolve_SeqSBAIJ_N_inplace()
37     xk = t + k * bs;                          /* Dk*xk = k-th block of x */  in MatSolve_SeqSBAIJ_N_inplace()
38     PetscCall(PetscArraycpy(xk_tmp, xk, bs)); /* xk_tmp <- xk */  in MatSolve_SeqSBAIJ_N_inplace()
51     PetscKernel_w_gets_A_times_v(bs, xk_tmp, diag, xk); /* xk <- diag * xk */  in MatSolve_SeqSBAIJ_N_inplace()
57     xk = t + k * bs; /* xk */  in MatSolve_SeqSBAIJ_N_inplace()
63       PetscKernel_v_gets_v_plus_A_times_w(bs, xk, v, xj); /* xk <- xk + v*xj */  in MatSolve_SeqSBAIJ_N_inplace()
69     for (k1 = 0; k1 < bs; k1++) x[idx + k1] = *xk++;  in MatSolve_SeqSBAIJ_N_inplace()
97   PetscScalar     *xk, *xj, *xk_tmp;  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()  local
104     xk = x + k * bs;                          /* Dk*xk = k-th block of x */  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
105     PetscCall(PetscArraycpy(xk_tmp, xk, bs)); /* xk_tmp <- xk */  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
118     PetscKernel_w_gets_A_times_v(bs, xk_tmp, diag, xk); /* xk <- diag * xk */  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
129   PetscScalar     *xk, *xj;  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()  local
134     xk = x + k * bs; /* xk */  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
140       PetscKernel_v_gets_v_plus_A_times_w(bs, xk, v, xj); /* xk <- xk + v*xj */  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
1681   PetscScalar       *x, xk, *t;  in MatSolve_SeqSBAIJ_1()  local
1695     xk = t[k];  in MatSolve_SeqSBAIJ_1()
1697     for (j = 0; j < nz; j++) t[vj[j]] += v[j] * xk;  in MatSolve_SeqSBAIJ_1()
1698     t[k] = xk * v[nz]; /* v[nz] = 1/D(k) */  in MatSolve_SeqSBAIJ_1()
1723   PetscScalar       *x, xk, *t;  in MatSolve_SeqSBAIJ_1_inplace()  local
1738     xk = t[k];  in MatSolve_SeqSBAIJ_1_inplace()
1740     while (nz--) t[*vj++] += (*v++) * xk;  in MatSolve_SeqSBAIJ_1_inplace()
1741     t[k] = xk * aa[ai[k]]; /* aa[k] = 1/D(k) */  in MatSolve_SeqSBAIJ_1_inplace()
1767   PetscScalar       *x, xk;  in MatForwardSolve_SeqSBAIJ_1()  local
1781     xk = x[k];  in MatForwardSolve_SeqSBAIJ_1()
1783     while (nz--) x[*vj++] += (*v++) * xk;  in MatForwardSolve_SeqSBAIJ_1()
1787     x[k] = xk * PetscSqrtReal(diagk);  in MatForwardSolve_SeqSBAIJ_1()
1803   PetscScalar       *x, xk;  in MatForwardSolve_SeqSBAIJ_1_inplace()  local
1817     xk = x[k];  in MatForwardSolve_SeqSBAIJ_1_inplace()
1819     while (nz--) x[*vj++] += (*v++) * xk;  in MatForwardSolve_SeqSBAIJ_1_inplace()
1823     x[k] = xk * PetscSqrtReal(diagk);  in MatForwardSolve_SeqSBAIJ_1_inplace()
2123   PetscScalar       *x, xk;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2136     xk = x[k];  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2138     while (nz--) x[*vj++] += (*v++) * xk;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2139     x[k] = xk * aa[ai[k]]; /* note: aa[diag[k]] = 1/D(k) */  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2146     xk = x[k];  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2148     while (nz--) xk += (*v++) * x[*vj++];  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2149     x[k] = xk;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()