Lines Matching refs:nz
169 b->maxnz = b->nz = bdiag[0] + 1; in MatLUFactorSymbolic_SeqAIJ()
225 PetscInt i, j, k, nz, nzL, row, *pj; in MatLUFactorNumeric_SeqAIJ() local
250 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ()
251 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqAIJ()
270 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ()
272 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ()
275 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ()
277 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ()
280 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
283 for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqAIJ()
300 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ()
302 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatLUFactorNumeric_SeqAIJ()
303 PetscCall(PetscLogFlops(1 + 2.0 * nz)); in MatLUFactorNumeric_SeqAIJ()
313 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ()
314 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqAIJ()
322 nz = bdiag[i] - bdiag[i + 1] - 1; in MatLUFactorNumeric_SeqAIJ()
323 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqAIJ()
403 PetscInt nz, row, i, j, n = A->rmap->n, diag; in MatLUFactorNumeric_SeqAIJ_inplace() local
429 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ_inplace()
430 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqAIJ_inplace()
447 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ_inplace()
449 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ_inplace()
452 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_inplace()
455 for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqAIJ_inplace()
466 nz = bi[row + 1] - ddiag[row] - 1; in MatLUFactorNumeric_SeqAIJ_inplace()
467 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatLUFactorNumeric_SeqAIJ_inplace()
468 PetscCall(PetscLogFlops(1 + 2.0 * nz)); in MatLUFactorNumeric_SeqAIJ_inplace()
475 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ_inplace()
478 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqAIJ_inplace()
561 PetscInt *ajtmp, nz, row; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm() local
588 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
589 for (j = 0; j < nz; j++) rs += PetscAbsScalar(vtmp[j]); in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
618 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
619 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
637 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
641 for (j = 0; j < nz; j++) ajtmp[j] = ics[ajtmp[j]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
642 PetscCall(PetscSortIntWithScalarArray(nz, ajtmp, v)); in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
645 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
660 nz = ai[r[row] + 1] - diag[r[row]] - 1; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
661 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
662 PetscCall(PetscLogFlops(1 + 2.0 * nz)); in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
669 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
673 for (j = 0; j < nz; j++) { in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
746 PetscInt nz; in MatSolve_SeqAIJ_inplace() local
773 nz = adiag[i] - ai[i]; in MatSolve_SeqAIJ_inplace()
775 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatSolve_SeqAIJ_inplace()
783 nz = ai[i + 1] - adiag[i] - 1; in MatSolve_SeqAIJ_inplace()
785 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatSolve_SeqAIJ_inplace()
793 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolve_SeqAIJ_inplace()
802 PetscInt nz, neq, ldb, ldx; in MatMatSolve_SeqAIJ_inplace() local
834 nz = adiag[i] - ai[i]; in MatMatSolve_SeqAIJ_inplace()
836 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatMatSolve_SeqAIJ_inplace()
843 nz = ai[i + 1] - adiag[i] - 1; in MatMatSolve_SeqAIJ_inplace()
845 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatMatSolve_SeqAIJ_inplace()
856 PetscCall(PetscLogFlops(B->cmap->n * (2.0 * a->nz - n))); in MatMatSolve_SeqAIJ_inplace()
866 PetscInt nz, neq, ldb, ldx; in MatMatSolve_SeqAIJ() local
896 nz = ai[i + 1] - ai[i]; in MatMatSolve_SeqAIJ()
898 PetscSparseDenseMinusDot(sum, tmp, v, vi, nz); in MatMatSolve_SeqAIJ()
900 v += nz; in MatMatSolve_SeqAIJ()
901 vi += nz; in MatMatSolve_SeqAIJ()
907 nz = adiag[i] - adiag[i + 1] - 1; in MatMatSolve_SeqAIJ()
909 PetscSparseDenseMinusDot(sum, tmp, v, vi, nz); in MatMatSolve_SeqAIJ()
910 x[c[i]] = tmp[i] = sum * v[nz]; /* v[nz] = aa[adiag[i]] */ in MatMatSolve_SeqAIJ()
920 PetscCall(PetscLogFlops(B->cmap->n * (2.0 * a->nz - n))); in MatMatSolve_SeqAIJ()
930 PetscInt nz, neq, ldb, ldx; in MatMatSolveTranspose_SeqAIJ() local
962 nz = adiag[i] - adiag[i + 1] - 1; in MatMatSolveTranspose_SeqAIJ()
964 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatMatSolveTranspose_SeqAIJ()
965 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatMatSolveTranspose_SeqAIJ()
973 nz = ai[i + 1] - ai[i]; in MatMatSolveTranspose_SeqAIJ()
975 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatMatSolveTranspose_SeqAIJ()
988 PetscCall(PetscLogFlops(B->cmap->n * (2.0 * a->nz - n))); in MatMatSolveTranspose_SeqAIJ()
998 PetscInt nz, row; in MatSolve_SeqAIJ_InplaceWithPerm() local
1024 nz = adiag[i] - ai[i]; in MatSolve_SeqAIJ_InplaceWithPerm()
1026 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatSolve_SeqAIJ_InplaceWithPerm()
1035 nz = ai[i + 1] - adiag[i] - 1; in MatSolve_SeqAIJ_InplaceWithPerm()
1037 PetscSparseDenseMinusDot(sum, tmps, v, vi, nz); in MatSolve_SeqAIJ_InplaceWithPerm()
1045 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolve_SeqAIJ_InplaceWithPerm()
1059 PetscInt adiag_i, i, nz, ai_i; in MatSolve_SeqAIJ_NaturalOrdering_inplace() local
1081 nz = adiag[i] - ai_i; in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1083 PetscSparseDenseMinusDot(sum, x, v, vi, nz); in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1092 nz = ai[i + 1] - adiag_i - 1; in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1094 PetscSparseDenseMinusDot(sum, x, v, vi, nz); in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1098 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1110 PetscInt nz; in MatSolveAdd_SeqAIJ_inplace() local
1135 nz = adiag[i] - ai[i]; in MatSolveAdd_SeqAIJ_inplace()
1137 for (j = 0; j < nz; j++) sum -= v[j] * tmp[vi[j]]; in MatSolveAdd_SeqAIJ_inplace()
1145 nz = ai[i + 1] - adiag[i] - 1; in MatSolveAdd_SeqAIJ_inplace()
1147 for (j = 0; j < nz; j++) sum -= v[j] * tmp[vi[j]]; in MatSolveAdd_SeqAIJ_inplace()
1157 PetscCall(PetscLogFlops(2.0 * a->nz)); in MatSolveAdd_SeqAIJ_inplace()
1166 PetscInt nz; in MatSolveAdd_SeqAIJ() local
1191 nz = ai[i + 1] - ai[i]; in MatSolveAdd_SeqAIJ()
1193 for (j = 0; j < nz; j++) sum -= v[j] * tmp[vi[j]]; in MatSolveAdd_SeqAIJ()
1195 v += nz; in MatSolveAdd_SeqAIJ()
1196 vi += nz; in MatSolveAdd_SeqAIJ()
1203 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveAdd_SeqAIJ()
1205 for (j = 0; j < nz; j++) sum -= v[j] * tmp[vi[j]]; in MatSolveAdd_SeqAIJ()
1206 tmp[i] = sum * v[nz]; in MatSolveAdd_SeqAIJ()
1208 v += nz + 1; in MatSolveAdd_SeqAIJ()
1209 vi += nz + 1; in MatSolveAdd_SeqAIJ()
1217 PetscCall(PetscLogFlops(2.0 * a->nz)); in MatSolveAdd_SeqAIJ()
1227 PetscInt nz; in MatSolveTranspose_SeqAIJ_inplace() local
1250 nz = ai[i + 1] - diag[i] - 1; in MatSolveTranspose_SeqAIJ_inplace()
1253 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqAIJ_inplace()
1261 nz = diag[i] - ai[i]; in MatSolveTranspose_SeqAIJ_inplace()
1263 for (j = 0; j > -nz; j--) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqAIJ_inplace()
1275 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolveTranspose_SeqAIJ_inplace()
1285 PetscInt nz; in MatSolveTranspose_SeqAIJ() local
1309 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveTranspose_SeqAIJ()
1311 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatSolveTranspose_SeqAIJ()
1312 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqAIJ()
1320 nz = ai[i + 1] - ai[i]; in MatSolveTranspose_SeqAIJ()
1322 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqAIJ()
1334 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolveTranspose_SeqAIJ()
1344 PetscInt nz; in MatSolveTransposeAdd_SeqAIJ_inplace() local
1368 nz = ai[i + 1] - diag[i] - 1; in MatSolveTransposeAdd_SeqAIJ_inplace()
1371 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTransposeAdd_SeqAIJ_inplace()
1379 nz = diag[i] - ai[i]; in MatSolveTransposeAdd_SeqAIJ_inplace()
1381 for (j = 0; j > -nz; j--) tmp[vi[j]] -= s1 * v[j]; in MatSolveTransposeAdd_SeqAIJ_inplace()
1393 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolveTransposeAdd_SeqAIJ_inplace()
1403 PetscInt nz; in MatSolveTransposeAdd_SeqAIJ() local
1428 nz = adiag[i] - adiag[i + 1] - 1; in MatSolveTransposeAdd_SeqAIJ()
1430 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatSolveTransposeAdd_SeqAIJ()
1431 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTransposeAdd_SeqAIJ()
1439 nz = ai[i + 1] - ai[i]; in MatSolveTransposeAdd_SeqAIJ()
1441 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTransposeAdd_SeqAIJ()
1453 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolveTransposeAdd_SeqAIJ()
1477 PetscInt i, j, k = 0, nz, *bi, *bj, *bdiag; in MatILUFactorSymbolic_SeqAIJ_ilu0() local
1504 nz = adiag[i] - ai[i]; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1505 bi[i + 1] = bi[i] + nz; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1507 for (j = 0; j < nz; j++) bj[k++] = aj[j]; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1513 nz = ai[i + 1] - adiag[i] - 1; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1515 for (j = 0; j < nz; j++) bj[k++] = aj[j]; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1518 bdiag[i] = bdiag[i + 1] + nz + 1; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1692 b->maxnz = b->nz = bdiag[0] + 1; in MatILUFactorSymbolic_SeqAIJ()
1712 PetscInt k, jmin, jmax, *c2r, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ() local
1734 nz = ai[i + 1] - ai[i]; in MatCholeskyFactorNumeric_SeqAIJ()
1735 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatCholeskyFactorNumeric_SeqAIJ()
1761 nz = bi[k + 1] - bi[k]; in MatCholeskyFactorNumeric_SeqAIJ()
1763 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatCholeskyFactorNumeric_SeqAIJ()
1879 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1901 nz = ai[i + 1] - ai[i]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1902 for (j = 0; j < nz; j++) rs += PetscAbsScalar(v[j]); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1925 nz = bi[k + 1] - bi[k]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1927 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1974 nz = bi[k + 1] - jmin; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1976 for (j = 0; j < nz; j++) rs += PetscAbsScalar(rtmp[bcol[j]]); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
2228 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqAIJ()
2395 b->maxnz = b->nz = ui[am]; in MatCholeskyFactorSymbolic_SeqAIJ()
2427 PetscInt i, nz; in MatSolve_SeqAIJ_NaturalOrdering() local
2441 nz = ai[i + 1] - ai[i]; in MatSolve_SeqAIJ_NaturalOrdering()
2443 PetscSparseDenseMinusDot(sum, x, v, vi, nz); in MatSolve_SeqAIJ_NaturalOrdering()
2444 v += nz; in MatSolve_SeqAIJ_NaturalOrdering()
2445 vi += nz; in MatSolve_SeqAIJ_NaturalOrdering()
2453 nz = adiag[i] - adiag[i + 1] - 1; in MatSolve_SeqAIJ_NaturalOrdering()
2455 PetscSparseDenseMinusDot(sum, x, v, vi, nz); in MatSolve_SeqAIJ_NaturalOrdering()
2456 x[i] = sum * v[nz]; /* x[i]=aa[adiag[i]]*sum; v++; */ in MatSolve_SeqAIJ_NaturalOrdering()
2459 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolve_SeqAIJ_NaturalOrdering()
2470 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag, nz; in MatSolve_SeqAIJ() local
2494 nz = ai[i + 1] - ai[i]; in MatSolve_SeqAIJ()
2496 PetscSparseDenseMinusDot(sum, tmp, v, vi, nz); in MatSolve_SeqAIJ()
2498 v += nz; in MatSolve_SeqAIJ()
2499 vi += nz; in MatSolve_SeqAIJ()
2506 nz = adiag[i] - adiag[i + 1] - 1; in MatSolve_SeqAIJ()
2508 PetscSparseDenseMinusDot(sum, tmp, v, vi, nz); in MatSolve_SeqAIJ()
2509 x[c[i]] = tmp[i] = sum * v[nz]; /* v[nz] = aa[adiag[i]] */ in MatSolve_SeqAIJ()
2517 PetscCall(PetscLogFlops(2.0 * a->nz - A->cmap->n)); in MatSolve_SeqAIJ()