QR had some issues.  Code now solves least squares (and square systems); new apply_q subroutine; qr and qr_blas pass all tests
qrinsertcol handles empty matrix (returns a column)