Home
last modified time | relevance | path

Searched refs:adapt (Results 1 – 25 of 58) sorted by relevance

123

/petsc/src/ts/adapt/interface/
H A Dtsadapt.c129 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()
[all …]
/petsc/src/ts/adapt/impls/history/
H A Dadapthist.c8 static PetscErrorCode TSAdaptChoose_History(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, P… in TSAdaptChoose_History() argument
11 TSAdapt_History *thadapt = (TSAdapt_History *)adapt->data; in TSAdaptChoose_History()
14 …PetscCheck(thadapt->hist, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ORDER, "Need to call TSAd… in TSAdaptChoose_History()
25 static PetscErrorCode TSAdaptReset_History(TSAdapt adapt) in TSAdaptReset_History() argument
27 TSAdapt_History *thadapt = (TSAdapt_History *)adapt->data; in TSAdaptReset_History()
34 static PetscErrorCode TSAdaptDestroy_History(TSAdapt adapt) in TSAdaptDestroy_History() argument
37 PetscCall(TSAdaptReset_History(adapt)); in TSAdaptDestroy_History()
38 PetscCall(PetscFree(adapt->data)); in TSAdaptDestroy_History()
43 PetscErrorCode TSAdaptHistorySetTSHistory(TSAdapt adapt, TSHistory hist, PetscBool backward) in TSAdaptHistorySetTSHistory() argument
50 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptHistorySetTSHistory()
[all …]
/petsc/src/ts/adapt/impls/basic/
H A Dadaptbasic.c4 static PetscErrorCode TSAdaptChoose_Basic(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, Pet… in TSAdaptChoose_Basic() argument
11 PetscReal safety = adapt->safety; in TSAdaptChoose_Basic()
20 PetscCall(TSEvaluateWLTE(ts, adapt->wnormtype, &order, &enorm)); in TSAdaptChoose_Basic()
21 …PetscCheck(enorm < 0 || order >= 1, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE,… in TSAdaptChoose_Basic()
23 …PetscCheck(adapt->candidates.n >= 1, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_WRONGSTATE… in TSAdaptChoose_Basic()
24 …adapt->candidates.inuse_set, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_WRONGSTATE, "The c… in TSAdaptChoose_Basic()
25 order = adapt->candidates.order[0]; in TSAdaptChoose_Basic()
29 PetscCall(TSErrorWeightedNorm(ts, ts->vec_sol, Y, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_Basic()
42 …if (!*accept) safety *= adapt->reject_safety; /* The last attempt also failed, shorten more aggres… in TSAdaptChoose_Basic()
43 if (h < (1 + PETSC_SQRT_MACHINE_EPSILON) * adapt->dt_min) { in TSAdaptChoose_Basic()
[all …]
/petsc/src/ts/adapt/impls/dsp/
H A Dadaptdsp.c43 static PetscErrorCode TSAdaptRestart_DSP(TSAdapt adapt) in TSAdaptRestart_DSP() argument
45 TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data; in TSAdaptRestart_DSP()
54 static PetscErrorCode TSAdaptRollBack_DSP(TSAdapt adapt) in TSAdaptRollBack_DSP() argument
56 TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data; in TSAdaptRollBack_DSP()
68 static PetscErrorCode TSAdaptChoose_DSP(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, Petsc… in TSAdaptChoose_DSP() argument
70 TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data; in TSAdaptChoose_DSP()
74 PetscReal safety = adapt->safety * (PetscReal)0.9; in TSAdaptChoose_DSP()
76 PetscReal hmin = adapt->dt_min * (1 + PETSC_SQRT_MACHINE_EPSILON); in TSAdaptChoose_DSP()
84 PetscCall(TSEvaluateWLTE(ts, adapt->wnormtype, &order, &enorm)); in TSAdaptChoose_DSP()
85 …PetscCheck(enorm < 0 || order >= 1, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE,… in TSAdaptChoose_DSP()
[all …]
/petsc/src/ts/adapt/impls/glee/
H A Dadaptglee.c8 static PetscErrorCode TSAdaptChoose_GLEE(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, Pets… in TSAdaptChoose_GLEE() argument
10 TSAdapt_GLEE *glee = (TSAdapt_GLEE *)adapt->data; in TSAdaptChoose_GLEE()
18 safety = adapt->safety; in TSAdaptChoose_GLEE()
20 order = adapt->candidates.order[0]; in TSAdaptChoose_GLEE()
26 if (!glee->Y && adapt->glee_use_local) { in TSAdaptChoose_GLEE()
34 …if (adapt->glee_use_local) PetscCall(VecAXPY(E, -1.0, glee->Y)); /* local error = current error - … in TSAdaptChoose_GLEE()
37 PetscCall(TSErrorWeightedENorm(ts, E, X, X, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_GLEE()
45 PetscCall(TSErrorWeightedNorm(ts, X, Y, adapt->wnormtype, &enorm, &enorma, &enormr)); in TSAdaptChoose_GLEE()
58 …if (!*accept) safety *= adapt->reject_safety; /* The last attempt also failed, shorten more aggres… in TSAdaptChoose_GLEE()
59 if (h < (1 + PETSC_SQRT_MACHINE_EPSILON) * adapt->dt_min) { in TSAdaptChoose_GLEE()
[all …]
/petsc/src/ts/impls/implicit/glle/
H A Dglleadapt.c118 PetscErrorCode TSGLLEAdaptSetType(TSGLLEAdapt adapt, TSGLLEAdaptType type) in TSGLLEAdaptSetType() argument
125 if (((PetscObject)adapt)->type_name) PetscUseTypeMethod(adapt, destroy); in TSGLLEAdaptSetType()
126 PetscCall((*r)(adapt)); in TSGLLEAdaptSetType()
127 PetscCall(PetscObjectChangeTypeName((PetscObject)adapt, type)); in TSGLLEAdaptSetType()
131 PetscErrorCode TSGLLEAdaptSetOptionsPrefix(TSGLLEAdapt adapt, const char prefix[]) in TSGLLEAdaptSetOptionsPrefix() argument
134 PetscCall(PetscObjectSetOptionsPrefix((PetscObject)adapt, prefix)); in TSGLLEAdaptSetOptionsPrefix()
138 PetscErrorCode TSGLLEAdaptView(TSGLLEAdapt adapt, PetscViewer viewer) in TSGLLEAdaptView() argument
145 PetscCall(PetscObjectPrintClassNamePrefixType((PetscObject)adapt, viewer)); in TSGLLEAdaptView()
146 if (adapt->ops->view) { in TSGLLEAdaptView()
148 PetscUseTypeMethod(adapt, view, viewer); in TSGLLEAdaptView()
[all …]
H A Dglle.c692 PetscErrorCode TSGLLEGetAdapt(TS ts, TSGLLEAdapt *adapt) in TSGLLEGetAdapt() argument
696 PetscAssertPointer(adapt, 2); in TSGLLEGetAdapt()
697 PetscUseMethod(ts, "TSGLLEGetAdapt_C", (TS, TSGLLEAdapt *), (ts, adapt)); in TSGLLEGetAdapt()
777 static PetscErrorCode TSGLLEGetAdapt_GLLE(TS ts, TSGLLEAdapt *adapt) in TSGLLEGetAdapt_GLLE() argument
782 if (!gl->adapt) { in TSGLLEGetAdapt_GLLE()
783 PetscCall(TSGLLEAdaptCreate(PetscObjectComm((PetscObject)ts), &gl->adapt)); in TSGLLEGetAdapt_GLLE()
784 PetscCall(PetscObjectIncrementTabLevel((PetscObject)gl->adapt, (PetscObject)ts, 1)); in TSGLLEGetAdapt_GLLE()
786 *adapt = gl->adapt; in TSGLLEGetAdapt_GLLE()
814 …PetscCall(TSGLLEAdaptChoose(gl->adapt, n, orders, errors, costs, cur, h, tleft, &next_sc, next_h, … in TSGLLEChooseNextScheme()
1036 if (gl->adapt) PetscCall(TSGLLEAdaptDestroy(&gl->adapt)); in TSDestroy_GLLE()
[all …]
/petsc/src/ts/adapt/impls/cfl/
H A Dadaptcfl.c3 static PetscErrorCode TSAdaptChoose_CFL(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, Petsc… in TSAdaptChoose_CFL() argument
11 PetscCall(TSAdaptCandidatesGet(adapt, &ncandidates, NULL, NULL, &ccflarray, NULL)); in TSAdaptChoose_CFL()
14 …PetscCheck(adapt->always_accept, PetscObjectComm((PetscObject)adapt), PETSC_ERR_SUP, "Step rejecti… in TSAdaptChoose_CFL()
19 if (adapt->always_accept) { in TSAdaptChoose_CFL()
20 …PetscCall(PetscInfo(adapt, "Step length %g with scheme of CFL coefficient %g did not satisfy user-… in TSAdaptChoose_CFL()
22 …PetscCall(PetscInfo(adapt, "Step length %g with scheme of CFL coefficient %g did not satisfy user-… in TSAdaptChoose_CFL()
28 hcfl = adapt->safety * cfltimestep * ccfl; in TSAdaptChoose_CFL()
29 if (hcfl < adapt->dt_min) { in TSAdaptChoose_CFL()
30 …adapt, "Cannot satisfy CFL constraint %g (with %g safety) at minimum time step %g with method coef… in TSAdaptChoose_CFL()
34 *next_h = PetscClipInterval(hcfl, adapt->dt_min, adapt->dt_max); in TSAdaptChoose_CFL()
[all …]
/petsc/src/ts/tests/
H A Dex14.c24 PetscErrorCode TestCheckStage(TSAdapt adapt, TS ts, PetscReal t, Vec X, PetscBool *accept) in TestCheckStage() argument
46 TSAdapt adapt; in TestExplicitTS() local
86 PetscCall(TSGetAdapt(ts, &adapt)); in TestExplicitTS()
87 PetscCall(TSAdaptSetCheckStage(adapt, TestCheckStage)); in TestExplicitTS()
93 PetscCall(TSAdaptSetCheckStage(adapt, NULL)); in TestExplicitTS()
105 TSAdapt adapt; in TestTSRK() local
115 PetscCall(TSGetAdapt(ts, &adapt)); in TestTSRK()
116 PetscCall(TSAdaptGetType(adapt, &adapttype)); in TestTSRK()
121 if (rk1 || rk3 || rk4) PetscCall(TSAdaptSetType(adapt, TSADAPTNONE)); in TestTSRK()
125 PetscCall(TSGetAdapt(ts, &adapt)); in TestTSRK()
[all …]
/petsc/src/dm/impls/plex/tutorials/
H A Dex10.c6 PetscBool adapt; /* Flag for adaptation of the surface mesh */ member
12 options->adapt = PETSC_FALSE; in ProcessOptions()
15 …apt", "Flag for adaptation of the surface mesh", "ex10.c", options->adapt, &options->adapt, NULL)); in ProcessOptions()
60 PetscBool adapt = PETSC_TRUE; in AdaptMesh() local
63 if (!ctx->adapt) PetscFunctionReturn(PETSC_SUCCESS); in AdaptMesh()
88 while (adapt) { in AdaptMesh()
93 adapt = PETSC_FALSE; in AdaptMesh()
137 adapt = PETSC_TRUE; in AdaptMesh()
/petsc/src/ts/adapt/impls/none/
H A Dadaptnone.c3 static PetscErrorCode TSAdaptChoose_None(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, Pets… in TSAdaptChoose_None() argument
22 PETSC_EXTERN PetscErrorCode TSAdaptCreate_None(TSAdapt adapt) in TSAdaptCreate_None() argument
25 adapt->ops->choose = TSAdaptChoose_None; in TSAdaptCreate_None()
/petsc/src/ts/tutorials/
H A Dex41.c168 TSAdapt adapt; in main() local
223 PetscCall(TSGetAdapt(ts, &adapt)); in main()
224 PetscCall(TSAdaptSetType(adapt, TSADAPTBASIC)); in main()
225 PetscCall(TSAdaptSetStepLimits(adapt, 0.0, 0.5)); in main()
248 PetscCall(TSGetAdapt(ts, &adapt)); in main()
249 PetscCall(TSAdaptSetType(adapt, TSADAPTHISTORY)); in main()
251 PetscCall(TSAdaptHistorySetTrajectory(adapt, tj, PETSC_FALSE)); in main()
252 PetscCall(TSAdaptHistoryGetStep(adapt, 0, &t0, &dt)); in main()
255 PetscCall(TSAdaptSetType(adapt, TSADAPTBASIC)); in main()
256 PetscCall(TSAdaptSetStepLimits(adapt, 0.0, 0.5)); in main()
H A Dex40.c168 TSAdapt adapt; in main() local
237 PetscCall(TSGetAdapt(ts, &adapt)); in main()
238 PetscCall(TSAdaptSetStepLimits(adapt, 0.0, 0.5)); in main()
263 PetscCall(TSGetAdapt(ts, &adapt)); in main()
264 PetscCall(TSAdaptSetType(adapt, TSADAPTHISTORY)); in main()
266 PetscCall(TSAdaptHistorySetTrajectory(adapt, tj, PETSC_FALSE)); in main()
267 PetscCall(TSAdaptHistoryGetStep(adapt, 0, &t0, &dt)); in main()
277 PetscCall(TSAdaptSetType(adapt, TSADAPTBASIC)); in main()
278 PetscCall(TSAdaptSetStepLimits(adapt, 0.0, 0.5)); in main()
H A Dex51.c62 TSAdapt adapt; in main() local
104 PetscCall(TSGetAdapt(ts, &adapt)); in main()
105 PetscCall(TSAdaptSetType(adapt, TSADAPTNONE)); in main()
H A Dex44.c127 TSAdapt adapt; in main() local
151 PetscCall(TSGetAdapt(ts, &adapt)); in main()
152 PetscCall(TSAdaptSetStepLimits(adapt, 0.0, 0.5)); in main()
/petsc/src/dm/impls/plex/tests/
H A Dex46.c127 DMLabel adapt; in main() local
131 PetscCall(DMLabelCreate(PETSC_COMM_SELF, "adapt", &adapt)); in main()
132 PetscCall(DMLabelSetDefaultValue(adapt, DM_ADAPT_COARSEN)); in main()
135 if (c % 2) PetscCall(DMLabelSetValue(adapt, c, DM_ADAPT_REFINE)); in main()
137 PetscCall(DMAdaptLabel(dm, adapt, &adm)); in main()
138 PetscCall(DMLabelDestroy(&adapt)); in main()
/petsc/src/dm/impls/forest/
H A Dforest.c211 PetscCall(DMDestroy(&forest->adapt)); in DMDestroy_Forest()
394 PetscErrorCode DMForestSetAdaptivityForest(DM dm, DM adapt) in DMForestSetAdaptivityForest() argument
402 if (adapt) PetscValidHeaderSpecific(adapt, DM_CLASSID, 2); in DMForestSetAdaptivityForest()
405 …PetscCheck(adapt == NULL || !dm->setupcalled, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_WRON… in DMForestSetAdaptivityForest()
408 adaptForest = (DM_Forest *)(adapt ? adapt->data : NULL); in DMForestSetAdaptivityForest()
417 PetscCall(PetscObjectReference((PetscObject)adapt)); in DMForestSetAdaptivityForest()
418 PetscCall(DMDestroy(&forest->adapt)); in DMForestSetAdaptivityForest()
419 forest->adapt = adapt; in DMForestSetAdaptivityForest()
422 PetscCall(DMSetCoarseDM(dm, adapt)); in DMForestSetAdaptivityForest()
426 PetscCall(DMSetFineDM(dm, adapt)); in DMForestSetAdaptivityForest()
[all …]
/petsc/src/ts/impls/explicit/euler/
H A Deuler.c22 PetscCall(TSAdaptCheckStage(ts->adapt, ts, ts->ptime, solution, &stageok)); in TSStep_Euler()
33 PetscCall(TSAdaptChoose(ts->adapt, ts, ts->time_step, NULL, &next_time_step, &accept)); in TSStep_Euler()
52 PetscCall(TSGetAdapt(ts, &ts->adapt)); in TSSetUp_Euler()
53 PetscCall(TSAdaptCandidatesClear(ts->adapt)); in TSSetUp_Euler()
/petsc/src/ts/impls/pseudo/
H A Dposindep.c88 TSAdapt adapt; in TSStep_Pseudo() local
103 PetscCall(TSGetAdapt(ts, &adapt)); in TSStep_Pseudo()
104 PetscCall(TSAdaptCheckStage(adapt, ts, ts->ptime + ts->time_step, ts->vec_sol, &accept)); in TSStep_Pseudo()
108 PetscCall(TSAdaptChoose(adapt, ts, ts->time_step, NULL, &next_time_step, &accept)); in TSStep_Pseudo()
354 static PetscErrorCode TSAdaptChoose_TSPseudo(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, … in TSAdaptChoose_TSPseudo() argument
370 static PetscErrorCode TSAdaptCheckStage_TSPseudo(TSAdapt adapt, TS ts, PetscReal t, Vec Y, PetscBoo… in TSAdaptCheckStage_TSPseudo() argument
394 static PetscErrorCode TSAdaptCreate_TSPseudo(TSAdapt adapt) in TSAdaptCreate_TSPseudo() argument
397 adapt->ops->choose = TSAdaptChoose_TSPseudo; in TSAdaptCreate_TSPseudo()
398 adapt->checkstage = TSAdaptCheckStage_TSPseudo; in TSAdaptCreate_TSPseudo()
/petsc/src/dm/impls/plex/
H A Dplexadapt.c168 PetscErrorCode (*adapt)(DM, Vec, DMLabel, DMLabel, DM *); in DMPlexRefine_Internal() local
196 adapt = fl->adapt; in DMPlexRefine_Internal()
206 adapt = fl->adapt; in DMPlexRefine_Internal()
219 if (adapt) { in DMPlexRefine_Internal()
220 PetscCall((*adapt)(dm, NULL, adaptLabel, NULL, dmRefined)); in DMPlexRefine_Internal()
/petsc/src/dm/impls/forest/tests/
H A Dex3.c14 PetscBool adapt = PETSC_FALSE, userSection = PETSC_FALSE; in main() local
19 PetscCall(PetscOptionsGetBool(NULL, NULL, "-adapt", &adapt, NULL)); in main()
36 if (adapt) { in main()
/petsc/src/ts/impls/bdf/
H A Dbdf.c235 PetscCall(TSAdaptCheckStage(ts->adapt, ts, bdf->time[0], bdf->work[0], accept)); in TSBDF_Restart()
276 PetscCall(TSAdaptCheckStage(ts->adapt, ts, bdf->time[0], bdf->work[0], &stageok)); in TSStep_BDF()
280 PetscCall(TSAdaptCandidatesClear(ts->adapt)); in TSStep_BDF()
281 …PetscCall(TSAdaptCandidateAdd(ts->adapt, BDF_SchemeName[bdf->k], bdf->k, 1, 1.0, 1.0, PETSC_TRUE)); in TSStep_BDF()
282 PetscCall(TSAdaptChoose(ts->adapt, ts, ts->time_step, NULL, &next_time_step, &accept)); in TSStep_BDF()
454 PetscCall(TSGetAdapt(ts, &ts->adapt)); in TSSetUp_BDF()
455 PetscCall(TSAdaptCandidatesClear(ts->adapt)); in TSSetUp_BDF()
456 PetscCall(TSAdaptGetClip(ts->adapt, &low, &high)); in TSSetUp_BDF()
457 PetscCall(TSAdaptSetClip(ts->adapt, low, PetscMin(high, two))); in TSSetUp_BDF()
/petsc/src/ts/impls/implicit/alpha/
H A Dalpha1.c140 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, X1, &stageok)); in TSAlpha_Restart()
150 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, X2, &stageok)); in TSAlpha_Restart()
207 PetscCall(TSAdaptCheckStage(ts->adapt, ts, th->stage_time, th->Xa, &stageok)); in TSStep_Alpha()
212 PetscCall(TSAdaptChoose(ts->adapt, ts, ts->time_step, NULL, &next_time_step, &accept)); in TSStep_Alpha()
359 PetscCall(TSGetAdapt(ts, &ts->adapt)); in TSSetUp_Alpha()
360 PetscCall(TSAdaptCandidatesClear(ts->adapt)); in TSSetUp_Alpha()
361 PetscCall(PetscObjectTypeCompare((PetscObject)ts->adapt, TSADAPTNONE, &match)); in TSSetUp_Alpha()
/petsc/src/dm/interface/
H A Ddmgenerate.c111 entry->adapt = alfnc; in DMGenerateRegister()
181 PetscCall((*fl->adapt)(dm, NULL, label, NULL, dmAdapt)); in DMAdaptLabel()
255 PetscCall((*fl->adapt)(dm, metric, bdLabel, rgLabel, dmAdapt)); in DMAdaptMetric()
/petsc/src/ts/interface/
H A Dts.c16 static PetscErrorCode TSAdaptSetDefaultType(TSAdapt adapt, TSAdaptType default_type) in TSAdaptSetDefaultType() argument
19 PetscValidHeaderSpecific(adapt, TSADAPT_CLASSID, 1); in TSAdaptSetDefaultType()
21 if (!((PetscObject)adapt)->type_name) PetscCall(TSAdaptSetType(adapt, default_type)); in TSAdaptSetDefaultType()
418 PetscCall(TSGetAdapt(ts, &ts->adapt)); in TSSetFromOptions()
419 PetscCall(TSAdaptSetDefaultType(ts->adapt, ts->default_adapt_type)); in TSSetFromOptions()
420 PetscCall(TSAdaptSetFromOptions(ts->adapt, PetscOptionsObject)); in TSSetFromOptions()
1948 PetscCall(TSAdaptView(ts->adapt, viewer)); in TSView()
1968 if (ts->adapt) PetscCall(TSAdaptView(ts->adapt, viewer)); in TSView()
1986 if (ts->adapt) PetscCall(TSAdaptView(ts->adapt, viewer)); in TSView()
2456 PetscCall(TSGetAdapt(ts, &ts->adapt)); in TSSetExactFinalTimeDefault()
[all …]

123