complete single-qubit

Author: sagnikpal2004

src/QuantumClifford.jl

@@ -53,7 +53,7 @@ export
     # Symbolic Clifford Ops
     AbstractSymbolicOperator, AbstractSingleQubitOperator, AbstractTwoQubitOperator,
     sHadamard, sPhase, sInvPhase, SingleQubitOperator, sId1, sX, sY, sZ,
-    sHadamardXY, sHadamardYZ, sSQRTX, sInvSQRTX, sSQRTY, sInvSQRTY, sCXYZ, sCZYX,
+    sHadamardXY, sHadamardYZ, sSQRTX, sInvSQRTX, sSQRTY, sInvSQRTY, sSQRTZ, sInvSQRTZ, sCXYZ, sCZYX,
     sCNOT, sCPHASE, sSWAP,
     sXCX, sXCY, sXCZ, sYCX, sYCY, sYCZ, sZCX, sZCY, sZCZ,
     sZCrY, sInvZCrY, sSWAPCX, sInvSWAPCX, sCZSWAP, sCXSWAP, sISWAP, sInvISWAP,

src/apply_right.jl

@@ -46,13 +46,16 @@ function apply_right!(l::CliffordOperator, r::sHadamardYZ)
     return l
 end
 
-# function apply_right!(l::CliffordOperator, r::sPhase)
-#     return l
-# end
+function apply_right!(l::CliffordOperator, r::sPhase)
+    l[r.q] = mul_right_log_i!(l[r.q], l[nqubits(l)+r.q])
+    apply_right!(l, sZ(r.q))
+    return l
+end
 
-# function apply_right!(l::CliffordOperator, r::sInvPhase)
-#     return l
-# end
+function apply_right!(l::CliffordOperator, r::sInvPhase)
+    l[r.q] = mul_right_log_i!(l[r.q], l[nqubits(l)+r.q])
+    return l
+end
 
 function apply_right!(l::CliffordOperator, r::sX)
     phases(tab(l))[nqubits(l)+r.q] ⊻= 0x02
@@ -104,13 +107,18 @@ function apply_right!(l::CliffordOperator, r::sInvSQRTZ)
     return l
 end
 
-# function apply_right!(l::CliffordOperator, r::sCXYZ)
-#     return l
-# end
+function apply_right!(l::CliffordOperator, r::sCXYZ)
+    rowswap!(tab(l), r.q, nqubits(l)+r.q)
+    l[r.q] = mul_right_log_i!(l[r.q], l[nqubits(l)+r.q])
+    return l
+end
 
-# function apply_right!(l::CliffordOperator, r::sCZYX)
-#     return l
-# end
+function apply_right!(l::CliffordOperator, r::sCZYX)
+    rowswap!(tab(l), r.q, nqubits(l)+r.q)
+    l[nqubits(l)+r.q] = mul_right_log_i!(l[nqubits(l)+r.q], l[r.q])
+    apply_right!(l, sX(r.q))
+    return l
+end
 
 function apply_right!(l::CliffordOperator, ::sId1)
     return l
@@ -128,6 +136,7 @@ function apply_right!(l::CliffordOperator, r::sSWAP)
 end
 
 # function apply_right!(l::CliffordOperator, r::sSWAPCX)
+#     rowswap!(tab(l), r.q1, r.q2)
 #     return l
 # end
 

test/test_apply_right.jl

@@ -47,8 +47,8 @@
             @test isequal(apply_right!(copy(l), sHadamard(q1)), apply_right_slow!(l, sHadamard(q1)))
             @test isequal(apply_right!(copy(l), sHadamardXY(q1)), apply_right_slow!(l, sHadamardXY(q1)))
             @test isequal(apply_right!(copy(l), sHadamardYZ(q1)), apply_right_slow!(l, sHadamardYZ(q1)))
-            # @test isequal(apply_right!(copy(l), sPhase(q1)), apply_right_slow!(l, sPhase(q1)))
-            # @test isequal(apply_right!(copy(l), sInvPhase(q1)), apply_right_slow!(l, sInvPhase(q1)))
+            @test isequal(apply_right!(copy(l), sPhase(q1)), apply_right_slow!(l, sPhase(q1)))
+            @test isequal(apply_right!(copy(l), sInvPhase(q1)), apply_right_slow!(l, sInvPhase(q1)))
             @test isequal(apply_right!(copy(l), sX(q1)), apply_right_slow!(l, sX(q1)))
             @test isequal(apply_right!(copy(l), sY(q1)), apply_right_slow!(l, sY(q1)))
             @test isequal(apply_right!(copy(l), sZ(q1)), apply_right_slow!(l, sZ(q1)))
@@ -58,8 +58,8 @@
             @test isequal(apply_right!(copy(l), sInvSQRTY(q1)), apply_right_slow!(l, sInvSQRTY(q1)))
             # @test isequal(apply_right!(copy(l), sSQRTZ(q1)), apply_right_slow!(l, sSQRTZ(q1)))
             # @test isequal(apply_right!(copy(l), sInvSQRTZ(q1)), apply_right_slow!(l, sInvSQRTZ(q1)))
-            # @test isequal(apply_right!(copy(l), sCXYZ(q1)), apply_right_slow!(l, sCXYZ(q1)))
-            # @test isequal(apply_right!(copy(l), sCZYX(q1)), apply_right_slow!(l, sCZYX(q1)))
+            @test isequal(apply_right!(copy(l), sCXYZ(q1)), apply_right_slow!(l, sCXYZ(q1)))
+            @test isequal(apply_right!(copy(l), sCZYX(q1)), apply_right_slow!(l, sCZYX(q1)))
             @test isequal(apply_right!(copy(l), sId1(q1)), apply_right_slow!(l, sId1(q1)))
         end
     end

test/test_inv.jl

@@ -5,7 +5,7 @@
         stabilizers = [S" I", S" X", S" Y", S" Z",
                        S"-I", S"-X", S"-Y", S"-Z",]
         for gate in subtypes(AbstractSingleQubitOperator)
-            gate == SingleQubitOperator && continue
+            gate ∈ [SingleQubitOperator, sSQRTZ, sInvSQRTZ] && continue
             for stab in stabilizers
                 @test apply_inv!(stab, gate(1)) == apply!(stab, inv(CliffordOperator(gate(1), 1)))
             end