Ofdft

I had to implement array_wrap in BaseField, otherwise the attributes that are not native to numpy.ndarray (such as rank) are screwed up by any function (such as np.sqrt) acting on a field.