We present an algorithm for variational inequalities VI(F, Y) that uses a p
rimal-dual version of the Analytic Center Cutting Plane Method. The point-t
o-set mapping F is assumed to be monotone, or pseudomonotone. Each computat
ion of a new analytic center requires at most four Newton iterations, in th
eory, and in practice one or sometimes two. Linear equalities that may be i
ncluded in the definition of the set Y are taken explicitly into account.
We report numerical experiments on several well-known variational inequalit
y problems as well as on one where the functional results from the solution
of large subproblems. The method is robust and competitive with algorithms
which use the same information as this one.