Suppose I have a triangulated mesh and I have solved a differential equation problem and have a solution phi. From phi I can construct a vector field, how would I go about computing the gradient of this vector field?
The mesh looks something like this. phi is solved for at the center of all of the cells in the picture.
I'm not sure how to handle the spatial discretization for a triangulated lattice. Any references would be appreciated.