1486   PetscInt           nz, k, k2, idx;  in MatSolve_SeqSBAIJ_2_inplace()  local
1506     k2 = k * 2;  in MatSolve_SeqSBAIJ_2_inplace()
1507     x0 = t[k2];  in MatSolve_SeqSBAIJ_2_inplace()
1508     x1 = t[k2 + 1];  in MatSolve_SeqSBAIJ_2_inplace()
1517     t[k2]     = diag[0] * x0 + diag[2] * x1;  in MatSolve_SeqSBAIJ_2_inplace()
1518     t[k2 + 1] = diag[1] * x0 + diag[3] * x1;  in MatSolve_SeqSBAIJ_2_inplace()
1525     k2 = k * 2;  in MatSolve_SeqSBAIJ_2_inplace()
1526     x0 = t[k2];  in MatSolve_SeqSBAIJ_2_inplace()
1527     x1 = t[k2 + 1];  in MatSolve_SeqSBAIJ_2_inplace()
1535     t[k2]      = x0;  in MatSolve_SeqSBAIJ_2_inplace()
1536     t[k2 + 1]  = x1;  in MatSolve_SeqSBAIJ_2_inplace()
1553   PetscInt         nz, k, k2;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()  local
1560     k2 = k * 2;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1561     x0 = x[k2];  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1562     x1 = x[k2 + 1]; /* Dk*xk = k-th block of x */  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1575     x[k2]     = diag[0] * x0 + diag[2] * x1;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1576     x[k2 + 1] = diag[1] * x0 + diag[3] * x1;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1585   PetscInt         nz, k, k2;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()  local
1592     k2 = k * 2;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1593     x0 = x[k2];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1594     x1 = x[k2 + 1]; /* xk */  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1605     x[k2]     = x0;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1606     x[k2 + 1] = x1;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()