Variational inequality (VI) is a fundamental mathematical framework for many classical problems. We present a path-following framework for finite-dimensional VIs with arbitrary continuous functions and compact convex domains. The approach first approximately reduces a general VI to a smooth VI on simplex. Its key innovation is to formulate the smooth VI on simplex on a fiber bundle called the fixed-point bundle. Exploiting this geometric structure, we systematically integrate starting point selection, path-following, and singularity avoidance. Without any assumptions such as monotonicity, the algorithm guarantees global linear convergence to nonsingular solutions. For singular solutions, it retains global linear reduction up to a fixed precision, after which convergence becomes sublinear as the required precision increases. In numerical experiments on 14400 randomly generated VIs of up to 800 dimensions, the algorithm succeeds in every instance, and iteration number increases only mildly with the dimension.