Modified the Lagrangian and primal infeasibility gradients to use a scaled matrix-vector product.
Hopefully, this improves the algorithm when the matrix-vector products are computed inexactly.