| /petsc/src/snes/impls/richardson/ |
| H A D | snesrichardson.c | 3 static PetscErrorCode SNESDestroy_NRichardson(SNES snes) in SNESDestroy_NRichardson() argument
6 PetscCall(PetscFree(snes->data)); in SNESDestroy_NRichardson()
10 static PetscErrorCode SNESSetUp_NRichardson(SNES snes) in SNESSetUp_NRichardson() argument
13 …PetscCheck(snes->npcside != PC_RIGHT, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONGSTATE, "NRichardson only… in SNESSetUp_NRichardson()
14 if (snes->functype == SNES_FUNCTION_DEFAULT) snes->functype = SNES_FUNCTION_UNPRECONDITIONED; in SNESSetUp_NRichardson()
18 static PetscErrorCode SNESSetFromOptions_NRichardson(SNES snes, PetscOptionItems PetscOptionsObject) in SNESSetFromOptions_NRichardson() argument
26 static PetscErrorCode SNESSolve_NRichardson(SNES snes) in SNESSolve_NRichardson() argument
34 …snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObject)snes), PE… in SNESSolve_NRichardson()
36 snes->reason = SNES_CONVERGED_ITERATING; in SNESSolve_NRichardson()
38 maxits = snes->max_its; /* maximum number of iterations */ in SNESSolve_NRichardson()
[all …]
|
| /petsc/src/snes/impls/ls/ |
| H A D | ls.c | 46 static PetscErrorCode SNESNEWTONLSCheckLocalMin_Private(SNES snes, Mat A, Vec F, PetscReal fnorm, P… in SNESNEWTONLSCheckLocalMin_Private() argument
55 PetscCall(SNESGetObjective(snes, &objective, NULL)); in SNESNEWTONLSCheckLocalMin_Private()
63 …PetscCall(PetscInfo(snes, "|| J^T F|| %14.12e near zero implies found a local minimum\n", (double)… in SNESNEWTONLSCheckLocalMin_Private()
77 …PetscCall(PetscInfo(snes, "(F^T J random)/(|| F ||*||random|| %14.12e near zero implies found a lo… in SNESNEWTONLSCheckLocalMin_Private()
88 static PetscErrorCode SNESNEWTONLSCheckResidual_Private(SNES snes, Mat A, Vec F, Vec X) in SNESNEWTONLSCheckResidual_Private() argument
96 PetscCall(SNESGetObjective(snes, &objective, NULL)); in SNESNEWTONLSCheckResidual_Private()
109 …if (a1 != 0.0) PetscCall(PetscInfo(snes, "||J^T(F-Ax)||/||F-AX|| %14.12e near zero implies inconsi… in SNESNEWTONLSCheckResidual_Private()
125 static PetscErrorCode SNESSolve_NEWTONLS(SNES snes) in SNESSolve_NEWTONLS() argument
139 …snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObject)snes), PE… in SNESSolve_NEWTONLS()
141 snes->numFailures = 0; in SNESSolve_NEWTONLS()
[all …]
|
| /petsc/src/snes/impls/ksponly/ |
| H A D | ksponly.c | 7 static PetscErrorCode SNESSolve_KSPONLY(SNES snes) in SNESSolve_KSPONLY() argument
9 SNES_KSPONLY *ksponly = (SNES_KSPONLY *)snes->data; in SNESSolve_KSPONLY()
15 …snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObject)snes), PE… in SNESSolve_KSPONLY()
17 snes->numFailures = 0; in SNESSolve_KSPONLY()
18 snes->numLinearSolveFailures = 0; in SNESSolve_KSPONLY()
19 snes->reason = SNES_CONVERGED_ITERATING; in SNESSolve_KSPONLY()
20 snes->iter = 0; in SNESSolve_KSPONLY()
21 snes->norm = 0.0; in SNESSolve_KSPONLY()
23 X = snes->vec_sol; in SNESSolve_KSPONLY()
24 F = snes->vec_func; in SNESSolve_KSPONLY()
[all …]
|
| /petsc/src/snes/interface/ |
| H A D | snes.c | 36 PetscErrorCode SNESSetErrorIfNotConverged(SNES snes, PetscBool flg) in SNESSetErrorIfNotConverged() argument
39 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESSetErrorIfNotConverged()
40 PetscValidLogicalCollectiveBool(snes, flg, 2); in SNESSetErrorIfNotConverged()
41 snes->errorifnotconverged = flg; in SNESSetErrorIfNotConverged()
60 PetscErrorCode SNESGetErrorIfNotConverged(SNES snes, PetscBool *flag) in SNESGetErrorIfNotConverged() argument
63 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESGetErrorIfNotConverged()
65 *flag = snes->errorifnotconverged; in SNESGetErrorIfNotConverged()
86 PetscErrorCode SNESSetAlwaysComputesFinalResidual(SNES snes, PetscBool flg) in SNESSetAlwaysComputesFinalResidual() argument
89 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESSetAlwaysComputesFinalResidual()
90 snes->alwayscomputesfinalresidual = flg; in SNESSetAlwaysComputesFinalResidual()
[all …]
|
| /petsc/src/snes/impls/gs/ |
| H A D | snesgs.c | 33 PetscErrorCode SNESNGSSetTolerances(SNES snes, PetscReal abstol, PetscReal rtol, PetscReal stol, Pe… in SNESNGSSetTolerances() argument
35 SNES_NGS *gs = (SNES_NGS *)snes->data; in SNESNGSSetTolerances()
38 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESNGSSetTolerances()
41 …PetscCheck(abstol >= 0.0, PetscObjectComm((PetscObject)snes), PETSC_ERR_ARG_OUTOFRANGE, "Absolute … in SNESNGSSetTolerances()
45 …PetscCheck(rtol >= 0.0 && rtol < 1.0, PetscObjectComm((PetscObject)snes), PETSC_ERR_ARG_OUTOFRANGE… in SNESNGSSetTolerances()
49 …PetscCheck(stol >= 0.0, PetscObjectComm((PetscObject)snes), PETSC_ERR_ARG_OUTOFRANGE, "Step tolera… in SNESNGSSetTolerances()
53 …PetscCheck(maxit >= 0, PetscObjectComm((PetscObject)snes), PETSC_ERR_ARG_OUTOFRANGE, "Maximum numb… in SNESNGSSetTolerances()
79 PetscErrorCode SNESNGSGetTolerances(SNES snes, PetscReal *atol, PetscReal *rtol, PetscReal *stol, P… in SNESNGSGetTolerances() argument
81 SNES_NGS *gs = (SNES_NGS *)snes->data; in SNESNGSGetTolerances()
84 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESNGSGetTolerances()
[all …]
|
| /petsc/src/snes/impls/qn/ |
| H A D | qn.c | 21 static PetscErrorCode SNESQNGetMatrix_Private(SNES snes, Mat *B) in SNESQNGetMatrix_Private() argument
23 SNES_QN *qn = (SNES_QN *)snes->data; in SNESQNGetMatrix_Private()
28 PetscCall(MatCreate(PetscObjectComm((PetscObject)snes), &qn->B)); in SNESQNGetMatrix_Private()
29 PetscCall(SNESGetOptionsPrefix(snes, &optionsprefix)); in SNESQNGetMatrix_Private()
62 static PetscErrorCode SNESSolve_QN(SNES snes) in SNESSolve_QN() argument
64 SNES_QN *qn = (SNES_QN *)snes->data; in SNESSolve_QN()
78 …snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObject)snes), PE… in SNESSolve_QN()
81 F = snes->vec_func; /* residual vector */ in SNESSolve_QN()
82 Y = snes->vec_sol_update; /* search direction generated by J^-1D*/ in SNESSolve_QN()
83 X = snes->vec_sol; /* solution vector */ in SNESSolve_QN()
[all …]
|
| /petsc/src/snes/impls/ncg/ |
| H A D | snesncg.c | 4 static PetscErrorCode SNESDestroy_NCG(SNES snes) in SNESDestroy_NCG() argument
7 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNCGSetType_C", NULL)); in SNESDestroy_NCG()
8 PetscCall(PetscFree(snes->data)); in SNESDestroy_NCG()
12 static PetscErrorCode SNESSetUp_NCG(SNES snes) in SNESSetUp_NCG() argument
15 PetscCall(SNESSetWorkVecs(snes, 2)); in SNESSetUp_NCG()
16 …PetscCheck(snes->npcside != PC_RIGHT, PetscObjectComm((PetscObject)snes), PETSC_ERR_ARG_WRONGSTATE… in SNESSetUp_NCG()
17 if (snes->functype == SNES_FUNCTION_DEFAULT) snes->functype = SNES_FUNCTION_UNPRECONDITIONED; in SNESSetUp_NCG()
25 SNES snes; in SNESLineSearchApply_NCGLinear() local
29 PetscCall(SNESLineSearchGetSNES(linesearch, &snes)); in SNESLineSearchApply_NCGLinear()
38 if (!snes->jacobian) PetscCall(SNESSetUpMatrices(snes)); in SNESLineSearchApply_NCGLinear()
[all …]
|
| /petsc/src/snes/impls/ngmres/ |
| H A D | snesngmres.c | 8 PetscErrorCode SNESReset_NGMRES(SNES snes) in SNESReset_NGMRES() argument
10 SNES_NGMRES *ngmres = (SNES_NGMRES *)snes->data; in SNESReset_NGMRES()
19 PetscErrorCode SNESDestroy_NGMRES(SNES snes) in SNESDestroy_NGMRES() argument
21 SNES_NGMRES *ngmres = (SNES_NGMRES *)snes->data; in SNESDestroy_NGMRES()
24 PetscCall(SNESReset_NGMRES(snes)); in SNESDestroy_NGMRES()
31 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNGMRESSetSelectType_C", NULL)); in SNESDestroy_NGMRES()
32 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNGMRESSetRestartType_C", NULL)); in SNESDestroy_NGMRES()
33 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNGMRESSetRestartFmRise_C", NULL)); in SNESDestroy_NGMRES()
34 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNGMRESGetRestartFmRise_C", NULL)); in SNESDestroy_NGMRES()
35 PetscCall(PetscFree(snes->data)); in SNESDestroy_NGMRES()
[all …]
|
| H A D | anderson.c | 3 static PetscErrorCode SNESSetFromOptions_Anderson(SNES snes, PetscOptionItems PetscOptionsObject) in SNESSetFromOptions_Anderson() argument
5 SNES_NGMRES *ngmres = (SNES_NGMRES *)snes->data; in SNESSetFromOptions_Anderson()
16 if (monitor) ngmres->monitor = PETSC_VIEWER_STDOUT_(PetscObjectComm((PetscObject)snes)); in SNESSetFromOptions_Anderson()
21 static PetscErrorCode SNESSolve_Anderson(SNES snes) in SNESSolve_Anderson() argument
23 SNES_NGMRES *ngmres = (SNES_NGMRES *)snes->data; in SNESSolve_Anderson()
39 …snes->xl && !snes->xu && !snes->ops->computevariablebounds, PetscObjectComm((PetscObject)snes), PE… in SNESSolve_Anderson()
43 snes->reason = SNES_CONVERGED_ITERATING; in SNESSolve_Anderson()
44 X = snes->vec_sol; in SNESSolve_Anderson()
45 F = snes->vec_func; in SNESSolve_Anderson()
46 B = snes->vec_rhs; in SNESSolve_Anderson()
[all …]
|
| /petsc/src/snes/impls/al/ |
| H A D | al.c | 40 static PetscErrorCode SNESNewtonALCheckArcLength(SNES snes, Vec XStep, PetscReal lambdaStep, PetscR… in SNESNewtonALCheckArcLength() argument
43 SNES_NEWTONAL *al = (SNES_NEWTONAL *)snes->data; in SNESNewtonALCheckArcLength()
50 …if (arcLengthError > PETSC_SQRT_MACHINE_EPSILON) PetscCall(PetscInfo(snes, "Arc length differs fro… in SNESNewtonALCheckArcLength()
64 static PetscErrorCode SNESNewtonALSetCorrectionType_NEWTONAL(SNES snes, SNESNewtonALCorrectionType … in SNESNewtonALSetCorrectionType_NEWTONAL() argument
66 SNES_NEWTONAL *al = (SNES_NEWTONAL *)snes->data; in SNESNewtonALSetCorrectionType_NEWTONAL()
89 PetscErrorCode SNESNewtonALSetCorrectionType(SNES snes, SNESNewtonALCorrectionType ctype) in SNESNewtonALSetCorrectionType() argument
92 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESNewtonALSetCorrectionType()
93 PetscValidLogicalCollectiveEnum(snes, ctype, 2); in SNESNewtonALSetCorrectionType()
94 …PetscTryMethod(snes, "SNESNewtonALSetCorrectionType_C", (SNES, SNESNewtonALCorrectionType), (snes,… in SNESNewtonALSetCorrectionType()
98 static PetscErrorCode SNESNewtonALSetFunction_NEWTONAL(SNES snes, SNESFunctionFn *func, PetscCtx ct… in SNESNewtonALSetFunction_NEWTONAL() argument
[all …]
|
| /petsc/src/snes/impls/composite/ |
| H A D | snescomposite.c | 11 SNES snes; member
47 static PetscErrorCode SNESCompositeApply_Multiplicative(SNES snes, Vec X, Vec B, Vec F, PetscReal *… in SNESCompositeApply_Multiplicative() argument
49 SNES_Composite *jac = (SNES_Composite *)snes->data; in SNESCompositeApply_Multiplicative()
55 …PetscCheck(next, PetscObjectComm((PetscObject)snes), PETSC_ERR_ARG_WRONGSTATE, "No composite SNESe… in SNESCompositeApply_Multiplicative()
56 if (snes->normschedule == SNES_NORM_ALWAYS) PetscCall(SNESSetInitialFunction(next->snes, F)); in SNESCompositeApply_Multiplicative()
57 PetscCall(SNESSolve(next->snes, B, X)); in SNESCompositeApply_Multiplicative()
58 PetscCall(SNESGetConvergedReason(next->snes, &reason)); in SNESCompositeApply_Multiplicative()
61 if (jac->innerFailures >= snes->maxFailures) { in SNESCompositeApply_Multiplicative()
62 snes->reason = SNES_DIVERGED_INNER; in SNESCompositeApply_Multiplicative()
69 if (next->snes->npcside == PC_RIGHT && next->snes->normschedule != SNES_NORM_NONE) { in SNESCompositeApply_Multiplicative()
[all …]
|
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 4 SNES snes; member
13 static PetscErrorCode SNESNewtonTRSetTolerances_TR(SNES snes, PetscReal delta_min, PetscReal delta_… in SNESNewtonTRSetTolerances_TR() argument
15 SNES_NEWTONTR *tr = (SNES_NEWTONTR *)snes->data; in SNESNewtonTRSetTolerances_TR()
27 static PetscErrorCode SNESNewtonTRGetTolerances_TR(SNES snes, PetscReal *delta_min, PetscReal *delt… in SNESNewtonTRGetTolerances_TR() argument
29 SNES_NEWTONTR *tr = (SNES_NEWTONTR *)snes->data; in SNESNewtonTRGetTolerances_TR()
38 static PetscErrorCode SNESComputeJacobian_MATLMVM(SNES snes, Vec X, Mat J, Mat B, void *dummy) in SNESComputeJacobian_MATLMVM() argument
42 PetscCall(MatLMVMUpdate(B, X, snes->vec_func)); in SNESComputeJacobian_MATLMVM()
47 PetscCall(MatLMVMUpdate(J, X, snes->vec_func)); in SNESComputeJacobian_MATLMVM()
57 SNES snes = ctx->snes; in SNESTR_KSPConverged_Private() local
58 SNES_NEWTONTR *neP = (SNES_NEWTONTR *)snes->data; in SNESTR_KSPConverged_Private()
[all …]
|
| /petsc/src/snes/impls/ms/ |
| H A D | ms.c | 236 static PetscErrorCode SNESMSStep_3Sstar(SNES snes, Vec X, Vec F) in SNESMSStep_3Sstar() argument
238 SNES_MS *ms = (SNES_MS *)snes->data; in SNESMSStep_3Sstar()
242 …Vec S1 = X, S2 = snes->work[1], S3 = snes->work[2], Y = snes->work[0], S1copy = snes-… in SNESMSStep_3Sstar()
252 if (i > 0) PetscCall(SNESComputeFunction(snes, S1, F)); in SNESMSStep_3Sstar()
253 PetscCall(KSPSolve(snes->ksp, F, Y)); in SNESMSStep_3Sstar()
264 static PetscErrorCode SNESMSStep_Basic(SNES snes, Vec X, Vec F) in SNESMSStep_Basic() argument
266 SNES_MS *ms = (SNES_MS *)snes->data; in SNESMSStep_Basic()
270 Vec X0 = snes->work[0]; in SNESMSStep_Basic()
275 if (i > 0) PetscCall(SNESComputeFunction(snes, X, F)); in SNESMSStep_Basic()
276 PetscCall(KSPSolve(snes->ksp, F, X)); in SNESMSStep_Basic()
[all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_snes.py | 12 def __call__(self, snes, x, f): argument
19 def __call__(self, snes, x, J, P): argument
30 def __call__(self, snes, x, f): argument
31 load = snes.getNewtonALLoadParameter()
38 def __call__(self, snes, x, f): argument
51 snes = PETSc.SNES()
52 snes.create(PETSc.COMM_SELF)
54 snes.setType(self.SNES_TYPE)
55 self.snes = snes
58 self.snes = None
[all …]
|
| /petsc/src/snes/impls/vi/ss/ |
| H A D | viss.c | 53 PetscErrorCode SNESVIComputeFunction(SNES snes, Vec X, Vec phi, void *functx) in SNESVIComputeFunction() argument
55 SNES_VINEWTONSSLS *vi = (SNES_VINEWTONSSLS *)snes->data; in SNESVIComputeFunction()
56 Vec Xl = snes->xl, Xu = snes->xu, F = snes->vec_func; in SNESVIComputeFunction()
62 PetscCall((*vi->computeuserfunction)(snes, X, F, functx)); in SNESVIComputeFunction()
96 static PetscErrorCode SNESVIComputeBsubdifferentialVectors(SNES snes, Vec X, Vec F, Mat jac, Vec Da… in SNESVIComputeBsubdifferentialVectors() argument
104 PetscCall(VecGetArray(snes->xl, &l)); in SNESVIComputeBsubdifferentialVectors()
105 PetscCall(VecGetArray(snes->xu, &u)); in SNESVIComputeBsubdifferentialVectors()
135 PetscCall(VecRestoreArray(snes->xl, &l)); in SNESVIComputeBsubdifferentialVectors()
136 PetscCall(VecRestoreArray(snes->xu, &u)); in SNESVIComputeBsubdifferentialVectors()
195 static PetscErrorCode SNESSolve_VINEWTONSSLS(SNES snes) in SNESSolve_VINEWTONSSLS() argument
[all …]
|
| /petsc/src/snes/impls/fas/ |
| H A D | fas.c | 6 static PetscErrorCode SNESReset_FAS(SNES snes) in SNESReset_FAS() argument
8 SNES_FAS *fas = (SNES_FAS *)snes->data; in SNESReset_FAS()
23 static PetscErrorCode SNESDestroy_FAS(SNES snes) in SNESDestroy_FAS() argument
25 SNES_FAS *fas = (SNES_FAS *)snes->data; in SNESDestroy_FAS()
29 PetscCall(SNESReset_FAS(snes)); in SNESDestroy_FAS()
35 static PetscErrorCode SNESFASSetUpLineSearch_Private(SNES snes, SNES smooth) in SNESFASSetUpLineSearch_Private() argument
44 if (!snes->linesearch) PetscFunctionReturn(PETSC_SUCCESS); in SNESFASSetUpLineSearch_Private()
45 PetscCall(SNESGetLineSearch(snes, &linesearch)); in SNESFASSetUpLineSearch_Private()
55 static PetscErrorCode SNESFASCycleSetUpSmoother_Private(SNES snes, SNES smooth) in SNESFASCycleSetUpSmoother_Private() argument
57 SNES_FAS *fas = (SNES_FAS *)snes->data; in SNESFASCycleSetUpSmoother_Private()
[all …]
|
| H A D | fasfunc.c | 16 PetscErrorCode SNESFASSetType(SNES snes, SNESFASType fastype) in SNESFASSetType() argument
21 PetscValidHeaderSpecificType(snes, SNES_CLASSID, 1, SNESFAS); in SNESFASSetType()
22 PetscValidLogicalCollectiveEnum(snes, fastype, 2); in SNESFASSetType()
23 fas = (SNES_FAS *)snes->data; in SNESFASSetType()
44 PetscErrorCode SNESFASGetType(SNES snes, SNESFASType *fastype) in SNESFASGetType() argument
49 PetscValidHeaderSpecificType(snes, SNES_CLASSID, 1, SNESFAS); in SNESFASGetType()
51 fas = (SNES_FAS *)snes->data; in SNESFASGetType()
74 PetscErrorCode SNESFASSetLevels(SNES snes, PetscInt levels, MPI_Comm *comms) in SNESFASSetLevels() argument
84 PetscValidHeaderSpecificType(snes, SNES_CLASSID, 1, SNESFAS); in SNESFASSetLevels()
85 fas = (SNES_FAS *)snes->data; in SNESFASSetLevels()
[all …]
|
| /petsc/src/snes/impls/vi/ |
| H A D | vi.c | 32 PetscErrorCode SNESVISetComputeVariableBounds(SNES snes, PetscErrorCode (*compute)(SNES snes, Vec l… in SNESVISetComputeVariableBounds() argument
37 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESVISetComputeVariableBounds()
38 PetscCall(PetscObjectQueryFunction((PetscObject)snes, "SNESVISetComputeVariableBounds_C", &f)); in SNESVISetComputeVariableBounds()
39 …if (f) PetscUseMethod(snes, "SNESVISetComputeVariableBounds_C", (SNES, PetscErrorCode (*)(SNES, Ve… in SNESVISetComputeVariableBounds()
40 else PetscCall(SNESVISetComputeVariableBounds_VI(snes, compute)); in SNESVISetComputeVariableBounds()
44 PetscErrorCode SNESVISetComputeVariableBounds_VI(SNES snes, SNESVIComputeVariableBoundsFn *compute) in SNESVISetComputeVariableBounds_VI() argument
47 snes->ops->computevariablebounds = compute; in SNESVISetComputeVariableBounds_VI()
51 static PetscErrorCode SNESVIMonitorResidual(SNES snes, PetscInt its, PetscReal fgnorm, PetscViewerA… in SNESVIMonitorResidual() argument
58 PetscCall(SNESGetFunction(snes, &F, NULL, NULL)); in SNESVIMonitorResidual()
59 PetscCall(SNESGetSolution(snes, &X)); in SNESVIMonitorResidual()
[all …]
|
| /petsc/src/snes/impls/shell/ |
| H A D | snesshell.c | 25 PetscErrorCode SNESShellSetSolve(SNES snes, PetscErrorCode (*solve)(SNES snes, Vec xout)) in SNESShellSetSolve() argument
28 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESShellSetSolve()
29 PetscTryMethod(snes, "SNESShellSetSolve_C", (SNES, PetscErrorCode (*)(SNES, Vec)), (snes, solve)); in SNESShellSetSolve()
33 static PetscErrorCode SNESDestroy_Shell(SNES snes) in SNESDestroy_Shell() argument
36 PetscCall(PetscFree(snes->data)); in SNESDestroy_Shell()
61 PetscErrorCode SNESShellGetContext(SNES snes, PetscCtxRt ctx) in SNESShellGetContext() argument
66 PetscValidHeaderSpecific(snes, SNES_CLASSID, 1); in SNESShellGetContext()
68 PetscCall(PetscObjectTypeCompare((PetscObject)snes, SNESSHELL, &flg)); in SNESShellGetContext()
70 else *(void **)ctx = ((SNES_Shell *)snes->data)->ctx; in SNESShellGetContext()
87 PetscErrorCode SNESShellSetContext(SNES snes, PetscCtx ctx) in SNESShellSetContext() argument
[all …]
|
| /petsc/src/snes/impls/ntrdc/ |
| H A D | ntrdc.c | 4 SNES snes; member
15 static PetscErrorCode SNESNewtonTRSetTolerances_TRDC(SNES snes, PetscReal delta_min, PetscReal delt… in SNESNewtonTRSetTolerances_TRDC() argument
17 SNES_NEWTONTRDC *tr = (SNES_NEWTONTRDC *)snes->data; in SNESNewtonTRSetTolerances_TRDC()
32 SNES snes = ctx->snes; in SNESTRDC_KSPConverged_Private() local
33 SNES_NEWTONTRDC *neP = (SNES_NEWTONTRDC *)snes->data; in SNESTRDC_KSPConverged_Private()
39 …if (*reason) PetscCall(PetscInfo(snes, "Default or user provided convergence test KSP iterations=%… in SNESTRDC_KSPConverged_Private()
44 …PetscCall(PetscInfo(snes, "Ending linear iteration early, delta=%g, length=%g\n", (double)neP->del… in SNESTRDC_KSPConverged_Private()
63 static PetscErrorCode SNESTRDC_Converged_Private(SNES snes, PetscInt it, PetscReal xnorm, PetscReal… in SNESTRDC_Converged_Private() argument
65 SNES_NEWTONTRDC *neP = (SNES_NEWTONTRDC *)snes->data; in SNESTRDC_Converged_Private()
70 …PetscCall(PetscInfo(snes, "Diverged due to too small a trust region %g<%g*%g\n", (double)neP->delt… in SNESTRDC_Converged_Private()
[all …]
|
| /petsc/src/snes/impls/multiblock/ |
| H A D | multiblock.c | 11 SNES snes; /* Solver for this block */ member
26 static PetscErrorCode SNESReset_Multiblock(SNES snes) in SNESReset_Multiblock() argument
28 SNES_Multiblock *mb = (SNES_Multiblock *)snes->data; in SNESReset_Multiblock()
33 PetscCall(SNESReset(blocks->snes)); in SNESReset_Multiblock()
53 static PetscErrorCode SNESDestroy_Multiblock(SNES snes) in SNESDestroy_Multiblock() argument
55 SNES_Multiblock *mb = (SNES_Multiblock *)snes->data; in SNESDestroy_Multiblock()
59 PetscCall(SNESReset_Multiblock(snes)); in SNESDestroy_Multiblock()
62 PetscCall(SNESDestroy(&blocks->snes)); in SNESDestroy_Multiblock()
68 PetscCall(PetscFree(snes->data)); in SNESDestroy_Multiblock()
73 static PetscErrorCode SNESMultiblockSetFieldsRuntime_Private(SNES snes) in SNESMultiblockSetFieldsRuntime_Private() argument
[all …]
|
| /petsc/include/petsc/private/ |
| H A D | snesimpl.h | 270 static inline PetscErrorCode SNESLogConvergenceHistory(SNES snes, PetscReal res, PetscInt its) in SNESLogConvergenceHistory() argument
273 PetscCall(PetscObjectSAWsTakeAccess((PetscObject)snes)); in SNESLogConvergenceHistory()
274 if (snes->conv_hist && snes->conv_hist_max > snes->conv_hist_len) { in SNESLogConvergenceHistory()
275 if (snes->conv_hist) snes->conv_hist[snes->conv_hist_len] = res; in SNESLogConvergenceHistory()
276 if (snes->conv_hist_its) snes->conv_hist_its[snes->conv_hist_len] = its; in SNESLogConvergenceHistory()
277 snes->conv_hist_len++; in SNESLogConvergenceHistory()
279 PetscCall(PetscObjectSAWsGrantAccess((PetscObject)snes)); in SNESLogConvergenceHistory()
351 #define SNESCheckFunctionDomainError(snes, fnorm) \ argument
354 …PetscCheck(!snes->errorifnotconverged, PetscObjectComm((PetscObject)snes), PETSC_ERR_NOT_CONVERGED…
357 …tscCallMPI(MPIU_Allreduce(&snes->functiondomainerror, &domainerror, 1, MPI_C_BOOL, MPI_LOR, PetscO…
[all …]
|
| /petsc/src/snes/interface/ftn-custom/ |
| H A D | zsnesf.c | 86 static PetscErrorCode ourtrprecheckfunction(SNES snes, Vec x, Vec y, PetscBool *changed_y, PetscCtx… in ourtrprecheckfunction() argument
90 …PetscCall(PetscObjectGetFortranCallback((PetscObject)snes, PETSC_FORTRAN_CALLBACK_CLASS, _cb.trpre… in ourtrprecheckfunction()
92 …seFortranCallback(snes, _cb.trprecheck, (SNES *, Vec *, Vec *, PetscBool *, void *, PetscErrorCode… in ourtrprecheckfunction()
95 PETSC_EXTERN void snesnewtontrsetprecheck_(SNES *snes, void (*func)(SNES, Vec, Vec, PetscBool *, vo… in snesnewtontrsetprecheck_() argument
97 …*ierr = PetscObjectSetFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, &_cb.trpre… in snesnewtontrsetprecheck_()
100 …*ierr = PetscObjectSetFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, &_cb.trpre… in snesnewtontrsetprecheck_()
103 *ierr = SNESNewtonTRSetPreCheck(*snes, ourtrprecheckfunction, NULL); in snesnewtontrsetprecheck_()
106 PETSC_EXTERN void snesnewtontrdcsetprecheck_(SNES *snes, void (*func)(SNES, Vec, Vec, PetscBool *, … in snesnewtontrdcsetprecheck_() argument
108 …*ierr = PetscObjectSetFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, &_cb.trpre… in snesnewtontrdcsetprecheck_()
111 …*ierr = PetscObjectSetFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, &_cb.trpre… in snesnewtontrdcsetprecheck_()
[all …]
|
| /petsc/src/snes/impls/nasm/ |
| H A D | nasm.c | 33 static PetscErrorCode SNESReset_NASM(SNES snes) in SNESReset_NASM() argument
35 SNES_NASM *nasm = (SNES_NASM *)snes->data; in SNESReset_NASM()
72 static PetscErrorCode SNESDestroy_NASM(SNES snes) in SNESDestroy_NASM() argument
75 PetscCall(SNESReset_NASM(snes)); in SNESDestroy_NASM()
76 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNASMSetType_C", NULL)); in SNESDestroy_NASM()
77 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNASMGetType_C", NULL)); in SNESDestroy_NASM()
78 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNASMSetSubdomains_C", NULL)); in SNESDestroy_NASM()
79 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNASMGetSubdomains_C", NULL)); in SNESDestroy_NASM()
80 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNASMSetDamping_C", NULL)); in SNESDestroy_NASM()
81 PetscCall(PetscObjectComposeFunction((PetscObject)snes, "SNESNASMGetDamping_C", NULL)); in SNESDestroy_NASM()
[all …]
|
| /petsc/src/snes/impls/vi/rs/ |
| H A D | virs.c | 19 PetscErrorCode SNESVIGetInactiveSet(SNES snes, IS *inact) in SNESVIGetInactiveSet() argument
21 SNES_VINEWTONRSLS *vi = (SNES_VINEWTONRSLS *)snes->data; in SNESVIGetInactiveSet()
249 static PetscErrorCode SNESCreateSubVectors_VINEWTONRSLS(SNES snes, PetscInt n, Vec *newv) in SNESCreateSubVectors_VINEWTONRSLS() argument
254 PetscCall(VecCreate(PetscObjectComm((PetscObject)snes), &v)); in SNESCreateSubVectors_VINEWTONRSLS()
262 static PetscErrorCode SNESVIResetPCandKSP(SNES snes, Mat Amat, Mat Pmat) in SNESVIResetPCandKSP() argument
267 PetscCall(SNESGetKSP(snes, &snesksp)); in SNESVIResetPCandKSP()
296 static PetscErrorCode SNESSolve_VINEWTONRSLS(SNES snes) in SNESSolve_VINEWTONRSLS() argument
298 SNES_VINEWTONRSLS *vi = (SNES_VINEWTONRSLS *)snes->data; in SNESSolve_VINEWTONRSLS()
309 PetscCall(SNESGetKSP(snes, &ksp)); in SNESSolve_VINEWTONRSLS()
313 snes->numFailures = 0; in SNESSolve_VINEWTONRSLS()
[all …]
|