thanks for a great C++ implementation! Everything is running smooth for me but I had some questions. Haven't gone in detail of the code yet (sitting with it next to the original matlab implementation) but is there a member function for obtaining the velocity field for any given set of points given by a 2D matrix?

Would be highly useful but have a hard time finding where in the code that would be done. Maybe a pointer to where I should start looking in the code.

As thanks I can contribute some more examples and tutorials. I have successfully written a wrapper of cpd to R for some basic functionality using Rcpp which might help the scientific community to adopt the C++ library faster.