Several methods have been recently proposed for estimating sparse Gaussian
graphical models using $\ell_{1}$ regularization on the inverse covariance
matrix. Despite recent advances, contemporary applications require methods that
are even faster in order to handle ill-conditioned high dimensional modern day
datasets. In this paper, we propose a new method, G