Skip to content

Commit de12355

Browse files
committed
test/decompositions: add test for FullPivHouseholderQR
1 parent 6c85efc commit de12355

File tree

1 file changed

+18
-0
lines changed

1 file changed

+18
-0
lines changed

unittest/python/test_QR.py

Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -20,3 +20,21 @@
2020

2121
Y = householder_qr_eye.solve(X)
2222
assert (X == Y).all()
23+
24+
# Test FullPivHouseholderQR decomposition
25+
fullpiv_householder_qr = eigenpy.FullPivHouseholderQR()
26+
fullpiv_householder_qr = eigenpy.FullPivHouseholderQR(rows, cols)
27+
fullpiv_householder_qr = eigenpy.FullPivHouseholderQR(A)
28+
29+
fullpiv_householder_qr = eigenpy.FullPivHouseholderQR(np.eye(rows, rows))
30+
X = rng.random((rows, 20))
31+
assert fullpiv_householder_qr.absDeterminant() == 1.0
32+
assert fullpiv_householder_qr.logAbsDeterminant() == 0.0
33+
34+
Y = fullpiv_householder_qr.solve(X)
35+
assert (X == Y).all()
36+
assert fullpiv_householder_qr.rank() == rows
37+
38+
fullpiv_householder_qr.setThreshold(1e-8)
39+
assert fullpiv_householder_qr.threshold() == 1e-8
40+
assert eigenpy.is_approx(np.eye(rows, rows), fullpiv_householder_qr.inverse())

0 commit comments

Comments
 (0)