Follow GeoMTK changes and use skeleton to calculate deformation matrix

Because the deformation matrix calculated through velocity gradient
tensor is not accurate, I use skeleton points (advected with centroid)
to fit the matrix, and the results look better.