search

uncomplicate / neanderthal

Fast Clojure Matrix Library
http://neanderthal.uncomplicate.org
Eclipse Public License 1.0
1.06k stars 56 forks source link

Full QR decomposition #94

Open atisharma opened 4 years ago

atisharma commented 4 years ago

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))).