Take a look at D language library Mir for N-Dimensional array and its comparison to Numpy [1].
Its native performance inside GLAS implementation is something to shout about and even comparable to OpenBLAS that is currently being used by Julia and Matlab [2][3].
[1]http://docs.algorithm.dlang.io/latest/mir_ndslice.html
[2]http://blog.mir.dlang.io/glas/benchmark/openblas/2016/09/23/...