qr, qr_blas, and form_qr can now handle size Nr>=Nc; doesn't yet take care of RHS of linear system