Robust and scalable PCA using Grassmann averages, in C++ with Matlab bindings - View it on GitHub