|
|
@ -17,8 +17,6 @@ |
|
|
|
from scipy import stats |
|
|
|
from scipy import stats |
|
|
|
import math |
|
|
|
import math |
|
|
|
|
|
|
|
|
|
|
|
_THRESHOLD = 1e-10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def scale(a, mul): |
|
|
|
def scale(a, mul): |
|
|
|
return [x * mul for x in a] |
|
|
|
return [x * mul for x in a] |
|
|
@ -28,18 +26,18 @@ def cmp(a, b): |
|
|
|
return stats.ttest_ind(a, b) |
|
|
|
return stats.ttest_ind(a, b) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def speedup(new, old): |
|
|
|
def speedup(new, old, threshold): |
|
|
|
if (len(set(new))) == 1 and new == old: return 0 |
|
|
|
if (len(set(new))) == 1 and new == old: return 0 |
|
|
|
s0, p0 = cmp(new, old) |
|
|
|
s0, p0 = cmp(new, old) |
|
|
|
if math.isnan(p0): return 0 |
|
|
|
if math.isnan(p0): return 0 |
|
|
|
if s0 == 0: return 0 |
|
|
|
if s0 == 0: return 0 |
|
|
|
if p0 > _THRESHOLD: return 0 |
|
|
|
if p0 > threshold: return 0 |
|
|
|
if s0 < 0: |
|
|
|
if s0 < 0: |
|
|
|
pct = 1 |
|
|
|
pct = 1 |
|
|
|
while pct < 101: |
|
|
|
while pct < 101: |
|
|
|
sp, pp = cmp(new, scale(old, 1 - pct / 100.0)) |
|
|
|
sp, pp = cmp(new, scale(old, 1 - pct / 100.0)) |
|
|
|
if sp > 0: break |
|
|
|
if sp > 0: break |
|
|
|
if pp > _THRESHOLD: break |
|
|
|
if pp > threshold: break |
|
|
|
pct += 1 |
|
|
|
pct += 1 |
|
|
|
return -(pct - 1) |
|
|
|
return -(pct - 1) |
|
|
|
else: |
|
|
|
else: |
|
|
@ -47,7 +45,7 @@ def speedup(new, old): |
|
|
|
while pct < 100000: |
|
|
|
while pct < 100000: |
|
|
|
sp, pp = cmp(new, scale(old, 1 + pct / 100.0)) |
|
|
|
sp, pp = cmp(new, scale(old, 1 + pct / 100.0)) |
|
|
|
if sp < 0: break |
|
|
|
if sp < 0: break |
|
|
|
if pp > _THRESHOLD: break |
|
|
|
if pp > threshold: break |
|
|
|
pct += 1 |
|
|
|
pct += 1 |
|
|
|
return pct - 1 |
|
|
|
return pct - 1 |
|
|
|
|
|
|
|
|
|
|
|