✨ Added clifford+t gateset

pyproject.toml

@@ -111,6 +111,8 @@ filterwarnings = [
     'ignore:.*Values in x.*:RuntimeWarning:',
     'ignore:.*divide by zero encountered in det.*:RuntimeWarning:',
     'ignore:.*invalid value encountered in det.*:RuntimeWarning:',
+    'ignore:.*invalid value encountered in divide.*:RuntimeWarning:',
+    'ignore:.*The iteration is not making good progress.*:RuntimeWarning:',
     # Qiskit 1.2 deprecations
     'ignore:.*The class ``qiskit.qobj.pulse_qobj.Pulse.*`` is deprecated as of qiskit 1.2.*:DeprecationWarning:qiskit.*',
     # Qiskit 1.3 deprecations
@@ -123,6 +125,7 @@ filterwarnings = [
     'ignore:.*Providing non-standard gates.*is deprecated for both ``transpile`` and ``generate_preset_pass_manager`` as of Qiskit 1.3.0.*:DeprecationWarning:qiskit.*',
     'ignore:.*is pending deprecation as of qiskit 1.3..*:PendingDeprecationWarning:',
     'ignore:.*The class ``qiskit.providers.exceptions.BackendPropertyError`` is deprecated as of qiskit 1.4.*:DeprecationWarning:',
+    'ignore:.*he SolovayKitaev algorithm relies *.:RuntimeWarning:qiskit.synthesis.discrete_basis.generate_basis_approximations',
 ]
 
 [tool.coverage]

src/mqt/bench/benchmark_generation.py

@@ -282,7 +282,31 @@ def get_native_gates_level(
     if file_precheck and path.is_file():
         return True
 
-    compiled = transpile(qc, basis_gates=gateset.gates, optimization_level=opt_level, seed_transpiler=10)
+    if gateset.name == "clifford+t":
+        from qiskit.transpiler import PassManager  # noqa: PLC0415
+        from qiskit.transpiler.passes.synthesis import SolovayKitaev  # noqa: PLC0415
+
+        # Transpile the circuit to single- and two-qubit gates including rotations
+        compiled_for_sk = transpile(
+            qc,
+            basis_gates=[*gateset.gates, "rx", "ry", "rz"],
+            optimization_level=opt_level,
+            seed_transpiler=10,
+        )
+        # Synthesize the rotations to Clifford+T gates
+        # Measurements are removed and added back after the synthesis to avoid errors in the Solovay-Kitaev pass
+        pm = PassManager(SolovayKitaev())
+        new_qc = pm.run(compiled_for_sk.remove_final_measurements(inplace=False))
+        new_qc.measure_all()
+        # Transpile once more to remove unnecessary gates and optimize the circuit
+        compiled = transpile(
+            new_qc,
+            basis_gates=gateset.gates,
+            optimization_level=opt_level,
+            seed_transpiler=10,
+        )
+    else:
+        compiled = transpile(qc, basis_gates=gateset.gates, optimization_level=opt_level, seed_transpiler=10)
 
     if return_qc:
         return compiled

src/mqt/bench/devices/__init__.py

@@ -88,6 +88,33 @@ def get_available_native_gatesets() -> list[Gateset]:
     for device in get_available_devices():
         if device.gateset not in available_gatesets:
             available_gatesets.append(device.gateset)
+    available_gatesets.append(
+        Gateset(
+            "clifford+t",
+            [
+                "i",
+                "x",
+                "y",
+                "z",
+                "h",
+                "s",
+                "sdg",
+                "t",
+                "tdg",
+                "sx",
+                "sxdg",
+                "cx",
+                "cy",
+                "cz",
+                "swap",
+                "iswap",
+                "dcx",
+                "ecr",
+                "measure",
+                "barrier",
+            ],
+        )
+    )
     return available_gatesets
 
 

tests/test_bench.py

@@ -128,7 +128,7 @@ def sample_filenames() -> list[str]:
         (vqesu2random, 3, True),
         (vqetwolocalrandom, 3, True),
         (wstate, 3, True),
-        # (shor, 3, False),
+        (shor, 3, False),
     ],
 )
 def test_quantumcircuit_alg_level(
@@ -568,6 +568,19 @@ def test_oqc_benchmarks() -> None:
     path.unlink()
 
 
+def test_clifford_t() -> None:
+    """Test the Clifford+T gateset."""
+    qc = get_benchmark(
+        benchmark_name="qft",
+        level="nativegates",
+        circuit_size=4,
+        gateset="clifford+t",
+    )
+
+    for gate_type in qc.count_ops():
+        assert gate_type in get_native_gateset_by_name("clifford+t").gates
+
+
 def test_get_module_for_benchmark() -> None:
     """Test the get_module_for_benchmark function."""
     for benchmark in get_supported_benchmarks():