[SciPy-User] qr decompostion gives negative q, r ?

Sturla Molden sturla@molden...
Wed Nov 21 08:15:15 CST 2012

On 21.11.2012 14:22, Alejandro Weinstein wrote:
> On Wed, Nov 21, 2012 at 1:34 AM, Virgil Stokes <vs@it.uu.se
> <mailto:vs@it.uu.se>> wrote:
>     But, again, it is an issue for the algorithm given in Table 3 (p.
>     2248 of paper).  Look at step 8. and equation (30).
> I see no Table  3, nor page 2248 in the paper you attached. Which paper
> are you talking about?

I think what he means is that in equation 22 (page 2591) S must be PD, 
and thus R = S' must be PD as well. Which (AFAIK) must be an error in 
the paper. Too bad referees didn't catch it.

They should have used Cholesky of P as S. So instead of the second 
equation in section IV which states that S = R', equation (6) should 
Cholesky factorize P to obtain S. And thus S would be PD for equation 
22. In conclusion, there is an error in the paper which was silenced by 
Matlab's QR.


More information about the SciPy-User mailing list