Add util function to do local LU decomposition
Fix TSQR qrQR to work with irregular matrices
This also fixes a couple of bugs in RowPartitionedMatrix's apply
function.
Add 2D-block QR using Householder reconstruction
Tested to work correctly in cases where
the row block size is a multiple of the column block size
Update TODO