| # SPDX-License-Identifier: GPL-2.0 |
| # |
| # Copyright (C) 2019, VMware Inc, Tzvetomir Stoyanov <[email protected]> |
| # Copyright (C) 2019, VMware Inc, Yordan Karadzhov <[email protected]> |
| |
| |
| import matplotlib.pyplot as plt |
| import matplotlib.lines as mlines |
| import numpy as np |
| import sys |
| |
| def newline(p1, p2): |
| ax = plt.gca() |
| xmin, xmax = ax.get_xbound() |
| |
| if(p2[0] == p1[0]): |
| xmin = xmax = p1[0] |
| ymin, ymax = ax.get_ybound() |
| else: |
| ymax = p1[1]+(p2[1]-p1[1])/(p2[0]-p1[0])*(xmax-p1[0]) |
| ymin = p1[1]+(p2[1]-p1[1])/(p2[0]-p1[0])*(xmin-p1[0]) |
| |
| l = mlines.Line2D([xmin,xmax], [ymin,ymax], color='red') |
| ax.add_line(l) |
| return l |
| |
| data = np.loadtxt(fname = sys.argv[1]) |
| x = data[:, 0] |
| y = data[:, 1] |
| |
| fig, ax = plt.subplots() |
| |
| ax.set_xlabel('samples (t)') |
| ax.set_ylabel('clock offset') |
| ax.set_title("$\delta$=%i ns" % (max(y) - min(y))) |
| |
| l = mlines.Line2D(x, y) |
| ax.add_line(l) |
| ax.set_xlim(min(x), max(x)) |
| ax.set_ylim(min(y), max(y) ) |
| |
| print(min(y), max(y), max(y) - min(y)) |
| |
| # Tweak spacing to prevent clipping of ylabel |
| fig.tight_layout() |
| plt.show() |
