Lines Matching refs:nz
32 PetscInt nzL, nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2() local
34 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
37 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_2()
39 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
42 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_2()
45 …for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, … in MatLUFactorNumeric_SeqBAIJ_2()
68 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_2()
69 for (PetscInt j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqBAIJ_2()
76 PetscCall(PetscLogFlops(16.0 * nz + 12)); /* flops = 2*bs^3*nz + 2*bs^3 - bs2) */ in MatLUFactorNumeric_SeqBAIJ_2()
84 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
85 … for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
101 nz = bdiag[i] - bdiag[i + 1] - 1; in MatLUFactorNumeric_SeqBAIJ_2()
102 … for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
121 PetscInt i, j, k, nz, nzL, row, *pj; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering() local
140 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
142 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
145 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
147 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
150 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
153 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ajtmp[j], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
174 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
175 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
182 PetscCall(PetscLogFlops(16.0 * nz + 12)); /* flops = 2*bs^3*nz + 2*bs^3 - bs2) */ in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
190 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
191 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
208 nz = bdiag[i] - bdiag[i + 1] - 1; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
209 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
230 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local
250 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
252 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_inplace()
258 nz = ai[idx + 1] - ai[idx]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
261 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_inplace()
287 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
289 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_inplace()
301 PetscCall(PetscLogFlops(16.0 * nz + 12.0)); in MatILUFactorNumeric_SeqBAIJ_2_inplace()
308 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
309 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_inplace()
341 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() local
358 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
360 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
365 nz = ai[i + 1] - ai[i]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
368 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
394 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
396 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
408 PetscCall(PetscLogFlops(16.0 * nz + 12.0)); in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
415 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
416 for (j = 0; j < nz; j++) { in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
456 PetscInt i, j, k, nz, nzL, row, *pj; in MatLUFactorNumeric_SeqBAIJ_1() local
478 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqBAIJ_1()
479 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqBAIJ_1()
498 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_1()
500 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqBAIJ_1()
503 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_1()
505 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqBAIJ_1()
508 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_1()
511 for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqBAIJ_1()
529 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqBAIJ_1()
530 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatLUFactorNumeric_SeqBAIJ_1()
531 PetscCall(PetscLogFlops(2.0 * nz)); in MatLUFactorNumeric_SeqBAIJ_1()
541 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_1()
542 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqBAIJ_1()
550 nz = bdiag[i] - bdiag[i + 1] - 1; in MatLUFactorNumeric_SeqBAIJ_1()
551 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqBAIJ_1()
619 PetscInt *ajtmpold, *ajtmp, nz, row, *ai = a->i, *aj = a->j; in MatILUFactorNumeric_SeqBAIJ_1_inplace() local
636 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
638 for (j = 0; j < nz; j++) rtmp[ajtmp[j]] = 0.0; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
641 nz = ai[r[i] + 1] - ai[r[i]]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
644 for (j = 0; j < nz; j++) rtmp[ic[ajtmpold[j]]] = v[j]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
654 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
655 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
656 PetscCall(PetscLogFlops(1.0 + 2.0 * nz)); in MatILUFactorNumeric_SeqBAIJ_1_inplace()
663 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
664 for (j = 0; j < nz; j++) pv[j] = rtmp[pj[j]]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
755 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
835 nz = bi[k + 1] - jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N()
836 if (nz) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
838 while (nz--) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
892 PetscInt k, jmin, *jl, *il, nexti, ili, *acol, *bcol, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering() local
913 nz = ai[k + 1] - ai[k]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
917 while (nz--) { in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
944 nz = bi[i + 1] - jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
945 if (nz > 0) { in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
948 while (nz--) rtmp[*bcol++] += uikdi * (*bval++); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
962 nz = bi[k + 1] - jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
963 if (nz) { in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
965 while (nz--) { in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
980 nz = bi[k + 1] - jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
981 if (nz) { in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
984 while (nz--) { in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
1215 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqBAIJ()
1395 b->maxnz = b->nz = ui[mbs]; in MatCholeskyFactorSymbolic_SeqBAIJ()
1428 PetscInt nz, bs = A->rmap->bs, bs2 = a->bs2; in MatSolve_SeqBAIJ_N_NaturalOrdering() local
1445 nz = ai[i + 1] - ai[i]; in MatSolve_SeqBAIJ_N_NaturalOrdering()
1448 for (k = 0; k < nz; k++) { in MatSolve_SeqBAIJ_N_NaturalOrdering()
1459 nz = adiag[i] - adiag[i + 1] - 1; in MatSolve_SeqBAIJ_N_NaturalOrdering()
1461 for (k = 0; k < nz; k++) { in MatSolve_SeqBAIJ_N_NaturalOrdering()
1471 PetscCall(PetscLogFlops(2.0 * (a->bs2) * (a->nz) - A->rmap->bs * A->cmap->n)); in MatSolve_SeqBAIJ_N_NaturalOrdering()
1481 PetscInt nz, bs = A->rmap->bs, bs2 = a->bs2; in MatSolve_SeqBAIJ_N() local
1502 nz = ai[i + 1] - ai[i]; in MatSolve_SeqBAIJ_N()
1505 for (m = 0; m < nz; m++) { in MatSolve_SeqBAIJ_N()
1516 nz = adiag[i] - adiag[i + 1] - 1; in MatSolve_SeqBAIJ_N()
1518 for (m = 0; m < nz; m++) { in MatSolve_SeqBAIJ_N()
1529 PetscCall(PetscLogFlops(2.0 * (a->bs2) * (a->nz) - A->rmap->bs * A->cmap->n)); in MatSolve_SeqBAIJ_N()