Questions about QR blocked algorithm described in # 191

Open forum for general discussions relating to PLASMA.

Questions about QR blocked algorithm described in # 191

Postby linyufly » Tue Aug 24, 2010 11:36 pm

As far as I am concerned, given m >= n, the result of QR(A_m,n) would be V_m,n T_n,n and R_m,n.

In LAPACK Working Note # 191, the DTSQRT is described below:

( Rkk) (Rkk)
( ) -----> (Vik, Tik, Rkk) = QR( )
( Aik ) (Aik )

Since the coupling of Rkk and Aik is a (2b)X(b) block, I think m = 2b and n = b here. Therefore I think Vik and Rkk are no longer bXb blocks but (2b)X(b). How can these blocks fit into original matrix? For Rkk, maybe we can just ignore the lower bXb zero block. But I don't think we can do the same thing for Vik.

Thank you!
linyufly
 
Posts: 1
Joined: Tue Aug 24, 2010 11:33 pm

Re: Questions about QR blocked algorithm described in # 191

Postby ltaief » Tue Sep 07, 2010 12:00 am

The size of Vik is indeed (2b)x(b). The lower part of size bxb contains the
computed reflectors and the upper part of size bxb is the Identity matrix
which does not need, therefore, to be stored.
Hope this helps!
Hatem
ltaief
 
Posts: 3
Joined: Fri Dec 11, 2009 11:32 am


Return to User discussion

Who is online

Users browsing this forum: No registered users and 1 guest