Merge branch 'surface_twist_fluctuation' of github.com:CFT-HY/HILA in…

…to surface_twist_fluctuation

.github/workflows/generate-doxygen-for-pages.yml

@@ -6,6 +6,7 @@ on:
   push:
     branches: 
       - documentation
+      - master
 
   # Allows you to run this workflow manually from the Actions tab
   workflow_dispatch:

.gitignore

@@ -16,4 +16,7 @@ docs/latex/
 docs/html/*
 !docs/html/images/
 latex/
-html/
+html/
+.vs/slnx.sqlite
+*config*
+*run_status

applications/suN_gauge/src/suN_gauge.cpp

@@ -12,12 +12,13 @@
 #include "gauge/stout_smear.h"
 #include "gauge/sun_heatbath.h"
 #include "gauge/sun_overrelax.h"
+#include "tools/checkpoint.h"
+
 
 #include <fftw3.h>
 
 // local includes
 #include "parameters.h"
-#include "checkpoint.h"
 
 /**
  * @brief Helper function to get valid z-coordinate index
@@ -71,12 +72,16 @@ void measure_stuff(const GaugeField<group> &U, const parameters &p) {
 template <typename group>
 void update(GaugeField<group> &U, const parameters &p, bool relax) {
 
-    foralldir(d) {
-        for (Parity par : {EVEN, ODD}) {
+    // go through dirs in random order
 
-            update_parity_dir(U, p, par, d, relax);
-        }
+    for (auto &dp : hila::shuffle_directions_and_parities()) {
+
+        update_parity_dir(U, p, dp.parity, dp.direction, relax);
     }
+
+    //   for (Parity par : {EVEN, ODD}) foralldir(d) {
+    //         update_parity_dir(U, p, par , d, relax);
+    //     }
 }
 
 /**
@@ -210,7 +215,7 @@ int main(int argc, char **argv) {
     hila::timer update_timer("Updates");
     hila::timer measure_timer("Measurements");
 
-    restore_checkpoint(U, start_traj, p);
+    restore_checkpoint(U, p.config_file, p.n_trajectories, start_traj);
 
     // We need random number here
     if (!hila::is_rng_seeded())
@@ -244,9 +249,9 @@ int main(int argc, char **argv) {
         }
 
         if (p.n_save > 0 && (trajectory + 1) % p.n_save == 0) {
-            checkpoint(U, trajectory, p);
+            checkpoint(U, p.config_file, p.n_trajectories, trajectory);
         }
     }
 
     hila::finishrun();
-}
+}

applications/suN_gauge_gf/Makefile

@@ -0,0 +1,61 @@
+
+# Following line(s) are printed with "make help".  Use columns 8 and 30
+#%     make [suN_gauge_gf]        - build suN_gauge_gf program
+#% suN_gauge_gf make options 
+#%     NCOL=<N>              - SU(N) gauge simulation program (default: 3)
+#%     SUN_OVERRELAX_dFJ=1   - use deForcrand-Jahn full overrelax (default: su2 subgroups)
+#%     GFLOWS=<K>            - Gauge action for gradient flow:
+#%                                                   K=0: "WILSON"
+#%                                                   K=1: "Bulk-Prevention"
+#%                                                   K=2: "Luscher-Weisz"
+#%                                                   K=3: "IWASAKI"         (default)
+#%                                                   K=4: "DBW2"
+
+# Give the location of the top level distribution directory wrt. this location.
+# Can be absolute or relative
+HILA_DIR := ../..
+
+# Number of colors
+NCOL := 3
+
+# Gauge action for gradient flow update: 
+# "WILSON": 0, "BP": 1, "LW": 2, "IWASAKI": 3, "DBW2": 4
+GFLOWS := 3
+
+
+ifndef ARCH
+ARCH := vanilla
+endif
+
+# Set default goal and arch
+.DEFAULT_GOAL := suN_gauge_gf
+
+# Read in the main makefile contents, incl. platforms
+include $(HILA_DIR)/libraries/main.mk
+
+#HILAPP_OPTS += -insert-includes
+#HILAPP_OPTS += -comment-pragmas
+HILAPP_OPTS += -check-init
+
+APP_OPTS += -DNDIM=4 -DNCOLOR=${NCOL} -DGFLOWACTION=${GFLOWS}
+
+ifdef SUN_OVERRELAX_dFJ
+APP_OPTS += -DSUN_OVERRELAX_dFJ
+endif
+
+# With multiple targets we want to use "make target", not "make build/target".
+# This is needed to carry the dependencies to build-subdir
+
+suN_gauge_gf: build/su${NCOL}_gauge_gf_fs${GFLOWS}_${ARCH} ; @:
+
+# Now the linking step for each target executable
+build/su${NCOL}_gauge_gf_fs${GFLOWS}_${ARCH}: Makefile build/suN_gauge_gf.o $(HILA_OBJECTS) $(HEADERS)
+	$(LD) -o $@ build/suN_gauge_gf.o $(HILA_OBJECTS) $(LDFLAGS) $(LDLIBS) 
+
+
+
+
+
+
+
+

applications/suN_gauge_gf/parameters

@@ -0,0 +1,16 @@
+# This is a parameter file for SUN_gauge_gf
+
+lattice size           12,12,12,12
+beta                   6.
+number of trajectories 302
+updates in trajectory  1
+overrelax steps        4
+thermalization trajs   50
+gflow freq             100
+gflow max lambda       12.
+gflow lambda step      0.1
+gflow abs. accuracy    1.0e-5
+gflow rel. accuracy    1.0e-5
+random seed            0
+trajs/saved            200
+config name            config