| /petsc/src/ts/impls/implicit/glle/ |
| H A D | glleadapt.c | 185 …etscInt cur, PetscReal h, PetscReal tleft, PetscInt *next_sc, PetscReal *next_h, PetscBool *finish) in TSGLLEAdaptChoose() argument
193 PetscAssertPointer(next_h, 10); in TSGLLEAdaptChoose()
195 …PetscUseTypeMethod(adapt, choose, n, orders, errors, cost, cur, h, tleft, next_sc, next_h, finish); in TSGLLEAdaptChoose()
226 …etscInt cur, PetscReal h, PetscReal tleft, PetscInt *next_sc, PetscReal *next_h, PetscBool *finish) in TSGLLEAdaptChoose_None() argument
230 *next_h = h; in TSGLLEAdaptChoose_None()
231 if (*next_h > tleft) { in TSGLLEAdaptChoose_None()
233 *next_h = tleft; in TSGLLEAdaptChoose_None()
254 …etscInt cur, PetscReal h, PetscReal tleft, PetscInt *next_sc, PetscReal *next_h, PetscBool *finish) in TSGLLEAdaptChoose_Size() argument
268 if (last_desired_h > 1e-14) *next_h = PetscSqrtReal(last_desired_h * sz->desired_h); in TSGLLEAdaptChoose_Size()
269 else *next_h = sz->desired_h; in TSGLLEAdaptChoose_Size()
[all …]
|
| H A D | glle.c | 376 …ep_Rescale(TSGLLEScheme sc, PetscReal h, TSGLLEScheme next_sc, PetscReal next_h, Vec Ydot[], Vec X… in TSGLLECompleteStep_Rescale() argument
395 …eAndModify(TSGLLEScheme sc, PetscReal h, TSGLLEScheme next_sc, PetscReal next_h, Vec Ydot[], Vec X… in TSGLLECompleteStep_RescaleAndModify() argument
406 ratio = next_h / h; in TSGLLECompleteStep_RescaleAndModify()
790 …PetscReal h, const PetscReal hmnorm[], PetscInt *next_scheme, PetscReal *next_h, PetscBool *finish) in TSGLLEChooseNextScheme() argument
814 …PetscCall(TSGLLEAdaptChoose(gl->adapt, n, orders, errors, costs, cur, h, tleft, &next_sc, next_h, … in TSGLLEChooseNextScheme()
817 …gl->schemes[*next_scheme]->r, gl->schemes[*next_scheme]->s, (double)*next_h, PetscBools[*finish])); in TSGLLEChooseNextScheme()
872 PetscReal h, hmnorm[4], enorm[3], next_h; in TSSolve_GLLE() local
976 PetscCall(TSGLLEChooseNextScheme(ts, h, hmnorm, &next_scheme, &next_h, &final_step)); in TSSolve_GLLE()
980 next_h = h; in TSSolve_GLLE()
987 …PetscCall((*gl->CompleteStep)(scheme, h, gl->schemes[next_scheme], next_h, Ydot, gl->Xold, gl->X)); in TSSolve_GLLE()
[all …]
|
| /petsc/src/ts/adapt/impls/none/ |
| H A D | adaptnone.c | 3 …se_None(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_None() argument
8 *next_h = h; /* Reuse the old step */ in TSAdaptChoose_None()
|
| /petsc/src/ts/adapt/interface/ |
| H A D | tsadapt.c | 914 …tChoose(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept) in TSAdaptChoose() argument
926 PetscAssertPointer(next_h, 5); in TSAdaptChoose()
932 *next_h = h; in TSAdaptChoose()
953 PetscUseTypeMethod(adapt, choose, ts, h, &scheme, next_h, accept, &wlte, &wltea, &wlter); in TSAdaptChoose()
955 …tscCheck(*next_h >= 0, PetscObjectComm((PetscObject)adapt), PETSC_ERR_ARG_OUTOFRANGE, "Computed st… in TSAdaptChoose()
979 …if (h <= *next_h && *next_h <= adapt->dt_eval_times_cached) *next_h = adapt->dt_eval_times_cached;… in TSAdaptChoose()
980 h1 = *next_h; in TSAdaptChoose()
994 if (t < tmax && tend > tmax) *next_h = hmax; in TSAdaptChoose()
995 if (t < tmax && tend < tmax && h1 * b > hmax) *next_h = hmax / 2; in TSAdaptChoose()
996 if (t < tmax && tend < tmax && h1 * a > hmax) *next_h = hmax; in TSAdaptChoose()
[all …]
|
| /petsc/src/ts/adapt/impls/basic/ |
| H A D | adaptbasic.c | 4 …e_Basic(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_Basic() argument
35 *next_h = h; /* Reuse the old step */ in TSAdaptChoose_Basic()
67 *next_h = PetscClipInterval(h_lte, adapt->dt_min, adapt->dt_max); in TSAdaptChoose_Basic()
|
| /petsc/src/ts/adapt/impls/glee/ |
| H A D | adaptglee.c | 8 …se_GLEE(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_GLEE() argument
50 *next_h = h; /* Reuse the old step */ in TSAdaptChoose_GLEE()
93 *next_h = PetscClipInterval(h_lte, adapt->dt_min, adapt->dt_max); in TSAdaptChoose_GLEE()
107 *next_h = PetscClipInterval(h_lte, adapt->dt_min, adapt->dt_max); in TSAdaptChoose_GLEE()
|
| /petsc/src/ts/adapt/impls/cfl/ |
| H A D | adaptcfl.c | 3 …ose_CFL(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_CFL() argument
34 *next_h = PetscClipInterval(hcfl, adapt->dt_min, adapt->dt_max); in TSAdaptChoose_CFL()
|
| /petsc/src/ts/adapt/impls/history/ |
| H A D | adapthist.c | 8 …History(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_History() argument
16 PetscCall(TSHistoryGetTimeStep(thadapt->hist, thadapt->bw, step + 1, next_h)); in TSAdaptChoose_History()
|
| /petsc/src/ts/adapt/impls/dsp/ |
| H A D | adaptdsp.c | 68 …ose_DSP(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_DSP() argument
102 *next_h = h; /* Reuse the old step size */ in TSAdaptChoose_DSP()
165 *next_h = PetscClipInterval(hnew, adapt->dt_min, adapt->dt_max); in TSAdaptChoose_DSP()
|
| /petsc/src/ts/impls/pseudo/ |
| H A D | posindep.c | 354 …SPseudo(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_TSPseudo() argument
360 PetscCallBack("TSPSEUDO callback time step", (*pseudo->dt)(ts, next_h, pseudo->dtctx)); in TSAdaptChoose_TSPseudo()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex53.c | 2085 …erzaghi(TSAdapt adapt, TS ts, PetscReal h, PetscInt *next_sc, PetscReal *next_h, PetscBool *accept… in TSAdaptChoose_Terzaghi() argument
2101 *next_h = dtInitial; in TSAdaptChoose_Terzaghi()
2105 *next_h = dtTarget < 0.0 ? dtInitial : dtTarget; in TSAdaptChoose_Terzaghi()
2110 *next_h = h; in TSAdaptChoose_Terzaghi()
|