15   PetscInt           nz, *vj, k, idx, k1;  in MatSolve_SeqSBAIJ_N_inplace()  local
40     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
41     xj = t + (*vj) * bs; /* *vj-th block of x, *vj>k */  in MatSolve_SeqSBAIJ_N_inplace()
45       vj++;  in MatSolve_SeqSBAIJ_N_inplace()
46       xj = t + (*vj) * bs;  in MatSolve_SeqSBAIJ_N_inplace()
59     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_N_inplace()
60     xj = t + (*vj) * bs;  in MatSolve_SeqSBAIJ_N_inplace()
64       vj++;  in MatSolve_SeqSBAIJ_N_inplace()
66       xj = t + (*vj) * bs;  in MatSolve_SeqSBAIJ_N_inplace()
95   const PetscInt  *vj, bs2 = bs * bs;  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()  local
107     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
108     xj = x + (size_t)(*vj) * bs; /* *vj-th block of x, *vj>k */  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
112       vj++;  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
113       xj = x + (size_t)(*vj) * bs;  in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
127   const PetscInt  *vj, bs2 = bs * bs;  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()  local
136     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
137     xj = x + (size_t)(*vj) * bs;  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
141       vj++;  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
143       xj = x + (size_t)(*vj) * bs;  in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
219   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *r, *vj;  in MatSolve_SeqSBAIJ_7_inplace()  local
247     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
257     tp = t + (*vj) * 7;  in MatSolve_SeqSBAIJ_7_inplace()
266       vj++;  in MatSolve_SeqSBAIJ_7_inplace()
267       tp = t + (*vj) * 7;  in MatSolve_SeqSBAIJ_7_inplace()
286     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_7_inplace()
297     tp = t + (*vj) * 7;  in MatSolve_SeqSBAIJ_7_inplace()
307       vj++;  in MatSolve_SeqSBAIJ_7_inplace()
308       tp = t + (*vj) * 7;  in MatSolve_SeqSBAIJ_7_inplace()
341   const PetscInt  *vj;  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()  local
355     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
356 …PetscPrefetchBlock(vj + nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
359       xp = x + (*vj) * 7;  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
368       vj++;  in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
390   const PetscInt  *vj;  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()  local
404     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
405 …PetscPrefetchBlock(vj - nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
408       xp = x + (*vj) * 7;  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
417       vj++;  in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
499   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *r, *vj;  in MatSolve_SeqSBAIJ_6_inplace()  local
526     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
535     tp = t + (*vj) * 6;  in MatSolve_SeqSBAIJ_6_inplace()
543       vj++;  in MatSolve_SeqSBAIJ_6_inplace()
544       tp = t + (*vj) * 6;  in MatSolve_SeqSBAIJ_6_inplace()
562     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_6_inplace()
572     tp = t + (*vj) * 6;  in MatSolve_SeqSBAIJ_6_inplace()
581       vj++;  in MatSolve_SeqSBAIJ_6_inplace()
582       tp = t + (*vj) * 6;  in MatSolve_SeqSBAIJ_6_inplace()
613   const PetscInt  *vj;  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()  local
626     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
627 …PetscPrefetchBlock(vj + nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
630       xp = x + (*vj) * 6;  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
638       vj++;  in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
658   const PetscInt  *vj;  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()  local
671     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
672 …PetscPrefetchBlock(vj - nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
675       xp = x + (*vj) * 6;  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
683       vj++;  in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
765   const PetscInt    *r, *vj;  in MatSolve_SeqSBAIJ_5_inplace()  local
791     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
800     tp = t + (*vj) * 5;  in MatSolve_SeqSBAIJ_5_inplace()
807       vj++;  in MatSolve_SeqSBAIJ_5_inplace()
808       tp = t + (*vj) * 5;  in MatSolve_SeqSBAIJ_5_inplace()
825     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_5_inplace()
834     tp = t + (*vj) * 5;  in MatSolve_SeqSBAIJ_5_inplace()
842       vj++;  in MatSolve_SeqSBAIJ_5_inplace()
843       tp = t + (*vj) * 5;  in MatSolve_SeqSBAIJ_5_inplace()
872   const PetscInt  *vj;  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()  local
884     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
885 …PetscPrefetchBlock(vj + nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
888       xp = x + (*vj) * 5;  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
895       vj++;  in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
915   const PetscInt  *vj;  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()  local
927     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
928 …PetscPrefetchBlock(vj - nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
931       xp = x + (*vj) * 5;  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
938       vj++;  in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
1019   const PetscInt    *r, *vj;  in MatSolve_SeqSBAIJ_4_inplace()  local
1044     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1052     tp = t + (*vj) * 4;  in MatSolve_SeqSBAIJ_4_inplace()
1058       vj++;  in MatSolve_SeqSBAIJ_4_inplace()
1059       tp = t + (*vj) * 4;  in MatSolve_SeqSBAIJ_4_inplace()
1075     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_4_inplace()
1083     tp = t + (*vj) * 4;  in MatSolve_SeqSBAIJ_4_inplace()
1090       vj++;  in MatSolve_SeqSBAIJ_4_inplace()
1091       tp = t + (*vj) * 4;  in MatSolve_SeqSBAIJ_4_inplace()
1118   const PetscInt  *vj;  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()  local
1129     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1130 …PetscPrefetchBlock(vj + nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1133       xp = x + (*vj) * 4;  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1139       vj++;  in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1158   const PetscInt  *vj;  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()  local
1169     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1170 …PetscPrefetchBlock(vj - nz, nz, 0, PETSC_PREFETCH_HINT_NTA);          /* Indices for the next row …  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1173       xp = x + (*vj) * 4;  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1179       vj++;  in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1260   const PetscInt    *vj;  in MatSolve_SeqSBAIJ_3_inplace()  local
1283     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1290     tp = t + (*vj) * 3;  in MatSolve_SeqSBAIJ_3_inplace()
1295       vj++;  in MatSolve_SeqSBAIJ_3_inplace()
1296       tp = t + (*vj) * 3;  in MatSolve_SeqSBAIJ_3_inplace()
1311     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_3_inplace()
1318     tp = t + (*vj) * 3;  in MatSolve_SeqSBAIJ_3_inplace()
1324       vj++;  in MatSolve_SeqSBAIJ_3_inplace()
1325       tp = t + (*vj) * 3;  in MatSolve_SeqSBAIJ_3_inplace()
1350   const PetscInt  *vj;  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()  local
1360     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1361 …PetscPrefetchBlock(vj + nz, nz, 0, PETSC_PREFETCH_HINT_NTA);        /* Indices for the next row (a…  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1364       xp = x + (*vj) * 3;  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1369       vj++;  in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1387   const PetscInt  *vj;  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()  local
1397     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1398 …PetscPrefetchBlock(vj - nz, nz, 0, PETSC_PREFETCH_HINT_NTA);        /* Indices for the next row (a…  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1401       xp = x + (*vj) * 3;  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1406       vj++;  in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1485   const PetscInt    *r, *vj;  in MatSolve_SeqSBAIJ_2_inplace()  local
1505     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1511       t[(*vj) * 2] += v[0] * x0 + v[1] * x1;  in MatSolve_SeqSBAIJ_2_inplace()
1512       t[(*vj) * 2 + 1] += v[2] * x0 + v[3] * x1;  in MatSolve_SeqSBAIJ_2_inplace()
1513       vj++;  in MatSolve_SeqSBAIJ_2_inplace()
1524     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_2_inplace()
1530       x0 += v[0] * t[(*vj) * 2] + v[2] * t[(*vj) * 2 + 1];  in MatSolve_SeqSBAIJ_2_inplace()
1531       x1 += v[1] * t[(*vj) * 2] + v[3] * t[(*vj) * 2 + 1];  in MatSolve_SeqSBAIJ_2_inplace()
1532       vj++;  in MatSolve_SeqSBAIJ_2_inplace()
1554   const PetscInt  *vj;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()  local
1559     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1564 …PetscPrefetchBlock(vj + nz, nz, 0, PETSC_PREFETCH_HINT_NTA);        /* Indices for the next row (a…  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1568       x[(*vj) * 2] += v[0] * x0 + v[1] * x1;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1569       x[(*vj) * 2 + 1] += v[2] * x0 + v[3] * x1;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1570       vj++;  in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1586   const PetscInt  *vj;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()  local
1591     vj = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1596 …PetscPrefetchBlock(vj - nz, nz, 0, PETSC_PREFETCH_HINT_NTA);        /* Indices for the next row (a…  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1600       x0 += v[0] * x[(*vj) * 2] + v[2] * x[(*vj) * 2 + 1];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1601       x1 += v[1] * x[(*vj) * 2] + v[3] * x[(*vj) * 2 + 1];  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1602       vj++;  in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1678   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj, *adiag = a->diag;  in MatSolve_SeqSBAIJ_1()  local
1694     vj = aj + ai[k];  in MatSolve_SeqSBAIJ_1()
1697     for (j = 0; j < nz; j++) t[vj[j]] += v[j] * xk;  in MatSolve_SeqSBAIJ_1()
1704     vj = aj + adiag[k] - 1;  in MatSolve_SeqSBAIJ_1()
1706     for (j = 0; j < nz; j++) t[k] += v[-j] * t[vj[-j]];  in MatSolve_SeqSBAIJ_1()
1721   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj;  in MatSolve_SeqSBAIJ_1_inplace()  local
1737     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1740     while (nz--) t[*vj++] += (*v++) * xk;  in MatSolve_SeqSBAIJ_1_inplace()
1747     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_inplace()
1749     while (nz--) t[k] += (*v++) * t[*vj++];  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
1780     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_1()
1783     while (nz--) x[*vj++] += (*v++) * xk;  in MatForwardSolve_SeqSBAIJ_1()
1800   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj;  in MatForwardSolve_SeqSBAIJ_1_inplace()  local
1816     vj = aj + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_inplace()
1819     while (nz--) x[*vj++] += (*v++) * xk;  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
1852     vj    = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_1()
1857     while (nz--) t[k] += (*v++) * t[*vj++];  in MatBackwardSolve_SeqSBAIJ_1()
1871   const PetscInt     mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj;  in MatBackwardSolve_SeqSBAIJ_1_inplace()  local
1887     vj    = aj + ai[k] + 1;  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1892     while (nz--) t[k] += (*v++) * t[*vj++];  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1911     const PetscInt    *vj, mbs = a->mbs, *ai = a->i, *aj = a->j, *rp;  in MatSolves_SeqSBAIJ_1()  local
1935       vj = aj + ai[k];  in MatSolves_SeqSBAIJ_1()
1938         for (i = 0; i < n; i++) t[n * (*vj) + i] += (*v) * t[n * k + i];  in MatSolves_SeqSBAIJ_1()
1940         vj++;  in MatSolves_SeqSBAIJ_1()
1948       vj = aj + ai[k] - 1;  in MatSolves_SeqSBAIJ_1()
1951         for (i = 0; i < n; i++) t[n * k + i] += (*v) * t[n * (*vj) + i];  in MatSolves_SeqSBAIJ_1()
1953         vj++;  in MatSolves_SeqSBAIJ_1()
1975     const PetscInt    *vj, mbs = a->mbs, *ai = a->i, *aj = a->j, *rp;  in MatSolves_SeqSBAIJ_1_inplace()  local
1999       vj = aj + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2002         for (i = 0; i < n; i++) t[n * (*vj) + i] += (*v) * t[n * k + i];  in MatSolves_SeqSBAIJ_1_inplace()
2004         vj++;  in MatSolves_SeqSBAIJ_1_inplace()
2012       vj = aj + ai[k];  in MatSolves_SeqSBAIJ_1_inplace()
2015         for (i = 0; i < n; i++) t[n * k + i] += (*v) * t[n * (*vj) + i];  in MatSolves_SeqSBAIJ_1_inplace()
2017         vj++;  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
2046     vj = aj + ai[i];  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2049     for (j = 0; j < nz; j++) x[vj[j]] += v[j] * xi;  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2056     vj = aj + adiag[i] - 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2058     for (j = 0; j < nz; j++) xi += v[-j] * x[vj[-j]];  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
2094       vj = aj + ai[i];  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2097       for (j = 0; j < nz; j++) x[vj[j]] += v[j] * xi;  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2104       vj = aj + adiag[i] - 1;  in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2106       for (j = 0; j < nz; j++) xi += v[-j] * x[vj[-j]];  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
2135     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2138     while (nz--) x[*vj++] += (*v++) * xk;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2145     vj = aj + ai[k] + 1;  in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2148     while (nz--) xk += (*v++) * x[*vj++];  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
2175     vj = aj + ai[k];  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2177     while (nz--) x[*vj++] += (*v++) * x[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
2205     vj = aj + ai[k] + 1;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2207     while (nz--) x[*vj++] += (*v++) * x[k];  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
2235     vj    = aj + ai[k];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2240     while (nz--) x[k] += (*v++) * x[*vj++];  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
2265     vj    = aj + ai[k] + 1;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2270     while (nz--) x[k] += (*v++) * x[*vj++];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2285   PetscInt       *prowl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()  local
2326       vj = rip[aj[j++]];  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2327       if (vj > k) {  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2332         } while (qm < vj);  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2333         PetscCheck(qm != vj, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Duplicate entry in A");  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2335         q[m]       = vj;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2336         q[vj]      = qm;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2337         levtmp[vj] = 0; /* initialize lev for nonzero element */  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2354         vj = ju[j];  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2358         } while (qm < vj);  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2359         if (qm != vj) { /* a fill */  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2361           q[m]       = vj;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2362           q[vj]      = qm;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2363           qm         = vj;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2364           levtmp[vj] = incrlev;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2365         } else if (levtmp[vj] > incrlev) levtmp[vj] = incrlev;  in MatICCFactorSymbolic_SeqSBAIJ_MSR()