Geometric Numerical Integration Via Auxiliary Variables