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.