205   TS_BasicSymplectic   *bsymp    = (TS_BasicSymplectic *)ts->data;  in TSStep_BasicSymplectic()  local
206   BasicSymplecticScheme scheme   = bsymp->scheme;  in TSStep_BasicSymplectic()
207   Vec                   solution = ts->vec_sol, update = bsymp->update, q, p, q_update, p_update;  in TSStep_BasicSymplectic()
208   IS                    is_q = bsymp->is_q, is_p = bsymp->is_p;  in TSStep_BasicSymplectic()
209   TS                    subts_q = bsymp->subts_q, subts_p = bsymp->subts_p;  in TSStep_BasicSymplectic()
281   TS_BasicSymplectic *bsymp = (TS_BasicSymplectic *)ts->data;  in TSSetUp_BasicSymplectic()  local
285   PetscCall(TSRHSSplitGetIS(ts, "position", &bsymp->is_q));  in TSSetUp_BasicSymplectic()
286   PetscCall(TSRHSSplitGetIS(ts, "momentum", &bsymp->is_p));  in TSSetUp_BasicSymplectic()
287 …PetscCheck(bsymp->is_q && bsymp->is_p, PetscObjectComm((PetscObject)ts), PETSC_ERR_USER, "Must set…  in TSSetUp_BasicSymplectic()
288   PetscCall(TSRHSSplitGetSubTS(ts, "position", &bsymp->subts_q));  in TSSetUp_BasicSymplectic()
289   PetscCall(TSRHSSplitGetSubTS(ts, "momentum", &bsymp->subts_p));  in TSSetUp_BasicSymplectic()
290 …PetscCheck(bsymp->subts_q && bsymp->subts_p, PetscObjectComm((PetscObject)ts), PETSC_ERR_USER, "Mu…  in TSSetUp_BasicSymplectic()
292   PetscCall(VecDuplicate(ts->vec_sol, &bsymp->update));  in TSSetUp_BasicSymplectic()
306   TS_BasicSymplectic *bsymp = (TS_BasicSymplectic *)ts->data;  in TSReset_BasicSymplectic()  local
309   PetscCall(VecDestroy(&bsymp->update));  in TSReset_BasicSymplectic()
325   TS_BasicSymplectic *bsymp = (TS_BasicSymplectic *)ts->data;  in TSSetFromOptions_BasicSymplectic()  local
338 …ethod", "TSBasicSymplecticSetType", (const char *const *)namelist, count, bsymp->scheme->name, &ch…  in TSSetFromOptions_BasicSymplectic()
348   TS_BasicSymplectic *bsymp  = (TS_BasicSymplectic *)ts->data;  in TSInterpolate_BasicSymplectic()  local
349   Vec                 update = bsymp->update;  in TSInterpolate_BasicSymplectic()
418   TS_BasicSymplectic       *bsymp = (TS_BasicSymplectic *)ts->data;  in TSBasicSymplecticSetType_BasicSymplectic()  local
423   if (bsymp->scheme) {  in TSBasicSymplecticSetType_BasicSymplectic()
424     PetscCall(PetscStrcmp(bsymp->scheme->name, bsymptype, &match));  in TSBasicSymplecticSetType_BasicSymplectic()
430       bsymp->scheme = &link->sch;  in TSBasicSymplecticSetType_BasicSymplectic()
439   TS_BasicSymplectic *bsymp = (TS_BasicSymplectic *)ts->data;  in TSBasicSymplecticGetType_BasicSymplectic()  local
442   *bsymptype = bsymp->scheme->name;  in TSBasicSymplecticGetType_BasicSymplectic()
493   TS_BasicSymplectic *bsymp;  in TSCreate_BasicSymplectic()  local
497   PetscCall(PetscNew(&bsymp));  in TSCreate_BasicSymplectic()
498   ts->data = (void *)bsymp;  in TSCreate_BasicSymplectic()