https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/DORGQR-First-dimension-of-matrix-A/m-p/817901#M4442 I do not understand why you want to form Q explicitly at the same time that you are trying to avoid forming B in full.<BR /><BR />Let B = [b1 b2 b3 ... ], Y = [y1 y2 y3 ...], e = [e1 e2 e3 ... ]<BR /><BR />where b1, b2, b3... are individual columns of B or, at your choosing, groups of columns that you are willing to process together.<BR /><BR />You call DGEQRF once to perfom the QR decomposition of X, with Q stored implicitly in the result.<BR /><BR />You then call DORMQR to solve for b1, which you then use to compute e1 = X b1 - y1, and similarly for e2, e3, ... For this to work, you need to save a copy of X before the call to DGEQRF, which overwrites X. Or, it may be just as easy to regenerate the necessary portions of X.<BR /><BR />Would this approach satisfy your needs?<BR /><BR />Minor note: you have inconsistent signs for the error matrix e in your post. Tue, 24 May 2011 14:32:31 GMT mecej4 2011-05-24T14:32:31Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/DORGQR-First-dimension-of-matrix-A/m-p/817900#M4441 <DIV><DIV>Hi All.</DIV><DIV>I have multiple regression model X * B = Y + e, where X - Txm, B - mxn, Y and e - Txn and I want to get residuals "e", without actualy calculation B. So I want to use QR decomposition of matrix X, and do it as following: e = Y - Q *(Q'*Y); So I want to use DORGQR twice. Assume, we use all elementary reflectors, and Q is Txm matrix.</DIV><DIV>Everything ok with the first call, I get Q'Y which is mxn matrix. But there is a problem with second call, in accordance with the article:</DIV><DIV><A href="http://software.intel.com/sites/products/documentation/hpc/mkl/mklman/index.htm#lse/functn_ormqr.htm" target="_blank">http://software.intel.com/sites/products/documentation/hpc/mkl/mklman/index.htm#lse/functn_ormqr.htm</A></DIV><DIV>I don't have any way to specify the number of rows in matrix Q. In section "Fortran 95 Interface Notes" (it's not very clear for me, why here) I see the note:</DIV><DIV>a:</DIV><DIV>Holds the matrix A of size (r,k).</DIV><DIV>r = m if side = 'L'.</DIV><DIV>r = n if side = 'R'</DIV><DIV>So I just want to multiply Q (Txm) and Q'R (mxn), but routine automatically assumes that Q is mxk which is definitely not what I want. May be somebody can give me a hint? What I'm doing wrong? Extraction Q matrix explicitly doesn't seem to be a good solution for me, I want to reduce memory usage.</DIV></DIV> Mon, 23 May 2011 21:43:05 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/DORGQR-First-dimension-of-matrix-A/m-p/817900#M4441 wisdoms_wind 2011-05-23T21:43:05Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/DORGQR-First-dimension-of-matrix-A/m-p/817901#M4442 I do not understand why you want to form Q explicitly at the same time that you are trying to avoid forming B in full.<BR /><BR />Let B = [b1 b2 b3 ... ], Y = [y1 y2 y3 ...], e = [e1 e2 e3 ... ]<BR /><BR />where b1, b2, b3... are individual columns of B or, at your choosing, groups of columns that you are willing to process together.<BR /><BR />You call DGEQRF once to perfom the QR decomposition of X, with Q stored implicitly in the result.<BR /><BR />You then call DORMQR to solve for b1, which you then use to compute e1 = X b1 - y1, and similarly for e2, e3, ... For this to work, you need to save a copy of X before the call to DGEQRF, which overwrites X. Or, it may be just as easy to regenerate the necessary portions of X.<BR /><BR />Would this approach satisfy your needs?<BR /><BR />Minor note: you have inconsistent signs for the error matrix e in your post. Tue, 24 May 2011 14:32:31 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/DORGQR-First-dimension-of-matrix-A/m-p/817901#M4442 mecej4 2011-05-24T14:32:31Z