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)
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)
