def METHOD_NAME(self):
    dp = u.Quantity([[[1, 0, 0]]], unit=u.N / u.m**3)
    B = u.Quantity([[[0, -1, 0]]], unit=u.T)
    result = drifts.diamagnetic_drift(2 * dp, 3 * B, self.n, self.q)
    assert_quantity_allclose(
        result, (2 / 3) * u.Quantity([[[0, 0, 1]]], unit=u.m / u.s)
    )