Find bases of matrix without multiplying This question is related to a solved problem in Gilbert Strang's 'Introduction to Linear Algebra'(Chapter 3,Question 3.6A, Page 190).

Q) Find bases and dimensions for all four fundamental subspaces of A if you know that
$A = \begin{bmatrix}1 & 0 & 0\\2 & 1 & 0 \\ 5 & 0 & 1\end{bmatrix}\begin{bmatrix}1 & 3 & 0 & 5\\0 & 0 & 1 & 6\\ 0 & 0 & 0 & 0\end{bmatrix} = LU = E^{-1}R$

Answer given in the text:

This matrix has pivots in columns $1$ and $3$. Its rank is $r=2$.
Column space : Basis $(1,2,5)$ and $(0,1,0)$ from $E^{-1}$

Why does he choose the first two columns of $E^{-1}$ as a basis ? If anything, the pivot columns of $R$ are an obvious choice for basis.

If you look closely at the right matrix, you'll see that


*

*It conflates the 2-dimensional subspace generated by $\{(1,0,0,0),(0,1,0,0)\}$ into the
1-dimensional subspace generated by $\{(1,0,0)\}$, because it's first two columns are linearly dependent. In particular, it maps vectors of the form $(3x,-x,0,0)$ to the zero vector, i.e. $\ker A \supset \{(3,-1,0,0)\}$. Columns $1,3,4$ are also obviously linearly dependent, and that yields $\ker A \supset \{5,0,6,-1\}$. 

*It's range (i.e. what you call column space I think) is generated by $\{(1,0,0),(0,1,0)\}$, because it's bottom row is all zeros. 
The left matrix is invertible (lower triangular matrix, non-zero diagonal), and it follows that


*

*$\textrm{rng } A = \textrm{span} \{(1,2,5),(0,1,0)\}$. This is simply the image of the range of the right matrix under the left matrix. This is easy to find because the right matrix's image is created by two canonical basis vectors, so we just had to pick the corresponding columns of the left matrix.

*$\textrm{ker } A \supset \textrm{span} \{(3,-1,0,0),(5,0,6,-1)\}$. And because of the dimension theorem, we can conclude from (1) that in fact $\textrm{ker } A = \textrm{span} \{(3,-1,0,0),(5,0,6,-1)\}$.
