def test_pauli_invert(self):
self.log.info("===== p3 =====")
self.log.info(self.p3)
self.assertEqual(str(self.p3), 'v = 1.0\t0.0\t1.0\t\nw = 0.0\t1.0\t1.0\t')
self.log.info("\tIn label form:")
self.log.info(self.p3.to_label())
self.assertEqual(self.p3.to_label(), 'ZXY')
self.log.info("\tIn matrix form:")
self.log.info(self.p3.to_matrix())
m = np.array([
[0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. - 1.j, 0. + 0.j],
[0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 1.j],
[0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. - 1.j, 0. + 0.j, 0. + 0.j, 0. + 0.j],
[0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 1.j, 0. + 0.j, 0. + 0.j],
[0. + 0.j, 0. + 0.j, 0. + 1.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j],
[0. + 0.j, 0. - 0.j, 0. + 0.j, 0. - 1.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j],
[0. + 1.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j],
[0. + 0.j, 0. - 1.j, 0. + 0.j, 0. - 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j, 0. + 0.j]])
self.assertTrue((self.p3.to_matrix() == m).all())
self.log.info("===== r =====")
r = inverse_pauli(self.p3)
self.assertEqual(str(r), 'v = 1.0\t0.0\t1.0\t\nw = 0.0\t1.0\t1.0\t')
self.log.info("In label form:")
self.log.info(r.to_label())
self.assertEqual(r.to_label(), 'ZXY')
self.log.info("\tIn matrix form:")
self.assertTrue((r.to_matrix() == m).all())
def test_pauli(self):
v = np.zeros(3)
w = np.zeros(3)
v[0] = 1
w[1] = 1
v[2] = 1
w[2] = 1
p = Pauli(v, w)
self.log.info(p)
self.log.info("In label form:")
self.log.info(p.to_label())
self.log.info("In matrix form:")
self.log.info(p.to_matrix())
q = random_pauli(2)
self.log.info(q)
r = inverse_pauli(p)
self.log.info("In label form:")
self.log.info(r.to_label())
self.log.info("Group in tensor order:")
grp = pauli_group(3, case=1)
for j in grp:
self.log.info(j.to_label())
self.log.info("Group in weight order:")
grp = pauli_group(3)
for j in grp:
self.log.info(j.to_label())
self.log.info("sign product:")
p1 = Pauli(np.array([0]), np.array([1]))
p2 = Pauli(np.array([1]), np.array([1]))
p3, sgn = sgn_prod(p1, p2)
self.log.info(p1.to_label())
self.log.info(p2.to_label())
self.log.info(p3.to_label())
self.log.info(sgn)
self.log.info("sign product reverse:")
p3, sgn = sgn_prod(p2, p1)
self.log.info(p2.to_label())
self.log.info(p1.to_label())
self.log.info(p3.to_label())
self.log.info(sgn)
def test_pauli(self):
v = np.zeros(3)
w = np.zeros(3)
v[0] = 1
w[1] = 1
v[2] = 1
w[2] = 1
p = Pauli(v, w)
self.log.info(p)
self.log.info("In label form:")
self.log.info(p.to_label())
self.log.info("In matrix form:")
self.log.info(p.to_matrix())
q = random_pauli(2)
self.log.info(q)
r = inverse_pauli(p)
self.log.info("In label form:")
self.log.info(r.to_label())
self.log.info("Group in tensor order:")
grp = pauli_group(3, case=1)
for j in grp:
self.log.info(j.to_label())
self.log.info("Group in weight order:")
grp = pauli_group(3)
for j in grp:
self.log.info(j.to_label())
self.log.info("sign product:")
p1 = Pauli(np.array([0]), np.array([1]))
p2 = Pauli(np.array([1]), np.array([1]))
p3, sgn = sgn_prod(p1, p2)
self.log.info(p1.to_label())
self.log.info(p2.to_label())
self.log.info(p3.to_label())
self.log.info(sgn)
self.log.info("sign product reverse:")
p3, sgn = sgn_prod(p2, p1)
self.log.info(p2.to_label())
self.log.info(p1.to_label())
self.log.info(p3.to_label())
self.log.info(sgn)
