Support save_statevector for QuantumCircuit from QASM3 string

qiskit_aer/library/default_qubits.py

@@ -48,8 +48,10 @@ def default_qubits(circuit, qubits=None):
         for tuple_element in tuples:
             if isinstance(tuple_element, QuantumRegister):
                 for j in range(tuple_element.size):
-                    qubits.append(tuple_element[j])
+                    if tuple_element[j] not in qubits:
+                        qubits.append(tuple_element[j])
             else:
-                qubits.append(tuple_element)
+                if tuple_element not in qubits:
+                    qubits.append(tuple_element)
 
     return qubits

releasenotes/notes/save_statevector_for_qasm3_circ-642ade99af3ff0d2.yaml

@@ -0,0 +1,7 @@
+---
+fixes:
+  - |
+    :meth:``QuantumCircuit.save_statevector()`` does not work if the circuit
+    is generated from OpenQASM3 text because its quantum registers have duplicated
+    qubit instances. With this commit, :meth:``QuantumCircuit.save_statevector()``
+    retrieves qubits from quantum registers without duplication.

test/terra/backends/aer_simulator/test_save_statevector.py

@@ -17,6 +17,7 @@
 import qiskit.quantum_info as qi
 from qiskit import QuantumCircuit, transpile
 from test.terra.backends.simulator_test_case import SimulatorTestCase, supported_methods
+from qiskit.qasm3 import dumps, loads
 
 
 @ddt
@@ -201,3 +202,40 @@ def test_save_statevector_cache_blocking(self, method, device):
         self.assertIn(label, simdata)
         value = simdata[label]
         self.assertEqual(value, target)
+
+    @supported_methods(
+        [
+            "automatic",
+            "statevector",
+            "matrix_product_state",
+            "extended_stabilizer",
+            "tensor_network",
+        ]
+    )
+    def test_save_statevector_for_qasm3_circuit(self, method, device):
+        """Test save statevector instruction"""
+        backend = self.backend(method=method, device=device)
+
+        # Stabilizer test circuit
+        circ = QuantumCircuit(3)
+        circ.h(0)
+        circ.sdg(0)
+        circ.cx(0, 1)
+        circ.cx(0, 2)
+
+        # Target statevector
+        target = qi.Statevector(circ)
+
+        circ = loads(dumps(circ))
+
+        # Add save to circuit
+        label = "state"
+        circ.save_statevector(label=label)
+
+        # Run
+        result = backend.run(transpile(circ, backend, optimization_level=0), shots=1).result()
+        self.assertTrue(result.success)
+        simdata = result.data(0)
+        self.assertIn(label, simdata)
+        value = simdata[label]
+        self.assertEqual(value, target)