Fix sign error in log density ratio for Langevin proposal (#41)

UCL · Oct 23, 2024 · 4e6a527 · 4e6a527
1 parent 697bc44
commit 4e6a527
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/langevin.R b/R/langevin.R
@@ -44,7 +44,7 @@ log_density_ratio_langevin <- function(
     proposed_state,
     target_distribution,
     scale_and_shape) {
-  sum(proposed_state$momentum()^2 - state$momentum()^2) / 2
+  sum(state$momentum()^2 - proposed_state$momentum()^2) / 2
 }
 
 #' Create a new Langevin proposal object.