Lines Matching refs:mbs
13 PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_N_inplace() local
31 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_N_inplace()
35 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_N_inplace()
55 for (k = mbs - 1; k >= 0; k--) { in MatSolve_SeqSBAIJ_N_inplace()
76 PetscCall(PetscLogFlops(4.0 * bs2 * a->nz - (bs + 2.0 * bs2) * mbs)); in MatSolve_SeqSBAIJ_N_inplace()
92 …ring(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscInt bs, Petsc… in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering() argument
102 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
124 …ring(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscInt bs, Petsc… in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering() argument
132 for (k = mbs - 1; k >= 0; 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
163 PetscCall(PetscArraycpy(x, b, bs * mbs)); /* x <- b */ in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
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()
171 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (bs + 2.0 * a->bs2) * mbs)); 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
187 PetscCall(PetscArraycpy(x, b, bs * mbs)); /* x <- b */ in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
188 PetscCall(MatForwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x)); in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
191 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - bs * mbs)); 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
207 PetscCall(PetscArraycpy(x, b, bs * mbs)); in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
208 PetscCall(MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x)); in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
211 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); 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
233 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_7_inplace()
245 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_7_inplace()
284 for (k = mbs - 1; k >= 0; k--) { 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()
336 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering() argument
344 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering()
385 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering() argument
393 for (k = mbs - 1; k >= 0; 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
445 PetscCall(PetscArraycpy(x, b, 7 * mbs)); /* x <- b */ in MatSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
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()
453 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); 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
468 PetscCall(PetscArraycpy(x, b, 7 * mbs)); in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
469 PetscCall(MatForwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x)); in MatForwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
472 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); 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
487 PetscCall(PetscArraycpy(x, b, 7 * mbs)); in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
488 PetscCall(MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering(ai, aj, aa, mbs, x)); in MatBackwardSolve_SeqSBAIJ_7_NaturalOrdering_inplace()
491 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); 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
513 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_6_inplace()
524 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_6_inplace()
560 for (k = mbs - 1; k >= 0; k--) { 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()
608 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering() argument
616 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering()
653 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering() argument
661 for (k = mbs - 1; k >= 0; 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
710 PetscCall(PetscArraycpy(x, b, 6 * mbs)); /* x <- b */ in MatSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
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()
718 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); 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
733 PetscCall(PetscArraycpy(x, b, 6 * mbs)); /* x <- b */ in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
734 PetscCall(MatForwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x)); in MatForwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
737 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); 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
752 PetscCall(PetscArraycpy(x, b, 6 * mbs)); /* x <- b */ in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
753 PetscCall(MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering(ai, aj, aa, mbs, x)); in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
756 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_6_NaturalOrdering_inplace()
764 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_5_inplace() local
779 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_5_inplace()
789 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_5_inplace()
823 for (k = mbs - 1; k >= 0; k--) { 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()
867 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering() argument
875 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering()
910 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering() argument
918 for (k = mbs - 1; k >= 0; 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
964 PetscCall(PetscArraycpy(x, b, 5 * mbs)); /* x <- b */ in MatSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
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()
972 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); 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
987 PetscCall(PetscArraycpy(x, b, 5 * mbs)); /* x <- b */ in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
988 PetscCall(MatForwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x)); in MatForwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
991 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); 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
1006 PetscCall(PetscArraycpy(x, b, 5 * mbs)); in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1007 PetscCall(MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering(ai, aj, aa, mbs, x)); in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1010 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_5_NaturalOrdering_inplace()
1018 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_4_inplace() local
1033 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_4_inplace()
1042 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_4_inplace()
1073 for (k = mbs - 1; k >= 0; k--) { 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()
1113 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering() argument
1121 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering()
1153 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering() argument
1161 for (k = mbs - 1; k >= 0; 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
1204 PetscCall(PetscArraycpy(x, b, 4 * mbs)); /* x <- b */ in MatSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
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()
1211 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); 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
1226 PetscCall(PetscArraycpy(x, b, 4 * mbs)); /* x <- b */ in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1227 PetscCall(MatForwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x)); in MatForwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1230 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); 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
1245 PetscCall(PetscArraycpy(x, b, 4 * mbs)); in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1246 PetscCall(MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering(ai, aj, aa, mbs, x)); in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1249 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_4_NaturalOrdering_inplace()
1257 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_3_inplace() local
1273 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_3_inplace()
1281 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_3_inplace()
1309 for (k = mbs - 1; k >= 0; k--) { 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()
1345 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering() argument
1353 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering()
1382 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering() argument
1390 for (k = mbs - 1; k >= 0; 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
1430 PetscCall(PetscArraycpy(x, b, 3 * mbs)); in MatSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
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()
1438 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); 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
1453 PetscCall(PetscArraycpy(x, b, 3 * mbs)); in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1454 PetscCall(MatForwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x)); in MatForwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1457 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); 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
1472 PetscCall(PetscArraycpy(x, b, 3 * mbs)); in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1473 PetscCall(MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering(ai, aj, aa, mbs, x)); in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1476 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); in MatBackwardSolve_SeqSBAIJ_3_NaturalOrdering_inplace()
1484 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_2_inplace() local
1498 for (k = 0; k < mbs; k++) { /* t <- perm(b) */ in MatSolve_SeqSBAIJ_2_inplace()
1503 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_2_inplace()
1522 for (k = mbs - 1; k >= 0; k--) { 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()
1549 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering() argument
1557 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering()
1581 …Ordering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa, PetscInt mbs, PetscScalar *x) in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering() argument
1589 for (k = mbs - 1; k >= 0; 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
1624 PetscCall(PetscArraycpy(x, b, 2 * mbs)); in MatSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
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()
1632 PetscCall(PetscLogFlops(4.0 * a->bs2 * a->nz - (A->rmap->bs + 2.0 * a->bs2) * mbs)); 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
1647 PetscCall(PetscArraycpy(x, b, 2 * mbs)); in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1648 PetscCall(MatForwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x)); in MatForwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1651 PetscCall(PetscLogFlops(2.0 * a->bs2 * a->nz - A->rmap->bs * mbs)); 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
1666 PetscCall(PetscArraycpy(x, b, 2 * mbs)); in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1667 PetscCall(MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering(ai, aj, aa, mbs, x)); in MatBackwardSolve_SeqSBAIJ_2_NaturalOrdering_inplace()
1670 PetscCall(PetscLogFlops(2.0 * a->bs2 * (a->nz - mbs))); 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
1691 for (k = 0; k < mbs; k++) t[k] = b[rp[k]]; in MatSolve_SeqSBAIJ_1()
1692 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_1()
1702 for (k = mbs - 1; k >= 0; k--) { in MatSolve_SeqSBAIJ_1()
1713 PetscCall(PetscLogFlops(4.0 * a->nz - 3.0 * mbs)); in MatSolve_SeqSBAIJ_1()
1721 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj; in MatSolve_SeqSBAIJ_1_inplace() local
1734 for (k = 0; k < mbs; k++) t[k] = b[rp[k]]; in MatSolve_SeqSBAIJ_1_inplace()
1735 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_1_inplace()
1745 for (k = mbs - 1; k >= 0; k--) { in MatSolve_SeqSBAIJ_1_inplace()
1756 PetscCall(PetscLogFlops(4.0 * a->nz - 3 * mbs)); 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
1777 for (k = 0; k < mbs; k++) x[k] = b[rp[k]]; in MatForwardSolve_SeqSBAIJ_1()
1778 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_1()
1792 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatForwardSolve_SeqSBAIJ_1()
1800 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj; in MatForwardSolve_SeqSBAIJ_1_inplace() local
1813 for (k = 0; k < mbs; k++) x[k] = b[rp[k]]; in MatForwardSolve_SeqSBAIJ_1_inplace()
1814 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_1_inplace()
1828 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); 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
1850 for (k = mbs - 1; k >= 0; k--) { in MatBackwardSolve_SeqSBAIJ_1()
1863 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1()
1871 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j, *rp, *vj; in MatBackwardSolve_SeqSBAIJ_1_inplace() local
1885 for (k = mbs - 1; k >= 0; k--) { in MatBackwardSolve_SeqSBAIJ_1_inplace()
1898 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1_inplace()
1911 const PetscInt *vj, mbs = a->mbs, *ai = a->i, *aj = a->j, *rp; in MatSolves_SeqSBAIJ_1() local
1930 for (k = 0; k < mbs; k++) { in MatSolves_SeqSBAIJ_1()
1931 … for (i = 0; i < n; i++) t[n * k + i] = b[rp[k] + i * mbs]; /* values are stored interlaced in t */ in MatSolves_SeqSBAIJ_1()
1933 for (k = 0; k < mbs; k++) { in MatSolves_SeqSBAIJ_1()
1946 for (k = mbs - 1; k >= 0; k--) { in MatSolves_SeqSBAIJ_1()
1955 for (i = 0; i < n; i++) x[rp[k] + i * mbs] = t[n * k + i]; in MatSolves_SeqSBAIJ_1()
1961 PetscCall(PetscLogFlops(bb->n * (4.0 * a->nz - 3.0 * mbs))); in MatSolves_SeqSBAIJ_1()
1975 const PetscInt *vj, mbs = a->mbs, *ai = a->i, *aj = a->j, *rp; in MatSolves_SeqSBAIJ_1_inplace() local
1994 for (k = 0; k < mbs; k++) { in MatSolves_SeqSBAIJ_1_inplace()
1995 … for (i = 0; i < n; i++) t[n * k + i] = b[rp[k] + i * mbs]; /* values are stored interlaced in t */ in MatSolves_SeqSBAIJ_1_inplace()
1997 for (k = 0; k < mbs; k++) { in MatSolves_SeqSBAIJ_1_inplace()
2010 for (k = mbs - 1; k >= 0; k--) { in MatSolves_SeqSBAIJ_1_inplace()
2019 for (i = 0; i < n; i++) x[rp[k] + i * mbs] = t[n * k + i]; in MatSolves_SeqSBAIJ_1_inplace()
2025 PetscCall(PetscLogFlops(bb->n * (4.0 * a->nz - 3.0 * mbs))); 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
2043 PetscCall(PetscArraycpy(x, b, mbs)); in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2044 for (i = 0; i < mbs; i++) { in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2053 for (i = mbs - 2; i >= 0; i--) { in MatSolve_SeqSBAIJ_1_NaturalOrdering()
2063 PetscCall(PetscLogFlops(4.0 * a->nz - 3 * mbs)); 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
2078 if (!mbs) PetscFunctionReturn(PETSC_SUCCESS); in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2091 PetscCall(PetscArraycpy(x, b, mbs)); in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2092 for (i = 0; i < mbs; i++) { in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2101 for (i = mbs - 2; i >= 0; i--) { in MatMatSolve_SeqSBAIJ_1_NaturalOrdering()
2114 PetscCall(PetscLogFlops(B->cmap->n * (4.0 * a->nz - 3 * mbs))); 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
2132 PetscCall(PetscArraycpy(x, b, mbs)); in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2133 for (k = 0; k < mbs; k++) { in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2143 for (k = mbs - 2; k >= 0; k--) { in MatSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2154 PetscCall(PetscLogFlops(4.0 * a->nz - 3 * mbs)); 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
2172 PetscCall(PetscArraycpy(x, b, mbs)); in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2173 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2184 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); 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
2202 PetscCall(PetscArraycpy(x, b, mbs)); in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2203 for (k = 0; k < mbs; k++) { in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2214 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); 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
2233 for (k = mbs - 1; k >= 0; k--) { in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2244 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); 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
2263 for (k = mbs - 1; k >= 0; k--) { in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2274 PetscCall(PetscLogFlops(2.0 * a->nz - mbs)); in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2282 const PetscInt *rip, mbs = a->mbs, *ai, *aj; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2301 umax = (PetscInt)(f * ai[mbs] + 1); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2303 umax += mbs + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2304 shift = mbs + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2305 PetscCall(PetscMalloc1(mbs + 1, &iu)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2307 iu[0] = mbs + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2308 juidx = mbs + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2310 PetscCall(PetscMalloc3(mbs, &prowl, mbs, &q, mbs, &levtmp)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2313 for (i = 0; i < mbs; i++) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2314 prowl[i] = mbs; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2319 for (k = 0; k < mbs; k++) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2321 q[k] = mbs; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2384 if (maxadd < nzk) maxadd = (mbs - k) * (nzk + 1) / 2; 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()
2434 PetscCall(PetscShmgetAllocateArray((iu[mbs] + 1) * a->bs2, sizeof(PetscScalar), (void **)&b->a)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2449 PetscCall(PetscMalloc1(bs * mbs + bs, &b->solve_work)); in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2452 b->maxnz = b->nz = iu[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()
2676 PetscInt bs = A->rmap->bs, am = a->mbs; in MatICCFactorSymbolic_SeqSBAIJ_inplace()