WIP: add option for automated regularization to mne.cov.compute_covariance

Todo:

- [x] show basic sensor examples / update existing examples
- [x] consolidate API including default values
- [x] deal with projs
why you need to do is:

compute the data_mean for all type of epochs. Subtract the mean of each type of epochs to each epoch in this Epochs. Then concatenate the data from all epochs to run cv on it.