[SciPy-Dev] QR-decomposition with Q but not Q
Thu Aug 11 02:02:21 CDT 2011
A QR decomposition is done in two steps: firstly,
R is calculated and an intermediate form of Q
(the so called "elementary reflectors"), and in a second
step, Q is actually calculated.
The second step, however, can be costly and
thus is not done if you are not interested in Q. The
QR decomposition code already allows for not
Very often, however, one is interested in Q only in
order to multiply it with a vector c. This can be done
without ever calculating Q.
I just wrote a patch and submitted it to github, it's here:
which does just this. I'm not sure if I should write both
on scipy-dev and github, scipy-dev just seems to
have faster responses...
+49 30 6392 1234
More information about the SciPy-Dev