QVNT

Advanced quantum computation simulator, written in Rust

Features

Ability to simulate up to 64 qubits. Common machine with 4-16 Gb of RAM is able to simulate 26-28 qubits, which is enough for several study cases;
Set of 1- or 2-qubits operations to build your own quantum circuits;
Quantum operations are tested and debugged to be safe in use;
Circuit execution is accelerated using multithreading Rayon library;
Complex quantum registers manipulations: tensor product of two registers and aliases for qubit to simplify interaction with register.

Usage

Add this lines to your Cargo.toml file to use QVNT crate:

[dependencies]
qvnt = { version = "0.4.4", features = ["multi-thread"] }

Quantum register and operators are controlled by bitmasks. Each bit in it will act on a specific qubit.

use qvnt::prelude::*;

//  Create quantum register with 10 qubits
let mut q_reg = QReg::new(10);
//  or with initial state, where 5th, 6th and 7th qubits are already in state |1>.
let mut q_reg = QReg::with_state(10, 0b0011100000);

//  Create qft (Quantum Fourier Transform) operation, acting on first 5 qubits in q_reg.
let op = op::qft(0b0000011111);

//  Apply created operation
q_reg.apply(&op);

//  Measure and write first 3 qubit, which leads to collapse of q_reg wave function.
//  Measured variable will contain one of the following values:
//  0b000, 0b001, 0b010, 0b011, 0b100, 0b101, 0b110, 0b111
let measured = q_reg.measure_mask(0b0000000111);

You're able to use VReg to simplify operations definition:

use qvnt::prelude::*;

let mut q_reg = QReg::new(10);
let q = q_reg.get_vreg();

//  Crate Hadamard operator, that act on odd qubits.
let op = op::h(q[1] | q[3] | q[5] | q[7] | q[9]);
//  This is equivalent to op::h(0b0101010101);

Implemented operations

Pauli's X, Y & Z operators;
Square and fourth root of Z - S & T operators;
Phase shift operator - phi;
1-qubit rotation operators - rx, ry & rz;
2-qubits rotation operators, aka Ising coupling gates, - rxx, ryy & rzz;
SWAP, iSWAP operators and square rooted ones;
Quantum Fourier and Hadamard Transform;
Universal U1, U2 and U3 operators;

ALL operators have inverse versions, accessing by .dgr() method:

use qvnt::prelude::*;

let usual_op = op::s(0b1);
//  Inverse S operator
let inverse_op = op::s(0b1).dgr();

Also, ALL these operators could be turned into controlled ones, using .c(...) method:

use qvnt::prelude::*;

let usual_op = op::x(0b001);
//  NOT gate, controlled by 2 qubits, aka CCNOT gate, aka Toffoli gate
let controlled_op = op::x(0b001).c(0b110).unwrap();

Controlled operation has to be unwrapped, since it could be None if its mask overlaps with the mask of operator. For example, this code will panic:

use qvnt::prelude::*;
let _ = op::x(0b001).c(0b001).unwrap();

License