a number of smaller changes

claudiopica · Sep 29, 2024 · 61a0f2a · 61a0f2a
1 parent 3160273
commit 61a0f2a
Show file tree

Hide file tree

Showing 8 changed files with 59 additions and 334 deletions.
diff --git a/Include/Utils/wilsonflow.h b/Include/Utils/wilsonflow.h
@@ -35,12 +35,6 @@ void WilsonFlow1(suNg_field *V, const double epsilon);
 void WilsonFlow3(suNg_field *V, const double epsilon);
 int WilsonFlow3_adaptative(suNg_field *V, double *epsilon, double *epsilon_new, double *delta);
 
-double WF_E(suNg_field *V);
-double WF_Esym(suNg_field *V);
-double WF_topo(suNg_field *V);
-
-void WF_E_T(double *E, suNg_field *V);
-void WF_Esym_T(double *Esym, suNg_field *V);
 data_storage_array *WF_update_and_measure(WF_integrator_type wft, suNg_field *V, double *tmax, double *eps, double *delta,
                                           int nmeas, storage_switch swc);
 

diff --git a/LibHR/Observables/loop_tools.c b/LibHR/Observables/loop_tools.c
@@ -14,8 +14,7 @@
 #include "io.h"
 #include "utils.h"
 #include "memory.h"
-#include "Update/representation.h"
-#include "Update/avr_plaquette.h"
+#include "update.h"
 #include "inverters.h"
 
 #if (!defined(BC_T_SF_ROTATED) && !defined(BC_T_SF) && !defined(FERMION_THETA))

diff --git a/LibHR/Update/mon_pg.c b/LibHR/Update/mon_pg.c
@@ -5,7 +5,6 @@
 
 #include "update.h"
 #include "libhr_core.h"
-#include "Update/avr_plaquette.h"
 #include "memory.h"
 #include "inverters.h"
 
@@ -30,7 +29,7 @@ static const spinor_field *pg_pseudofermion(monomial const *m) {
 static void pg_add_local_action(monomial const *m, scalar_field *loc_action) {
     mon_pg_par *par = (mon_pg_par *)(m->data.par);
     if (local_plaq_tmp == NULL) { local_plaq_tmp = alloc(local_plaq_tmp, 1, &glattice); }
-    local_plaquette(local_plaq_tmp);
+    local_plaquette(u_gauge, local_plaq_tmp);
     mul_add_assign(loc_action, -(par->beta / ((double)NG)), local_plaq_tmp);
 }
 

diff --git a/LibHR/Utils/background_field.c b/LibHR/Utils/background_field.c
@@ -82,7 +82,7 @@ void apply_background_field_zdir(suNg_field *V, double Q, int n) {
                     index = ipt(c[0], c[1], c[2], c[3]);
 
                     // calculate new plaquettes in direction 3,4
-                    cplaq(&tmp, index, 3, 0);
+                    cplaq(&tmp, u_gauge, index, 3, 0);
 
                     // now divide by e^{iQE}
                     phase = cos(Q * E) + I * sin(Q * E);
@@ -94,7 +94,7 @@ void apply_background_field_zdir(suNg_field *V, double Q, int n) {
                     start_sendrecv_suNg_field(V);
                     complete_sendrecv_suNg_field(V);
 
-                    cplaq(&tmp, index, 3, 0);
+                    cplaq(&tmp, u_gauge, index, 3, 0);
 
                     // compute absolute value of the difference between old and new for the real and imaginary parts independantly
                     diff_re += fabs(creal(tmp2 - tmp));