BUG : np.sum silently drops keepdims for sub-classes of ndarray

change test from `type(a) is not mu.ndarray` to `not isinstance(a, mu.ndarray)`

This bug can be demonstrated by:

```
class test_array(ndarray):
    pass

base_ = np.arange(30).reshape(5, 6)
tt = base_.view(test_array)

def sum_test(a):
    print(np.sum(a))
    print(np.sum(a, keepdims=True))

sum_test(base_)
sum_test(tt)
```
I don't quite see how this is supposed to work, is the default value still `_NotGiven`, so we are passing a class object around by default, or should I make an object for each function and change the test to `isinstance(keepdims, _NotGiven)`?