| /petsc/src/ts/trajectory/interface/ |
| H A D | traj.c | 52 PetscErrorCode TSTrajectorySet(TSTrajectory tj, TS ts, PetscInt stepnum, PetscReal time, Vec X) in TSTrajectorySet() argument
55 if (!tj) PetscFunctionReturn(PETSC_SUCCESS); in TSTrajectorySet()
56 PetscValidHeaderSpecific(tj, TSTRAJECTORY_CLASSID, 1); in TSTrajectorySet()
58 PetscValidLogicalCollectiveInt(tj, stepnum, 3); in TSTrajectorySet()
59 PetscValidLogicalCollectiveReal(tj, time, 4); in TSTrajectorySet()
61 …PetscCheck(tj->setupcalled, PetscObjectComm((PetscObject)tj), PETSC_ERR_ORDER, "TSTrajectorySetUp … in TSTrajectorySet()
62 …tj->monitor) PetscCall(PetscViewerASCIIPrintf(tj->monitor, "TSTrajectorySet: stepnum %" PetscInt_F… in TSTrajectorySet()
63 PetscCall(PetscLogEventBegin(TSTrajectory_Set, tj, ts, 0, 0)); in TSTrajectorySet()
64 PetscUseTypeMethod(tj, set, ts, stepnum, time, X); in TSTrajectorySet()
65 PetscCall(PetscLogEventEnd(TSTrajectory_Set, tj, ts, 0, 0)); in TSTrajectorySet()
[all …]
|
| /petsc/src/ts/trajectory/utils/ |
| H A D | reconstruct.c | 49 PetscErrorCode TSTrajectoryReconstruct_Private(TSTrajectory tj, TS ts, PetscReal t, Vec U, Vec Udot) in TSTrajectoryReconstruct_Private() argument
51 TSHistory tsh = tj->tsh; in TSTrajectoryReconstruct_Private()
62 …SETERRQ(PetscObjectComm((PetscObject)tj), PETSC_ERR_PLIB, "Requested time %g is outside the histor… in TSTrajectoryReconstruct_Private()
64 …if (tj->monitor) PetscCall(PetscViewerASCIIPrintf(tj->monitor, "Reconstructing at time %g, order %… in TSTrajectoryReconstruct_Private()
65 if (!tj->lag.T) { in TSTrajectoryReconstruct_Private()
66 PetscInt o = tj->lag.order + 1; in TSTrajectoryReconstruct_Private()
67 …PetscCall(PetscMalloc5(o, &tj->lag.L, o, &tj->lag.T, o, &tj->lag.WW, 2 * o, &tj->lag.TT, o, &tj->l… in TSTrajectoryReconstruct_Private()
68 for (i = 0; i < o; i++) tj->lag.T[i] = PETSC_MAX_REAL; in TSTrajectoryReconstruct_Private()
69 PetscCall(VecDuplicateVecs(U ? U : Udot, o, &tj->lag.W)); in TSTrajectoryReconstruct_Private()
72 PetscCall(PetscArrayzero(tj->lag.TT, 2 * (tj->lag.order + 1))); in TSTrajectoryReconstruct_Private()
[all …]
|
| /petsc/src/ts/trajectory/impls/basic/ |
| H A D | trajbasic.c | 11 static PetscErrorCode TSTrajectorySet_Basic(TSTrajectory tj, TS ts, PetscInt stepnum, PetscReal tim… in TSTrajectorySet_Basic() argument
13 TSTrajectory_Basic *tjbasic = (TSTrajectory_Basic *)tj->data; in TSTrajectorySet_Basic()
18 PetscCall(PetscSNPrintf(filename, sizeof(filename), tj->dirfiletemplate, stepnum)); in TSTrajectorySet_Basic()
23 if (stepnum && !tj->solution_only) { in TSTrajectorySet_Basic()
49 static PetscErrorCode TSTrajectorySetFromOptions_Basic(TSTrajectory tj, PetscOptionItems PetscOptio… in TSTrajectorySetFromOptions_Basic() argument
57 static PetscErrorCode TSTrajectoryGet_Basic(TSTrajectory tj, TS ts, PetscInt stepnum, PetscReal *t) in TSTrajectoryGet_Basic() argument
65 PetscCall(PetscSNPrintf(filename, sizeof(filename), tj->dirfiletemplate, stepnum)); in TSTrajectoryGet_Basic()
66 …PetscCall(PetscViewerBinaryOpen(PetscObjectComm((PetscObject)tj), filename, FILE_MODE_READ, &viewe… in TSTrajectoryGet_Basic()
71 if (stepnum && !tj->solution_only) { in TSTrajectoryGet_Basic()
81 if (tj->adjoint_solve_mode) PetscCall(TSSetTimeStep(ts, -(*t) + timepre)); in TSTrajectoryGet_Basic()
[all …]
|
| /petsc/src/ts/trajectory/impls/memory/ |
| H A D | trajmemory.c | 337 static PetscErrorCode StackDumpAll(TSTrajectory tj, TS ts, Stack *stack, PetscInt id) in StackDumpAll() argument
342 TJScheduler *tjsch = (TJScheduler *)tj->data; in StackDumpAll()
348 if (tj->monitor) { in StackDumpAll()
349 PetscCall(PetscViewerASCIIPushTab(tj->monitor)); in StackDumpAll()
350 PetscCall(PetscViewerASCIIPrintf(tj->monitor, "Dump stack id %" PetscInt_FMT " to file\n", id)); in StackDumpAll()
351 PetscCall(PetscViewerASCIIPopTab(tj->monitor)); in StackDumpAll()
353 …PetscCall(PetscSNPrintf(filename, sizeof(filename), "%s/TS-STACK%06" PetscInt_FMT ".bin", tj->dirn… in StackDumpAll()
362 PetscCall(PetscLogEventBegin(TSTrajectory_DiskWrite, tj, ts, 0, 0)); in StackDumpAll()
364 PetscCall(PetscLogEventEnd(TSTrajectory_DiskWrite, tj, ts, 0, 0)); in StackDumpAll()
370 PetscCall(PetscLogEventBegin(TSTrajectory_DiskWrite, tj, ts, 0, 0)); in StackDumpAll()
[all …]
|
| /petsc/src/ts/trajectory/impls/singlefile/ |
| H A D | singlefile.c | 7 static PetscErrorCode TSTrajectorySet_Singlefile(TSTrajectory tj, TS ts, PetscInt stepnum, PetscRea… in TSTrajectorySet_Singlefile() argument
9 TSTrajectory_Singlefile *sf = (TSTrajectory_Singlefile *)tj->data; in TSTrajectorySet_Singlefile()
17 PetscCall(PetscObjectGetName((PetscObject)tj, &filename)); in TSTrajectorySet_Singlefile()
25 static PetscErrorCode TSTrajectoryDestroy_Singlefile(TSTrajectory tj) in TSTrajectoryDestroy_Singlefile() argument
27 TSTrajectory_Singlefile *sf = (TSTrajectory_Singlefile *)tj->data; in TSTrajectoryDestroy_Singlefile()
43 PETSC_EXTERN PetscErrorCode TSTrajectoryCreate_Singlefile(TSTrajectory tj, TS ts) in TSTrajectoryCreate_Singlefile() argument
49 tj->data = sf; in TSTrajectoryCreate_Singlefile()
50 tj->ops->set = TSTrajectorySet_Singlefile; in TSTrajectoryCreate_Singlefile()
51 tj->ops->get = NULL; in TSTrajectoryCreate_Singlefile()
52 tj->ops->destroy = TSTrajectoryDestroy_Singlefile; in TSTrajectoryCreate_Singlefile()
|
| /petsc/src/ts/trajectory/impls/visualization/ |
| H A D | trajvisualization.c | 13 static PetscErrorCode TSTrajectorySet_Visualization(TSTrajectory tj, TS ts, PetscInt stepnum, Petsc… in TSTrajectorySet_Visualization() argument
29 if (tj->names) { in TSTrajectorySet_Visualization()
32 PetscCall(PetscViewerBinaryWriteStringArray(bnames, (const char *const *)tj->names)); in TSTrajectorySet_Visualization()
37 if (!tj->transform) { in TSTrajectorySet_Visualization()
41 PetscCall((*tj->transform)(tj->transformctx, X, &XX)); in TSTrajectorySet_Visualization()
51 if (!tj->transform) { in TSTrajectorySet_Visualization()
55 PetscCall((*tj->transform)(tj->transformctx, X, &XX)); in TSTrajectorySet_Visualization()
84 PETSC_EXTERN PetscErrorCode TSTrajectoryCreate_Visualization(TSTrajectory tj, TS ts) in TSTrajectoryCreate_Visualization() argument
87 tj->ops->set = TSTrajectorySet_Visualization; in TSTrajectoryCreate_Visualization()
88 tj->setupcalled = PETSC_TRUE; in TSTrajectoryCreate_Visualization()
|
| /petsc/src/ts/tests/ |
| H A D | ex13.c | 22 TSTrajectory tj; in main() local
41 PetscCall(TSGetTrajectory(ts, &tj)); in main()
42 PetscCall(TSTrajectorySetType(tj, ts, TSTRAJECTORYBASIC)); in main()
43 PetscCall(TSTrajectorySetFromOptions(tj, ts)); in main()
44 PetscCall(TSTrajectorySetSolutionOnly(tj, PETSC_TRUE)); in main()
45 PetscCall(TSTrajectorySetUp(tj, ts)); in main()
61 PetscCall(TSTrajectorySet(tj, ts, II[i], TT[i], W)); in main()
67 PetscCall(TSTrajectoryGetVecs(tj, ts, PETSC_DECIDE, &testtime, W, Wdot)); in main()
83 PetscCall(TSTrajectoryGetVecs(tj, ts, PETSC_DECIDE, &testtime, W, NULL)); in main()
94 PetscCall(TSTrajectoryGetVecs(tj, ts, PETSC_DECIDE, &testtime, NULL, Wdot)); in main()
[all …]
|
| /petsc/src/ts/adapt/impls/history/ |
| H A D | adapthist.c | 152 PetscErrorCode TSAdaptHistorySetTrajectory(TSAdapt adapt, TSTrajectory tj, PetscBool backward) in TSAdaptHistorySetTrajectory() argument
158 PetscValidHeaderSpecific(tj, TSTRAJECTORY_CLASSID, 2); in TSAdaptHistorySetTrajectory()
162 PetscCall(TSAdaptHistorySetTSHistory(adapt, tj->tsh, backward)); in TSAdaptHistorySetTrajectory()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex41.c | 238 TSTrajectory tj; in main() local
250 PetscCall(TSGetTrajectory(ts, &tj)); in main()
251 PetscCall(TSAdaptHistorySetTrajectory(adapt, tj, PETSC_FALSE)); in main()
|
| H A D | extchem.c | 81 TSTrajectory tj; in main() local
171 PetscCall(TSGetTrajectory(ts, &tj)); in main()
172 if (tj) { in main()
173 PetscCall(TSTrajectorySetVariableNames(tj, (const char *const *)user.snames)); in main()
174 …PetscCall(TSTrajectorySetTransform(tj, (PetscErrorCode (*)(void *, Vec, Vec *))MassFractionToMoleF… in main()
|
| H A D | ex40.c | 253 TSTrajectory tj; in main() local
265 PetscCall(TSGetTrajectory(ts, &tj)); in main()
266 PetscCall(TSAdaptHistorySetTrajectory(adapt, tj, PETSC_FALSE)); in main()
|
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 1073 TSTrajectory tj; in TSAdjointSetUp() local
1081 PetscCall(TSGetTrajectory(ts, &tj)); in TSAdjointSetUp()
1082 PetscCall(PetscObjectTypeCompare((PetscObject)tj, TSTRAJECTORYBASIC, &match)); in TSAdjointSetUp()
1085 PetscCall(TSTrajectoryGetSolutionOnly(tj, &solution_only)); in TSAdjointSetUp()
1088 PetscCall(TSTrajectorySetUseHistory(tj, PETSC_FALSE)); /* not use TSHistory */ in TSAdjointSetUp()
|