I found the construction of high-dimension array very problematic. It is not friendly if you want to create a high-dimensional (dim>3) array by hand. Also, the column major will lead to extra complications here, for example:
julia> a=reshape([1,2,3,4,5,6,7,8],(2,2,2))
2×2×2 Array{Int64,3}:
[:, :, 1] =
1 3
2 4
[:, :, 2] =
5 7
6 8
which is a little conter-intuitive as we expect
2×2×2 Array{Int64,3}:
[:, :, 1] =
1 2
3 4
[:, :, 2] =
5 6
7 8
Also, can’t we create a high-dimension array from nested list?
In conclusion, my questions are how to construct high-dimension matrices and how to understand the column major problem in high-dimension arrays.