Solving Constrained Variational Inequalities via an Interior Point Method