121   PetscDS        ds;  in SetupProblem()  local
128   PetscCall(DMGetDS(dm, &ds));  in SetupProblem()
130   PetscCall(PetscDSSetResidual(ds, 0, f0_v, f1_v));  in SetupProblem()
131   PetscCall(PetscDSSetResidual(ds, 1, f0_q_linear, NULL));  in SetupProblem()
132   PetscCall(PetscDSSetJacobian(ds, 0, 0, g0_vu, NULL, NULL, NULL));  in SetupProblem()
133   PetscCall(PetscDSSetJacobian(ds, 0, 1, NULL, NULL, g2_vp, NULL));  in SetupProblem()
134   PetscCall(PetscDSSetJacobian(ds, 1, 0, NULL, g1_qu, NULL, NULL));  in SetupProblem()
137 …PetscCall(PetscDSAddBoundary(ds, DM_BC_NATURAL, "Boundary Integral", label, 1, &id, 0, 0, NULL, (P…  in SetupProblem()
138 …PetscCall(PetscDSGetBoundary(ds, bd, &wf, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NU…  in SetupProblem()
141   PetscCall(PetscDSSetExactSolution(ds, 0, linear_u, NULL));  in SetupProblem()
142   PetscCall(PetscDSSetExactSolution(ds, 1, linear_p, NULL));  in SetupProblem()
150   PetscDS   ds;  in SetupDiscretization()  local
162     PetscCall(DMGetDS(mesh_sum, &ds));  in SetupDiscretization()
163     PetscCall(PetscDSGetForceQuad(ds, &force));  in SetupDiscretization()
164     if (force) PetscCall(PetscDSSetForceQuad(ds, PETSC_FALSE));  in SetupDiscretization()