& QUIC: Sparse Inverse Covariance Estimation for a Million Variables

Neural Information Processing Systems 

However, it requires solving a difficult non-smooth log-determinant program with number of parameters scaling quadratically with the number of Gaussian variables. State-of-the-art methods thus do not scale to problems with more than 20, 000 variables.