Full QR decomposition
atisharma opened this issue · 0 comments
atisharma commented
From the documentation, it seems org returns Q in the reduced QR decomposition (e.g. (org (qrf M)) ).
Whereas, MKL's ?orgqr allows the call to specify the number of columns to be returned for Q, allowing the full QR decomposition to be found. This is important to calculate the kernel of a matrix, for example.
Is there a way to specify the number of columns (or ask for the full Q) in org?
As a workaround, I can generate the full Q by right-padding the M with zeros, as per (org (qrf (hcat M Z))).