In the documentation for the QR decomposition with pivoting there is very little mention on how to obtain the Q and R matrices that result from the application of the function.

The one -somewhat suspicious- guideline, is to use the ogqr (real matrices) function.

However, this function is ignorant of the permutation vector.

Typically for the regular QR the input matrix is overwritten and the upper triangular part contains R and the lower contains info that allow for the extraction of Q - using ogqr.

So, how does it work now ? How does the pivoting info (i.e. column permutation) surfaces here ?

Thank you very much in advance,


PS: I am using MKL 10.3, Win7 and visual studio 2010 (i.e. c interface)




3 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

Apologies , for a moment I thought that orgqpr was mkl provided and not part of mkl.
I think I am good now.
Thank you,

..and not part of LAPACK ! ;-)

Leave a Comment

Please sign in to add a comment. Not a member? Join today