17   const MatScalar   *aa = a->a, *v, *diag;  in MatSolve_SeqSBAIJ_N_inplace()  local
36     v  = aa + bs2 * ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
50     diag = aa + k * bs2;                                /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_N_inplace()
56     v  = aa + bs2 * ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
92 …N_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()  argument
103     v  = aa + bs2 * ai[k];  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
117     diag = aa + k * bs2;                                /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
124 …N_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()  argument
133     v  = aa + bs2 * ai[k];  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
154   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()  local
164   PetscCall(MatForwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x));  in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
167   PetscCall(MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x));  in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
180   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()  local
188   PetscCall(MatForwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x));  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
200   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()  local
208   PetscCall(MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x));  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
221   const MatScalar   *aa = a->a, *v, *d;  in MatSolve_SeqSBAIJ_7_inplace()  local
246     v  = aa + 49 * ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
272     d     = aa + k * 49; /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_7_inplace()
285     v  = aa + 49 * ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
336 …7_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()  argument
345     v  = aa + 49 * ai[k];  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
372     d     = aa + k * 49; /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
385 …7_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()  argument
394     v  = aa + 49 * ai[k];  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
436   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_7_NaturalOrdering_inplace()  local
446   PetscCall(MatForwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
449   PetscCall(MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
461   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()  local
469   PetscCall(MatForwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
480   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()  local
488   PetscCall(MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
501   const MatScalar   *aa = a->a, *v, *d;  in MatSolve_SeqSBAIJ_6_inplace()  local
525     v  = aa + 36 * ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
549     d     = aa + k * 36; /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_6_inplace()
561     v  = aa + 36 * ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
608 …6_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()  argument
617     v  = aa + 36 * ai[k];  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
642     d     = aa + k * 36; /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
653 …6_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()  argument
662     v  = aa + 36 * ai[k];  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
701   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_6_NaturalOrdering_inplace()  local
711   PetscCall(MatForwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
714   PetscCall(MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
726   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()  local
734   PetscCall(MatForwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
745   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()  local
753   PetscCall(MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
767   const MatScalar   *aa = a->a, *v, *diag;  in MatSolve_SeqSBAIJ_5_inplace()  local
790     v  = aa + 25 * ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
813     diag  = aa + k * 25; /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_5_inplace()
824     v  = aa + 25 * ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
867 …5_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()  argument
876     v  = aa + 25 * ai[k];  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
899     diag  = aa + k * 25; /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
910 …5_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()  argument
919     v  = aa + 25 * ai[k];  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
955   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_5_NaturalOrdering_inplace()  local
965   PetscCall(MatForwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
968   PetscCall(MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
980   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()  local
988   PetscCall(MatForwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
999   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()  local
1007   PetscCall(MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1021   const MatScalar   *aa = a->a, *v, *diag;  in MatSolve_SeqSBAIJ_4_inplace()  local
1043     v  = aa + 16 * ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1064     diag  = aa + k * 16; /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_4_inplace()
1074     v  = aa + 16 * ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1113 …4_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()  argument
1122     v  = aa + 16 * ai[k];  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1143     diag  = aa + k * 16; /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1153 …4_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()  argument
1162     v  = aa + 16 * ai[k];  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1195   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_4_NaturalOrdering_inplace()  local
1205   PetscCall(MatForwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1208   PetscCall(MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1219   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()  local
1227   PetscCall(MatForwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1238   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()  local
1246   PetscCall(MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1261   const MatScalar   *aa = a->a, *v, *diag;  in MatSolve_SeqSBAIJ_3_inplace()  local
1282     v  = aa + 9 * ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1301     diag  = aa + k * 9; /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_3_inplace()
1310     v  = aa + 9 * ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1345 …3_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()  argument
1354     v  = aa + 9 * ai[k];  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1373     diag  = aa + k * 9; /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1382 …3_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()  argument
1391     v  = aa + 9 * ai[k];  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1421   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_3_NaturalOrdering_inplace()  local
1431   PetscCall(MatForwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1434   PetscCall(MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1446   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()  local
1454   PetscCall(MatForwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1465   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()  local
1473   PetscCall(MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1487   const MatScalar   *aa = a->a, *v, *diag;  in MatSolve_SeqSBAIJ_2_inplace()  local
1504     v  = aa + 4 * ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1516     diag      = aa + k * 4; /* ptr to inv(Dk) */  in MatSolve_SeqSBAIJ_2_inplace()
1523     v  = aa + 4 * ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1549 …2_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()  argument
1558     v  = aa + 4 * ai[k];  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1574     diag      = aa + k * 4; /* ptr to inv(Dk) */  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1581 …2_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, Petsc…  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()  argument
1590     v  = aa + 4 * ai[k];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1615   const MatScalar   *aa = a->a;  in MatSolve_SeqSBAIJ_2_NaturalOrdering_inplace()  local
1625   PetscCall(MatForwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1628   PetscCall(MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x));  in MatSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1640   const MatScalar   *aa = a->a;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()  local
1648   PetscCall(MatForwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1659   const MatScalar   *aa = a->a;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()  local
1667   PetscCall(MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1679   const MatScalar   *aa = a->a, *v;  in MatSolve_SeqSBAIJ_1()  local
1693     v  = aa + ai[k];  in MatSolve_SeqSBAIJ_1()
1703     v  = aa + adiag[k] - 1;  in MatSolve_SeqSBAIJ_1()
1722   const MatScalar   *aa = a->a, *v;  in MatSolve_SeqSBAIJ_1_inplace()  local
1736     v  = aa + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1741     t[k] = xk * aa[ai[k]]; /* aa[k] = 1/D(k) */  in MatSolve_SeqSBAIJ_1_inplace()
1746     v  = aa + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1765   const MatScalar   *aa = a->a, *v;  in MatForwardSolve_SeqSBAIJ_1()  local
1779     v  = aa + ai[k];  in MatForwardSolve_SeqSBAIJ_1()
1785     diagk = PetscRealPart(aa[adiag[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1()
1786 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatForwardSolve_SeqSBAIJ_1()
1801   const MatScalar   *aa = a->a, *v;  in MatForwardSolve_SeqSBAIJ_1_inplace()  local
1815     v  = aa + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_inplace()
1821     diagk = PetscRealPart(aa[ai[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1_inplace()
1822 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatForwardSolve_SeqSBAIJ_1_inplace()
1837   const MatScalar   *aa = a->a, *v;  in MatBackwardSolve_SeqSBAIJ_1()  local
1851     v     = aa + ai[k];  in MatBackwardSolve_SeqSBAIJ_1()
1853     diagk = PetscRealPart(aa[adiag[k]]);  in MatBackwardSolve_SeqSBAIJ_1()
1854 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatBackwardSolve_SeqSBAIJ_1()
1872   const MatScalar   *aa = a->a, *v;  in MatBackwardSolve_SeqSBAIJ_1_inplace()  local
1886     v     = aa + ai[k] + 1;  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1888     diagk = PetscRealPart(aa[ai[k]]);  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1889 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1912     const MatScalar   *aa = a->a, *v;  in MatSolves_SeqSBAIJ_1()  local
1934       v  = aa + ai[k];  in MatSolves_SeqSBAIJ_1()
1942       for (i = 0; i < n; i++) t[n * k + i] *= aa[nz]; /* note: aa[nz] = 1/D(k) */  in MatSolves_SeqSBAIJ_1()
1947       v  = aa + ai[k] - 1;  in MatSolves_SeqSBAIJ_1()
1976     const MatScalar   *aa = a->a, *v;  in MatSolves_SeqSBAIJ_1_inplace()  local
1998       v  = aa + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2006       for (i = 0; i < n; i++) t[n * k + i] *= aa[k]; /* note: aa[k] = 1/D(k) */  in MatSolves_SeqSBAIJ_1_inplace()
2011       v  = aa + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2034   const MatScalar   *aa = a->a, *v;  in MatSolve_SeqSBAIJ_1_NaturalOrdering()  local
2045     v  = aa + ai[i];  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2055     v  = aa + adiag[i] - 1; /* end of row i, excluding diag */  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2071   const MatScalar   *aa = a->a, *v;  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()  local
2093       v  = aa + ai[i];  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2103       v  = aa + adiag[i] - 1; /* end of row i, excluding diag */  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2122   const MatScalar   *aa = a->a, *v;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2134     v  = aa + ai[k] + 1;  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()
2144     v  = aa + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2162   const MatScalar   *aa = a->a, *v;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()  local
2174     v  = aa + ai[k];  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2178     diagk = PetscRealPart(aa[adiag[k]]); /* note: aa[adiag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2179 …aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal (%g,%g) must be real and no…  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2192   const MatScalar   *aa = a->a, *v;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2204     v  = aa + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2208     diagk = PetscRealPart(aa[ai[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2209 …aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal (%g,%g) must be real and nonne…  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2222   const MatScalar   *aa = a->a, *v;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()  local
2234     v     = aa + ai[k];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2236     diagk = PetscRealPart(aa[adiag[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2237 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2252   const MatScalar   *aa = a->a, *v;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2264     v     = aa + ai[k] + 1;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2266     diagk = PetscRealPart(aa[ai[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2267 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()