quantrs2-core 0.2.0

Core types and traits for the QuantRS2 quantum computing framework
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134

//! SIMD Performance Benchmarks
//!
//! This benchmark suite measures the performance of SIMD-accelerated
//! quantum gate operations.

use criterion::{black_box, criterion_group, criterion_main, BenchmarkId, Criterion};
use quantrs2_core::simd_enhanced::{RotationAxis, SimdGateEngine};
use scirs2_core::Complex64;

fn benchmark_rx_gates(c: &mut Criterion) {
    let engine = SimdGateEngine::new();

    let mut group = c.benchmark_group("rx_gate");

    for num_qubits in &[4, 6, 8, 10, 12] {
        let size = 1 << num_qubits;
        let mut state: Vec<Complex64> = vec![Complex64::new(0.0, 0.0); size];
        state[0] = Complex64::new(1.0, 0.0);

        group.bench_with_input(
            BenchmarkId::from_parameter(num_qubits),
            num_qubits,
            |b, &num_qubits| {
                b.iter(|| {
                    let mut state_copy = state.clone();
                    engine
                        .apply_rotation_gate(
                            black_box(&mut state_copy),
                            black_box(0),
                            black_box(RotationAxis::X),
                            black_box(std::f64::consts::PI / 4.0),
                        )
                        .unwrap();
                });
            },
        );
    }

    group.finish();
}

fn benchmark_cnot_gates(c: &mut Criterion) {
    let engine = SimdGateEngine::new();

    let mut group = c.benchmark_group("cnot_gate");

    for num_qubits in &[4, 6, 8, 10, 12] {
        let size = 1 << num_qubits;
        let mut state: Vec<Complex64> = vec![Complex64::new(0.0, 0.0); size];
        state[0] = Complex64::new(1.0, 0.0);

        group.bench_with_input(
            BenchmarkId::from_parameter(num_qubits),
            num_qubits,
            |b, &num_qubits| {
                b.iter(|| {
                    let mut state_copy = state.clone();
                    engine
                        .apply_cnot(black_box(&mut state_copy), black_box(0), black_box(1))
                        .unwrap();
                });
            },
        );
    }

    group.finish();
}

fn benchmark_batch_gates(c: &mut Criterion) {
    let engine = SimdGateEngine::new();

    let mut group = c.benchmark_group("batch_gates");

    for num_qubits in &[4, 6, 8] {
        let size = 1 << num_qubits;
        let mut state: Vec<Complex64> = vec![Complex64::new(0.0, 0.0); size];
        state[0] = Complex64::new(1.0, 0.0);

        // Create a batch of gates
        let gates: Vec<_> = (0..*num_qubits)
            .map(|i| (i, RotationAxis::X, std::f64::consts::PI / 4.0))
            .collect();

        group.bench_with_input(
            BenchmarkId::from_parameter(num_qubits),
            num_qubits,
            |b, _| {
                b.iter(|| {
                    let mut state_copy = state.clone();
                    engine
                        .batch_apply_single_qubit(black_box(&mut state_copy), black_box(&gates))
                        .unwrap();
                });
            },
        );
    }

    group.finish();
}

fn benchmark_fidelity(c: &mut Criterion) {
    let engine = SimdGateEngine::new();

    let mut group = c.benchmark_group("fidelity");

    for num_qubits in &[4, 6, 8, 10, 12] {
        let size = 1 << num_qubits;
        let state1: Vec<Complex64> = vec![Complex64::new(1.0 / (size as f64).sqrt(), 0.0); size];
        let state2: Vec<Complex64> = vec![Complex64::new(1.0 / (size as f64).sqrt(), 0.0); size];

        group.bench_with_input(
            BenchmarkId::from_parameter(num_qubits),
            num_qubits,
            |b, _| {
                b.iter(|| {
                    engine
                        .fidelity(black_box(&state1), black_box(&state2))
                        .unwrap()
                });
            },
        );
    }

    group.finish();
}

criterion_group!(
    benches,
    benchmark_rx_gates,
    benchmark_cnot_gates,
    benchmark_batch_gates,
    benchmark_fidelity
);
criterion_main!(benches);