96   DMLabel        label;  in SetupPrimalProblem()  local
100   PetscCall(DMGetRegionNumDS(dm, 0, &label, NULL, &ds, NULL));  in SetupPrimalProblem()
102   PetscCall(PetscWeakFormSetIndexResidual(wf, label, 1, 0, 0, 0, f0_quad_u, 0, f1_u));  in SetupPrimalProblem()
103 …PetscCall(PetscWeakFormSetIndexJacobian(wf, label, 1, 0, 0, 0, 0, NULL, 0, NULL, 0, NULL, 0, g3_uu…  in SetupPrimalProblem()
105   PetscCall(DMGetRegionNumDS(dm, 1, &label, NULL, &ds, NULL));  in SetupPrimalProblem()
107   PetscCall(PetscWeakFormSetIndexResidual(wf, label, 1, 0, 0, 0, f0_quad_u, 0, f1_u));  in SetupPrimalProblem()
108 …PetscCall(PetscWeakFormSetIndexJacobian(wf, label, 1, 0, 0, 0, 0, NULL, 0, NULL, 0, NULL, 0, g3_uu…  in SetupPrimalProblem()
109   PetscCall(PetscWeakFormSetIndexResidual(wf, label, 1, 1, 0, 0, f0_quad_p, 0, NULL));  in SetupPrimalProblem()
110 …PetscCall(PetscWeakFormSetIndexJacobian(wf, label, 1, 1, 1, 0, 0, g0_pp, 0, NULL, 0, NULL, 0, NULL…  in SetupPrimalProblem()
113   PetscCall(DMGetLabel(dm, "marker", &label));  in SetupPrimalProblem()
114 …PetscCall(DMAddBoundary(dm, DM_BC_ESSENTIAL, "wall", label, 1, &id, 0, 0, NULL, (PetscVoidFn *)qua…  in SetupPrimalProblem()