I was looking for a solution that minimizes memory allocation because I’m solving several systems in a loop. Since my systems have full rank, the cheaper LQ factorization can achieve the same result as pivoted QR.
I was wondering if there is any function available which can use the LQ factorization itself, without constructing the L and Q matrices explicitly.