'Regression test for https://github.com/astropy/astropy/issues/2003'

def test_slice_to_quantity(self):

    a = np.random.uniform(size=(10, 8))
    (x, y, z) = ((a[:, 1:4].T * u.km) / u.s)
    total = np.sum((((a[:, 1] * u.km) / u.s) - x))
    assert isinstance(total, u.Quantity)
    assert (total == ((0.0 * u.km) / u.s))
