Code Development

Build & Test

# Build
make               # -> build/dqmc
make clean         # remove build/

# Test
cd test && make && ./test_greens

Smoke Test

After any code change:

dqmc-util gen Nx=4 Ny=4 U=4 dt=0.1 L=20 n_sweep_warm=50 n_sweep_meas=100
build/dqmc sim_0.h5
dqmc-util summary sim_0.h5

Adding New Measurements

1. Add to X-macro list in

   src/rc/sim_types.h

   :

   // Add to MEAS_EQLT_LIST (or MEAS_UNEQLT_LIST for unequal-time)
   #define MEAS_EQLT_LIST \
       /* ... */ \
       X(my_meas, num_ij)  // pick an appropriate size expression
   

2. Implement measurement kernel in

   src/rc/meas.c

3. Add analysis function in

   dqmc_util/analyze_hub.py

   (and update reshaping/symmetrization if needed):

   @observable(
       description="my new observable",
       requires=("meas_eqlt/n_sample", "meas_eqlt/sign", "meas_eqlt/my_meas"),
   )
   def _my_meas(ctx):
       return core.jackknife_noniid(
           ctx.data["meas_eqlt/n_sample"],
           ctx.data["meas_eqlt/sign"],
           ctx.data["meas_eqlt/my_meas"],
       )
   

Key Code Patterns

Dual compilation: Files in

src/rc/

compile twice. Use

RC(name)

macro for real/complex dispatch.

Memory: Use

my_calloc()

for 64-byte aligned allocations.

Profiling: Wrap with

profile_begin(name)

/

profile_end(name)

.

Definition of Done

• make

  succeeds without warnings

• cd test && make && ./test_greens

  passes
• Smoke run completes
• If touching measurements: verify via

  dqmc-util summary