From cd2b5e58b6b5965bdefc6cdc2305fd32aa8e4b3e Mon Sep 17 00:00:00 2001
From: Joshua S <joshuatsdev@gmail.com>
Date: Thu, 12 Dec 2024 14:07:36 -0500
Subject: [PATCH] Moved V_proj to H_full timeslice

---
 applications/false-vacuum-decay/hmc.py | 2 +-
 qlat/qlat/include/qlat/qm-action.h     | 8 ++++----
 2 files changed, 5 insertions(+), 5 deletions(-)

diff --git a/applications/false-vacuum-decay/hmc.py b/applications/false-vacuum-decay/hmc.py
index ab8ac2eb..69229be2 100755
--- a/applications/false-vacuum-decay/hmc.py
+++ b/applications/false-vacuum-decay/hmc.py
@@ -458,7 +458,7 @@ def main():
     # The number of trajectories to calculate
     n_traj = 50000
     #
-    version = "8-1"
+    version = "8-2"
     date = datetime.datetime.now().date()
     # The number of steps to take in a single trajectory
     steps = 10
diff --git a/qlat/qlat/include/qlat/qm-action.h b/qlat/qlat/include/qlat/qm-action.h
index 8146be5e..99e0dc26 100644
--- a/qlat/qlat/include/qlat/qm-action.h
+++ b/qlat/qlat/include/qlat/qm-action.h
@@ -84,7 +84,7 @@ struct QMAction {
     else if(t<t_TV_start+t_full1)
         return V_full(x);
     else if(t==t_TV_start+t_full1)
-      return V_proj(x);
+      return V_full(x) + V_proj(x);
     else if(t<t_TV_start+t_full1+t_FV_out)
       return V_FV_out(x);
     else if(t<t_TV_start+t_full1+t_FV_out+t_FV_mid)
@@ -92,7 +92,7 @@ struct QMAction {
     else if(t<t_TV_start+t_full1+2*t_FV_out+t_FV_mid-1)
       return V_FV_out(x);
     else if(t==t_TV_start+t_full1+2*t_FV_out+t_FV_mid-1)
-      return V_proj(x);
+      return V_full(x) + V_proj(x);
     else if(t<t_TV_start+t_full1+2*t_FV_out+t_FV_mid+t_full2)
       return V_full(x);
     else {
@@ -179,13 +179,13 @@ struct QMAction {
   
   inline double V_proj(const double x)
   {
-    return -P*log(1-exp(-(V_FV_out(x) - V_full(x) + epsilon)*dt));
+    return -P*log(1-exp(-(V_FV_out(x) - V_full(x) + epsilon + epsilon*x)*dt)) / dt;
   }
   
   inline double dV_proj(const double x)
   {
     double Vbar = V_FV_out(x) - V_full(x);
-    return -P*((dV_FV_out(x) - dV_full(x))*dt*exp(-(Vbar + epsilon)*dt))/(1-exp(-(Vbar + epsilon)*dt));
+    return -P*((dV_FV_out(x) - dV_full(x) + epsilon)*exp(-(Vbar + epsilon + epsilon*x)*dt))/(1-exp(-(Vbar + epsilon + epsilon*x)*dt));
   }
 
   inline double V_FV_mid(const double x)