Follow statistical convention in ndarray order of axes, e.g. storing `num_vec` of `k`-dimensional vectors should be done in ndarray of shape `(num_vec, k)`, and not `(k, num_vec)` as it's done now.

[jacobnzw]