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)
```
This is a bit deceptive with `_NotGiven`. I think the default needs better explanation. Also in the docstrings of all the other functions changed by this PR.