CrysQuant::GateMatrix

Base_00 = State_0.⊗(State_0.⊤)

Base_01 = State_0.⊗(State_1.⊤)

Base_10 = State_1.⊗(State_0.⊤)

Base_11 = State_1.⊗(State_1.⊤)

Bases = [[Base_00, Base_01], [Base_10, Base_11]]

CCNOT = (Base_00.⊗(Unity4)) + (Base_11.⊗(CNOT))

CNOT = (Base_00.⊗(Unity)) + (Base_11.⊗(Pauli_X))

CSWAP = (Base_00.⊗(Unity4)) + (Base_11.⊗(SWAP))

Hadamard = Matrix[[1, 1], [1, -1]] / (Math.sqrt(2))

Helper_A = Matrix[[1, 0], [0, Complex.new(1 / 2, 1 / 2)]]

Helper_B = Matrix[[Complex.new(1 / 2, 1 / 2), 0], [0, 1]]

Pauli_X = Matrix[[0, 1], [1, 0]]

Pauli_Y = Matrix[[0, -I], [I, 0]]

Pauli_Z = Matrix[[1, 0], [0, -1]]

Rotation_Phase_45 = Matrix[[1, 0], [0, (Complex.new(1, 1)) / (Math.sqrt(2))]]

Rotation_Phase_90 = Matrix[[1, 0], [0, I]]

Rotation_X_90 = Matrix[[1, -I], [-I, 1]] / (Math.sqrt(2))

Rotation_Y_90 = Matrix[[1, -1], [1, 1]] / (Math.sqrt(2))

Sqrt_NOT = Matrix[[Complex.new(1, 1), Complex.new(1, -1)], [Complex.new(1, -1), Complex.new(1, 1)]] / 2

Sqrt_SWAP = ((Helper_B.⊗(Base_11)) + (Helper_A.⊗(Base_00))) + (((Base_01.⊗(Base_10)) + (Base_10.⊗(Base_01))) * (Complex.new(1 / 2, -1 / 2)))

SWAP = (((Base_00.⊗(Base_00)) + (Base_01.⊗(Base_10))) + (Base_10.⊗(Base_01))) + (Base_11.⊗(Base_11))

Unity = Matrix.identity(2)

Unity4 = Matrix.identity(4)

Zero = Matrix[[0, 0], [0, 0]]