228   PetscDS        ds;  in SetupPrimalProblem()  local
236   PetscCall(DMGetDS(dm, &ds));  in SetupPrimalProblem()
237   PetscCall(PetscDSSetResidual(ds, 0, f0_q, f1_q));  in SetupPrimalProblem()
238   PetscCall(PetscDSSetJacobian(ds, 0, 0, g0_qq, NULL, NULL, NULL));  in SetupPrimalProblem()
239   PetscCall(PetscDSSetJacobian(ds, 0, 1, NULL, NULL, g2_qu, NULL));  in SetupPrimalProblem()
240   PetscCall(PetscDSSetJacobian(ds, 1, 0, NULL, g1_uq, NULL, NULL));  in SetupPrimalProblem()
243     PetscCall(PetscDSSetResidual(ds, 1, f0_linear_u, NULL));  in SetupPrimalProblem()
245 …PetscCall(PetscDSGetBoundary(ds, bd, &wf, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NU…  in SetupPrimalProblem()
247     PetscCall(PetscDSSetExactSolution(ds, 0, linear_q, user));  in SetupPrimalProblem()
248     PetscCall(PetscDSSetExactSolution(ds, 1, linear_u, user));  in SetupPrimalProblem()
251     PetscCall(PetscDSSetResidual(ds, 1, f0_quadratic_u, NULL));  in SetupPrimalProblem()
253 …PetscCall(PetscDSGetBoundary(ds, bd, &wf, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NU…  in SetupPrimalProblem()
255     PetscCall(PetscDSSetExactSolution(ds, 0, quadratic_q, user));  in SetupPrimalProblem()
256     PetscCall(PetscDSSetExactSolution(ds, 1, quadratic_u, user));  in SetupPrimalProblem()
259     PetscCall(PetscDSSetResidual(ds, 1, f0_quartic_u, NULL));  in SetupPrimalProblem()
260     PetscCall(PetscDSSetExactSolution(ds, 0, quartic_q, user));  in SetupPrimalProblem()
261     PetscCall(PetscDSSetExactSolution(ds, 1, quartic_u, user));  in SetupPrimalProblem()
264     PetscCall(PetscDSSetResidual(ds, 1, f0_quartic_u, NULL));  in SetupPrimalProblem()
265     PetscCall(PetscDSSetExactSolution(ds, 0, quartic_q, user));  in SetupPrimalProblem()
266     PetscCall(PetscDSSetExactSolution(ds, 1, quartic_u, user));  in SetupPrimalProblem()