pca: use CUBLAS scale function instead of axpy
pca: add selector between float or double implementation of gemm function
pca: implement patch mean multiplication
pca: compute covariance matrix