From e2b8d999c072ce16610f63448cc2bee14bd55888 Mon Sep 17 00:00:00 2001 From: Jacob <902786+higj@users.noreply.github.com> Date: Thu, 14 Mar 2024 10:23:31 +0200 Subject: [PATCH] Account for the fact that beta is actually beta_np in the bosonic primitive kinetic energy estimator --- src/REPLICA/bosonic_exchange.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/REPLICA/bosonic_exchange.cpp b/src/REPLICA/bosonic_exchange.cpp index 44aae36d65..794dfc6520 100644 --- 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; } /* ---------------------------------------------------------------------- */