e = ExponentialModelFitter(sim)
params = e.fit()
sim_tau = abs(params[0])
- assert (sim_tau - tau)/tau, 'simulation tau = %g != %g = tau' % (
- sim_tau, tau)
+ for s,t,n in [(sim_tau, tau, 'tau')]:
+ assert (s - t)/t < 0.1, 'simulation %s = %g != %g = %s' % (n,s,t,n)
return sim.residual(theory)