A coordinate-free theory of virtual holonomic constraints