Why not simply add **kwargs to this function since it is a wrapper, extending the `assert_array_almost_equal` in `numpy.testing`?  (i.e. not adding the decimal keyword explicitly)?

I don't really care either way, but it was already performing a check with a default precision of 6 by default and control would be given to all keyword arguments of the numpy function if passing **kwargs.