QR Decomposition With Pivoting

Open discussion for MAGMA library (Matrix Algebra on GPU and Multicore Architectures)
Post Reply
aada3
Posts: 1
Joined: Wed Jun 29, 2016 4:30 pm

QR Decomposition With Pivoting

Post by aada3 » Wed Jun 29, 2016 4:44 pm

I am trying to perform the QR Decomposition of a square, complex matrix with pivoting. Can I use cqeqp3_gpu to get R and the pivoting vector, followed by cungqr_gpu to output Q? The cungqr_gpu description mentions that it uses the outputs of cgeqrf_gpu, but does not mention cgeqp3_gpu.

mgates3
Posts: 918
Joined: Fri Jan 06, 2012 2:13 pm

Re: QR Decomposition With Pivoting

Post by mgates3 » Mon Jul 11, 2016 10:19 am

Yes, it appears that cungqr works with the output of cgeqp3. You can check by computing || A*P - Q*R || / (m * ||A||), where P is the permutation defined by jpvt. [See LAPACK working note 41 for error checks.]

-mark

Post Reply