def test_unitary_qsd_5qubits(self):
u = unitary_group.rvs(16)
gate = unitary(u, 'qsd')
state = get_state(gate)
transpiled_circuit = qiskit.transpile(gate, basis_gates=['u', 'cx'])
n_cx = transpiled_circuit.count_ops()['cx']
self.assertTrue(n_cx <= 120)
self.assertTrue(np.allclose(u[:, 0], state))
def test_initialize(self):
a = np.random.rand(16) + np.random.rand(16) * 1j
a = a / np.linalg.norm(a)
circ = initialize(a)
state = get_state(circ)
self.assertTrue(np.allclose(a, state))
def test_unitary_qsd_4qubits(self):
u = unitary_group.rvs(16)
gate = unitary(u, 'qsd')
state = get_state(gate)
self.assertTrue(np.allclose(u[:, 0], state))
circuit = qiskit.QuantumCircuit(4)
circuit.x(0)
circuit.append(gate, circuit.qubits)
state = get_state(circuit)
self.assertTrue(np.allclose(u[:, 1], state))
circuit = qiskit.QuantumCircuit(4)
circuit.x(1)
circuit.append(gate, circuit.qubits)
state = get_state(circuit)
self.assertTrue(np.allclose(u[:, 2], state))
circuit = qiskit.QuantumCircuit(4)
circuit.x([0, 1, 2, 3])
circuit.append(gate, circuit.qubits)
state = get_state(circuit)
self.assertTrue(np.allclose(u[:, 15], state))
def test_pivoting_3nonzero(self):
vector = [
-1 / np.sqrt(4), 0, 0, 0, 0, 0, 0, 1 / np.sqrt(4), 0, 0, 0,
1 / np.sqrt(4), 0, 0, 0, 1 / np.sqrt(4)
]
vector2 = {}
for k, value in enumerate(vector):
if value != 0:
index = format(k, '04b')
vector2[index] = value
circ = initialize(vector)
n_qubits = int(np.log2(len(vector)))
txt = "{0:0" + str(n_qubits) + "b}"
index_zero = txt.format(1)
index_nonzero = txt.format(7)
circuit, next_state = _pivoting(index_zero,
index_nonzero,
2,
state=vector2)
circ.compose(circuit.reverse_bits(), circ.qubits, inplace=True)
index_zero = txt.format(2)
index_nonzero = txt.format(9)
circuit, next_state = _pivoting(index_zero,
index_nonzero,
2,
state=next_state)
circ.compose(circuit.reverse_bits(), circ.qubits, inplace=True)
vector2 = get_state(circ)
self.assertTrue(
np.allclose(vector2[:3], [vector[0], vector[7], vector[11]]))
def test_sparse_initialize(self):
s = 3
n = 8
vector = np.zeros(2**n)
for k in range(2**s):
index = np.random.randint(0, 2**n)
while vector[index] != 0.0:
index = np.random.randint(0, 2**n)
vector[index] = np.random.rand() # + np.random.rand() * 1j
vector = vector / np.linalg.norm(vector)
vector2 = {}
for index, value in enumerate(vector):
if not np.isclose(value, 0.0):
txt = '{0:0' + str(n) + 'b}'
index_txt = txt.format(index)
vector2[index_txt] = vector[index]
circ = sparse_initialize(vector2)
calc_vector = get_state(circ)
self.assertTrue(np.allclose(vector, calc_vector))
def test_unitary_csd_2qubits(self):
u = unitary_group.rvs(4)
gate = unitary(u)
state = get_state(gate)
self.assertTrue(np.allclose(u[:, 0], state))