1. Download it from milc_download.
2. Goto hvy_qpot , copy Makefile from main directory using cd hvy_qpot/ && cp ../Makefile .
3. Edit the Makefile (choose double precison and comment out WANTQIO = true using WANTQIO= # true)
4. Make the executable su3_hqp_hyp using make su3_hqp_hyp
5. Remember that the parameter files are in /test/ directory of each applications of MILC. This directory also contains the sample input and sample output, checked them recently and they agrees upon over all observables.
6. The binary files are in milc_qcd-7.8.1/binary_samples
7. Note that the each executable applications has control.c file which contains the main code of applications.

Gauge configurations:
/root/project1/build_SIMULATeQCD/work/rishabh_test/l328f21b6315m00281m0759Nb1_s1_v/nersc.1000
Header Info:
BEGIN_HEADER
DATATYPE = 4D_SU3_GAUGE
DIMENSION_1 = 32
DIMENSION_2 = 32
DIMENSION_3 = 32
DIMENSION_4 = 8
CHECKSUM = 6a3f7dcd
LINK_TRACE = -0.0003083589643
PLAQUETTE = 0.5369923862
FLOATING_POINT = IEEE32BIG
END_HEADER

Pure gauge SU3 Wilson loops
MIMD version 6
Machine = Scalar processor, with 1 nodes
start: Mon Apr 7 09:47:52 2025

type 0 for no prompts, 1 for prompts, or 2 for proofreading
nx 32
ny 32
nz 32
nt 8
LAYOUT = Hypercubes, options = with fixed node_geometry
ON EACH NODE 32 x 32 x 32 x 8
Mallocing 306.2 MBytes per node for lattice

no_smear_level 0
smear_num[0] 0
smear_num[1] 0
smear_num[2] 0
smear_num[3] 0
smear_num[4] 0
hyp_alpha2 0.6
hyp_alpha3 0.3
off_axis_flag 0
max_t 6
reload_serial ../../binary_samples/rishabh_test/l328f21b6315m00281m0759Nb1_s1_v/nersc.1000
forget
reading as archive format with byte reversal
reading Archive header:
DATATYPE = 4D_SU3_GAUGE
DIMENSION_1 = 32
DIMENSION_2 = 32
DIMENSION_3 = 32
DIMENSION_4 = 8
CHECKSUM = 6a3f7dcd
LINK_TRACE = -0.0003083589643
PLAQUETTE = 0.5369923862
FLOATING_POINT = IEEE32BIG
END_HEADER
Restored archive gauge configuration serially from file ../../binary_samples/rishabh_test/l328f21b6315m00281m0759Nb1_s1_v/nersc.1000
Archive style checksum = 6a3f7dcd OK
Time to reload gauge configuration = 1.332240e-01
CHECK PLAQ: 1.6103238379248912e+00 1.6116304769281531e+00
CHECK NERSC LINKTR: -3.0835895704682276e-04 CKSUM: 6a3f7dcd
Reunitarized for double precision. Max deviation 1.56e-07 changed to 5.55e-16
Time to check unitarity = 1.323130e-01
RUNNING COMPLETED
Time = 2.311590e-01 seconds

get_i(0): EOF on input.
exit: Mon Apr 7 09:47:53 2025

[2025-04-07 09:53:47] DEBUG: [Rank 0] MPI File [00] ( 32 32 32 8 ) ( 32 32 32 8 ) ( 0 0 0 0 )
[2025-04-07 09:53:47] INFO: Polyakov Loop (Real) =0.00101486
[2025-04-07 09:53:47] INFO: plaquette=0.536992
[2025-04-07 09:53:47] INFO: rectangle=0.277423
[2025-04-07 09:53:47] INFO: Polyakov Loop =(0.00101486,0.000356387)
Node:0,Smear_count:0,Length=1,Height=1,Wilson Loop:5.369923861926766e-01
Node:0,Smear_count:0,Length=1,Height=2,Wilson Loop:2.705562740443333e-01
Node:0,Smear_count:0,Length=2,Height=1,Wilson Loop:2.842903503747579e-01
Node:0,Smear_count:0,Length=2,Height=2,Wilson Loop:8.919544525263672e-02
[2025-04-07 09:53:47] INFO: Finalizing MPI

[2025-04-07 10:03:57] DEBUG: [Rank 0] MPI File [00] ( 32 32 32 8 ) ( 32 32 32 8 ) ( 0 0 0 0 )
[2025-04-07 10:03:58] INFO: Polyakov Loop (Real) =0.00101486
[2025-04-07 10:03:58] INFO: plaquette=0.536992
[2025-04-07 10:03:58] INFO: rectangle=0.277423
[2025-04-07 10:03:58] INFO: Polyakov Loop =(0.00101486,0.000356387)
Node:0,Smear_count:0,Length=1,Height=1,Wilson Loop:5.367746130250461e-01
Node:0,Smear_count:0,Length=1,Height=2,Wilson Loop:2.841646226670829e-01
Node:0,Smear_count:0,Length=2,Height=1,Wilson Loop:2.843749465270896e-01
Node:0,Smear_count:0,Length=2,Height=2,Wilson Loop:9.350367800802106e-02
[2025-04-07 10:03:58] INFO: Finalizing MPI

Observable	MILC	SIMULATeQCD
spatial plaquette	1.6103238379248912e+00 /3 = 0.53677461264163	0.5367746130250461
space-time plaquette	1.6116304769281531e+00 /3 = 0.53721015897605	0.5369923861926766

We can also try the same with many such gauge configurations for zero temperature;

Beta=6.354, configurations are synced from bielefeld cluster(/work/DATA/petreczk).
These are the values for 400 measurements for spatial plaquette by MILC : plaquette_ss_milc
These are the values for 400 measurements for space time plaquette by MILC : plaquette_st_milc
These are the values for 400 measurements for spatial plaquette by SIMULATeQCD : plaquette_ss_simqcd
These are the values for 400 measurements for space time plaquette by SIMULATeQCD : plaquette_st_simqcd

Comparing first 10 of these for spatial plaquette:

conf no.	MILC results	SIMULATeQCD
1	0.541731084762594	0.541731082059970
2	0.541756428886104	0.541756426158047
3	0.541847370086132	0.541847367337050
4	0.541573897412149	0.541573894781463
5	0.541936493962628	0.541936491284987
6	0.541708205165146	0.541708202516068
7	0.541393033181526	0.541393030479058
8	0.541404096591367	0.541404093873545
9	0.541676734249788	0.541676731502657
10	0.541651900195796	0.541651897543551

Comparing first 10 of these for space time plaquette (there are only 3 space time plaquette in 3 space and one time case):

conf no.	MILC results	SIMULATeQCD
1	0.541531052446446	0.541531049746394
2	0.541662986685118	0.541662983967113
3	0.541826254419194	0.541826251746843
4	0.541713703185204	0.541713700614260
5	0.541759608436956	0.541759605770859
6	0.541780611686284	0.541780609010254
7	0.541439430890162	0.541439428111445
8	0.541614248623674	0.541614245881417
9	0.541758876290544	0.541758873584825
10	0.541515946280264	0.541515943552465

The space space plaquette comparision plot:
The space time plaquette comparision plot: