def setUp(self): self.N = 17 # Number of rows to create -- esssentially arbitrary # Random complex matrix of unit-length vectors self.psi = random_array((self.N, 2), iscomplex=True) self.psi = normalize(self.psi) # Random real matrix of unit length vectors self.X = random_array((self.N, 3)) self.X = normalize(self.X)
def test_su_to_SU(self): """Test that the image of generic Cayley map ends up in SU(2).""" N = 16 a = random_array((N, 3)) A = hatmap(a) U = cayley(A) self.assertTrue(is_SU(U).all())
def testCayleyKlein(self): """Test whether generic Cayley and Cayley-Klein agree on su(2).""" N = 16 a = random_array((N, 3)) A = hatmap(a) U1 = cayley(A) U2 = cayley_klein(a) self.assertTrue(is_float_equal(U1, U2))
def setUp(self): self.N = 5 self.a = random_array((self.N, 3)) self.A = hatmap(self.a)