sim_alpha = params[0]
sim_rho = abs(params[1])
for s,t,n in [(sim_alpha, alpha, 'alpha'), (sim_rho, rho, 'rho')]:
- assert (s - t)/w < 3, (
+ assert abs(s - t)/w < 3, (
'simulation %s = %g != %g = %s (bin width = %g)' % (n,s,t,n,w))
return sim.residual(theory)