quant-iron 2.0.0

A high-performance, hardware-accelerated modular quantum computing library with a focus on physical applications. Quant-Iron provides tools to represent quantum states, apply standard quantum gates, perform measurements, build quantum circuits, and implement quantum algorithms.
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
40
41
42
43
44
45

__kernel void rotate_z_kernel(
    __global float2 *state_vector,
    const int num_qubits,
    const int target_qubit,
    __global const int *control_qubits,
    const int num_control_qubits,
    const float cos_half_angle,
    const float sin_half_angle
) {
    int global_id = get_global_id(0); // Corresponds to the basis state index i

    // Check control qubits
    bool controls_met = true;
    if (num_control_qubits > 0) {
        for (int i = 0; i < num_control_qubits; ++i) {
            if (!((global_id >> control_qubits[i]) & 1)) {
                controls_met = false;
                break;
            }
        }
    }

    if (controls_met) {
        float2 amp = state_vector[global_id];
        float phase_re;
        float phase_im;

        // Determine phase based on target qubit state
        if (!((global_id >> target_qubit) & 1)) { // Target qubit is |0>
            phase_re = cos_half_angle;
            phase_im = -sin_half_angle; // exp(-i * angle / 2)
        } else { // Target qubit is |1>
            phase_re = cos_half_angle;
            phase_im = sin_half_angle;  // exp(i * angle / 2)
        }

        // Apply phase: (a+bi) * (phase_re + i*phase_im)
        // = (a*phase_re - b*phase_im) + i*(a*phase_im + b*phase_re)
        float original_real = amp.x;
        float original_imag = amp.y;
        amp.x = original_real * phase_re - original_imag * phase_im;
        amp.y = original_real * phase_im + original_imag * phase_re;
        state_vector[global_id] = amp;
    }
}