Variational Discretizations of Gauge Field Theories Using Group-Equivariant Interpolation

We describe a systematic mathematical approach to the geometric discretization of gauge field theories that is based on dirac and multi-dirac mechanics and geometry, which provide a unified mathematical framework for describing lagrangian and hamiltonian mechanics and field theories, as well as degenerate, interconnected, and nonholonomic systems. variational integrators yield geometric structure-preserving numerical methods that automatically preserve the symplectic form and momentum maps, and exhibit excellent long-time energy stability. the construction of momentum-preserving variational integrators relies on the use of group-equivariant function spaces, and we describe a general construction for functions taking values in symmetric spaces. this is motivated by the geometric discretization of general relativity, which is a second-order covariant gauge field theory on the symmetric space of lorentzian metrics.