What's the difference to decompose_epochs? Only that covariances are passed rather than computed on the fly?
I think if we should mingle this in decompose_epochs.