In legacy DOLFIN you would have to project the quantity (grad(u)) into an appropriate function space instead of interpolating it as in DOLFINx. Please note that we advice all new users to use DOLFINx: The new DOLFINx solver is now recommended over DOLFIN