Fixed bug in Hasenbusch

Author: ctpeterson

src/stagg_pv_hmc/staghmc_spv.nim

@@ -614,7 +614,8 @@ proc generate_pseudoferms(s: Staggered) =
                var hphit = newOneOf(hphi[fld_ind])
                s.Ddag(phi[fld_ind], psi, hmass)
                s.Ddag(hphit, hpsi, masses[0])
-               s.solve(hphi[fld_ind], hphit, -masses[1], spa[fld_ind])
+               #s.solve(hphi[fld_ind], hphit, -masses[1], spa[fld_ind]) Jan 21, 2026 bug
+               s.solve(hphi[fld_ind], hphit, -hmass, spa[fld_ind])
             of "regress": s.Ddag(phi[fld_ind], psi, masses[0])
             else: qexError "Invalid option for pseudofermion preconditioner"
       else:

src/stagg_pv_hmc/staghmc_spv_init_and_io.nim

@@ -2,6 +2,7 @@
 
 import qex # Import qex
 import gauge/gaugeUtils # For handling gauge fields
+import physics/[stagSolve] # CG initialization
 import layout # For initializing lattice layout
 import rng # For handing RNG
 import parseopt # For parsing command line arguments
@@ -12,7 +13,6 @@ import times # Timing for generating seeds
 import os # For checking if files exist
 import staghmc_spv_rng # For certain RNG operations
 import options # For optional IO behavior with fields
-import physics/[stagSolve]
 
 #[ ~~~~ Command line and XML inputs ~~~~ ]#