#include "petscds.h" PetscErrorCode PetscDSSetDiscretization(PetscDS prob, PetscInt f, PetscObject disc)