Is there a reason why PetscDSSetDiscretization takes PetscObject instead of PetscFE? The whole FEM API changed since I last looked at it, and having untyped routines makes it somewhat difficult to track the changes. Geoffrey