circus 0.1.1

Quantum circuit simulator
Documentation 
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

use super::Gate;
use crate::{State, PW};

#[derive(Clone, Copy, Debug)]
pub struct CNotGate {
    pub target: usize,
    pub control: usize,
}

impl Gate for CNotGate {
    fn apply(&self, state: &mut State) {
        let b5 = self.target >> 5;
        let c5 = self.control >> 5;
        let pwb = PW[self.target & 31];
        let pwc = PW[self.control & 31];
        for i in 0..2 * state.n {
            if state.x[i][b5] & pwb > 0 {
                state.x[i][c5] ^= pwc;
            }
            if state.z[i][c5] & pwc > 0 {
                state.z[i][b5] ^= pwb;
            }
            if (state.x[i][b5] & pwb > 0)
                && (state.z[i][c5] & pwc > 0)
                && (state.x[i][c5] & pwc > 0)
                && (state.z[i][b5] & pwb > 0)
            {
                state.r[i] = (state.r[i] + 2) % 4;
            }
            if (state.x[i][b5] & pwb > 0)
                && (state.z[i][c5] & pwc > 0)
                && !(state.x[i][c5] & pwc > 0)
                && !(state.z[i][b5] & pwb > 0)
            {
                state.r[i] = (state.r[i] + 2) % 4;
            }
        }
    }
}