void applyU(GaugeField &force, GaugeField &U)
double computeMomAction(const GaugeField &mom)
Compute and return global the momentum action 1/2 mom^2.
bool forceMonitor()
Whether we are monitoring the force or not.
void flushForceMonitor()
Flush any outstanding force monitoring information.
void updateMomentum(GaugeField &mom, double coeff, GaugeField &force, const char *fname)