Took me a moment to figure out what 'atol' and 'rtol' are: they're not explained until the np.allclose call below.