From 4f1e66b04413f922311104886f8901d200cb2fae Mon Sep 17 00:00:00 2001 From: "Henrique B.R" Date: Sat, 4 Jul 2020 03:53:06 +0100 Subject: [PATCH] Fixed HMC SU(N) integrator which was causing fields to leave Lie Algebra manifold for N>2 --- Grid/qcd/action/scalar/ScalarImpl.h | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/Grid/qcd/action/scalar/ScalarImpl.h b/Grid/qcd/action/scalar/ScalarImpl.h index 203e1824..d33cdc79 100644 --- a/Grid/qcd/action/scalar/ScalarImpl.h +++ b/Grid/qcd/action/scalar/ScalarImpl.h @@ -38,7 +38,6 @@ public: static inline void update_field(Field& P, Field& U, double ep) { U += P*ep; - std::cout << "Field updated. Epsilon = " << std::setprecision(10) << ep << std::endl; } static inline RealD FieldSquareNorm(Field& U) { @@ -175,14 +174,13 @@ public: P *= scale; } - static inline Field projectForce(Field& P) {return P;} + static inline Field projectForce(Field& P) {return Ta(P);} static inline void update_field(Field &P, Field &U, double ep) { #ifndef USE_FFT_ACCELERATION double t0=usecond(); U += P*ep; - std::cout << "Field updated. Epsilon = " << std::setprecision(10) << ep << std::endl; double t1=usecond(); double total_time = (t1-t0)/1e6; std::cout << GridLogIntegrator << "Total time for updating field (s) : " << total_time << std::endl;