@karlnapf sorry forgot to mention it before. I tested it against [Matlab implementation](http://www.mathworks.in/help/stats/princomp.html). I also validated the results using the brml code that you pointed to. BRML code uses SVD always but results match with Matlab implementation and our Shogun implementation.