> depends on how you interpret the rows and columns of a matrix
I assume you know everything I am writing below, but maybe it's a clearer version of what I was trying to say in my other post.
Computationally, the definition of matrix multiplication AB = C is identical to multiplying A by each column of B, and concatenating the results as columns of C. The definition works both ways as you said, but this specific interpretation is one consequence of that definition.
Ax is a linear combination of the columns of A, with the values of x as coefficients. This is a direct translation of plugging the x values into a system of linear equations where the coefficients of each equation in the system form the rows of A:
a11 * x1 + a12 * x2 = y1
a21 * x1 + a22 * x2 = y2
[ a11 a12 ; a21 a22 ] [ x1 ; x2 ] = [ y1 ; y2 ]
Moreover, in a matrix that corresponds to a system of non-redundant linear equations, the columns correspond to a basis for the image of the linear transformation represented by A. This follows immediately from the definition of a basis.
You can go on like this, but the overall idea is that if you define "a matrix" such that each row contains the coefficients of one equation in a system of linear equations, then you are setting up row vectors to correspond to transformations (as in, y = a1x1 + a2x2), and then it becomes very ergonomic and natural to represent vectors as column vectors.
In short: if transformations go on the left and vectors go on the right, then transformations are rows and vectors are columns.