コード例 #1
0
 def test_determinant_of_4x4(self):
     a = Matrix([[-2, -8, 3, 5], [-3, 1, 7, 3], [1, 2, -9, 6],
                 [-6, 7, 7, -9]])
     assert a.cofactor(0, 0) == 690
     assert a.cofactor(0, 1) == 447
     assert a.cofactor(0, 2) == 210
     assert a.cofactor(0, 3) == 51
     assert a.determinant() == -4071
コード例 #2
0
 def test_inverse(self):
     a = Matrix([[-5, 2, 6, -8], [1, -5, 1, 8], [7, 7, -6, -7],
                 [1, -3, 7, 4]])
     b = a.inverse()
     assert a.determinant() == 532
     assert a.cofactor(2, 3) == -160
     assert b[3, 2] == -160 / 532
     assert a.cofactor(3, 2) == 105
     assert b[2, 3] == 105 / 532
     assert b == Matrix([[0.21805, 0.45113, 0.24060, -0.04511],
                         [-0.80827, -1.45677, -0.44361, 0.52068],
                         [-0.07895, -0.22368, -0.05263, 0.19737],
                         [-0.52256, -0.81391, -0.30075, 0.30639]])
コード例 #3
0
 def test_determinant_of_3x3(self):
     a = Matrix([[1, 2, 6], [-5, 8, -4], [2, 6, 4]])
     assert a.cofactor(0, 0) == 56
     assert a.cofactor(0, 1) == 12
     assert a.cofactor(0, 2) == -46
     assert a.determinant() == -196
コード例 #4
0
 def test_determinant_2x2_matrix(self):
     a = Matrix([[1, 5], [-3, 2]])
     assert a.determinant() == 17
コード例 #5
0
 def test_noninvertible_matrix_for_invertibility(self):
     a = Matrix([[-4, 2, -2, 3], [9, 6, 2, 6], [0, -5, 1, -5], [0, 0, 0,
                                                                0]])
     assert a.determinant() == 0
     assert not a.invertible
コード例 #6
0
 def test_invertible_matrix_for_invertibility(self):
     a = Matrix([[6, 4, 4, 4], [5, 5, 7, 6], [4, -9, 3, -7], [9, 1, 7, -6]])
     assert a.determinant() == -2120
     assert a.invertible