338   PetscDS        prob;  in SetupProblem()  local
343   PetscCall(DMGetDS(dm, &prob));  in SetupProblem()
345   PetscCall(PetscDSSetResidual(prob, 1, f0_q, NULL));  in SetupProblem()
346   PetscCall(PetscDSSetResidual(prob, 2, f0_w, NULL));  in SetupProblem()
347   PetscCall(PetscDSSetJacobian(prob, 0, 0, g0_vu, NULL, NULL, NULL));  in SetupProblem()
348   PetscCall(PetscDSSetJacobian(prob, 1, 0, NULL, g1_qu, NULL, NULL));  in SetupProblem()
349   PetscCall(PetscDSSetJacobian(prob, 1, 1, g0_qp, NULL, NULL, NULL));  in SetupProblem()
350   PetscCall(PetscDSSetJacobian(prob, 2, 0, NULL, g1_wu, NULL, NULL));  in SetupProblem()
351   PetscCall(PetscDSSetJacobian(prob, 2, 1, g0_wp, NULL, NULL, NULL));  in SetupProblem()
352   PetscCall(PetscDSSetJacobian(prob, 2, 2, g0_wd, NULL, NULL, NULL));  in SetupProblem()
356   PetscCall(PetscDSSetExactSolution(prob, 2, zero_func, NULL));  in SetupProblem()
360     PetscCall(PetscDSSetResidual(prob, 0, f0_v_linear, NULL));  in SetupProblem()
361     PetscCall(PetscDSSetBdResidual(prob, 0, f0_bd_u_linear, NULL));  in SetupProblem()
362     PetscCall(PetscDSSetExactSolution(prob, 0, linear_u, NULL));  in SetupProblem()
363     PetscCall(PetscDSSetExactSolution(prob, 1, linear_p, NULL));  in SetupProblem()
366     PetscCall(PetscDSSetResidual(prob, 0, f0_v_sinusoid, NULL));  in SetupProblem()
367     PetscCall(PetscDSSetBdResidual(prob, 0, f0_bd_u_sinusoid, NULL));  in SetupProblem()
368     PetscCall(PetscDSSetExactSolution(prob, 0, sinusoid_u, NULL));  in SetupProblem()
369     PetscCall(PetscDSSetExactSolution(prob, 1, sinusoid_p, NULL));  in SetupProblem()
376 …PetscCall(PetscDSAddBoundary(prob, DM_BC_NATURAL, "Boundary Integral", label, 1, &id, 0, 0, NULL, …  in SetupProblem()