Gradient of a scalar field

Is there an function/class that computes the gradient of a scalar field and the divergence of a vector field. For example, if a extract the density from the MultiBlockLattice3D using computeRho() and store it in MultiScalarField3D, is there an option where I can obtain the gradient of the density field?
Similarly, given a velocity field via a MultiTensorField3D, is there an in-built option to compute the divergence of the velocity.