""" Inputs proposed state, returns fwd/bak probabilities.""" # TODO: Not currently implemented... fwdprob = bakprob = 0 return fwdprob, bakprob ## if __name__ == "__main__": import cPickle as pickle import matplotlib.pyplot as plt import matplotlib.cm as cm import time LightBase.destroy_windows() def plot_scores(sampler): scores = sampler.store_info["score"] proposal_scores = sampler.store_info["proposal_score"] accepted = sampler.store_info["accepted"] accepted_idx = np.flatnonzero(accepted) N = scores.size # Plot plt.figure(20) plt.clf() plt.subplot(2, 1, 1) plt.plot(np.vstack((proposal_scores, scores)).T) plt.plot(accepted_idx, scores[accepted_idx], "go") plt.subplot(2, 1, 2) plt.plot(np.zeros(N), "k")
# TODO: Not currently implemented... fwdprob = bakprob = 0 return fwdprob, bakprob ## if __name__ == "__main__": import cPickle as pickle import matplotlib.pyplot as plt import matplotlib.cm as cm import time LightBase.destroy_windows() def plot_scores(sampler): scores = sampler.store_info["score"] proposal_scores = sampler.store_info["proposal_score"] accepted = sampler.store_info["accepted"] accepted_idx = np.flatnonzero(accepted) N = scores.size # Plot plt.figure(20) plt.clf() plt.subplot(2, 1, 1) plt.plot(np.vstack((proposal_scores, scores)).T) plt.plot(accepted_idx, scores[accepted_idx], "go") plt.subplot(2, 1, 2) plt.plot(np.zeros(N), "k")