powered by
opt.covmat(xtx, xtr, bound, Lagrangian) tib.covmat(xtx, gen.inverse.diag, beta, Lagrangian)