for debugging: add a method to compute gradients of local geodesic fe functions by finite-difference approximation (only 2d) [[Imported from SVN: r4065]]
