3691   PetscInt Nf, Nfn, fn;  in PetscDSSelectDiscretizations()  local
3700   for (fn = 0; fn < numFields; ++fn) {  in PetscDSSelectDiscretizations()
3701     const PetscInt f = fields ? fields[fn] : fn;  in PetscDSSelectDiscretizations()
3712       PetscCall(PetscDSSetDiscretization(newprob, fn, (PetscObject)fe));  in PetscDSSelectDiscretizations()
3715       PetscCall(PetscDSSetDiscretization(newprob, fn, disc));  in PetscDSSelectDiscretizations()
3740   PetscInt Nf, Nfn, fn, gn;  in PetscDSSelectEquations()  local
3749   for (fn = 0; fn < numFields; ++fn) {  in PetscDSSelectEquations()
3750     const PetscInt  f = fields ? fields[fn] : fn;  in PetscDSSelectEquations()
3761     PetscCall(PetscDSSetObjective(newprob, fn, obj));  in PetscDSSelectEquations()
3762     PetscCall(PetscDSSetResidual(newprob, fn, f0, f1));  in PetscDSSelectEquations()
3763     PetscCall(PetscDSSetBdResidual(newprob, fn, f0Bd, f1Bd));  in PetscDSSelectEquations()
3764     PetscCall(PetscDSSetRiemannSolver(newprob, fn, r));  in PetscDSSelectEquations()
3775       PetscCall(PetscDSSetJacobian(newprob, fn, gn, g0, g1, g2, g3));  in PetscDSSelectEquations()
3776       PetscCall(PetscDSSetJacobianPreconditioner(newprob, fn, gn, g0p, g1p, g2p, g3p));  in PetscDSSelectEquations()
3777       PetscCall(PetscDSSetBdJacobian(newprob, fn, gn, g0Bd, g1Bd, g2Bd, g3Bd));  in PetscDSSelectEquations()