7   PetscInt  nt;           /*  Number of time steps; as of now, must be divisible by 8 */  member
112   user.nt = 8;  in main()
113   PetscCall(PetscOptionsInt("-nt", "Number of time steps", "", user.nt, &user.nt, NULL));  in main()
128   user.n  = user.m * (user.nt + 1);      /*  number of variables */  in main()
129   user.ht = (PetscReal)1 / user.nt;  in main()
134   PetscCall(VecSetSizes(user.u, PETSC_DECIDE, user.n - user.m * user.nt));  in main()
135   PetscCall(VecSetSizes(user.y, PETSC_DECIDE, user.m * user.nt));  in main()
136   PetscCall(VecSetSizes(user.c, PETSC_DECIDE, user.m * user.nt));  in main()
233   PetscCall(Scatter_i(user->y, user->yi, user->yi_scatter, user->nt));  in FormFunction()
262   PetscCall(Scatter_i(user->y, user->yi, user->yi_scatter, user->nt));  in FormGradient()
271   PetscCall(Scatter_i(user->ywork, user->yiwork, user->yi_scatter, user->nt));  in FormGradient()
276   PetscCall(Gather_i(user->ywork, user->yiwork, user->yi_scatter, user->nt));  in FormGradient()
294   PetscCall(Scatter_i(user->y, user->yi, user->yi_scatter, user->nt));  in FormFunctionGradient()
304   PetscCall(Scatter_i(user->ywork, user->yiwork, user->yi_scatter, user->nt));  in FormFunctionGradient()
309   PetscCall(Gather_i(user->ywork, user->yiwork, user->yi_scatter, user->nt));  in FormFunctionGradient()
371   PetscCall(Scatter_i(X, user->yi, user->yi_scatter, user->nt));  in StateMatMult()
373   for (i = 1; i < user->nt; i++) {  in StateMatMult()
377   PetscCall(Gather_i(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatMult()
388   PetscCall(Scatter_i(X, user->yi, user->yi_scatter, user->nt));  in StateMatMultTranspose()
389   for (i = 0; i < user->nt - 1; i++) {  in StateMatMultTranspose()
393   i = user->nt - 1;  in StateMatMultTranspose()
395   PetscCall(Gather_i(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatMultTranspose()
437   PetscCall(Scatter_i(user->y, user->yi, user->yi_scatter, user->nt));  in DesignMatMult()
438   for (i = 0; i < user->nt; i++) {  in DesignMatMult()
447   PetscCall(Gather_i(Y, user->yiwork, user->yi_scatter, user->nt));  in DesignMatMult()
468   PetscCall(Scatter_i(user->y, user->yi, user->yi_scatter, user->nt));  in DesignMatMultTranspose()
469   PetscCall(Scatter_i(X, user->yiwork, user->yi_scatter, user->nt));  in DesignMatMultTranspose()
470   for (i = 0; i < user->nt; i++) {  in DesignMatMultTranspose()
520   PetscCall(Scatter_i(X, user->yi, user->yi_scatter, user->nt));  in StateMatInvMult()
525   for (i = 1; i < user->nt; i++) {  in StateMatInvMult()
531   PetscCall(Gather_i(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatInvMult()
543   PetscCall(Scatter_i(X, user->yi, user->yi_scatter, user->nt));  in StateMatInvTransposeMult()
545   i = user->nt - 1;  in StateMatInvTransposeMult()
551   for (i = user->nt - 2; i >= 0; i--) {  in StateMatInvTransposeMult()
559   PetscCall(Gather_i(Y, user->yiwork, user->yi_scatter, user->nt));  in StateMatInvTransposeMult()
601   PetscCall(Scatter_i(user->y, user->yi, user->yi_scatter, user->nt));  in FormConstraints()
603   for (i = 1; i < user->nt; i++) {  in FormConstraints()
607   PetscCall(Gather_i(C, user->yiwork, user->yi_scatter, user->nt));  in FormConstraints()
622 PetscErrorCode Scatter_i(Vec y, Vec *yi, VecScatter *scat, PetscInt nt)  in Scatter_i()  argument
627   for (i = 0; i < nt; i++) {  in Scatter_i()
644 PetscErrorCode Gather_i(Vec y, Vec *yi, VecScatter *scat, PetscInt nt)  in Gather_i()  argument
649   for (i = 0; i < nt; i++) {  in Gather_i()
694   stime = (PetscReal)user->nt / user->ns;  in ParabolicInitialize()
701   PetscCall(VecSetSizes(user->q, PETSC_DECIDE, n * user->nt));  in ParabolicInitialize()
904   PetscCall(VecSetSizes(user->d, PETSC_DECIDE, user->ndata * user->nt));  in ParabolicInitialize()
918 …SC_COMM_WORLD, PETSC_DETERMINE, PETSC_DETERMINE, user->m * user->nt, user->m * user->nt, user, &us…  in ParabolicInitialize()
941 …PetscCall(MatCreateShell(PETSC_COMM_WORLD, PETSC_DETERMINE, PETSC_DETERMINE, user->m * user->nt, u…  in ParabolicInitialize()
946 …SC_COMM_WORLD, PETSC_DETERMINE, PETSC_DETERMINE, user->m * user->nt, user->m * user->nt, user, &us…  in ParabolicInitialize()
965   PetscCall(PetscMalloc1(user->nt * user->m, &user->yi_scatter));  in ParabolicInitialize()
969   PetscCall(VecDuplicateVecs(yi, user->nt, &user->yi));  in ParabolicInitialize()
970   PetscCall(VecDuplicateVecs(yi, user->nt, &user->yiwork));  in ParabolicInitialize()
974   for (i = 0; i < user->nt; i++) {  in ParabolicInitialize()
1006   for (i = 1; i < user->nt; i++) PetscCall(VecSet(user->yiwork[i], 0.0));  in ParabolicInitialize()
1007   PetscCall(Gather_i(user->q, user->yiwork, user->yi_scatter, user->nt));  in ParabolicInitialize()
1012   PetscCall(VecSetSizes(user->ytrue, PETSC_DECIDE, n * user->nt));  in ParabolicInitialize()
1171   PetscCall(Scatter_i(user->ywork, user->yiwork, user->yi_scatter, user->nt));  in ParabolicInitialize()
1214   PetscCall(VecDestroyVecs(user->nt, &user->yi));  in ParabolicDestroy()
1215   PetscCall(VecDestroyVecs(user->nt, &user->yiwork));  in ParabolicDestroy()
1237   for (i = 0; i < user->nt; i++) PetscCall(VecScatterDestroy(&user->yi_scatter[i]));  in ParabolicDestroy()