7   PetscInt    nt;    /*  Number of time steps */  member
110   user.nt = 16;  in main()
111   PetscCall(PetscOptionsInt("-nt", "Number of time steps", "", user.nt, &user.nt, NULL));  in main()
121   user.m     = user.mx * user.mx * user.nt;     /*  number of constraints */  in main()
122   user.n     = user.mx * user.mx * 3 * user.nt; /*  number of variables */  in main()
123   user.ht    = user.T / user.nt;                /*  Time step */  in main()
306   PetscCall(Scatter_yi(user->u, user->ui, user->ui_scatter, user->nt));  in FormJacobianState()
307 …ll(Scatter_uxi_uyi(user->u, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in FormJacobianState()
308   for (i = 0; i < user->nt; i++) {  in FormJacobianState()
339   PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt));  in StateMatMult()
343   for (i = 1; i < user->nt; i++) {  in StateMatMult()
349   PetscCall(Gather_yi(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatMult()
360   PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt));  in StateMatMultTranspose()
362   for (i = 0; i < user->nt - 1; i++) {  in StateMatMultTranspose()
369   i                 = user->nt - 1;  in StateMatMultTranspose()
372   PetscCall(Gather_yi(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatMultTranspose()
416   PetscCall(Scatter_yi(user->y, user->yi, user->yi_scatter, user->nt));  in DesignMatMult()
417 …(Scatter_uxi_uyi(X, user->uxiwork, user->uxi_scatter, user->uyiwork, user->uyi_scatter, user->nt));  in DesignMatMult()
418   for (i = 0; i < user->nt; i++) {  in DesignMatMult()
425   PetscCall(Gather_yi(Y, user->ziwork, user->yi_scatter, user->nt));  in DesignMatMult()
436   PetscCall(Scatter_yi(user->y, user->yi, user->yi_scatter, user->nt));  in DesignMatMultTranspose()
437   PetscCall(Scatter_yi(X, user->yiwork, user->yi_scatter, user->nt));  in DesignMatMultTranspose()
438   for (i = 0; i < user->nt; i++) {  in DesignMatMultTranspose()
446   PetscCall(Gather_yi(Y, user->uiwork, user->ui_scatter, user->nt));  in DesignMatMultTranspose()
491   PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt));  in StateMatInvMult()
492   PetscCall(Scatter_yi(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatInvMult()
493 …ll(Scatter_uxi_uyi(user->u, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in StateMatInvMult()
500   for (i = 1; i < user->nt; i++) {  in StateMatInvMult()
509   PetscCall(Gather_yi(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatInvMult()
521   PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt));  in StateMatInvTransposeMult()
522   PetscCall(Scatter_yi(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatInvTransposeMult()
523 …ll(Scatter_uxi_uyi(user->u, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in StateMatInvTransposeMult()
525   i                 = user->nt - 1;  in StateMatInvTransposeMult()
532   for (i = user->nt - 2; i >= 0; i--) {  in StateMatInvTransposeMult()
541   PetscCall(Gather_yi(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatInvTransposeMult()
583   PetscCall(Scatter_yi(user->y, user->yi, user->yi_scatter, user->nt));  in FormConstraints()
584 …ll(Scatter_uxi_uyi(user->u, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in FormConstraints()
589   for (i = 1; i < user->nt; i++) {  in FormConstraints()
596   PetscCall(Gather_yi(C, user->yiwork, user->yi_scatter, user->nt));  in FormConstraints()
611 …rCode Scatter_uxi_uyi(Vec u, Vec *uxi, VecScatter *scatx, Vec *uyi, VecScatter *scaty, PetscInt nt)  in Scatter_uxi_uyi()  argument
616   for (i = 0; i < nt; i++) {  in Scatter_uxi_uyi()
635 …orCode Gather_uxi_uyi(Vec u, Vec *uxi, VecScatter *scatx, Vec *uyi, VecScatter *scaty, PetscInt nt)  in Gather_uxi_uyi()  argument
640   for (i = 0; i < nt; i++) {  in Gather_uxi_uyi()
649 PetscErrorCode Scatter_yi(Vec y, Vec *yi, VecScatter *scat, PetscInt nt)  in Scatter_yi()  argument
654   for (i = 0; i < nt; i++) {  in Scatter_yi()
661 PetscErrorCode Gather_yi(Vec y, Vec *yi, VecScatter *scat, PetscInt nt)  in Gather_yi()  argument
666   for (i = 0; i < nt; i++) {  in Gather_yi()
793   PetscCall(VecSetSizes(user->q, PETSC_DECIDE, n * user->nt));  in HyperbolicInitialize()
875   PetscCall(PetscMalloc1(user->nt * user->mx * user->mx, &user->yi_scatter));  in HyperbolicInitialize()
879   PetscCall(VecDuplicateVecs(yi, user->nt, &user->yi));  in HyperbolicInitialize()
880   PetscCall(VecDuplicateVecs(yi, user->nt, &user->yiwork));  in HyperbolicInitialize()
881   PetscCall(VecDuplicateVecs(yi, user->nt, &user->ziwork));  in HyperbolicInitialize()
882   for (i = 0; i < user->nt; i++) {  in HyperbolicInitialize()
893   PetscCall(PetscMalloc1(user->nt * user->mx * user->mx, &user->uxi_scatter));  in HyperbolicInitialize()
894   PetscCall(PetscMalloc1(user->nt * user->mx * user->mx, &user->uyi_scatter));  in HyperbolicInitialize()
895   PetscCall(PetscMalloc1(user->nt * user->mx * user->mx, &user->ux_scatter));  in HyperbolicInitialize()
896   PetscCall(PetscMalloc1(user->nt * user->mx * user->mx, &user->uy_scatter));  in HyperbolicInitialize()
897   PetscCall(PetscMalloc1(2 * user->nt * user->mx * user->mx, &user->ui_scatter));  in HyperbolicInitialize()
904   PetscCall(VecDuplicateVecs(uxi, user->nt, &user->uxi));  in HyperbolicInitialize()
905   PetscCall(VecDuplicateVecs(uxi, user->nt, &user->uyi));  in HyperbolicInitialize()
906   PetscCall(VecDuplicateVecs(uxi, user->nt, &user->uxiwork));  in HyperbolicInitialize()
907   PetscCall(VecDuplicateVecs(uxi, user->nt, &user->uyiwork));  in HyperbolicInitialize()
908   PetscCall(VecDuplicateVecs(ui, user->nt, &user->ui));  in HyperbolicInitialize()
909   PetscCall(VecDuplicateVecs(ui, user->nt, &user->uiwork));  in HyperbolicInitialize()
910   for (i = 0; i < user->nt; i++) {  in HyperbolicInitialize()
954   for (i = 1; i < user->nt; i++) PetscCall(VecSet(user->yiwork[i], 0.0));  in HyperbolicInitialize()
955   PetscCall(Gather_yi(user->q, user->yiwork, user->yi_scatter, user->nt));  in HyperbolicInitialize()
959   for (i = 0; i < user->nt; i++) {  in HyperbolicInitialize()
966 …Gather_uxi_uyi(user->utrue, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in HyperbolicInitialize()
970   for (i = 0; i < user->nt; i++) {  in HyperbolicInitialize()
976 …ll(Gather_uxi_uyi(user->ur, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in HyperbolicInitialize()
980   PetscCall(MatSetSizes(user->LT, PETSC_DECIDE, PETSC_DECIDE, 2 * n * user->nt, n * user->nt));  in HyperbolicInitialize()
1041 …ll(Scatter_uxi_uyi(user->u, user->uxi, user->uxi_scatter, user->uyi, user->uyi_scatter, user->nt));  in HyperbolicInitialize()
1044   for (i = 0; i < user->nt; i++) {  in HyperbolicInitialize()
1070   PetscCall(VecSetSizes(user->ytrue, PETSC_DECIDE, n * user->nt));  in HyperbolicInitialize()
1137   for (i = 0; i < user->nt; i++) {  in HyperbolicDestroy()
1150   PetscCall(VecDestroyVecs(user->nt, &user->yi));  in HyperbolicDestroy()
1151   PetscCall(VecDestroyVecs(user->nt, &user->yiwork));  in HyperbolicDestroy()
1152   PetscCall(VecDestroyVecs(user->nt, &user->ziwork));  in HyperbolicDestroy()
1153   PetscCall(VecDestroyVecs(user->nt, &user->uxi));  in HyperbolicDestroy()
1154   PetscCall(VecDestroyVecs(user->nt, &user->uyi));  in HyperbolicDestroy()
1155   PetscCall(VecDestroyVecs(user->nt, &user->uxiwork));  in HyperbolicDestroy()
1156   PetscCall(VecDestroyVecs(user->nt, &user->uyiwork));  in HyperbolicDestroy()
1157   PetscCall(VecDestroyVecs(user->nt, &user->ui));  in HyperbolicDestroy()
1158   PetscCall(VecDestroyVecs(user->nt, &user->uiwork));  in HyperbolicDestroy()
1171   for (i = 0; i < user->nt; i++) {  in HyperbolicDestroy()