added fast calculation of residuals to OMP. At the moment I recompute the inverse of AtA at each step after exiting c++ so I do not take advantage of the cholesky updating. Speed could be improved by doing leave p out cross validation inside c++.