Include MATPOWER hessian function for Lagrangian. This should be possible by setting the lambda and entries corresponding to the copy nodes to zero.