13   PetscInt           mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatSolve_SeqSBAIJ_N_inplace()  local
36     v  = aa + bs2 * ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
39     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
40     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
56     v  = aa + bs2 * ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
58     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
59     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
92 static PetscErrorCode MatForwardSolve_SeqSBAIJ_N_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()  argument
103     v  = aa + bs2 * ai[k];  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
106     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
107     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
124 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()  argument
133     v  = aa + bs2 * ai[k];  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
135     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
136     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
152   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
178   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
198   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
219   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *r, *vj;  in MatSolve_SeqSBAIJ_7_inplace()  local
246     v  = aa + 49 * ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
247     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
256     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
285     v  = aa + 49 * ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
286     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
295     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
336 static PetscErrorCode MatForwardSolve_SeqSBAIJ_7_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()  argument
345     v  = aa + 49 * ai[k];  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
354     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
355     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
385 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()  argument
394     v  = aa + 49 * ai[k];  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
403     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
404     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
435   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
460   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()  local
469   PetscCall(MatForwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
479   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()  local
488   PetscCall(MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
499   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *r, *vj;  in MatSolve_SeqSBAIJ_6_inplace()  local
525     v  = aa + 36 * ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
526     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
534     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
561     v  = aa + 36 * ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
562     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
570     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
608 static PetscErrorCode MatForwardSolve_SeqSBAIJ_6_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()  argument
617     v  = aa + 36 * ai[k];  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
625     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
626     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
653 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()  argument
662     v  = aa + 36 * ai[k];  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
670     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
671     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
700   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
725   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()  local
734   PetscCall(MatForwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
744   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()  local
753   PetscCall(MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
764   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatSolve_SeqSBAIJ_5_inplace()  local
790     v  = aa + 25 * ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
791     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
798     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
824     v  = aa + 25 * ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
825     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
832     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
867 static PetscErrorCode MatForwardSolve_SeqSBAIJ_5_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()  argument
876     v  = aa + 25 * ai[k];  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
883     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
884     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
910 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()  argument
919     v  = aa + 25 * ai[k];  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
926     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
927     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
954   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
979   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()  local
988   PetscCall(MatForwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
998   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()  local
1007   PetscCall(MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1018   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatSolve_SeqSBAIJ_4_inplace()  local
1043     v  = aa + 16 * ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1044     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1050     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1074     v  = aa + 16 * ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1075     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1081     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1113 static PetscErrorCode MatForwardSolve_SeqSBAIJ_4_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()  argument
1122     v  = aa + 16 * ai[k];  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1128     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1129     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1153 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()  argument
1162     v  = aa + 16 * ai[k];  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1168     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1169     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1194   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
1218   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()  local
1227   PetscCall(MatForwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1237   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()  local
1246   PetscCall(MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1257   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatSolve_SeqSBAIJ_3_inplace()  local
1282     v  = aa + 9 * ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1283     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1288     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1310     v  = aa + 9 * ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1311     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1316     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1345 static PetscErrorCode MatForwardSolve_SeqSBAIJ_3_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()  argument
1354     v  = aa + 9 * ai[k];  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1359     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1360     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1382 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()  argument
1391     v  = aa + 9 * ai[k];  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1396     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1397     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1420   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
1445   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()  local
1454   PetscCall(MatForwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1464   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()  local
1473   PetscCall(MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1484   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatSolve_SeqSBAIJ_2_inplace()  local
1504     v  = aa + 4 * ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1505     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1509     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1523     v  = aa + 4 * ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1524     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1528     nz = ai[k + 1] - ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1549 static PetscErrorCode MatForwardSolve_SeqSBAIJ_2_NaturalOrdering(const PetscInt *ai, const PetscInt…  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()  argument
1558     v  = aa + 4 * ai[k];  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1559     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1563     nz = ai[k + 1] - ai[k];  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1581 static PetscErrorCode MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering(const PetscInt *ai, const PetscIn…  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()  argument
1590     v  = aa + 4 * ai[k];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1591     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1595     nz = ai[k + 1] - ai[k];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1614   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  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()
1639   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()  local
1648   PetscCall(MatForwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x));  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1658   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()  local
1667   PetscCall(MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x));  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1678   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj, *adiag = a->diag;  in MatSolve_SeqSBAIJ_1()  local
1693     v  = aa + ai[k];  in MatSolve_SeqSBAIJ_1()
1694     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_1()
1696     nz = ai[k + 1] - ai[k] - 1;  in MatSolve_SeqSBAIJ_1()
1705     nz = ai[k + 1] - ai[k] - 1;  in MatSolve_SeqSBAIJ_1()
1721   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj;  in MatSolve_SeqSBAIJ_1_inplace()  local
1736     v  = aa + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1737     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1739     nz = ai[k + 1] - 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()
1747     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1748     nz = ai[k + 1] - ai[k] - 1;  in MatSolve_SeqSBAIJ_1_inplace()
1764   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj, *adiag = a->diag;  in MatForwardSolve_SeqSBAIJ_1()  local
1779     v  = aa + ai[k];  in MatForwardSolve_SeqSBAIJ_1()
1780     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_1()
1782     nz = ai[k + 1] - ai[k] - 1;  in MatForwardSolve_SeqSBAIJ_1()
1800   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj;  in MatForwardSolve_SeqSBAIJ_1_inplace()  local
1815     v  = aa + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_inplace()
1816     vj = aj + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_inplace()
1818     nz = ai[k + 1] - 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()
1836   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj, *adiag = a->diag;  in MatBackwardSolve_SeqSBAIJ_1()  local
1851     v     = aa + ai[k];  in MatBackwardSolve_SeqSBAIJ_1()
1852     vj    = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_1()
1856     nz   = ai[k + 1] - ai[k] - 1;  in MatBackwardSolve_SeqSBAIJ_1()
1871   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj;  in MatBackwardSolve_SeqSBAIJ_1_inplace()  local
1886     v     = aa + ai[k] + 1;  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1887     vj    = aj + 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()
1891     nz   = ai[k + 1] - ai[k] - 1;  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1911     const PetscInt    *vj, mbs = a->mbs, *ai = a->i, *aj = a->j, *rp;  in MatSolves_SeqSBAIJ_1()  local
1934       v  = aa + ai[k];  in MatSolves_SeqSBAIJ_1()
1935       vj = aj + ai[k];  in MatSolves_SeqSBAIJ_1()
1936       nz = ai[k + 1] - ai[k] - 1;  in MatSolves_SeqSBAIJ_1()
1947       v  = aa + ai[k] - 1;  in MatSolves_SeqSBAIJ_1()
1948       vj = aj + ai[k] - 1;  in MatSolves_SeqSBAIJ_1()
1949       nz = ai[k + 1] - ai[k] - 1;  in MatSolves_SeqSBAIJ_1()
1975     const PetscInt    *vj, mbs = a->mbs, *ai = a->i, *aj = a->j, *rp;  in MatSolves_SeqSBAIJ_1_inplace()  local
1998       v  = aa + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
1999       vj = aj + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2000       nz = ai[k + 1] - ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2011       v  = aa + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2012       vj = aj + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2013       nz = ai[k + 1] - ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2033   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *vj, *adiag = a->diag;  in MatSolve_SeqSBAIJ_1_NaturalOrdering()  local
2045     v  = aa + ai[i];  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2046     vj = aj + ai[i];  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2048     nz = ai[i + 1] - ai[i] - 1; /* exclude diag[i] */  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2057     nz = ai[i + 1] - ai[i] - 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2070   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *vj, *adiag = a->diag;  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()  local
2093       v  = aa + ai[i];  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2094       vj = aj + ai[i];  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2096       nz = ai[i + 1] - ai[i] - 1; /* exclude diag[i] */  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2105       nz = ai[i + 1] - ai[i] - 1;  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2121   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *vj;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2134     v  = aa + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2135     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2137     nz = ai[k + 1] - ai[k] - 1; /* exclude diag[k] */  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()
2145     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2147     nz = ai[k + 1] - ai[k] - 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2161   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *adiag = a->diag, *vj;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()  local
2174     v  = aa + ai[k];  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2175     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2176     nz = ai[k + 1] - ai[k] - 1; /* exclude diag[k] */  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2191   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *vj;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2204     v  = aa + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2205     vj = aj + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2206     nz = ai[k + 1] - ai[k] - 1; /* exclude diag[k] */  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 …ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal (%g,%g) must be real and nonnegat…  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2221   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *adiag = a->diag, *vj;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()  local
2234     v     = aa + ai[k];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2235     vj    = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2239     nz   = ai[k + 1] - ai[k] - 1;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2251   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *vj;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2264     v     = aa + ai[k] + 1;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2265     vj    = aj + 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()
2269     nz   = ai[k + 1] - ai[k] - 1;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2282   const PetscInt *rip, mbs    = a->mbs, *ai, *aj;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()  local
2296   ai         = a->i;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2301   umax = (PetscInt)(f * ai[mbs] + 1);  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2323     nz = ai[rip[k] + 1] - ai[rip[k]];  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2324     j  = ai[rip[k]];  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2411   if (ai[mbs] != 0) {  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2412     PetscReal af = ((PetscReal)iu[mbs]) / ((PetscReal)ai[mbs]);  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2456   if (ai[mbs] != 0) {  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2457     B->info.fill_ratio_needed = ((PetscReal)iu[mbs]) / ((PetscReal)ai[mbs]);  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2474   PetscInt           bs = A->rmap->bs, am = a->mbs, *ai = a->i, *aj = a->j;  in MatICCFactorSymbolic_SeqSBAIJ()  local
2506       ncols     = ai[i + 1] - ai[i];  in MatICCFactorSymbolic_SeqSBAIJ()
2513       aj    = a->j + ai[i] + 1; /* 1st entry of U(i,:) without diagonal */  in MatICCFactorSymbolic_SeqSBAIJ()
2514       ncols = ai[i + 1] - ai[i] - 1;  in MatICCFactorSymbolic_SeqSBAIJ()
2535     PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space));  in MatICCFactorSymbolic_SeqSBAIJ()
2539     PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space_lvl));  in MatICCFactorSymbolic_SeqSBAIJ()
2546       ncols = ai[k + 1] - ai[k];  in MatICCFactorSymbolic_SeqSBAIJ()
2548       cols = aj + ai[k];  in MatICCFactorSymbolic_SeqSBAIJ()
2652   if (ai[am] != 0) {  in MatICCFactorSymbolic_SeqSBAIJ()
2653     fact->info.fill_ratio_needed = ((PetscReal)ui[am]) / ai[am];  in MatICCFactorSymbolic_SeqSBAIJ()
2658   if (ai[am] != 0) {  in MatICCFactorSymbolic_SeqSBAIJ()
2677   const PetscInt    *cols, *rip, *ai = a->i, *aj = a->j;  in MatICCFactorSymbolic_SeqSBAIJ_inplace()  local
2733     PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space));  in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2737     PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space_lvl));  in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2744       ncols = ai[rip[k] + 1] - ai[rip[k]];  in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2745       cols  = aj + ai[rip[k]];  in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2816     if (ai[am] != 0) {  in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2817       ratio_needed = ((PetscReal)ui[am]) / ((PetscReal)ai[am]);  in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2850   if (ai[am] != 0) {  in MatICCFactorSymbolic_SeqSBAIJ_inplace()