84       PetscInt  tid = LagrangeGetId(t, tj->lag.order + 1, tj->lag.T, tj->lag.TT);  in TSTrajectoryReconstruct_Private()  local
85       if (tid < 0) continue;  in TSTrajectoryReconstruct_Private()
86 …, "Reusing snapshot %" PetscInt_FMT ", step %" PetscInt_FMT ", time %g\n", tid, tshhist_id[s], (do…  in TSTrajectoryReconstruct_Private()
87       tj->lag.TT[tid]                         = PETSC_TRUE;  in TSTrajectoryReconstruct_Private()
88       tj->lag.WW[cnt]                         = tj->lag.W[tid];  in TSTrajectoryReconstruct_Private()
97       PetscInt  tid;  in TSTrajectoryReconstruct_Private()  local
100       tid = LagrangeGetId(t, tj->lag.order + 1, tj->lag.T, tj->lag.TT);  in TSTrajectoryReconstruct_Private()
101 …   PetscCheck(tid < 0, PetscObjectComm((PetscObject)tj), PETSC_ERR_PLIB, "This should not happen");  in TSTrajectoryReconstruct_Private()
102       tid = -tid - 1;  in TSTrajectoryReconstruct_Private()
104         if (tj->lag.T[tid] < PETSC_MAX_REAL) {  in TSTrajectoryReconstruct_Private()
105 …f(tj->monitor, "Discarding snapshot %" PetscInt_FMT " at time %g\n", tid, (double)tj->lag.T[tid]));  in TSTrajectoryReconstruct_Private()
107           PetscCall(PetscViewerASCIIPrintf(tj->monitor, "New snapshot %" PetscInt_FMT "\n", tid));  in TSTrajectoryReconstruct_Private()
111       PetscCall(TSTrajectoryGetVecs(tj, ts, tshhist_id[s], &t, tj->lag.W[tid], NULL));  in TSTrajectoryReconstruct_Private()
112       tj->lag.T[tid] = t;  in TSTrajectoryReconstruct_Private()
114       tj->lag.TT[tid]                         = PETSC_TRUE;  in TSTrajectoryReconstruct_Private()
115       tj->lag.WW[cnt]                         = tj->lag.W[tid];  in TSTrajectoryReconstruct_Private()
124     PetscInt tid = LagrangeGetId(t, tj->lag.order + 1, tj->lag.T, tj->lag.TT);  in TSTrajectoryReconstruct_Private()  local
129     if (tid < 0) {  in TSTrajectoryReconstruct_Private()
130       tid = -tid - 1;  in TSTrajectoryReconstruct_Private()
132         if (tj->lag.T[tid] < PETSC_MAX_REAL) {  in TSTrajectoryReconstruct_Private()
133 …f(tj->monitor, "Discarding snapshot %" PetscInt_FMT " at time %g\n", tid, (double)tj->lag.T[tid]));  in TSTrajectoryReconstruct_Private()
135           PetscCall(PetscViewerASCIIPrintf(tj->monitor, "New snapshot %" PetscInt_FMT "\n", tid));  in TSTrajectoryReconstruct_Private()
139       PetscCall(TSTrajectoryGetVecs(tj, ts, tshhist_id[id], &t, tj->lag.W[tid], NULL));  in TSTrajectoryReconstruct_Private()
141       tj->lag.T[tid] = t;  in TSTrajectoryReconstruct_Private()
143 …r, "Reusing snapshot %" PetscInt_FMT " step %" PetscInt_FMT ", time %g\n", tid, tshhist_id[id], (d…  in TSTrajectoryReconstruct_Private()
145     PetscCall(VecCopy(tj->lag.W[tid], U));  in TSTrajectoryReconstruct_Private()