Add FractionalISwap gate

bsc-quantic · Mar 18, 2024 · 24233b9 · 24233b9
1 parent 0bde6bd
commit 24233b9
Show file tree

Hide file tree

Showing 2 changed files with 22 additions and 0 deletions.
diff --git a/src/Array.jl b/src/Array.jl
@@ -109,6 +109,15 @@ Matrix{T}(::Swap) where {T} = Matrix{T}([1 0 0 0; 0 0 1 0; 0 1 0 0; 0 0 0 1])
 
 Matrix{T}(::ISwap) where {T} = Matrix{T}([1 0 0 0; 0 0 1im 0; 0 1im 0 0; 0 0 0 1])
 
+Matrix{T}(op::FractionalISwap) where {T} = Matrix{T}(
+    [
+        1 0 0 0
+        0 cospi(op.exponent / 2) 1im*sinpi(op.exponent / 2)*cispi(2 * op.phase_exponent) 0
+        0 1im*sinpi(op.exponent / 2)*cispi(-2 * op.phase_exponent) cospi(op.exponent / 2) 0
+        0 0 0 1
+    ],
+)
+
 Matrix{T}(::FSwap) where {T} = Matrix{T}([1 0 0 0; 0 0 1 0; 0 1 0 0; 0 0 0 -1])
 
 Matrix{T}(op::FSim) where {T} = Matrix{T}([
@@ -165,6 +174,13 @@ Array{T,4}(::Swap) where {T} = Array{T}([1; 0;; 0; 0;;; 0; 0;; 1; 0;;;; 0; 1;; 0
 
 Array{T,4}(::ISwap) where {T} = Array{T}([1; 0;; 0; 0;;; 0; 0;; 1im; 0;;;; 0; 1im;; 0; 0;;; 0; 0;; 0; 1])
 
+Array{T,4}(op::FractionalISwap) where {T} = Array{T}(
+    [
+        1; 0;; 0; 0;;; 0; cospi(op.exponent / 2);; 1im*sinpi(op.exponent / 2)*cispi(-2 * op.phase_exponent); 0;;;;
+        0; 1im*sinpi(op.exponent / 2)*cispi(-2 * op.phase_exponent);; cospi(op.exponent / 2); 0;;; 0; 0;; 0; 1
+    ],
+)
+
 Array{T,4}(::Swap) where {T} = Array{T}([1; 0;; 0; 0;;; 0; 0;; 1; 0;;;; 0; 1;; 0; 0;;; 0; 0;; 0; -1])
 
 Array{T,4}(op::FSim) where {T} = Array{T}(

diff --git a/src/Gate.jl b/src/Gate.jl
@@ -334,6 +334,12 @@ The ISWAP gate.
 """
 @gatedecl ISwap n = 2
 
+# TODO move `exponent`, `phase_exponent` to `ISwap` operator?
+@gatedecl FractionalISwap n = 2 begin
+    exponent::Float64 = 1
+    phase_exponent::Float64 = 0
+end
+
 """
     FSwap(lane1, lane2)