From ab0b1f3cefedb35327a19ec1b6afd560bfdf802d Mon Sep 17 00:00:00 2001 From: Thibaut Horel Date: Mon, 14 Sep 2015 23:08:02 -0400 Subject: Import supplements and repo reorganization --- experiments/plot3d.py | 23 ----------------------- 1 file changed, 23 deletions(-) delete mode 100644 experiments/plot3d.py (limited to 'experiments/plot3d.py') diff --git a/experiments/plot3d.py b/experiments/plot3d.py deleted file mode 100644 index b1b337f..0000000 --- a/experiments/plot3d.py +++ /dev/null @@ -1,23 +0,0 @@ -from mpl_toolkits.mplot3d import Axes3D -from matplotlib import cm -import matplotlib.pyplot as plt -import numpy as np - -with open("out.log") as fh: - values = [map(float, line.strip().split()) for line in fh] - #values = [(b, a, l) for (b, a, l) in values if b >= 0.04] - am = max(values, key=lambda x: x[4]) - # am[0] = 1./am[0] - print am - alpha, delta, beta, r , l = zip(*values) - # alpha = 1./np.array(alpha) - - fig = plt.figure(figsize=(12, 8)) - ax = fig.gca(projection='3d') - - ax.plot_trisurf(alpha, delta, l, cmap=cm.jet, linewidth=0.001) - plt.xlabel("alpha") - plt.ylabel("delta") - ax.set_zlabel('Likelihood') - #plt.savefig("ll.pdf") - plt.show() -- cgit v1.2.3-70-g09d2