Account for the fact that beta is actually beta_np in the bosonic primitive kinetic energy estimator

2024-03-14 10:23:31 +02:00
parent ba8d1587d4
commit e2b8d999c0
1 changed files with 2 additions and 2 deletions
--- a/src/REPLICA/bosonic_exchange.cpp
+++ b/src/REPLICA/bosonic_exchange.cpp
@ -377,8 +377,8 @@ double BosonicExchange::prim_estimator()

    prim_est[m] = sig / sig_denom_m;
  }
-  
-  return 0.5 * np * domain->dimension * nbosons / beta + prim_est[nbosons];
+
+  return 0.5 * domain->dimension * nbosons / beta + prim_est[nbosons] / np;
 }

 /* ---------------------------------------------------------------------- */