Source file `gate.ml`

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
(*
 * Copyright 2025 Elias GAUTHIER <elias.gauthier@etu.u-bordeaux.fr>
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *)
open Qubit

let x q =
  let new_alpha = Complex.cadd
    (Complex.cmul Complex.zero q.alpha)
    (Complex.cmul Complex.one q.beta) in
  let new_beta = Complex.cadd
    (Complex.cmul Complex.one q.alpha)
    (Complex.cmul Complex.zero q.beta) in
  q.alpha <- new_alpha;
  q.beta <- new_beta

let y q =
  let i = { Complex.re = 0.0; im = 1.0 } in
  let minus_i = { Complex.re = 0.0; im = -1.0 } in
  let new_alpha = Complex.cmul minus_i q.beta in
  let new_beta = Complex.cmul i q.alpha in
  q.alpha <- new_alpha;
  q.beta <- new_beta

let z q =
  let new_alpha = Complex.cadd
    (Complex.cmul Complex.one q.alpha)
    (Complex.cmul Complex.zero q.beta) in
  let new_beta = Complex.cadd
    (Complex.cmul Complex.zero q.alpha)
    (Complex.cmul Complex.minus_one q.beta) in
  q.alpha <- new_alpha;
  q.beta <- new_beta

let h q =
  let sqrt2_inv = 1.0 /. sqrt 2.0 in
  let factor = { Complex.re = sqrt2_inv; im = 0.0 } in
  let new_alpha = Complex.cmul factor (Complex.cadd q.alpha q.beta) in
  let new_beta = Complex.cmul factor (Complex.csub q.alpha q.beta) in
  q.alpha <- new_alpha;
  q.beta <- new_beta