def testDiagonalization(self): vals = jacobi_math.diagonalize_2by2(array([2, 1, 1, 2]).reshape(2, 2)) assert_array_almost_equal(vals[1], array([3, 0, 0, 1]).reshape(2, 2), 10) assert_array_almost_equal( vals[0], ((1 / math.sqrt(2)) * array([1, 1, 1, -1])).reshape(2, 2))
def testDiagonalization(self): vals = jacobi_math.diagonalize_2by2(array([2,1,1,2]).reshape(2,2)) assert_array_almost_equal(vals[1], array([3,0,0,1]).reshape(2,2), 10) assert_array_almost_equal(vals[0], ((1/math.sqrt(2))*array([1,1,1,-1])).reshape(2,2))
def testSanity(self): for matr in self.matrices: vals = jacobi_math.diagonalize_2by2(matr) res = dot(dot(vals[0], vals[1]), vals[2]) assert_array_almost_equal(matr, dot(dot(vals[0], vals[1]), vals[2]), 10)