aboutsummaryrefslogtreecommitdiff
path: root/libpsi/src/core/component.rs
diff options
context:
space:
mode:
authorhachem <im@hachem.wtf>2024-09-24 21:33:51 +0200
committerhachem <im@hachem.wtf>2024-09-24 21:33:51 +0200
commitc07ab7911bd50d0060841690a863333a9a3f226d (patch)
tree5deb12fc8234d553318e4615ae21ddfcc6db4ab4 /libpsi/src/core/component.rs
parent26010345eedc716c70e8a4219005db19b266c960 (diff)
rename project
Diffstat (limited to 'libpsi/src/core/component.rs')
-rw-r--r--libpsi/src/core/component.rs141
1 files changed, 0 insertions, 141 deletions
diff --git a/libpsi/src/core/component.rs b/libpsi/src/core/component.rs
deleted file mode 100644
index 41b7cc7..0000000
--- a/libpsi/src/core/component.rs
+++ /dev/null
@@ -1,141 +0,0 @@
-use crate::{complex, ColumnVector, Complex, Matrix, Vector, VectorMatrix};
-use core::{fmt, ops};
-
-pub type QuantumBit = ColumnVector<Complex<f64>>;
-pub type QuantumGate = Matrix<Complex<f64>>;
-
-#[macro_export]
-macro_rules! count {
- () => { 0 };
- ($head:expr $(,$tail:expr)*) => { 1 + count!($( $tail ),*) };
-}
-
-#[macro_export]
-macro_rules! qubit {
- ($(($re:expr, $im:expr)),*) => {
- {
- let mut vector = Vec::new();
- $(
- vector.push(complex!($re, $im));
- )*
- QuantumBit::new(vector)
- }
- };
-}
-
-#[macro_export]
-macro_rules! quantum_register {
- ($($bit:expr),*) => {
- {
- const N: usize = count!($($bit),*);
- let mut bits: [QuantumBit; N] = [$($bit),*];
- QuantumRegister::from(&mut bits)
- }
- };
-}
-
-pub struct ClassicalRegister {
- bits: Vec<i32>,
-}
-
-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;
-
- if alpha_norm > beta_norm {
- 0
- } else {
- 1
- }
- }
-
- pub fn state_0() -> QuantumBit {
- QuantumBit::new(vec![complex!(1.0, 0.0), complex!(0.0, 0.0)])
- }
-
- pub fn state_1() -> QuantumBit {
- QuantumBit::new(vec![complex!(0.0, 0.0), complex!(1.0, 0.0)])
- }
-}
-
-impl ClassicalRegister {
- pub fn new(count: usize) -> ClassicalRegister {
- ClassicalRegister {
- bits: Vec::with_capacity(count),
- }
- }
-
- pub fn set_bits(&mut self, bits: Vec<i32>) {
- self.bits = bits;
- }
-}
-
-impl QuantumRegister {
- fn update(&mut self) {
- let matrices: Vec<Matrix<Complex<f64>>> =
- self.qubits.iter().map(|qubit| qubit.to_matrix()).collect();
- let mut new_result = matrices[0].clone();
- for matrix in &matrices[1..] {
- new_result = new_result.kronecker(matrix);
- }
-
- self.state = ColumnVector::from_matrix(&new_result);
- }
-
- pub fn from(bits: &mut [QuantumBit]) -> QuantumRegister {
- let mut register = QuantumRegister {
- qubits: bits.to_vec(),
- state: ColumnVector::new(vec![]),
- };
-
- register.update();
- register
- }
-
- pub fn measure(&self, classical_register: &mut ClassicalRegister) {
- classical_register.set_bits(self.qubits.iter().map(|qubit| qubit.measure()).collect());
- }
-
- pub fn apply(&mut self, gate: &QuantumGate, index: usize) {
- let result: ColumnVector<Complex<f64>> = self.state.mul_matrix(gate).unwrap();
- self.qubits[index] = result;
- self.update();
- }
-}
-
-impl ops::Index<usize> for QuantumRegister {
- type Output = QuantumBit;
-
- fn index(&self, index: usize) -> &Self::Output {
- &self.qubits[index]
- }
-}
-
-impl ops::IndexMut<usize> for QuantumRegister {
- fn index_mut(&mut self, index: usize) -> &mut Self::Output {
- &mut self.qubits[index]
- }
-}
-
-impl ops::Index<usize> for ClassicalRegister {
- type Output = i32;
-
- fn index(&self, index: usize) -> &Self::Output {
- &self.bits[index]
- }
-}
-
-impl ops::IndexMut<usize> for ClassicalRegister {
- fn index_mut(&mut self, index: usize) -> &mut Self::Output {
- &mut self.bits[index]
- }
-}