Tests and Bug fixes

src/PauliOperators.jl

@@ -17,7 +17,7 @@ include("function_base.jl")
 include("function_convert.jl")
 include("function_commute.jl")
 include("function_otimes.jl")
-
+include("function_osum.jl")
 
 # Exports
 export Pauli

src/function_mul.jl

@@ -204,7 +204,7 @@
 Multiply two `Pauli`'s together
 """
 function Base.:*(p1::Pauli{N}, p2::Pauli{N}) where {N}
-    θ = (p1.θ + p2.θ + phase(p1.pauli) + phase(p2.pauli) + phase(p1.pauli, p2.pauli)) % 4
+    θ = (p1.θ + p2.θ + phase(p1.pauli, p2.pauli)) % 4
     return Pauli{N}(θ, p1.pauli*p2.pauli)
 end
 
@@ -215,7 +215,7 @@
     # θ = (p1.θ + phase(p2)) % 4
     # θ += (2*count_ones(p1.x & p2.z)) % 4
     # return Pauli{N}(θ, z, x)
-    θ = (p1.θ + phase(p1.pauli) + phase(p2) + phase(p1.pauli, p2)) % 4
+    θ = (p1.θ + phase(p1.pauli, p2)) % 4
     return Pauli{N}(θ, p1.pauli*p2)
 end
 
@@ -225,7 +225,7 @@
     # θ = (phase(p1) + p2.θ) % 4
     # θ += (2*count_ones(p1.x & p2.z)) % 4
     # return Pauli{N}(θ, z, x)
-    θ = (p2.θ + phase(p1) + phase(p2.pauli) + phase(p1, p2.pauli)) % 4
+    θ = (p2.θ + phase(p1, p2.pauli)) % 4
     return Pauli{N}(θ,p1*p2.pauli)
 end
 
@@ -277,4 +277,4 @@
     tmp = p.pauli.x ⊻ ψ.v
     sign = count_ones(p.pauli.z & tmp) % 2
     return p.coeff*(-1)^sign, KetBitString{N}(tmp)
-end
+end

src/function_otimes.jl

@@ -1,7 +1,7 @@
 otimes(p1::Pauli{N}, p2::Pauli{M}) where {N,M} = Pauli{N+M}((p1.θ + p2.θ)%4, p1.pauli ⊗ p2.pauli)
 otimes(p1::FixedPhasePauli{N}, p2::FixedPhasePauli{M}) where {N,M} = FixedPhasePauli{N+M}(p1.z | p2.z << N, p1.x | p2.x << N)
-otimes(p::Pauli{N}, fpp::FixedPhasePauli{M}) where {N,M} = otimes(p, Pauli{M}(0,fpp.z, fpp.x))
-otimes(fpp::FixedPhasePauli{M}, p::Pauli{N}) where {N,M} = otimes(Pauli{M}(0,fpp.z, fpp.x), p)
+otimes(p::Pauli{N}, fpp::FixedPhasePauli{M}) where {N,M} = otimes(p, Pauli{M}(0,fpp))
+otimes(fpp::FixedPhasePauli{M}, p::Pauli{N}) where {N,M} = otimes(Pauli{M}(0,fpp), p)
 
 function otimes(ps::PauliSum{N}, p::Pauli{M}) where {N,M}
     out = PauliSum(N+M)
@@ -27,4 +27,4 @@
     end
     return out 
 end
-const ⊗ = otimes
+const ⊗ = otimes

src/operations.jl

@@ -138,7 +138,7 @@
 Check if operator is diagonal in the computational (z) basis. E.g., does this operator consist of only I and/or Z?
 """
 function is_diagonal(p::Pauli)
-    return count_ones(p.x) == 0
+    return count_ones(p.pauli.x) == 0
 end
 
 
@@ -158,6 +158,6 @@
 
     println(p)
 
-    count_ones(p.z & ket.v) % 2 == 0 || return -(1im)^p.θ
+    count_ones(p.pauli.z & ket.v) % 2 == 0 || return -(1im)^p.θ
     return (1im)^p.θ 
 end

test/tests.jl

@@ -91,6 +91,20 @@ using Random
 
 
     println("Test Multiply")
+    for i in 1:10
+        N = 8
+        a = random_Pauli(N)
+        b = random_Pauli(N)
+
+        display(a * b)
+        # println()
+        # display(s2)
+
+        @test all(Matrix(a) * Matrix(b) - Matrix(a * b) .≈ 0)
+    end
+
+
+    println("Test Addition")
     for i in 1:10
         N = 8
         a = random_Pauli(N)
@@ -100,14 +114,12 @@ using Random
 
         s1 = a + b + d
         s2 = a + c + d
-        display(s1 * s2)
-        # println()
-        # display(s2)
 
-        @test all(Matrix(s1) * Matrix(s2) - Matrix(s1 * s2) .≈ 0)
+        @test all(Matrix(s1) + Matrix(s2) - Matrix(s1 + s2) .≈ 0)
     end
+
 
-
+    
     a = random_Pauli(6)
     b = random_Pauli(6)
     s = a + b