Comment by Schiphol

The usual matrix notation is a_ij for the element in the ith row and the jth column. If you order axes in this way you get the numpy convention.