def test_get_pauli_2_qubit_proj_3_1_2(self): num_qubits = 3 qubit_pro1 = 1 qubit_pro2 = 2 paulis = get_pauli_n_qubit_projection(num_qubits, [qubit_pro1, qubit_pro2]) self.assertTrue(len(paulis) == 16) X = np.array([[0, 1], [1, 0]], dtype=np.complex128) Y = np.array([[0, -1j], [1j, 0]], dtype=np.complex128) Z = np.array([[1, 0], [0, -1]], dtype=np.complex128) I = np.array([[1, 0], [0, 1]], dtype=np.complex128) self.assertTrue(self.in_array(np.kron(np.kron(I, X), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Y), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Z), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, I), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, X), X), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Y), X), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Z), X), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, I), X), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, X), Y), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Y), Y), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Z), Y), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, I), Y), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, X), Z), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Y), Z), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, Z), Z), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, I), Z), paulis))
def test_fixedgate_get_unitary(self): reset_tensor_cache() fg = FixedGate("Test", 4, 2, (0, 1)) with tf.Session() as sess: sess.run(tf.global_variables_initializer()) gate = fg.get_unitary(sess) paulis = get_pauli_n_qubit_projection(4, (0, 1)) H = pauli_dot_product([0.25] * 16, paulis) U = la.expm(1j * H) self.assertTrue(hilbert_schmidt_distance(gate, U) <= 1e-16)
def test_fixedgate_get_herm ( self ): reset_tensor_cache() fg = FixedGate( "Test", 4, 2, (0, 1) ) herm = fg.get_herm() with tf.Session() as sess: sess.run( tf.global_variables_initializer() ) herm = herm.eval() paulis = get_pauli_n_qubit_projection( 4, (0, 1) ) H = pauli_dot_product( [ 0.25 ] * 16, paulis ) self.assertTrue( np.array_equal( herm, H ) )
def test_genericgate_get_unitary(self): reset_tensor_cache() lm = LocationModel(4, 2) gg = GenericGate("Test", 4, 2, lm, loc_vals=[1, 0, 0, 0, 0, 0]) with tf.Session() as sess: sess.run(tf.global_variables_initializer()) gate = gg.get_unitary(sess) paulis = get_pauli_n_qubit_projection(4, (0, 1)) H = pauli_dot_product([0.25] * 16, paulis) U = la.expm(1j * H) self.assertTrue(hilbert_schmidt_distance(gate, U) <= 1e-16)
def test_genericgate_get_herm(self): reset_tensor_cache() lm = LocationModel(4, 2) gg = GenericGate("Test", 4, 2, lm, loc_vals=[1, 0, 0, 0, 0, 0]) herm = gg.get_herm() with tf.Session() as sess: sess.run(tf.global_variables_initializer()) herm = herm.eval() paulis = get_pauli_n_qubit_projection(4, list(lm.locations)[0]) H = pauli_dot_product([0.25] * 16, paulis) self.assertTrue(np.allclose(herm, H, rtol=0, atol=1e-15))
def test_get_pauli_1_qubit_proj_3_0(self): num_qubits = 3 qubit_proj = 0 paulis = get_pauli_n_qubit_projection(num_qubits, [qubit_proj]) self.assertTrue(len(paulis) == 4) X = np.array([[0, 1], [1, 0]], dtype=np.complex128) Y = np.array([[0, -1j], [1j, 0]], dtype=np.complex128) Z = np.array([[1, 0], [0, -1]], dtype=np.complex128) I = np.array([[1, 0], [0, 1]], dtype=np.complex128) self.assertTrue(self.in_array(np.kron(np.kron(X, I), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(Y, I), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(Z, I), I), paulis)) self.assertTrue(self.in_array(np.kron(np.kron(I, I), I), paulis))