Nice, good to have an expert

BTW: I switch to using ``numpy.testing.assert_allclose`` since we require numpy 1.5