Lines Matching refs:nz
15 PetscInt nz, *vj, k, idx, k1; in MatSolve_SeqSBAIJ_N_inplace() local
39 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_N_inplace()
42 while (nz--) { in MatSolve_SeqSBAIJ_N_inplace()
58 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_N_inplace()
61 while (nz--) { in MatSolve_SeqSBAIJ_N_inplace()
76 PetscCall(PetscLogFlops(4.0 * bs2 * a->nz - (bs + 2.0 * bs2) * mbs)); in MatSolve_SeqSBAIJ_N_inplace()
94 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering() local
106 nz = ai[k + 1] - ai[k]; in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
109 while (nz--) { in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
126 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering() local
135 nz = ai[k + 1] - ai[k]; in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
138 while (nz--) { in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
171 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
191 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - bs * mbs)); in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
211 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
220 PetscInt nz, k, idx; in MatSolve_SeqSBAIJ_7_inplace() local
256 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_7_inplace()
258 while (nz--) { in MatSolve_SeqSBAIJ_7_inplace()
295 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_7_inplace()
298 while (nz--) { in MatSolve_SeqSBAIJ_7_inplace()
332 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_7_inplace()
340 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering() local
354 nz = ai[k + 1] - 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()
357 …PetscPrefetchBlock(v + 49 * nz, 49 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
358 while (nz--) { in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
389 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering() local
403 nz = ai[k + 1] - 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()
406 …PetscPrefetchBlock(v - 49 * nz, 49 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
407 while (nz--) { in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering()
453 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
472 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
491 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
500 PetscInt nz, k, idx; in MatSolve_SeqSBAIJ_6_inplace() local
534 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_6_inplace()
536 while (nz--) { in MatSolve_SeqSBAIJ_6_inplace()
570 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_6_inplace()
573 while (nz--) { in MatSolve_SeqSBAIJ_6_inplace()
604 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_6_inplace()
612 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering() local
625 nz = ai[k + 1] - 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()
628 …PetscPrefetchBlock(v + 36 * nz, 36 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
629 while (nz--) { in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
657 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering() local
670 nz = ai[k + 1] - 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()
673 …PetscPrefetchBlock(v - 36 * nz, 36 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
674 while (nz--) { in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering()
718 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
737 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
756 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
766 PetscInt nz, k, idx; in MatSolve_SeqSBAIJ_5_inplace() local
798 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_5_inplace()
801 while (nz--) { in MatSolve_SeqSBAIJ_5_inplace()
832 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_5_inplace()
835 while (nz--) { in MatSolve_SeqSBAIJ_5_inplace()
863 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_5_inplace()
871 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering() local
883 nz = ai[k + 1] - 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()
886 …PetscPrefetchBlock(v + 25 * nz, 25 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
887 while (nz--) { in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
914 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering() local
926 nz = ai[k + 1] - 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()
929 …PetscPrefetchBlock(v - 25 * nz, 25 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
930 while (nz--) { in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering()
972 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
991 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1010 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1020 PetscInt nz, k, idx; in MatSolve_SeqSBAIJ_4_inplace() local
1050 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_4_inplace()
1053 while (nz--) { in MatSolve_SeqSBAIJ_4_inplace()
1081 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_4_inplace()
1084 while (nz--) { in MatSolve_SeqSBAIJ_4_inplace()
1109 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_4_inplace()
1117 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering() local
1128 nz = ai[k + 1] - 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()
1131 …PetscPrefetchBlock(v + 16 * nz, 16 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1132 while (nz--) { in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1157 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering() local
1168 nz = ai[k + 1] - 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()
1171 …PetscPrefetchBlock(v - 16 * nz, 16 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row … in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1172 while (nz--) { in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering()
1211 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1230 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1249 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1259 PetscInt nz, k, idx; in MatSolve_SeqSBAIJ_3_inplace() local
1288 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_3_inplace()
1291 while (nz--) { in MatSolve_SeqSBAIJ_3_inplace()
1316 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_3_inplace()
1319 while (nz--) { in MatSolve_SeqSBAIJ_3_inplace()
1341 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_3_inplace()
1349 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering() local
1359 nz = ai[k + 1] - 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()
1362 … PetscPrefetchBlock(v + 9 * nz, 9 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1363 while (nz--) { in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1386 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering() local
1396 nz = ai[k + 1] - 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()
1399 … PetscPrefetchBlock(v - 9 * nz, 9 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1400 while (nz--) { in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering()
1438 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1457 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1476 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1486 PetscInt nz, k, k2, idx; in MatSolve_SeqSBAIJ_2_inplace() local
1509 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_2_inplace()
1510 while (nz--) { in MatSolve_SeqSBAIJ_2_inplace()
1528 nz = ai[k + 1] - ai[k]; in MatSolve_SeqSBAIJ_2_inplace()
1529 while (nz--) { in MatSolve_SeqSBAIJ_2_inplace()
1545 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_2_inplace()
1553 PetscInt nz, k, k2; in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering() local
1563 nz = ai[k + 1] - 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()
1565 … PetscPrefetchBlock(v + 4 * nz, 4 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1566 while (nz--) { in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1585 PetscInt nz, k, k2; in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering() local
1595 nz = ai[k + 1] - 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()
1597 … PetscPrefetchBlock(v - 4 * nz, 4 * nz, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1598 while (nz--) { in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering()
1632 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); in MatSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1651 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1670 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1682 PetscInt nz, k, j; in MatSolve_SeqSBAIJ_1() local
1696 nz = ai[k + 1] - ai[k] - 1; 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()
1705 nz = ai[k + 1] - ai[k] - 1; in MatSolve_SeqSBAIJ_1()
1706 for (j = 0; j < nz; j++) t[k] += v[-j] * t[vj[-j]]; in MatSolve_SeqSBAIJ_1()
1713 PetscCall(PetscLogFlops(4.0 * a->nz - 3.0 * mbs)); in MatSolve_SeqSBAIJ_1()
1725 PetscInt nz, k; in MatSolve_SeqSBAIJ_1_inplace() local
1739 nz = ai[k + 1] - ai[k] - 1; in MatSolve_SeqSBAIJ_1_inplace()
1740 while (nz--) t[*vj++] += (*v++) * xk; in MatSolve_SeqSBAIJ_1_inplace()
1748 nz = ai[k + 1] - ai[k] - 1; in MatSolve_SeqSBAIJ_1_inplace()
1749 while (nz--) t[k] += (*v++) * t[*vj++]; in MatSolve_SeqSBAIJ_1_inplace()
1756 PetscCall(PetscLogFlops(4.0 * a->nz - 3 * mbs)); in MatSolve_SeqSBAIJ_1_inplace()
1769 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_1() local
1782 nz = ai[k + 1] - ai[k] - 1; in MatForwardSolve_SeqSBAIJ_1()
1783 while (nz--) x[*vj++] += (*v++) * xk; in MatForwardSolve_SeqSBAIJ_1()
1792 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatForwardSolve_SeqSBAIJ_1()
1805 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_1_inplace() local
1818 nz = ai[k + 1] - ai[k] - 1; in MatForwardSolve_SeqSBAIJ_1_inplace()
1819 while (nz--) x[*vj++] += (*v++) * xk; in MatForwardSolve_SeqSBAIJ_1_inplace()
1828 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatForwardSolve_SeqSBAIJ_1_inplace()
1841 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_1() local
1856 nz = ai[k + 1] - ai[k] - 1; in MatBackwardSolve_SeqSBAIJ_1()
1857 while (nz--) t[k] += (*v++) * t[*vj++]; in MatBackwardSolve_SeqSBAIJ_1()
1863 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1()
1876 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_1_inplace() local
1891 nz = ai[k + 1] - ai[k] - 1; in MatBackwardSolve_SeqSBAIJ_1_inplace()
1892 while (nz--) t[k] += (*v++) * t[*vj++]; in MatBackwardSolve_SeqSBAIJ_1_inplace()
1898 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1_inplace()
1915 PetscInt nz, k, n, i, j; in MatSolves_SeqSBAIJ_1() local
1936 nz = ai[k + 1] - ai[k] - 1; in MatSolves_SeqSBAIJ_1()
1937 for (j = 0; j < nz; j++) { 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()
1949 nz = ai[k + 1] - ai[k] - 1; in MatSolves_SeqSBAIJ_1()
1950 for (j = 0; j < nz; j++) { in MatSolves_SeqSBAIJ_1()
1961 PetscCall(PetscLogFlops(bb->n * (4.0 * a->nz - 3.0 * mbs))); in MatSolves_SeqSBAIJ_1()
1979 PetscInt nz, k, n, i; in MatSolves_SeqSBAIJ_1_inplace() local
2000 nz = ai[k + 1] - ai[k]; in MatSolves_SeqSBAIJ_1_inplace()
2001 while (nz--) { in MatSolves_SeqSBAIJ_1_inplace()
2013 nz = ai[k + 1] - ai[k]; in MatSolves_SeqSBAIJ_1_inplace()
2014 while (nz--) { in MatSolves_SeqSBAIJ_1_inplace()
2025 PetscCall(PetscLogFlops(bb->n * (4.0 * a->nz - 3.0 * mbs))); in MatSolves_SeqSBAIJ_1_inplace()
2037 PetscInt nz, i, j; in MatSolve_SeqSBAIJ_1_NaturalOrdering() local
2048 nz = ai[i + 1] - ai[i] - 1; /* exclude diag[i] */ in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2049 for (j = 0; j < nz; j++) x[vj[j]] += v[j] * xi; in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2050 x[i] = xi * v[nz]; /* v[nz] = aa[diag[i]] = 1/D(i) */ in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2057 nz = ai[i + 1] - ai[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()
2063 PetscCall(PetscLogFlops(4.0 * a->nz - 3 * mbs)); in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2074 PetscInt nz, i, j, neq, ldb, ldx; in MatMatSolve_SeqSBAIJ_1_NaturalOrdering() local
2096 nz = ai[i + 1] - ai[i] - 1; /* exclude diag[i] */ in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2097 for (j = 0; j < nz; j++) x[vj[j]] += v[j] * xi; in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2098 x[i] = xi * v[nz]; /* v[nz] = aa[diag[i]] = 1/D(i) */ in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2105 nz = ai[i + 1] - ai[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()
2114 PetscCall(PetscLogFlops(B->cmap->n * (4.0 * a->nz - 3 * mbs))); in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2125 PetscInt nz, k; in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace() local
2137 nz = ai[k + 1] - ai[k] - 1; /* exclude diag[k] */ in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2138 while (nz--) x[*vj++] += (*v++) * xk; in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2147 nz = ai[k + 1] - ai[k] - 1; in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2148 while (nz--) xk += (*v++) * x[*vj++]; in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2154 PetscCall(PetscLogFlops(4.0 * a->nz - 3 * mbs)); in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2166 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering() local
2176 nz = ai[k + 1] - ai[k] - 1; /* exclude diag[k] */ in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2177 while (nz--) x[*vj++] += (*v++) * x[k]; in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2184 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2196 PetscInt nz, k; in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace() local
2206 nz = ai[k + 1] - ai[k] - 1; /* exclude diag[k] */ in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2207 while (nz--) x[*vj++] += (*v++) * x[k]; in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2214 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2226 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering() local
2239 nz = ai[k + 1] - ai[k] - 1; in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2240 while (nz--) x[k] += (*v++) * x[*vj++]; in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2244 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2256 PetscInt nz, k; in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace() local
2269 nz = ai[k + 1] - ai[k] - 1; in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2270 while (nz--) x[k] += (*v++) * x[*vj++]; in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2274 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2286 PetscInt incrlev, *lev, shift, prow, nz; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2323 nz = ai[rip[k] + 1] - ai[rip[k]]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2325 while (nz--) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2452 b->maxnz = b->nz = iu[mbs]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2648 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqSBAIJ()
2844 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqSBAIJ_inplace()