コード例 #1
0
        """ 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")
コード例 #2
0
ファイル: main.py プロジェクト: acs1903/genthor
        
        # 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")