Biegamqr factorization wnumpy który zwraca listęndarrays, a mianowicieQiR:

>>> [q,r] = np.linalg.qr(np.array([1,0,0,0,1,1,1,1,1]).reshape(3,3))

R jest tablicą dwuwymiarową, mającą obrócone linie zerowe na dole (sprawdzone nawet dla wszystkich przykładów w moim zestawie testowym):

>>> print r
[[ 1.41421356  0.70710678  0.70710678]
 [ 0.          1.22474487  1.22474487]
 [ 0.          0.          0.        ]]

. Teraz chcę się podzielićR w dwóch macierzachR_~:

[[ 1.41421356  0.70710678  0.70710678]
 [ 0.          1.22474487  1.22474487]]

iR_0:

[[ 0.          0.          0.        ]]

(wyodrębnianie wszystkich linii zerowych). Wydaje się być blisko tego rozwiązania:usuwanie wierszy w tablicy numpy.

EDYTOWAĆ:
Jeszcze bardziej interesujące:np.linalg.qr() zwraca an x n-matryca. Nie, czego bym się spodziewał:

A := n x m
Q := n x m
R := n x m