Add Rustiq-based synthesis for PauliEvolutionGate

crates/accelerate/Cargo.toml

@@ -26,6 +26,7 @@ qiskit-circuit.workspace = true
 thiserror.workspace = true
 ndarray_einsum_beta = "0.7"
 once_cell = "1.20.2"
+rustiq-core = "0.0.8"
 bytemuck.workspace = true
 
 [dependencies.smallvec]

crates/accelerate/src/synthesis/clifford/greedy_synthesis.rs

@@ -15,8 +15,8 @@ use indexmap::IndexSet;
 use ndarray::{s, ArrayView2};
 use smallvec::smallvec;
 
-use crate::synthesis::clifford::utils::CliffordGatesVec;
 use crate::synthesis::clifford::utils::{adjust_final_pauli_gates, SymplecticMatrix};
+use crate::synthesis::clifford::utils::{Clifford, CliffordGatesVec};
 use qiskit_circuit::operations::StandardGate;
 use qiskit_circuit::Qubit;
 
@@ -437,3 +437,14 @@ impl GreedyCliffordSynthesis<'_> {
         Ok((self.num_qubits, clifford_gates))
     }
 }
+
+/// Resynthesizes a clifford circuit using the greedy Clifford synthesis algorithm.
+pub fn resynthesize_clifford_circuit(
+    num_qubits: usize,
+    gates: &CliffordGatesVec,
+) -> Result<CliffordGatesVec, String> {
+    let sim_clifford = Clifford::from_gate_sequence(gates, num_qubits)?;
+    let mut synthesis = GreedyCliffordSynthesis::new(sim_clifford.tableau.view())?;
+    let (_, new_gates) = synthesis.run()?;
+    Ok(new_gates)
+}

crates/accelerate/src/synthesis/clifford/mod.rs

@@ -11,9 +11,9 @@
 // that they have been altered from the originals.
 
 mod bm_synthesis;
-mod greedy_synthesis;
+pub(crate) mod greedy_synthesis;
 mod random_clifford;
-mod utils;
+pub(crate) mod utils;
 
 use crate::synthesis::clifford::bm_synthesis::synth_clifford_bm_inner;
 use crate::synthesis::clifford::greedy_synthesis::GreedyCliffordSynthesis;

crates/accelerate/src/synthesis/clifford/utils.rs

@@ -149,6 +149,30 @@ impl Clifford {
         azip!((z in &mut z, &x in &x) *z ^= x);
     }
 
+    /// Modifies the tableau in-place by appending SX-gate
+    pub fn append_sx(&mut self, qubit: usize) {
+        let (mut x, z, mut p) = self.tableau.multi_slice_mut((
+            s![.., qubit],
+            s![.., self.num_qubits + qubit],
+            s![.., 2 * self.num_qubits],
+        ));
+
+        azip!((p in &mut p, &x in &x, &z in &z)  *p ^= !x & z);
+        azip!((&z in &z, x in &mut x) *x ^= z);
+    }
+
+    /// Modifies the tableau in-place by appending SXDG-gate
+    pub fn append_sxdg(&mut self, qubit: usize) {
+        let (mut x, z, mut p) = self.tableau.multi_slice_mut((
+            s![.., qubit],
+            s![.., self.num_qubits + qubit],
+            s![.., 2 * self.num_qubits],
+        ));
+
+        azip!((p in &mut p, &x in &x, &z in &z)  *p ^= x & z);
+        azip!((&z in &z, x in &mut x) *x ^= z);
+    }
+
     /// Modifies the tableau in-place by appending H-gate
     pub fn append_h(&mut self, qubit: usize) {
         let (mut x, mut z, mut p) = self.tableau.multi_slice_mut((
@@ -227,6 +251,18 @@ impl Clifford {
                     clifford.append_s(qubits[0].index());
                     Ok(())
                 }
+                StandardGate::SdgGate => {
+                    clifford.append_sdg(qubits[0].0 as usize);
+                    Ok(())
+                }
+                StandardGate::SXGate => {
+                    clifford.append_sx(qubits[0].0 as usize);
+                    Ok(())
+                }
+                StandardGate::SXdgGate => {
+                    clifford.append_sxdg(qubits[0].0 as usize);
+                    Ok(())
+                }
                 StandardGate::HGate => {
                     clifford.append_h(qubits[0].index());
                     Ok(())
@@ -283,21 +319,18 @@ pub fn adjust_final_pauli_gates(
     // add pauli gates
     for qubit in 0..num_qubits {
         if delta_phase_pre[qubit] && delta_phase_pre[qubit + num_qubits] {
-            // println!("=> Adding Y-gate on {}", qubit);
             gate_seq.push((
                 StandardGate::YGate,
                 smallvec![],
                 smallvec![Qubit::new(qubit)],
             ));
         } else if delta_phase_pre[qubit] {
-            // println!("=> Adding Z-gate on {}", qubit);
             gate_seq.push((
                 StandardGate::ZGate,
                 smallvec![],
                 smallvec![Qubit::new(qubit)],
             ));
         } else if delta_phase_pre[qubit + num_qubits] {
-            // println!("=> Adding X-gate on {}", qubit);
             gate_seq.push((
                 StandardGate::XGate,
                 smallvec![],