Add Hoelbling operator to docs

.github/workflows/Tests.yml

@@ -10,7 +10,7 @@ on:
 
 jobs:
   test:
-    runs-on: [ubuntu-latest, windows-latest]
+    runs-on: ubuntu-latest
     timeout-minutes: 20
     strategy:
       fail-fast: false

docs/src/dirac.md

@@ -24,3 +24,7 @@ StaggeredDiracOperator
 ```@docs
 StaggeredEOPreDiracOperator
 ```
+
+```@docs
+StaggeredHoelblingDiracOperator
+```

src/diracoperators/staggered_hoelbling.jl

@@ -1,14 +1,16 @@
 """
-    StaggeredHoelblingDiracOperator(f::AbstractField, mass; bc_str="antiperiodic")
+    StaggeredHoelblingDiracOperator{MT}(f::AbstractField, mass; bc_str="antiperiodic")
     StaggeredHoelblingDiracOperator(D::StaggeredHoelblingDiracOperator, U::Gaugefield)
 
-Create a free Hölbling mass split Staggered Dirac Operator with mass `mass`.
+Create a free Hölbling mass split Staggered Dirac Operator (arXiv:1009.5362) with mass
+`mass`.
+The type-parameter `MT` determines the kind of operator that is used:
+- `MT = 1`: M12 + M34
+- `MT = 2`: M13 + M24
 
 `bc_str` can either be `"periodic"` or `"antiperiodic"` and specifies the boundary
 condition in the time direction.
 
-If `csw ≠ 0`, a clover term is included. 
-
 This object cannot be directly applied to a fermion vector, since it lacks a gauge
 background.
 A Wilson Dirac operator with gauge background is created by applying it to a `Gaugefield`

src/main/runbuild.jl

@@ -142,7 +142,8 @@ function metabuild!(
                     therm=true,
                 )
             end
-            # @level1("|  Elapsed time:\t$(updatetime) [s] @ $(current_time())\n-")
+
+            @level1("|  Elapsed time:\t$(updatetime) [s] @ $(string(current_time()))")
         end
     end
 
@@ -168,7 +169,7 @@ function metabuild!(
                 numaccepts += accepted
             end
 
-            @level1("|  Elapsed time:\t$(updatetime) [s] @ $(current_time())\n")
+            @level1("|  Elapsed time:\t$(updatetime) [s] @ $(string(current_time()))")
             # all procs send their CVs to all other procs and update their copy of the bias
             CVs = mpi_allgather(U.CV::Float64, comm)
             accepteds = mpi_allgather(accepted::Bool, comm)
@@ -185,13 +186,8 @@ function metabuild!(
         end
     end
 
-    # @level1("└\nTotal elapsed time:\t$(convert_seconds(runtime_all))\n@ $(current_time())")
+    print_total_time(runtime_all)
     flush(stdout)
-    # close all the I/O streams
-    close(updatemethod)
-    close(measurements)
-    close(measurements_with_flow)
-    close(bias)
     close(MetaIO.__GlobalLogger[])
     isinteractive() && set_global_logger!(1) # Reset logger if run from REPL
     return nothing

src/main/runsim.jl

@@ -261,6 +261,8 @@ function metaqcd!(
     @level1("└ Total elapsed time:\t$(runtime_therm) [s]\n")
     recalc_CV!(U, bias) # need to recalc cv since it was not updated during therm
 
+    mpi_barrier()
+
     @level1("┌ Production:")
     _, runtime_all = @timed begin
         numaccepts = 0.0