import numpy as np
import matplotlib.pyplot as plt
import sys
import re
def fhex(s):
return float.fromhex(s)
def parse(f):
xs = []
gs = []
ys = []
es = []
r = re.compile(r"[^ (]+\(([^ )]*)\) got ([^ ]+) want ([^ ]+) [^ ]+ ulp err ([^ ]+)")
for line in f:
m = r.match(line)
if m:
x = fhex(m.group(1))
g = fhex(m.group(2))
y = fhex(m.group(3))
e = float(m.group(4))
xs.append(x)
gs.append(g)
ys.append(y)
es.append(e)
elif line.startswith("PASS") or line.startswith("FAIL"):
print(line)
return xs, gs, ys, es
def plot(xs, gs, ys, es):
if len(xs) < 2:
print("not enough samples")
return
a = min(xs)
b = max(xs)
fig, (ax0, ax1) = plt.subplots(nrows=2)
es = np.abs(es) emax = max(es)
ax0.text(a + (b - a) * 0.7, emax * 0.8, "%s\n%g" % (emax.hex(), emax))
ax0.plot(xs, es, "r.")
ax0.grid()
ax1.plot(xs, ys, "r.", label="want")
ax1.plot(xs, gs, "b.", label="got")
ax1.grid()
ax1.legend()
plt.show()
xs, gs, ys, es = parse(sys.stdin)
plot(xs, gs, ys, es)