54   UserCtx  user;  in main()  local
60   PetscCall(DMCompositeCreate(PETSC_COMM_WORLD, &user.packer));  in main()
61   PetscCall(DMRedundantCreate(PETSC_COMM_WORLD, 0, 1, &user.red1));  in main()
62   PetscCall(DMCompositeAddDM(user.packer, user.red1));  in main()
63   PetscCall(DMDACreate1d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, 5, 1, 1, NULL, &user.da1));  in main()
64   PetscCall(DMSetFromOptions(user.da1));  in main()
65   PetscCall(DMSetUp(user.da1));  in main()
66   PetscCall(DMCompositeAddDM(user.packer, user.da1));  in main()
67   PetscCall(DMDACreate1d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, 5, 1, 1, NULL, &user.da2));  in main()
68   PetscCall(DMSetFromOptions(user.da2));  in main()
69   PetscCall(DMSetUp(user.da2));  in main()
70   PetscCall(DMDASetFieldName(user.da1, 0, "u"));  in main()
71   PetscCall(DMDASetFieldName(user.da2, 0, "lambda"));  in main()
72   PetscCall(DMCompositeAddDM(user.packer, user.da2));  in main()
73   PetscCall(DMCreateGlobalVector(user.packer, &U));  in main()
77 …PetscCall(PetscViewerDrawOpen(PETSC_COMM_WORLD, 0, "u - state variables", -1, -1, -1, -1, &user.u_…  in main()
78 …awOpen(PETSC_COMM_WORLD, 0, "lambda - Lagrange multipliers", -1, -1, -1, -1, &user.lambda_viewer));  in main()
79 …n(PETSC_COMM_WORLD, 0, "fu - derivative w.r.t. state variables", -1, -1, -1, -1, &user.fu_viewer));  in main()
80 …RLD, 0, "flambda - derivative w.r.t. Lagrange multipliers", -1, -1, -1, -1, &user.flambda_viewer));  in main()
84   PetscCall(SNESSetFunction(snes, FU, FormFunction, &user));  in main()
86   PetscCall(SNESMonitorSet(snes, Monitor, &user, 0));  in main()
91   PetscCall(DMDestroy(&user.red1));  in main()
92   PetscCall(DMDestroy(&user.da1));  in main()
93   PetscCall(DMDestroy(&user.da2));  in main()
94   PetscCall(DMDestroy(&user.packer));  in main()
97   PetscCall(PetscViewerDestroy(&user.u_viewer));  in main()
98   PetscCall(PetscViewerDestroy(&user.lambda_viewer));  in main()
99   PetscCall(PetscViewerDestroy(&user.fu_viewer));  in main()
100   PetscCall(PetscViewerDestroy(&user.flambda_viewer));  in main()
111   UserCtx     *user = (UserCtx *)dummy;  in FormFunction()  local
117   PetscCall(DMCompositeGetLocalVectors(user->packer, &vw, &vu, &vlambda));  in FormFunction()
118   PetscCall(DMCompositeGetLocalVectors(user->packer, &vfw, &vfu, &vflambda));  in FormFunction()
119   PetscCall(DMCompositeScatter(user->packer, U, vw, vu, vlambda));  in FormFunction()
121   PetscCall(DMDAGetCorners(user->da1, &xs, NULL, NULL, &xm, NULL, NULL));  in FormFunction()
122   PetscCall(DMDAGetInfo(user->da1, 0, &N, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0));  in FormFunction()
125   PetscCall(DMDAVecGetArray(user->da1, vu, &u));  in FormFunction()
126   PetscCall(DMDAVecGetArray(user->da1, vfu, &fu));  in FormFunction()
127   PetscCall(DMDAVecGetArray(user->da1, vlambda, &lambda));  in FormFunction()
128   PetscCall(DMDAVecGetArray(user->da1, vflambda, &flambda));  in FormFunction()
155   PetscCall(DMDAVecRestoreArray(user->da1, vu, &u));  in FormFunction()
156   PetscCall(DMDAVecRestoreArray(user->da1, vfu, &fu));  in FormFunction()
157   PetscCall(DMDAVecRestoreArray(user->da1, vlambda, &lambda));  in FormFunction()
158   PetscCall(DMDAVecRestoreArray(user->da1, vflambda, &flambda));  in FormFunction()
160   PetscCall(DMCompositeGather(user->packer, INSERT_VALUES, FU, vfw, vfu, vflambda));  in FormFunction()
161   PetscCall(DMCompositeRestoreLocalVectors(user->packer, &vw, &vu, &vlambda));  in FormFunction()
162   PetscCall(DMCompositeRestoreLocalVectors(user->packer, &vfw, &vfu, &vflambda));  in FormFunction()
168   UserCtx *user = (UserCtx *)dummy;  in Monitor()  local
173   PetscCall(DMCompositeGetAccess(user->packer, U, &w, &u, &lambda));  in Monitor()
174   PetscCall(VecView(u, user->u_viewer));  in Monitor()
175   PetscCall(VecView(lambda, user->lambda_viewer));  in Monitor()
176   PetscCall(DMCompositeRestoreAccess(user->packer, U, &w, &u, &lambda));  in Monitor()
179   PetscCall(DMCompositeGetAccess(user->packer, F, &w, &u, &lambda));  in Monitor()
180   PetscCall(VecView(u, user->fu_viewer));  in Monitor()
181   PetscCall(VecView(lambda, user->flambda_viewer));  in Monitor()
182   PetscCall(DMCompositeRestoreAccess(user->packer, F, &w, &u, &lambda));  in Monitor()