Horizontal concatenation (hcat) for Tableaux and Stabilizers

CHANGELOG.md

@@ -5,6 +5,10 @@
 
 # News
 
+## v0.9.5 - dev
+
+- Addition of a horizontal concatenation for Tableaus -- `hcat`.
+
 ## v0.9.4 - 2024-06-28
 
 - Addition of a constructor for concatenated quantum codes -- `Concat`.

src/QuantumClifford.jl

@@ -922,6 +922,20 @@
 
 Base.vcat(stabs::Stabilizer...) = Stabilizer(vcat((tab(s) for s in stabs)...))
 
+function Base.hcat(tabs::Tableau...)
+    # Ensure all Tableaus have the same number of rows
+    nqubits = tabs[1].nqubits
+    if !all(t -> t.nqubits == nqubits, tabs)
+        throw(ArgumentError("All Tableaus must have the same number of qubits for horizontal concatenation."))
+    end
+    Tableau(
+    tabs[1].phases, #does not make sense, better method needed for phases concatenation
+    sum(s.nqubits for s in tabs),
+    hcat((s.xzs for s in tabs)...))
+end
+
+Base.hcat(stabs::Stabilizer...) = Stabilizer(hcat((tab(s) for s in stabs)...))
+
 ##############################
 # Unitary Clifford Operations
 ##############################