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)/t < 0.1, 'simulation %s = %g != %g = %s' % (n,s,t,n))
+ assert (s - t)/w < 2, (
+ 'simulation %s = %g != %g = %s (bin width = %g)' % (n,s,t,n,w))
return sim.residual(theory)