33 void applyU(GaugeField &force, GaugeField &U);
double computeMomAction(const GaugeField &mom)
Compute and return global the momentum action 1/2 mom^2.
void applyU(GaugeField &force, GaugeField &U)
void updateMomentum(GaugeField &mom, double coeff, GaugeField &force)