[WIP] Prediction variance in Ridge

I started working on returning prediction variances in Ridge (for the Cholesky solver only). The n_samples > n_features (primal) case works. The n_features > n_samples (dual) case is a work in progress.