implement sSQRTXX, sInvSQRTXX, sSQRTYY, sInvSQRTYY

QuantumSavory · Sep 24, 2024 · f574c34 · f574c34
1 parent 05491d9
commit f574c34
Show file tree

Hide file tree

Showing 3 changed files with 20 additions and 1 deletion.
diff --git a/src/QuantumClifford.jl b/src/QuantumClifford.jl
@@ -51,7 +51,7 @@ export
     sHadamardXY, sHadamardYZ, sSQRTX, sInvSQRTX, sSQRTY, sInvSQRTY, sCXYZ, sCZYX,
     sCNOT, sCPHASE, sSWAP,
     sXCX, sXCY, sXCZ, sYCX, sYCY, sYCZ, sZCX, sZCY, sZCZ,
-    sZCrY, sInvZCrY,
+    sZCrY, sInvZCrY, sSQRTXX, sInvSQRTXX, sSQRTYY, sInvSQRTYY,
     # Misc Ops
     SparseGate,
     sMX, sMY, sMZ, PauliMeasurement, Reset, sMRX, sMRY, sMRZ,

diff --git a/src/symbolic_cliffords.jl b/src/symbolic_cliffords.jl
@@ -326,6 +326,12 @@ end
 @qubitop2 ZCrY  (x1, x1⊻z1⊻x2⊻z2, x1⊻x2, x1⊻z2, ~iszero((x1 & ~z1 & x2) | (x1 & ~z1 & ~z2) | (x1 & x2 & ~z2)))
 @qubitop2 InvZCrY (x1, x1⊻z1⊻x2⊻z2, x1⊻x2, x1⊻z2, ~iszero((x1 & z1 & ~x2 & ~z2) | (x1 & ~z1 & ~x2 & z2) | (x1 & z1 & ~x2 & z2) | (x1 & z1 & x2 & z2)))
 
+@qubitop2 SQRTXX      (z1⊻z2⊻x1, z1      , z1⊻x2⊻z2, z2      , ~iszero((~x1 & z1 &~z2) | ( ~z1 &~x2 & z2)))
+@qubitop2 InvSQRTXX   (z1⊻z2⊻x1, z1      , z1⊻x2⊻z2, z2      , ~iszero((x1 & z1 & ~z2) | ( ~z1 & x2 & z2)))
+
+@qubitop2 SQRTYY     (z1⊻x2⊻z2, x1⊻z2⊻x2, x1⊻z1⊻z2, x1⊻x2⊻z1, ~iszero((~x1 &~z1 & x2 &~z2) | ( x1 &~z1 &~x2 &~z2) | ( x1 &~z1 & x2 & z2) | ( x1 & z1 & x2 &~z2)))
+@qubitop2 InvSQRTYY  (z1⊻x2⊻z2, x1⊻z2⊻x2, x1⊻z1⊻z2, x1⊻x2⊻z1, ~iszero(( x1 & z1 &~x2 & z2) | (~x1 & z1 & x2 & z2) | (~x1 & z1 &~x2 &~z2) | (~x1 &~z1 &~x2 & z2)))
+
 #=
 To get the boolean formulas for the phase, it is easiest to first write down the truth table for the phase:
 for i in 0:15
@@ -384,6 +390,10 @@ LinearAlgebra.inv(op::sYCY) = sYCY(op.q1, op.q2)
 LinearAlgebra.inv(op::sYCZ) = sYCZ(op.q1, op.q2)
 LinearAlgebra.inv(op::sZCrY) = sInvZCrY(op.q1, op.q2)
 LinearAlgebra.inv(op::sInvZCrY) = sZCrY(op.q1, op.q2)
+LinearAlgebra.inv(op::sSQRTXX) = sInvSQRTXX(op.q1, op.q2)
+LinearAlgebra.inv(op::sInvSQRTXX) = sSQRTXX(op.q1, op.q2)
+LinearAlgebra.inv(op::sSQRTYY) = sInvSQRTYY(op.q1, op.q2)
+LinearAlgebra.inv(op::sInvSQRTYY) = sSQRTYY(op.q1, op.q2)
 
 ##############################
 # Functions that perform direct application of common operators without needing an operator instance

diff --git a/test/test_symcliff.jl b/test/test_symcliff.jl
@@ -76,6 +76,7 @@
             @test CliffordOperator(inv(random_op), i) == inv(CliffordOperator(random_op, i))
             @test CliffordOperator(inv(SingleQubitOperator(random_op)), i) == inv(CliffordOperator(random_op, i))
         end
+    end
 
     @testset "Consistency checks with Stim" begin
        # see https://github.com/quantumlib/Stim/blob/main/doc/gates.md
@@ -101,4 +102,12 @@
             @test CliffordOperator(inv(sInvZCrY(n₁, n₂)), n₁) == inv(CliffordOperator(sInvZCrY(n₁, n₂), n₁))
         end
     end
+
+    @testset "Consistency check with STIM conventions" begin
+        # see https://github.com/quantumlib/Stim/blob/main/doc/gates.md
+        @test CliffordOperator(sSQRTXX)       == C"XI IX -YX -XY"
+        @test CliffordOperator(sInvSQRTXX)    == C"XI IX YX XY"
+        @test CliffordOperator(sSQRTYY)       == C"-ZY -YZ XY YX"
+        @test CliffordOperator(sInvSQRTYY)    == C"ZY YZ -XY -YX"
+    end
 end