From 12ed742471c479ae37e3431753244fec66f97fc7 Mon Sep 17 00:00:00 2001
From: "gregory.ashton" <gregory.ashton@ligo.org>
Date: Tue, 4 Oct 2016 15:23:01 +0200
Subject: [PATCH] Adds tswap acceptance fractions
---
pyfstat.py | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/pyfstat.py b/pyfstat.py
index 760ea89..14abc36 100755
--- a/pyfstat.py
+++ b/pyfstat.py
@@ -689,6 +689,9 @@ class MCMCSearch(BaseSearchClass):
sampler.run_mcmc(p0, n)
logging.info("Mean acceptance fraction: {0:.3f}"
.format(np.mean(sampler.acceptance_fraction)))
+ if self.ntemps > 1:
+ logging.info("Tswap acceptance fraction: {}"
+ .format(sampler.tswap_acceptance_fraction))
fig, axes = self.plot_walkers(sampler, symbols=self.theta_symbols)
fig.savefig('{}/{}_init_{}_walkers.png'.format(
self.outdir, self.label, j))
@@ -705,6 +708,9 @@ class MCMCSearch(BaseSearchClass):
sampler.run_mcmc(p0, nburn+nprod)
logging.info("Mean acceptance fraction: {0:.3f}"
.format(np.mean(sampler.acceptance_fraction)))
+ if self.ntemps > 1:
+ logging.info("Tswap acceptance fraction: {}"
+ .format(sampler.tswap_acceptance_fraction))
fig, axes = self.plot_walkers(sampler, symbols=self.theta_symbols)
fig.savefig('{}/{}_walkers.png'.format(self.outdir, self.label))
--
GitLab