install.packages('splitFeas')
wood_inv_solve
computes the inverse of the Hessian term of the majorization of the proximity function
using the Woodbury formula. The function mmqn_step
invokes wood_inv_solve
instead of ddg if the argument
woodbury=TRUE
. This should be used when p << n.