Home
last modified time | relevance | path

Searched refs:fact (Results 1 – 25 of 60) sorted by relevance

123

/petsc/src/ksp/pc/impls/factor/qr/
H A Dqr.c34 ((PC_Factor *)dir)->fact = pc->pmat; in PCSetUp_QR()
39 …!((PC_Factor *)dir)->fact) PetscCall(MatGetFactor(pc->pmat, ((PC_Factor *)dir)->solvertype, MAT_FA… in PCSetUp_QR()
40 …PetscCall(MatQRFactorSymbolic(((PC_Factor *)dir)->fact, pc->pmat, dir->col, &((PC_Factor *)dir)->i… in PCSetUp_QR()
41 PetscCall(MatGetInfo(((PC_Factor *)dir)->fact, MAT_LOCAL, &info)); in PCSetUp_QR()
44 …PetscCall(MatQRFactorSymbolic(((PC_Factor *)dir)->fact, pc->pmat, dir->col, &((PC_Factor *)dir)->i… in PCSetUp_QR()
45 PetscCall(MatGetInfo(((PC_Factor *)dir)->fact, MAT_LOCAL, &info)); in PCSetUp_QR()
48 PetscCall(MatFactorGetError(((PC_Factor *)dir)->fact, &err)); in PCSetUp_QR()
50 PetscCall(MatFactorGetError(((PC_Factor *)dir)->fact, &err)); in PCSetUp_QR()
56 PetscCall(MatQRFactorNumeric(((PC_Factor *)dir)->fact, pc->pmat, &((PC_Factor *)dir)->info)); in PCSetUp_QR()
57 PetscCall(MatFactorGetError(((PC_Factor *)dir)->fact, &err)); in PCSetUp_QR()
[all …]
/petsc/src/ksp/pc/impls/factor/
H A Dfactor.c16 PC_Factor *fact = (PC_Factor *)pc->data; in PCFactorSetDefaultOrdering_Factor() local
17 …PetscCall(MatSolverTypeGet(fact->solvertype, ((PetscObject)pc->pmat)->type_name, fact->factortype,… in PCFactorSetDefaultOrdering_Factor()
19 if (!fact->fact) { in PCFactorSetDefaultOrdering_Factor()
22 PetscCall(MatGetFactor(pc->pmat, fact->solvertype, fact->factortype, &fact->fact)); in PCFactorSetDefaultOrdering_Factor()
25 if (!fact->fact) PetscFunctionReturn(PETSC_SUCCESS); in PCFactorSetDefaultOrdering_Factor()
26 if (!fact->fact->assembled) { in PCFactorSetDefaultOrdering_Factor()
27 PetscCall(PetscStrcmp(fact->solvertype, fact->fact->solvertype, &flg)); in PCFactorSetDefaultOrdering_Factor()
30 PetscCall(MatGetFactor(pc->pmat, fact->solvertype, fact->factortype, &B)); in PCFactorSetDefaultOrdering_Factor()
31 PetscCall(MatHeaderReplace(fact->fact, &B)); in PCFactorSetDefaultOrdering_Factor()
34 if (!fact->ordering) { in PCFactorSetDefaultOrdering_Factor()
[all …]
H A Dfactimpl.c9 if (!icc->fact) PetscCall(MatGetFactor(pc->pmat, icc->solvertype, icc->factortype, &icc->fact)); in PCFactorSetUpMatSolverType_Factor()
10 …PetscCheck(icc->fact, PetscObjectComm((PetscObject)pc->pmat), PETSC_ERR_SUP, "MatFactor type %s no… in PCFactorSetUpMatSolverType_Factor()
168 …PetscCheck(ilu->fact, PetscObjectComm((PetscObject)pc), PETSC_ERR_ORDER, "Matrix not yet factored;… in PCFactorGetMatrix_Factor()
169 *mat = ilu->fact; in PCFactorGetMatrix_Factor()
181 if (lu->fact && lu->fact->assembled) { in PCFactorSetMatSolverType_Factor()
185 PetscCall(MatFactorGetSolverType(lu->fact, &ltype)); in PCFactorSetMatSolverType_Factor()
293 if (factor->fact) { in PCView_Factor()
294 PetscCall(MatFactorGetCanUseOrdering(factor->fact, &canuseordering)); in PCView_Factor()
298 if (!factor->fact->assembled) { in PCView_Factor()
299 … PetscCall(PetscViewerASCIIPrintf(viewer, " matrix solver type: %s\n", factor->fact->solvertype)); in PCView_Factor()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaijfact3.c10 PetscErrorCode MatSeqBAIJSetNumericFactorization(Mat fact, PetscBool natural) in MatSeqBAIJSetNumericFactorization() argument
14 switch (fact->rmap->bs) { in MatSeqBAIJSetNumericFactorization()
16 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_1; in MatSeqBAIJSetNumericFactorization()
19 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
22 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
25 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
28 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_5_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
31 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_6_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
34 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_7_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
38 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqBAIJ_9_NaturalOrdering; in MatSeqBAIJSetNumericFactorization()
[all …]
H A Dbaijfact2.c252 static PetscErrorCode MatILUFactorSymbolic_SeqBAIJ_ilu0(Mat fact, Mat A, IS isrow, IS iscol, const … in MatILUFactorSymbolic_SeqBAIJ_ilu0() argument
260 PetscCall(MatDuplicateNoCreate_SeqBAIJ(fact, A, MAT_DO_NOT_COPY_VALUES, PETSC_FALSE)); in MatILUFactorSymbolic_SeqBAIJ_ilu0()
261 b = (Mat_SeqBAIJ *)fact->data; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
269 fact->preallocated = PETSC_TRUE; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
270 fact->assembled = PETSC_TRUE; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
311 PetscErrorCode MatILUFactorSymbolic_SeqBAIJ(Mat fact, Mat A, IS isrow, IS iscol, const MatFactorInf… in MatILUFactorSymbolic_SeqBAIJ() argument
351 PetscCall(MatILUFactorSymbolic_SeqBAIJ_ilu0(fact, A, isrow, iscol, info)); in MatILUFactorSymbolic_SeqBAIJ()
352 PetscCall(MatSeqBAIJSetNumericFactorization(fact, both_identity)); in MatILUFactorSymbolic_SeqBAIJ()
354 fact->factortype = MAT_FACTOR_ILU; in MatILUFactorSymbolic_SeqBAIJ()
355 fact->info.factor_mallocs = 0; in MatILUFactorSymbolic_SeqBAIJ()
[all …]
/petsc/src/ksp/pc/impls/factor/icc/
H A Dicc.c22 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)icc)->fact, &canuseordering)); in PCSetUp_ICC()
27 …PetscCall(MatICCFactorSymbolic(((PC_Factor *)icc)->fact, pc->pmat, perm, &((PC_Factor *)icc)->info… in PCSetUp_ICC()
29 PetscCall(MatDestroy(&((PC_Factor *)icc)->fact)); in PCSetUp_ICC()
31 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)icc)->fact, &canuseordering)); in PCSetUp_ICC()
36 …PetscCall(MatICCFactorSymbolic(((PC_Factor *)icc)->fact, pc->pmat, perm, &((PC_Factor *)icc)->info… in PCSetUp_ICC()
38 PetscCall(MatGetInfo(((PC_Factor *)icc)->fact, MAT_LOCAL, &info)); in PCSetUp_ICC()
44 PetscCall(MatFactorGetError(((PC_Factor *)icc)->fact, &err)); in PCSetUp_ICC()
50 …PetscCall(MatCholeskyFactorNumeric(((PC_Factor *)icc)->fact, pc->pmat, &((PC_Factor *)icc)->info)); in PCSetUp_ICC()
51 PetscCall(MatFactorGetError(((PC_Factor *)icc)->fact, &err)); in PCSetUp_ICC()
60 PetscCall(MatFactorGetSolverType(((PC_Factor *)icc)->fact, &solverpackage)); in PCSetUp_ICC()
[all …]
/petsc/src/ksp/pc/impls/factor/cholesky/
H A Dcholesky.c58 ((PC_Factor *)dir)->fact = pc->pmat; in PCSetUp_Cholesky()
66 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)dir)->fact, &canuseordering)); in PCSetUp_Cholesky()
90 …PetscCall(MatCholeskyFactorSymbolic(((PC_Factor *)dir)->fact, pc->pmat, dir->row, &((PC_Factor *)d… in PCSetUp_Cholesky()
91 PetscCall(MatGetInfo(((PC_Factor *)dir)->fact, MAT_LOCAL, &info)); in PCSetUp_Cholesky()
94 PetscCall(MatDestroy(&((PC_Factor *)dir)->fact)); in PCSetUp_Cholesky()
99 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)dir)->fact, &canuseordering)); in PCSetUp_Cholesky()
120 …PetscCall(MatCholeskyFactorSymbolic(((PC_Factor *)dir)->fact, pc->pmat, dir->row, &((PC_Factor *)d… in PCSetUp_Cholesky()
121 PetscCall(MatGetInfo(((PC_Factor *)dir)->fact, MAT_LOCAL, &info)); in PCSetUp_Cholesky()
124 PetscCall(MatFactorGetError(((PC_Factor *)dir)->fact, &err)); in PCSetUp_Cholesky()
126 PetscCall(MatFactorClearError(((PC_Factor *)dir)->fact)); in PCSetUp_Cholesky()
[all …]
/petsc/src/ksp/pc/impls/factor/lu/
H A Dlu.c72 ((PC_Factor *)dir)->fact = pc->pmat; in PCSetUp_LU()
80 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)dir)->fact, &canuseordering)); in PCSetUp_LU()
91 …PetscCall(MatLUFactorSymbolic(((PC_Factor *)dir)->fact, pc->pmat, dir->row, dir->col, &((PC_Factor… in PCSetUp_LU()
92 PetscCall(MatGetInfo(((PC_Factor *)dir)->fact, MAT_LOCAL, &info)); in PCSetUp_LU()
97 PetscCall(MatDestroy(&((PC_Factor *)dir)->fact)); in PCSetUp_LU()
100 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)dir)->fact, &canuseordering)); in PCSetUp_LU()
114 …PetscCall(MatLUFactorSymbolic(((PC_Factor *)dir)->fact, pc->pmat, dir->row, dir->col, &((PC_Factor… in PCSetUp_LU()
115 PetscCall(MatGetInfo(((PC_Factor *)dir)->fact, MAT_LOCAL, &info)); in PCSetUp_LU()
118 PetscCall(MatFactorGetError(((PC_Factor *)dir)->fact, &err)); in PCSetUp_LU()
120 PetscCall(MatFactorClearError(((PC_Factor *)dir)->fact)); in PCSetUp_LU()
[all …]
/petsc/src/ksp/pc/impls/factor/ilu/
H A Dilu.c22 if (!ilu->hdr.inplace) PetscCall(MatDestroy(&((PC_Factor *)ilu)->fact)); in PCReset_ILU()
101 PetscCall(MatDestroy(&((PC_Factor *)ilu)->fact)); in PCSetUp_ILU()
117 ((PC_Factor *)ilu)->fact = pc->pmat; in PCSetUp_ILU()
126 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)ilu)->fact, &canuseordering)); in PCSetUp_ILU()
133 …PetscCall(MatILUFactorSymbolic(((PC_Factor *)ilu)->fact, pc->pmat, ilu->row, ilu->col, &((PC_Facto… in PCSetUp_ILU()
134 PetscCall(MatGetInfo(((PC_Factor *)ilu)->fact, MAT_LOCAL, &info)); in PCSetUp_ILU()
137 PetscCall(MatDestroy(&((PC_Factor *)ilu)->fact)); in PCSetUp_ILU()
142 PetscCall(MatFactorGetCanUseOrdering(((PC_Factor *)ilu)->fact, &canuseordering)); in PCSetUp_ILU()
153 …PetscCall(MatILUFactorSymbolic(((PC_Factor *)ilu)->fact, pc->pmat, ilu->row, ilu->col, &((PC_Facto… in PCSetUp_ILU()
154 PetscCall(MatGetInfo(((PC_Factor *)ilu)->fact, MAT_LOCAL, &info)); in PCSetUp_ILU()
[all …]
/petsc/src/mat/tests/
H A Dex126f.F9011 Mat A, fact
81 PetscCallA(MatGetFactor(A, MATSOLVERMUMPS, MAT_FACTOR_LU, fact, ierr))
84 PetscCallA(MatGetFactor(A, MATSOLVERPETSC, MAT_FACTOR_LU, fact, ierr))
86 PetscCallA(MatLUFactorSymbolic(fact, A, perm, iperm, info, ierr))
90 PetscCallA(MatLUFactorNumeric(fact, A, info, ierr))
91 PetscCallA(MatSolve(fact, b, x, ierr))
92 PetscCallA(MatDestroy(fact, ierr))
H A Dex117.c10 Mat mat, fact, B; in main() local
101 PetscCall(MatGetFactor(mat, MATSOLVERPETSC, MAT_FACTOR_CHOLESKY, &fact)); in main()
102 PetscCall(MatCholeskyFactorSymbolic(fact, mat, perm, &info)); in main()
103 PetscCall(MatCholeskyFactorNumeric(fact, mat, &info)); in main()
107 PetscCall(MatDestroy(&fact)); in main()
/petsc/src/mat/impls/aij/seq/bas/
H A Dbasfactor.c5 static PetscErrorCode MatICCFactorSymbolic_SeqAIJ_Bas(Mat fact, Mat A, IS perm, const MatFactorInfo… in MatICCFactorSymbolic_SeqAIJ_Bas() argument
62 b = (Mat_SeqSBAIJ *)fact->data; in MatICCFactorSymbolic_SeqAIJ_Bas()
83 fact->info.factor_mallocs = reallocs; in MatICCFactorSymbolic_SeqAIJ_Bas()
84 fact->info.fill_ratio_given = fill; in MatICCFactorSymbolic_SeqAIJ_Bas()
86 fact->info.fill_ratio_needed = (PetscReal)ui[am] / (PetscReal)ai[am]; in MatICCFactorSymbolic_SeqAIJ_Bas()
88 fact->info.fill_ratio_needed = 0.0; in MatICCFactorSymbolic_SeqAIJ_Bas()
/petsc/src/mat/impls/cdiagonal/
H A Dcdiagonal.c412 static PetscErrorCode MatFactorNumeric_ConstantDiagonal(Mat fact, Mat A, const MatFactorInfo *info) in MatFactorNumeric_ConstantDiagonal() argument
414 …nstantDiagonal *actx = (Mat_ConstantDiagonal *)A->data, *fctx = (Mat_ConstantDiagonal *)fact->data; in MatFactorNumeric_ConstantDiagonal()
417 if (actx->diag == 0.0) fact->factorerrortype = MAT_FACTOR_NUMERIC_ZEROPIVOT; in MatFactorNumeric_ConstantDiagonal()
418 else fact->factorerrortype = MAT_FACTOR_NOERROR; in MatFactorNumeric_ConstantDiagonal()
420 fact->ops->solve = MatMult_ConstantDiagonal; in MatFactorNumeric_ConstantDiagonal()
424 static PetscErrorCode MatFactorSymbolic_LU_ConstantDiagonal(Mat fact, Mat A, IS isrow, IS iscol, co… in MatFactorSymbolic_LU_ConstantDiagonal() argument
427 fact->ops->lufactornumeric = MatFactorNumeric_ConstantDiagonal; in MatFactorSymbolic_LU_ConstantDiagonal()
431 static PetscErrorCode MatFactorSymbolic_Cholesky_ConstantDiagonal(Mat fact, Mat A, IS isrow, const … in MatFactorSymbolic_Cholesky_ConstantDiagonal() argument
434 fact->ops->choleskyfactornumeric = MatFactorNumeric_ConstantDiagonal; in MatFactorSymbolic_Cholesky_ConstantDiagonal()
/petsc/src/ts/utils/dmplexlandau/
H A Dland_tensors.h132 PetscReal fact, u1u2, diff[3], udiff2, u12, u22, wsq, rsq, tt; in LandauTensor3DRelativistic() local
151 …fact = -rsq / (g1 * g2 * PetscSqrtReal(wsq)); /* flip sign. papers use du/dt = C, PETSc uses form … in LandauTensor3DRelativistic()
153 …for (j = 0; j < 3; ++j) U[i][j] = fact * (-diff[i] * diff[j] / wsq + (PetscSqrtReal(rsq) - 1.) * (… in LandauTensor3DRelativistic()
154 U[i][i] += fact; in LandauTensor3DRelativistic()
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact.c8 Mat_SeqSBAIJ *fact = (Mat_SeqSBAIJ *)F->data; in MatGetInertia_SeqSBAIJ() local
9 MatScalar *dd = fact->a; in MatGetInertia_SeqSBAIJ()
10 PetscInt mbs = fact->mbs, bs = F->rmap->bs, i, nneg_tmp, npos_tmp, *fi = fact->diag; in MatGetInertia_SeqSBAIJ()
25 if (PetscRealPart(dd[fact->i[i]]) > 0.0) npos_tmp++; in MatGetInertia_SeqSBAIJ()
26 else if (PetscRealPart(dd[fact->i[i]]) < 0.0) nneg_tmp++; in MatGetInertia_SeqSBAIJ()
220 PetscErrorCode MatCholeskyFactorSymbolic_SeqSBAIJ(Mat fact, Mat A, IS perm, const MatFactorInfo *in… in MatCholeskyFactorSymbolic_SeqSBAIJ() argument
239 PetscCall(MatCholeskyFactorSymbolic_SeqSBAIJ_inplace(fact, A, perm, info)); in MatCholeskyFactorSymbolic_SeqSBAIJ()
341 PetscCall(MatSeqSBAIJSetPreallocation(fact, bs, MAT_SKIP_ALLOCATION, NULL)); in MatCholeskyFactorSymbolic_SeqSBAIJ()
343 b = (Mat_SeqSBAIJ *)fact->data; in MatCholeskyFactorSymbolic_SeqSBAIJ()
364 fact->info.factor_mallocs = reallocs; in MatCholeskyFactorSymbolic_SeqSBAIJ()
[all …]
/petsc/src/tao/linesearch/impls/armijo/
H A Darmijo.c92 PetscReal fact, ref, gdx; in TaoLineSearchApply_Armijo() local
187 fact = armP->sigma; in TaoLineSearchApply_Armijo()
189 fact = armP->sigma * gdx; in TaoLineSearchApply_Armijo()
219 if (armP->nondescending && *f <= ref - ls->step * fact * ref) break; in TaoLineSearchApply_Armijo()
220 if (!armP->nondescending && *f <= ref + ls->step * fact) break; in TaoLineSearchApply_Armijo()
/petsc/src/ksp/pc/impls/telescope/
H A Dtelescope.c761 static PetscErrorCode PCTelescopeGetReductionFactor_Telescope(PC pc, PetscInt *fact) in PCTelescopeGetReductionFactor_Telescope() argument
766 if (fact) *fact = red->redfactor; in PCTelescopeGetReductionFactor_Telescope()
770 static PetscErrorCode PCTelescopeSetReductionFactor_Telescope(PC pc, PetscInt fact) in PCTelescopeSetReductionFactor_Telescope() argument
777 …heck(fact > 0, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONG, "Reduction factor of telesco… in PCTelescopeSetReductionFactor_Telescope()
778 …(fact <= size, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONG, "Reduction factor of telesco… in PCTelescopeSetReductionFactor_Telescope()
779 red->redfactor = fact; in PCTelescopeSetReductionFactor_Telescope()
884 PetscErrorCode PCTelescopeGetReductionFactor(PC pc, PetscInt *fact) in PCTelescopeGetReductionFactor() argument
887 PetscUseMethod(pc, "PCTelescopeGetReductionFactor_C", (PC, PetscInt *), (pc, fact)); in PCTelescopeGetReductionFactor()
907 PetscErrorCode PCTelescopeSetReductionFactor(PC pc, PetscInt fact) in PCTelescopeSetReductionFactor() argument
910 PetscTryMethod(pc, "PCTelescopeSetReductionFactor_C", (PC, PetscInt), (pc, fact)); in PCTelescopeSetReductionFactor()
/petsc/src/mat/impls/aij/seq/
H A Daijfact.c1473 PetscErrorCode MatILUFactorSymbolic_SeqAIJ_ilu0(Mat fact, Mat A, IS isrow, IS iscol, const MatFacto… in MatILUFactorSymbolic_SeqAIJ_ilu0() argument
1483 PetscCall(MatDuplicateNoCreate_SeqAIJ(fact, A, MAT_DO_NOT_COPY_VALUES, PETSC_FALSE)); in MatILUFactorSymbolic_SeqAIJ_ilu0()
1484 b = (Mat_SeqAIJ *)fact->data; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1521 fact->factortype = MAT_FACTOR_ILU; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1522 fact->info.factor_mallocs = 0; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1523 fact->info.fill_ratio_given = info->fill; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1524 fact->info.fill_ratio_needed = 1.0; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1525 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqAIJ; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1526 PetscCall(MatSeqAIJCheckInode_FactorLU(fact)); in MatILUFactorSymbolic_SeqAIJ_ilu0()
1528 b = (Mat_SeqAIJ *)fact->data; in MatILUFactorSymbolic_SeqAIJ_ilu0()
[all …]
/petsc/src/mat/interface/
H A Dmatrix.c3231 PetscErrorCode MatLUFactorSymbolic(Mat fact, Mat mat, IS row, IS col, const MatFactorInfo *info) in MatLUFactorSymbolic() argument
3236 PetscValidHeaderSpecific(fact, MAT_CLASSID, 1); in MatLUFactorSymbolic()
3241 PetscValidType(fact, 1); in MatLUFactorSymbolic()
3251 if (!fact->trivialsymbolic) PetscCall(PetscLogEventBegin(MAT_LUFactorSymbolic, mat, row, col, 0)); in MatLUFactorSymbolic()
3252 PetscUseTypeMethod(fact, lufactorsymbolic, mat, row, col, info); in MatLUFactorSymbolic()
3253 if (!fact->trivialsymbolic) PetscCall(PetscLogEventEnd(MAT_LUFactorSymbolic, mat, row, col, 0)); in MatLUFactorSymbolic()
3254 PetscCall(PetscObjectStateIncrease((PetscObject)fact)); in MatLUFactorSymbolic()
3284 PetscErrorCode MatLUFactorNumeric(Mat fact, Mat mat, const MatFactorInfo *info) in MatLUFactorNumeric() argument
3289 PetscValidHeaderSpecific(fact, MAT_CLASSID, 1); in MatLUFactorNumeric()
3291 PetscValidType(fact, 1); in MatLUFactorNumeric()
[all …]
/petsc/src/snes/impls/vi/
H A Dvi.c101 PetscInt i, n, act[2] = {0, 0}, fact[2], N; in SNESMonitorVI() local
133 PetscCallMPI(MPIU_Allreduce(act, fact, 2, MPIU_INT, MPI_SUM, PetscObjectComm((PetscObject)snes))); in SNESMonitorVI()
138 if (snes->ntruebounds) tmp = ((double)(fact[0] + fact[1])) / ((double)snes->ntruebounds); in SNESMonitorVI()
140 …t of bounded %g\n", its, (double)fnorm, fact[0], fact_bound[0], fact[1], fact_bound[1], ((double)(… in SNESMonitorVI()
/petsc/src/tao/linesearch/impls/owarmijo/
H A Dowarmijo.c116 PetscReal fact, ref, gdx; in TaoLineSearchApply_OWArmijo() local
126 fact = 0.0; in TaoLineSearchApply_OWArmijo()
204 if (armP->nondescending) fact = armP->sigma; in TaoLineSearchApply_OWArmijo()
245 if (armP->nondescending && *f <= ref - ls->step * fact * ref) break; in TaoLineSearchApply_OWArmijo()
/petsc/src/ts/tests/
H A Dex30.c313 …PetscReal v2 = 0, fact = (dim == 2) ? 2.0 * PETSC_PI * coords[p * dim + 0] : 1, w = fact * wp[p] *… in gridToParticles_private() local
333 …const PetscReal fact = dim == 2 ? 2.0 * PETSC_PI * coords[p * dim + 0] : 1, w = fact * wp[p] * ctx… in gridToParticles_private() local
408 …PetscReal v2 = 0, fact = (dim == 2) ? 2.0 * PETSC_PI * coords[pp * dim + 0] : 1, w = fact * wp[pp]… in PostStep() local
429 …PetscReal fact = (dim == 2) ? 2.0 * PETSC_PI * coords[pp * dim + 0] : 1, w = fact * wp[pp] * ctx->… in PostStep() local
635 …PetscReal v2 = 0, fact = dim == 2 ? 2.0 * PETSC_PI * x[0] : 1, w = fact * wp_t[grid][tid][pp] * ct… in go() local
675 …const PetscReal fact = dim == 2 ? 2.0 * PETSC_PI * xx_t[grid][tid][pp] : 1, w = fact * ctx->n_0 * … in go() local
791 …PetscReal v2 = 0, fact = (dim == 2) ? 2.0 * PETSC_PI * coords[pp * dim + 0] : 1, w = fact * wp[pp]… in go() local
812 …PetscReal fact = (dim == 2) ? 2.0 * PETSC_PI * coords[pp * dim + 0] : 1, w = fact * wp[pp] * ctx->… in go() local
/petsc/src/mat/impls/aij/seq/seqhipsparse/
H A Daijhipsparse.hip.cxx1247 static PetscErrorCode MatSolve_SeqAIJHIPSPARSE_ILU0(Mat fact, Vec b, Vec x) in MatSolve_SeqAIJHIPSPARSE_ILU0() argument
1249 Mat_SeqAIJHIPSPARSETriFactors *fs = (Mat_SeqAIJHIPSPARSETriFactors *)fact->spptr; in MatSolve_SeqAIJHIPSPARSE_ILU0()
1250 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)fact->data; in MatSolve_SeqAIJHIPSPARSE_ILU0()
1282 PetscCall(PetscLogGpuFlops(2.0 * aij->nz - fact->rmap->n)); in MatSolve_SeqAIJHIPSPARSE_ILU0()
1286 static PetscErrorCode MatSolveTranspose_SeqAIJHIPSPARSE_ILU0(Mat fact, Vec b, Vec x) in MatSolveTranspose_SeqAIJHIPSPARSE_ILU0() argument
1288 Mat_SeqAIJHIPSPARSETriFactors *fs = (Mat_SeqAIJHIPSPARSETriFactors *)fact->spptr; in MatSolveTranspose_SeqAIJHIPSPARSE_ILU0()
1289 Mat_SeqAIJ *aij = (Mat_SeqAIJ *)fact->data; in MatSolveTranspose_SeqAIJHIPSPARSE_ILU0()
1339 PetscCall(PetscLogGpuFlops(2.0 * aij->nz - fact->rmap->n)); in MatSolveTranspose_SeqAIJHIPSPARSE_ILU0()
1343 static PetscErrorCode MatILUFactorNumeric_SeqAIJHIPSPARSE_ILU0(Mat fact, Mat A, const MatFactorInfo… in MatILUFactorNumeric_SeqAIJHIPSPARSE_ILU0() argument
1345 Mat_SeqAIJHIPSPARSETriFactors *fs = (Mat_SeqAIJHIPSPARSETriFactors *)fact->spptr; in MatILUFactorNumeric_SeqAIJHIPSPARSE_ILU0()
[all …]
/petsc/src/mat/impls/dense/seq/
H A Ddense.c819 static PetscErrorCode MatLUFactorNumeric_SeqDense(Mat fact, Mat A, const MatFactorInfo *info) in MatLUFactorNumeric_SeqDense() argument
822 PetscCall(MatDuplicateNoCreate_SeqDense(fact, A, MAT_COPY_VALUES)); in MatLUFactorNumeric_SeqDense()
823 PetscUseTypeMethod(fact, lufactor, NULL, NULL, info); in MatLUFactorNumeric_SeqDense()
827 PetscErrorCode MatLUFactorSymbolic_SeqDense(Mat fact, Mat A, IS row, IS col, PETSC_UNUSED const Mat… in MatLUFactorSymbolic_SeqDense() argument
830 fact->preallocated = PETSC_TRUE; in MatLUFactorSymbolic_SeqDense()
831 fact->assembled = PETSC_TRUE; in MatLUFactorSymbolic_SeqDense()
832 fact->ops->lufactornumeric = MatLUFactorNumeric_SeqDense; in MatLUFactorSymbolic_SeqDense()
897 static PetscErrorCode MatCholeskyFactorNumeric_SeqDense(Mat fact, Mat A, const MatFactorInfo *info) in MatCholeskyFactorNumeric_SeqDense() argument
900 PetscCall(MatDuplicateNoCreate_SeqDense(fact, A, MAT_COPY_VALUES)); in MatCholeskyFactorNumeric_SeqDense()
901 PetscUseTypeMethod(fact, choleskyfactor, NULL, info); in MatCholeskyFactorNumeric_SeqDense()
[all …]
/petsc/src/ts/event/
H A Dtsevent.c402 PetscInt fact = 2; in TSEventRecorderResize() local
406 PetscCall(PetscMalloc1(fact * event->recsize, &time)); in TSEventRecorderResize()
407 PetscCall(PetscMalloc1(fact * event->recsize, &stepnum)); in TSEventRecorderResize()
408 PetscCall(PetscMalloc1(fact * event->recsize, &nevents)); in TSEventRecorderResize()
409 PetscCall(PetscMalloc1(fact * event->recsize, &eventidx)); in TSEventRecorderResize()
410 …for (PetscInt i = 0; i < fact * event->recsize; i++) PetscCall(PetscMalloc1(event->nevents, &event… in TSEventRecorderResize()
432 event->recsize *= fact; in TSEventRecorderResize()

123