The CreateSortedStridePerm does not work correctly when it comes to non C ordered arrays that include a 1 dimensional axis somewhere in the middle (and are not contiguous as this is typically optimized anyways). Its a corner case, but since I noticed thought I would post it.

For these cases it will fail to sort the array and this results unnecessarily slow operations for some things. (IE: e = np.empty(400,400,400)[::2,::2,::2]; e[:,0,:] = 5) is much faster then e[:,0:1,:] = 5. Also this means that e.copy('k') is not Fortran order, which maybe it never was, but seems not what one would expect.