From 4864635c327c2fd58bcddabd1668e8126cb71293 Mon Sep 17 00:00:00 2001 From: William Lamb <60004701+astrolamb@users.noreply.github.com> Date: Tue, 21 Oct 2025 13:55:06 -0500 Subject: [PATCH 1/2] attempting to fix swapping bug --- QuickCW/QuickMCMCUtils.py | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/QuickCW/QuickMCMCUtils.py b/QuickCW/QuickMCMCUtils.py index 367a4f5..c13cdd0 100644 --- a/QuickCW/QuickMCMCUtils.py +++ b/QuickCW/QuickMCMCUtils.py @@ -175,13 +175,17 @@ def do_pt_swap(n_chain, samples, itrb, Ts, a_yes, a_no, x0s, FLIs, log_likelihoo FLIs_new = [] x0s_new = [] fisher_diag_new = np.zeros_like(fisher_diag) + samples_new = np.zeros_like(samples) + log_likelihood_new = np.zeros_like(log_likelihood) for j in range(n_chain): - samples[j,itrb+1,:] = samples[swap_map[j],itrb,:] + samples_new[j,itrb+1,:] = samples[swap_map[j],itrb,:] fisher_diag_new[j,:] = fisher_diag[swap_map[j],:] - log_likelihood[j,itrb+1] = log_likelihood[swap_map[j],itrb] + log_likelihood_new[j,itrb+1] = log_likelihood[swap_map[j],itrb] FLIs_new.append(FLIs[swap_map[j]]) x0s_new.append(x0s[swap_map[j]]) + samples_new = samples + log_likelihood = log_likelihood_new fisher_diag[:] = fisher_diag_new FLIs[:] = List(FLIs_new) x0s[:] = List(x0s_new) From 69a84ab3c1a5ad63f9f1269894d11ddb95734df7 Mon Sep 17 00:00:00 2001 From: William Lamb <60004701+astrolamb@users.noreply.github.com> Date: Tue, 21 Oct 2025 13:59:21 -0500 Subject: [PATCH 2/2] minor bug fix --- QuickCW/QuickMCMCUtils.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/QuickCW/QuickMCMCUtils.py b/QuickCW/QuickMCMCUtils.py index c13cdd0..9dc59aa 100644 --- a/QuickCW/QuickMCMCUtils.py +++ b/QuickCW/QuickMCMCUtils.py @@ -184,7 +184,7 @@ def do_pt_swap(n_chain, samples, itrb, Ts, a_yes, a_no, x0s, FLIs, log_likelihoo FLIs_new.append(FLIs[swap_map[j]]) x0s_new.append(x0s[swap_map[j]]) - samples_new = samples + samples = samples_new log_likelihood = log_likelihood_new fisher_diag[:] = fisher_diag_new FLIs[:] = List(FLIs_new)