aboutsummaryrefslogtreecommitdiff
path: root/libpsi-core/src/core
diff options
context:
space:
mode:
Diffstat (limited to 'libpsi-core/src/core')
-rw-r--r--libpsi-core/src/core/circuit.rs25
-rw-r--r--libpsi-core/src/core/component.rs68
-rw-r--r--libpsi-core/src/core/gates.rs41
-rw-r--r--libpsi-core/src/core/runtime.rs19
4 files changed, 76 insertions, 77 deletions
diff --git a/libpsi-core/src/core/circuit.rs b/libpsi-core/src/core/circuit.rs
index 635066a..4d9d286 100644
--- a/libpsi-core/src/core/circuit.rs
+++ b/libpsi-core/src/core/circuit.rs
@@ -1,20 +1,19 @@
use super::{ClassicalRegister, QuantumGate, QuantumRegister};
#[derive(Clone)]
-pub struct QuantumCircuitInstruction {
- pub target_bit: usize,
- pub control_indices: Vec<usize>,
- pub gate: QuantumGate,
+pub enum Instruction<'a> {
+ ApplyGate(&'a QuantumGate, &'a [usize], usize),
+ Measure(usize),
}
-pub struct QuantumCircuit {
+pub struct QuantumCircuit<'a> {
quantum_register: QuantumRegister,
classical_register: ClassicalRegister,
- instructions: Vec<QuantumCircuitInstruction>,
+ instructions: Vec<Instruction<'a>>,
}
-impl QuantumCircuit {
- pub fn new(quantum_bit_count: usize, classical_bit_count: usize) -> QuantumCircuit {
+impl<'a> QuantumCircuit<'a> {
+ pub fn new(quantum_bit_count: usize, classical_bit_count: usize) -> QuantumCircuit<'a> {
QuantumCircuit {
quantum_register: QuantumRegister::new(quantum_bit_count),
classical_register: ClassicalRegister::new(classical_bit_count),
@@ -22,12 +21,8 @@ impl QuantumCircuit {
}
}
- pub fn apply(&mut self, gate: &QuantumGate, control_bits: &[usize], target_bit: usize) {
- self.instructions.push(QuantumCircuitInstruction {
- gate: gate.clone(),
- target_bit,
- control_indices: control_bits.to_vec(),
- });
+ pub fn execute(&mut self, instruction: Instruction<'a>) {
+ self.instructions.push(instruction);
}
pub fn get_quantum_register(&self) -> QuantumRegister {
@@ -38,7 +33,7 @@ impl QuantumCircuit {
self.classical_register.clone()
}
- pub fn get_instructions(&self) -> Vec<QuantumCircuitInstruction> {
+ pub fn get_instructions(&self) -> Vec<Instruction> {
self.instructions.clone()
}
}
diff --git a/libpsi-core/src/core/component.rs b/libpsi-core/src/core/component.rs
index d2dc781..7748ac5 100644
--- a/libpsi-core/src/core/component.rs
+++ b/libpsi-core/src/core/component.rs
@@ -1,11 +1,25 @@
-use rand::Rng;
-
use crate::{complex, ColumnVector, Complex, Matrix, Vector, VectorMatrix};
-use core::ops;
+use core::{fmt, ops};
+
+pub type QuantumState = ColumnVector<Complex<f64>>;
+pub type QuantumBit = QuantumState;
+
+#[derive(Clone)]
+pub struct QuantumGate {
+ pub name: &'static str,
+ pub matrix: Matrix<Complex<f64>>,
+}
-pub type QuantumBit = ColumnVector<Complex<f64>>;
-pub type QuantumState = Matrix<Complex<f64>>;
-pub type QuantumGate = (&'static str, QuantumState);
+#[derive(Clone)]
+pub struct ClassicalRegister {
+ bits: Vec<i32>,
+}
+
+#[derive(Clone)]
+pub struct QuantumRegister {
+ state_vector: QuantumState,
+ qubits: Vec<QuantumBit>,
+}
#[macro_export]
macro_rules! count {
@@ -37,35 +51,7 @@ macro_rules! quantum_register {
};
}
-#[derive(Clone)]
-pub struct ClassicalRegister {
- bits: Vec<i32>,
-}
-
-#[derive(Clone)]
-pub struct QuantumRegister {
- state: ColumnVector<Complex<f64>>,
- qubits: Vec<QuantumBit>,
-}
-
impl QuantumBit {
- pub fn measure(&self) -> i32 {
- let alpha_abs = self[0].abs();
- let beta_abs = self[1].abs();
-
- let alpha_norm = alpha_abs * alpha_abs;
- let beta_norm = beta_abs * beta_abs;
-
- let mut rng = rand::thread_rng();
- let random_value = rng.gen_range(0.0..(alpha_norm as f32 + beta_norm as f32));
-
- if random_value < alpha_norm as f32 {
- 0
- } else {
- 1
- }
- }
-
pub fn state_0() -> QuantumBit {
QuantumBit::new(vec![complex!(1.0, 0.0), complex!(0.0, 0.0)])
}
@@ -75,6 +61,12 @@ impl QuantumBit {
}
}
+impl fmt::Display for QuantumGate {
+ fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
+ write!(f, "{}", self.name)
+ }
+}
+
impl ClassicalRegister {
pub fn new(count: usize) -> ClassicalRegister {
ClassicalRegister {
@@ -100,7 +92,7 @@ impl QuantumRegister {
new_result = new_result.kronecker(matrix);
}
- self.state = ColumnVector::from_matrix(&new_result);
+ self.state_vector = ColumnVector::from_matrix(&new_result);
}
pub fn new(count: usize) -> QuantumRegister {
@@ -110,7 +102,7 @@ impl QuantumRegister {
pub fn from(bits: &mut [QuantumBit]) -> QuantumRegister {
let mut register = QuantumRegister {
qubits: bits.to_vec(),
- state: ColumnVector::new(vec![]),
+ state_vector: ColumnVector::new(vec![]),
};
register.update();
@@ -121,8 +113,8 @@ impl QuantumRegister {
self.qubits.clone()
}
- pub fn get_state(&self) -> ColumnVector<Complex<f64>> {
- self.state.clone()
+ pub fn get_state(&self) -> QuantumState {
+ self.state_vector.clone()
}
}
diff --git a/libpsi-core/src/core/gates.rs b/libpsi-core/src/core/gates.rs
index 094fec9..97254fc 100644
--- a/libpsi-core/src/core/gates.rs
+++ b/libpsi-core/src/core/gates.rs
@@ -2,21 +2,36 @@ use crate::{complex, matrix, QuantumGate};
#[rustfmt::skip]
lazy_static::lazy_static! {
- pub static ref HADAMARD: QuantumGate = ("H", matrix!([complex!(1.0, 0.0), complex!( 1.0, 0.0)];
- [complex!(1.0, 0.0), complex!(-1.0, 0.0)]) *
- complex!(1.0/2.0_f64.sqrt(), 0.0));
+ pub static ref HADAMARD: QuantumGate = QuantumGate {
+ name: "H",
+ matrix: matrix!([complex!(1.0, 0.0), complex!( 1.0, 0.0)];
+ [complex!(1.0, 0.0), complex!(-1.0, 0.0)]) *
+ complex!(1.0/2.0_f64.sqrt(), 0.0)
+ };
- pub static ref PAULI_X: QuantumGate = ("Pauli-X", matrix!([complex!(0.0, 0.0), complex!(1.0, 0.0)];
- [complex!(1.0, 0.0), complex!(0.0, 0.0)]));
+ pub static ref PAULI_X: QuantumGate = QuantumGate {
+ name: "Pauli-X",
+ matrix: matrix!([complex!(0.0, 0.0), complex!(1.0, 0.0)];
+ [complex!(1.0, 0.0), complex!(0.0, 0.0)])
+ };
- pub static ref PAULI_Y: QuantumGate = ("Pauli-Y", matrix!([complex!(0.0, 0.0), complex!(0.0, -1.0)];
- [complex!(0.0, 1.0), complex!(0.0, 0.0)]));
+ pub static ref PAULI_Y: QuantumGate = QuantumGate {
+ name: "Pauli-Y",
+ matrix: matrix!([complex!(0.0, 0.0), complex!(0.0, -1.0)];
+ [complex!(0.0, 1.0), complex!(0.0, 0.0)])
+ };
- pub static ref PAULI_Z: QuantumGate = ("Pauli-Z", matrix!([complex!(1.0, 0.0), complex!( 0.0, 0.0)];
- [complex!(0.0, 0.0), complex!(-1.0, 0.0)]));
+ pub static ref PAULI_Z: QuantumGate = QuantumGate {
+ name: "Pauli-Z",
+ matrix: matrix!([complex!(1.0, 0.0), complex!( 0.0, 0.0)];
+ [complex!(0.0, 0.0), complex!(-1.0, 0.0)])
+ };
- pub static ref CNOT: QuantumGate = ("CNOT", matrix!([complex!(1.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0)];
- [complex!(0.0, 0.0), complex!(1.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0)];
- [complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(1.0, 0.0)];
- [complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(1.0, 0.0), complex!(0.0, 0.0)]));
+ pub static ref CNOT: QuantumGate = QuantumGate {
+ name: "CNOT",
+ matrix: matrix!([complex!(1.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0)];
+ [complex!(0.0, 0.0), complex!(1.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0)];
+ [complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(1.0, 0.0)];
+ [complex!(0.0, 0.0), complex!(0.0, 0.0), complex!(1.0, 0.0), complex!(0.0, 0.0)])
+ };
}
diff --git a/libpsi-core/src/core/runtime.rs b/libpsi-core/src/core/runtime.rs
index 19ac191..f3e52ac 100644
--- a/libpsi-core/src/core/runtime.rs
+++ b/libpsi-core/src/core/runtime.rs
@@ -1,25 +1,22 @@
-use super::{ClassicalRegister, QuantumCircuit, QuantumRegister};
+use super::{ClassicalRegister, QuantumCircuit, QuantumRegister, QuantumState};
-#[allow(unused)]
-pub struct ExecutionData {
+pub struct RuntimeData {
pub quantum_register: QuantumRegister,
pub classical_states: ClassicalRegister,
}
-impl ExecutionData {
+impl RuntimeData {
pub fn new(quantum_register: QuantumRegister, classical_states: ClassicalRegister) -> Self {
- ExecutionData {
+ RuntimeData {
quantum_register,
classical_states,
}
}
}
-pub trait RuntimeBase {
- fn get_circuit(&self) -> &QuantumCircuit;
-}
-
-pub trait Runtime<'a, T: RuntimeBase = Self> {
+pub trait Runtime<'a> {
fn new(circuit: &'a QuantumCircuit) -> Self;
- fn execute(&self, repeat: usize) -> Vec<ExecutionData>;
+ fn execute(&self, repeat: usize) -> Vec<RuntimeData>;
+ fn get_state(&self) -> QuantumState;
+ fn get_circuit(&self) -> &QuantumCircuit;
}