Lines Matching refs:adapt
129 PetscErrorCode TSAdaptSetType(TSAdapt adapt, TSAdaptType type) in TSAdaptSetType() argument
135 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetType()
137 PetscCall(PetscObjectTypeCompare((PetscObject)adapt, type, &match)); in TSAdaptSetType()
140 …PetscCheck(r, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_UNKNOWN_TYPE, "Unknown TSAdapt ty… in TSAdaptSetType()
141 PetscTryTypeMethod(adapt, destroy); in TSAdaptSetType()
142 PetscCall(PetscMemzero(adapt->ops, sizeof(struct _TSAdaptOps))); in TSAdaptSetType()
143 PetscCall(PetscObjectChangeTypeName((PetscObject)adapt, type)); in TSAdaptSetType()
144 PetscCall((*r)(adapt)); in TSAdaptSetType()
163 PetscErrorCode TSAdaptGetType(TSAdapt adapt, TSAdaptType *type) in TSAdaptGetType() argument
166 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptGetType()
168 *type = ((PetscObject)adapt)->type_name; in TSAdaptGetType()
172 PetscErrorCode TSAdaptSetOptionsPrefix(TSAdapt adapt, const char prefix[]) in TSAdaptSetOptionsPrefix() argument
175 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetOptionsPrefix()
176 PetscCall(PetscObjectSetOptionsPrefix((PetscObject)adapt, prefix)); in TSAdaptSetOptionsPrefix()
198 PetscErrorCode TSAdaptLoad(TSAdapt adapt, PetscViewer viewer) in TSAdaptLoad() argument
204 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptLoad()
210 PetscCall(TSAdaptSetType(adapt, type)); in TSAdaptLoad()
211 PetscTryTypeMethod(adapt, load, viewer); in TSAdaptLoad()
246 PetscErrorCode TSAdaptView(TSAdapt adapt, PetscViewer viewer) in TSAdaptView() argument
251 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptView()
252 if (!viewer) PetscCall(PetscViewerASCIIGetStdout(PetscObjectComm((PetscObject)adapt), &viewer)); in TSAdaptView()
254 PetscCheckSameComm(adapt, 1, viewer, 2); in TSAdaptView()
258 PetscCall(PetscObjectPrintClassNamePrefixType((PetscObject)adapt, viewer)); in TSAdaptView()
259 PetscCall(PetscObjectTypeCompare((PetscObject)adapt, TSADAPTNONE, &isnone)); in TSAdaptView()
260 PetscCall(PetscObjectTypeCompare((PetscObject)adapt, TSADAPTGLEE, &isglee)); in TSAdaptView()
262 … if (adapt->always_accept) PetscCall(PetscViewerASCIIPrintf(viewer, " always accepting steps\n")); in TSAdaptView()
263 PetscCall(PetscViewerASCIIPrintf(viewer, " safety factor %g\n", (double)adapt->safety)); in TSAdaptView()
264 …IIPrintf(viewer, " extra safety factor after step rejection %g\n", (double)adapt->reject_safety)); in TSAdaptView()
265 … PetscCall(PetscViewerASCIIPrintf(viewer, " clip fastest increase %g\n", (double)adapt->clip[1])); in TSAdaptView()
266 … PetscCall(PetscViewerASCIIPrintf(viewer, " clip fastest decrease %g\n", (double)adapt->clip[0])); in TSAdaptView()
267 …PetscCall(PetscViewerASCIIPrintf(viewer, " maximum allowed timestep %g\n", (double)adapt->dt_max)… in TSAdaptView()
268 …PetscCall(PetscViewerASCIIPrintf(viewer, " minimum allowed timestep %g\n", (double)adapt->dt_min)… in TSAdaptView()
269 …Printf(viewer, " maximum solution absolute value to be ignored %g\n", (double)adapt->ignore_max)); in TSAdaptView()
272 if (adapt->glee_use_local) { in TSAdaptView()
279 PetscTryTypeMethod(adapt, view, viewer); in TSAdaptView()
285 PetscCall(PetscStrncpy(type, ((PetscObject)adapt)->type_name, 256)); in TSAdaptView()
287 } else PetscTryTypeMethod(adapt, view, viewer); in TSAdaptView()
303 PetscErrorCode TSAdaptReset(TSAdapt adapt) in TSAdaptReset() argument
306 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptReset()
307 PetscTryTypeMethod(adapt, reset); in TSAdaptReset()
311 PetscErrorCode TSAdaptDestroy(TSAdapt *adapt) in TSAdaptDestroy() argument
314 if (!*adapt) PetscFunctionReturn(PETSC_SUCCESS); in TSAdaptDestroy()
315 PetscValidHeaderSpecific(*adapt, TSADAPT_CLASSID, 1); in TSAdaptDestroy()
316 if (--((PetscObject)*adapt)->refct > 0) { in TSAdaptDestroy()
317 *adapt = NULL; in TSAdaptDestroy()
321 PetscCall(TSAdaptReset(*adapt)); in TSAdaptDestroy()
323 PetscTryTypeMethod(*adapt, destroy); in TSAdaptDestroy()
324 PetscCall(PetscViewerDestroy(&(*adapt)->monitor)); in TSAdaptDestroy()
325 PetscCall(PetscHeaderDestroy(adapt)); in TSAdaptDestroy()
345 PetscErrorCode TSAdaptSetMonitor(TSAdapt adapt, PetscBool flg) in TSAdaptSetMonitor() argument
348 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetMonitor()
349 PetscValidLogicalCollectiveBool(adapt, flg, 2); in TSAdaptSetMonitor()
351 …if (!adapt->monitor) PetscCall(PetscViewerASCIIOpen(PetscObjectComm((PetscObject)adapt), "stdout",… in TSAdaptSetMonitor()
353 PetscCall(PetscViewerDestroy(&adapt->monitor)); in TSAdaptSetMonitor()
378 PetscErrorCode TSAdaptSetCheckStage(TSAdapt adapt, PetscErrorCode (*func)(TSAdapt adapt, TS ts, Pet… in TSAdaptSetCheckStage() argument
381 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetCheckStage()
382 adapt->checkstage = func; in TSAdaptSetCheckStage()
403 PetscErrorCode TSAdaptSetAlwaysAccept(TSAdapt adapt, PetscBool flag) in TSAdaptSetAlwaysAccept() argument
406 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetAlwaysAccept()
407 PetscValidLogicalCollectiveBool(adapt, flag, 2); in TSAdaptSetAlwaysAccept()
408 adapt->always_accept = flag; in TSAdaptSetAlwaysAccept()
436 PetscErrorCode TSAdaptSetSafety(TSAdapt adapt, PetscReal safety, PetscReal reject_safety) in TSAdaptSetSafety() argument
439 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetSafety()
440 PetscValidLogicalCollectiveReal(adapt, safety, 2); in TSAdaptSetSafety()
441 PetscValidLogicalCollectiveReal(adapt, reject_safety, 3); in TSAdaptSetSafety()
446 if (safety != (PetscReal)PETSC_CURRENT) adapt->safety = safety; in TSAdaptSetSafety()
447 if (reject_safety != (PetscReal)PETSC_CURRENT) adapt->reject_safety = reject_safety; in TSAdaptSetSafety()
467 PetscErrorCode TSAdaptGetSafety(TSAdapt adapt, PetscReal *safety, PetscReal *reject_safety) in TSAdaptGetSafety() argument
470 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptGetSafety()
473 if (safety) *safety = adapt->safety; in TSAdaptGetSafety()
474 if (reject_safety) *reject_safety = adapt->reject_safety; in TSAdaptGetSafety()
495 PetscErrorCode TSAdaptSetMaxIgnore(TSAdapt adapt, PetscReal max_ignore) in TSAdaptSetMaxIgnore() argument
498 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetMaxIgnore()
499 PetscValidLogicalCollectiveReal(adapt, max_ignore, 2); in TSAdaptSetMaxIgnore()
500 adapt->ignore_max = max_ignore; in TSAdaptSetMaxIgnore()
520 PetscErrorCode TSAdaptGetMaxIgnore(TSAdapt adapt, PetscReal *max_ignore) in TSAdaptGetMaxIgnore() argument
523 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptGetMaxIgnore()
525 *max_ignore = adapt->ignore_max; in TSAdaptGetMaxIgnore()
552 PetscErrorCode TSAdaptSetClip(TSAdapt adapt, PetscReal low, PetscReal high) in TSAdaptSetClip() argument
555 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetClip()
556 PetscValidLogicalCollectiveReal(adapt, low, 2); in TSAdaptSetClip()
557 PetscValidLogicalCollectiveReal(adapt, high, 3); in TSAdaptSetClip()
561 if (low != (PetscReal)PETSC_CURRENT) adapt->clip[0] = low; in TSAdaptSetClip()
562 if (high != (PetscReal)PETSC_CURRENT) adapt->clip[1] = high; in TSAdaptSetClip()
582 PetscErrorCode TSAdaptGetClip(TSAdapt adapt, PetscReal *low, PetscReal *high) in TSAdaptGetClip() argument
585 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptGetClip()
588 if (low) *low = adapt->clip[0]; in TSAdaptGetClip()
589 if (high) *high = adapt->clip[1]; in TSAdaptGetClip()
609 PetscErrorCode TSAdaptSetScaleSolveFailed(TSAdapt adapt, PetscReal scale) in TSAdaptSetScaleSolveFailed() argument
612 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetScaleSolveFailed()
613 PetscValidLogicalCollectiveReal(adapt, scale, 2); in TSAdaptSetScaleSolveFailed()
616 adapt->scale_solve_failed = scale; in TSAdaptSetScaleSolveFailed()
635 PetscErrorCode TSAdaptGetScaleSolveFailed(TSAdapt adapt, PetscReal *scale) in TSAdaptGetScaleSolveFailed() argument
638 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptGetScaleSolveFailed()
640 if (scale) *scale = adapt->scale_solve_failed; in TSAdaptGetScaleSolveFailed()
668 PetscErrorCode TSAdaptSetStepLimits(TSAdapt adapt, PetscReal hmin, PetscReal hmax) in TSAdaptSetStepLimits() argument
671 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetStepLimits()
672 PetscValidLogicalCollectiveReal(adapt, hmin, 2); in TSAdaptSetStepLimits()
673 PetscValidLogicalCollectiveReal(adapt, hmax, 3); in TSAdaptSetStepLimits()
676 if (hmin != (PetscReal)PETSC_CURRENT) adapt->dt_min = hmin; in TSAdaptSetStepLimits()
677 if (hmax != (PetscReal)PETSC_CURRENT) adapt->dt_max = hmax; in TSAdaptSetStepLimits()
678 hmin = adapt->dt_min; in TSAdaptSetStepLimits()
679 hmax = adapt->dt_max; in TSAdaptSetStepLimits()
700 PetscErrorCode TSAdaptGetStepLimits(TSAdapt adapt, PetscReal *hmin, PetscReal *hmax) in TSAdaptGetStepLimits() argument
703 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptGetStepLimits()
706 if (hmin) *hmin = adapt->dt_min; in TSAdaptGetStepLimits()
707 if (hmax) *hmax = adapt->dt_max; in TSAdaptGetStepLimits()
740 PetscErrorCode TSAdaptSetFromOptions(TSAdapt adapt, PetscOptionItems PetscOptionsObject) in TSAdaptSetFromOptions() argument
748 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetFromOptions()
752 …daptivity", "TSAdaptSetType", TSAdaptList, ((PetscObject)adapt)->type_name ? ((PetscObject)adapt)-… in TSAdaptSetFromOptions()
753 if (flg || !((PetscObject)adapt)->type_name) PetscCall(TSAdaptSetType(adapt, type)); in TSAdaptSetFromOptions()
755 …adapt_always_accept", "Always accept the step", "TSAdaptSetAlwaysAccept", adapt->always_accept, &f… in TSAdaptSetFromOptions()
756 if (set) PetscCall(TSAdaptSetAlwaysAccept(adapt, flg)); in TSAdaptSetFromOptions()
758 safety = adapt->safety; in TSAdaptSetFromOptions()
759 reject_safety = adapt->reject_safety; in TSAdaptSetFromOptions()
762 if (set || flg) PetscCall(TSAdaptSetSafety(adapt, safety, reject_safety)); in TSAdaptSetFromOptions()
765 clip[0] = adapt->clip[0]; in TSAdaptSetFromOptions()
766 clip[1] = adapt->clip[1]; in TSAdaptSetFromOptions()
768 …PetscCheck(!set || (two == 2), PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE, "Mus… in TSAdaptSetFromOptions()
769 if (set) PetscCall(TSAdaptSetClip(adapt, clip[0], clip[1])); in TSAdaptSetFromOptions()
771 hmin = adapt->dt_min; in TSAdaptSetFromOptions()
772 hmax = adapt->dt_max; in TSAdaptSetFromOptions()
775 if (set || flg) PetscCall(TSAdaptSetStepLimits(adapt, hmin, hmax)); in TSAdaptSetFromOptions()
777 …gnores (absolute) solution values smaller than this value", "", adapt->ignore_max, &adapt->ignore_… in TSAdaptSetFromOptions()
778 …EE adaptor uses local error estimation for step control", "", adapt->glee_use_local, &adapt->glee_… in TSAdaptSetFromOptions()
780 …"Scale step by this factor if solve fails", "TSAdaptSetScaleSolveFailed", adapt->scale_solve_faile… in TSAdaptSetFromOptions()
781 if (set) PetscCall(TSAdaptSetScaleSolveFailed(adapt, scale)); in TSAdaptSetFromOptions()
783 …computed for error estimation", "", NormTypes, (PetscEnum)adapt->wnormtype, (PetscEnum *)&adapt->w… in TSAdaptSetFromOptions()
784 …PetscCheck(adapt->wnormtype == NORM_2 || adapt->wnormtype == NORM_INFINITY, PetscObjectComm((Petsc… in TSAdaptSetFromOptions()
786 … to failed solver", "TSAdaptSetTimeStepIncreaseDelay", adapt->timestepjustdecreased_delay, &adapt-… in TSAdaptSetFromOptions()
788 …nitor", "Print choices made by adaptive controller", "TSAdaptSetMonitor", adapt->monitor ? PETSC_T… in TSAdaptSetFromOptions()
789 if (set) PetscCall(TSAdaptSetMonitor(adapt, flg)); in TSAdaptSetFromOptions()
791 PetscTryTypeMethod(adapt, setfromoptions, PetscOptionsObject); in TSAdaptSetFromOptions()
808 PetscErrorCode TSAdaptCandidatesClear(TSAdapt adapt) in TSAdaptCandidatesClear() argument
811 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptCandidatesClear()
812 PetscCall(PetscMemzero(&adapt->candidates, sizeof(adapt->candidates))); in TSAdaptCandidatesClear()
834 PetscErrorCode TSAdaptCandidateAdd(TSAdapt adapt, const char name[], PetscInt order, PetscInt stage… in TSAdaptCandidateAdd() argument
839 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptCandidateAdd()
840 …PetscCheck(order >= 1, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE, "Classical o… in TSAdaptCandidateAdd()
842 …PetscCheck(!adapt->candidates.inuse_set, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_WRONGS… in TSAdaptCandidateAdd()
843 adapt->candidates.inuse_set = PETSC_TRUE; in TSAdaptCandidateAdd()
846 c = inuse ? 0 : !adapt->candidates.inuse_set + adapt->candidates.n; in TSAdaptCandidateAdd()
848 adapt->candidates.name[c] = name; in TSAdaptCandidateAdd()
849 adapt->candidates.order[c] = order; in TSAdaptCandidateAdd()
850 adapt->candidates.stageorder[c] = stageorder; in TSAdaptCandidateAdd()
851 adapt->candidates.ccfl[c] = ccfl; in TSAdaptCandidateAdd()
852 adapt->candidates.cost[c] = cost; in TSAdaptCandidateAdd()
853 adapt->candidates.n++; in TSAdaptCandidateAdd()
879 PetscErrorCode TSAdaptCandidatesGet(TSAdapt adapt, PetscInt *n, const PetscInt **order, const Petsc… in TSAdaptCandidatesGet() argument
882 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptCandidatesGet()
883 if (n) *n = adapt->candidates.n; in TSAdaptCandidatesGet()
884 if (order) *order = adapt->candidates.order; in TSAdaptCandidatesGet()
885 if (stageorder) *stageorder = adapt->candidates.stageorder; in TSAdaptCandidatesGet()
886 if (ccfl) *ccfl = adapt->candidates.ccfl; in TSAdaptCandidatesGet()
887 if (cost) *cost = adapt->candidates.cost; in TSAdaptCandidatesGet()
914 PetscErrorCode TSAdaptChoose(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_… in TSAdaptChoose() argument
916 PetscInt ncandidates = adapt->candidates.n; in TSAdaptChoose()
923 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptChoose()
934 if (adapt->monitor) { in TSAdaptChoose()
935 PetscCall(PetscViewerASCIIAddTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptChoose()
945 …PetscCall(PetscViewerASCIIPrintf(adapt->monitor, " TSAdapt does not interfere, step %3" PetscIn… in TSAdaptChoose()
946 …} else PetscCall(PetscViewerASCIIPrintf(adapt->monitor, " TSAdapt does not interfere, step %3" … in TSAdaptChoose()
948 PetscCall(PetscViewerASCIISubtractTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptChoose()
953 PetscUseTypeMethod(adapt, choose, ts, h, &scheme, next_h, accept, &wlte, &wltea, &wlter); in TSAdaptChoose()
954 … (ncandidates <= 0 || scheme < ncandidates), PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OU… in TSAdaptChoose()
955 …PetscCheck(*next_h >= 0, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE, "Computed … in TSAdaptChoose()
961 PetscReal a = (PetscReal)(1.0 + adapt->matchstepfac[0]); in TSAdaptChoose()
962 PetscReal b = adapt->matchstepfac[1]; in TSAdaptChoose()
979 …if (h <= *next_h && *next_h <= adapt->dt_eval_times_cached) *next_h = adapt->dt_eval_times_cached;… in TSAdaptChoose()
984 …PetscCheck(ts->eval_times->worktol == 0, PetscObjectComm((PetscObject)adapt), PETSC_ERR_PLIB, "Une… in TSAdaptChoose()
997 …if (ts->eval_times && h1 != *next_h && !adapt->dt_eval_times_cached) adapt->dt_eval_times_cached =… in TSAdaptChoose()
998 …if (ts->eval_times && h1 == *next_h && adapt->dt_eval_times_cached) adapt->dt_eval_times_cached = … in TSAdaptChoose()
1000 if (adapt->monitor) { in TSAdaptChoose()
1001 const char *sc_name = (scheme < ncandidates) ? adapt->candidates.name[scheme] : ""; in TSAdaptChoose()
1002 PetscCall(PetscViewerASCIIAddTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptChoose()
1004 …CIIPrintf(adapt->monitor, " TSAdapt %s %s %" PetscInt_FMT ":%s step %3" PetscInt_FMT " %s t=%-1… in TSAdaptChoose()
1007 …adapt->monitor, " TSAdapt %s %s %" PetscInt_FMT ":%s step %3" PetscInt_FMT " %s t=%-11g+%10.3e … in TSAdaptChoose()
1010 PetscCall(PetscViewerASCIISubtractTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptChoose()
1040 PetscErrorCode TSAdaptSetTimeStepIncreaseDelay(TSAdapt adapt, PetscInt cnt) in TSAdaptSetTimeStepIncreaseDelay() argument
1043 adapt->timestepjustdecreased_delay = cnt; in TSAdaptSetTimeStepIncreaseDelay()
1065 PetscErrorCode TSAdaptCheckStage(TSAdapt adapt, TS ts, PetscReal t, Vec Y, PetscBool *accept) in TSAdaptCheckStage() argument
1071 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptCheckStage()
1084 if (adapt->monitor) { in TSAdaptCheckStage()
1085 PetscCall(PetscViewerASCIIAddTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptCheckStage()
1086 …adapt->monitor, " TSAdapt %s step %3" PetscInt_FMT " stage rejected t=%-11g+%10.3e, nonlinear s… in TSAdaptCheckStage()
1088 PetscCall(PetscViewerASCIISubtractTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptCheckStage()
1093 if (*accept && adapt->checkstage) PetscCall((*adapt->checkstage)(adapt, ts, t, Y, accept)); in TSAdaptCheckStage()
1099 if (adapt->monitor) { in TSAdaptCheckStage()
1100 PetscCall(PetscViewerASCIIAddTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptCheckStage()
1101 …ll(PetscViewerASCIIPrintf(adapt->monitor, " TSAdapt %s step %3" PetscInt_FMT " stage rejected b… in TSAdaptCheckStage()
1102 PetscCall(PetscViewerASCIISubtractTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptCheckStage()
1110 new_dt = dt * adapt->scale_solve_failed; in TSAdaptCheckStage()
1112 adapt->timestepjustdecreased += adapt->timestepjustdecreased_delay; in TSAdaptCheckStage()
1113 if (adapt->monitor) { in TSAdaptCheckStage()
1114 PetscCall(PetscViewerASCIIAddTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptCheckStage()
1115 …(adapt->monitor, " TSAdapt %s step %3" PetscInt_FMT " stage rejected (SNES reason %s) t=%-11g+%… in TSAdaptCheckStage()
1117 PetscCall(PetscViewerASCIISubtractTab(adapt->monitor, ((PetscObject)adapt)->tablevel)); in TSAdaptCheckStage()
1143 TSAdapt adapt; in TSAdaptCreate() local
1149 …PetscCall(PetscHeaderCreate(adapt, TSADAPT_CLASSID, "TSAdapt", "Time stepping adaptivity", "TS", c… in TSAdaptCreate()
1150 adapt->always_accept = PETSC_FALSE; in TSAdaptCreate()
1151 adapt->safety = 0.9; in TSAdaptCreate()
1152 adapt->reject_safety = 0.5; in TSAdaptCreate()
1153 adapt->clip[0] = 0.1; in TSAdaptCreate()
1154 adapt->clip[1] = 10.; in TSAdaptCreate()
1155 adapt->dt_min = 1e-20; in TSAdaptCreate()
1156 adapt->dt_max = 1e+20; in TSAdaptCreate()
1157 adapt->ignore_max = -1.0; in TSAdaptCreate()
1158 adapt->glee_use_local = PETSC_TRUE; in TSAdaptCreate()
1159 adapt->scale_solve_failed = 0.25; in TSAdaptCreate()
1162 adapt->matchstepfac[0] = 0.01; /* allow 1% step size increase in the last step */ in TSAdaptCreate()
1163 …adapt->matchstepfac[1] = 2.0; /* halve last step if it is greater than what remains d… in TSAdaptCreate()
1164 adapt->wnormtype = NORM_2; in TSAdaptCreate()
1165 adapt->timestepjustdecreased_delay = 0; in TSAdaptCreate()
1166 *inadapt = adapt; in TSAdaptCreate()