Quipper.Circuit

Low-level quantum circuit implementation. This is our backend implementation of quantum circuits. Note: there is no run-time error checking at the moment.

At its heart, a circuit is a list of gates. All well-definedness checking (e.g. input arity, output arity, and checking that the intermediate gates are connected to legitimate wires) is done dynamically, at circuit generation time, and is not stored within the circuit itself. This allows circuits to be produced and consumed lazily.

Implementation note: this file is in the intermediate stage of a code refactoring, and should be considered "under renovation".

Quantum circuit data type

type Wire = Int

Wire identifier. Wires are currently identified by an integer, but the users of this interface should be oblivious to this.

data Wiretype

Wire type. A wire is either quantum or classical.

Constructors

Qbit

Quantum wire.

Cbit

Classical wire.

Instances

Eq Wiretype

Methods

(==) :: Wiretype -> Wiretype -> Bool

(/=) :: Wiretype -> Wiretype -> Bool

Show Wiretype

Methods

showsPrec :: Int -> Wiretype -> ShowS

show :: Wiretype -> String

showList :: [Wiretype] -> ShowS

type Arity = IntMap Wiretype

An arity, also known as a typing context, is a map from a finite set of wires to wire types.

data Signed a

A signed item of type a. Signed x True represents a positive item, and Signed x False represents a negative item.

When used with wires in a circuit, a positive sign is used to represent a positive control, i.e., a filled dot, and a negative sign is used to represent a negative control, i.e., an empty dot.

Constructors

Signed a Bool

Instances

Show a => Show (Signed a)

Methods

showsPrec :: Int -> Signed a -> ShowS

show :: Signed a -> String

showList :: [Signed a] -> ShowS

ControlSource (Signed Wire)

Methods

to_control :: Signed Wire -> ControlList

ControlSource (Signed Bit)

Methods

to_control :: Signed Bit -> ControlList

ControlSource (Signed Qubit)

Methods

to_control :: Signed Qubit -> ControlList

QCData a => QCData (Signed a)

Methods

qcdata_mapM :: Monad m => Signed a -> (q -> m q') -> (c -> m c') -> QCType q c (Signed a) -> m (QCType q' c' (Signed a))

qcdata_zip :: Signed a -> q -> c -> q' -> c' -> QCType q c (Signed a) -> QCType q' c' (Signed a) -> ErrMsg -> QCType (q, q') (c, c') (Signed a)

qcdata_promote :: BType (Signed a) -> Signed a -> ErrMsg -> BType (Signed a)

Labelable a String => Labelable (Signed a) String

Methods

label_rec :: Signed a -> String -> LabelMonad ()

Labelable a String => Labelable (Signed a) (Signed String)

Methods

label_rec :: Signed a -> Signed String -> LabelMonad ()

type QCType x y (Signed a)

type QCType x y (Signed a) = Signed (QCType x y a)

type QTypeB (Signed a)

type QTypeB (Signed a) = Signed (QTypeB a)

from_signed :: Signed a -> a

Extract the underlying item of a signed item.

get_sign :: Signed a -> Bool

Extract the sign of a signed item: True is positive, and False is negative.

type Controls = [Signed Wire]

A list of controlled wires, possibly empty.

type Timestep = Double

A time step is a small floating point number used as a parameter to certain gates, such as rotation gates or the [exp −iZt] gate.

type InverseFlag = Bool

A flag that, if True, indicates that the gate is inverted.

type NoControlFlag = Bool

A flag that, if True, indicates that the gate is controllable, but any further controls on the gate should be ignored. This is used, e.g., for circuits consisting of a basis change, some operation, and the inverse basis change. When controlling such a circuit, it is sufficient to control the middle operation, so the gates belonging to the basis change and its inverse will have the NoControlFlag set.

data ControllableFlag

A flag, to specify if the corresponding subroutine can be controlled. Either no control allowed, or all controls, or only classical.

Constructors

NoCtl

AllCtl

OnlyClassicalCtl

Instances

Eq ControllableFlag

Methods

(==) :: ControllableFlag -> ControllableFlag -> Bool

(/=) :: ControllableFlag -> ControllableFlag -> Bool

Ord ControllableFlag

Methods

compare :: ControllableFlag -> ControllableFlag -> Ordering

(<) :: ControllableFlag -> ControllableFlag -> Bool

(<=) :: ControllableFlag -> ControllableFlag -> Bool

(>) :: ControllableFlag -> ControllableFlag -> Bool

(>=) :: ControllableFlag -> ControllableFlag -> Bool

max :: ControllableFlag -> ControllableFlag -> ControllableFlag

min :: ControllableFlag -> ControllableFlag -> ControllableFlag

Show ControllableFlag

Methods

showsPrec :: Int -> ControllableFlag -> ShowS

show :: ControllableFlag -> String

showList :: [ControllableFlag] -> ShowS

data BoxId

An identifier for a subroutine. A boxed subroutine is currently identified by a pair of: the user-defined name of the subroutine; and a value uniquely identifying the type and shape of the argument.

For now, we represent the shape as a string, because this gives an easy total Ord instance, needed for Data.Map. However, in principle, one could also use a pair of a type representation and a shape term. The implementation of this may change later.

Constructors

BoxId String String

Instances

Eq BoxId

Methods

(==) :: BoxId -> BoxId -> Bool

(/=) :: BoxId -> BoxId -> Bool

Ord BoxId

Methods

compare :: BoxId -> BoxId -> Ordering

(<) :: BoxId -> BoxId -> Bool

(<=) :: BoxId -> BoxId -> Bool

(>) :: BoxId -> BoxId -> Bool

(>=) :: BoxId -> BoxId -> Bool

max :: BoxId -> BoxId -> BoxId

min :: BoxId -> BoxId -> BoxId

Show BoxId

Methods

showsPrec :: Int -> BoxId -> ShowS

show :: BoxId -> String

showList :: [BoxId] -> ShowS

data RepeatFlag

A flag that indicates how many times a particular subroutine should be repeated. If non-zero, it implies some constraints on the type of the subroutine.

Constructors

RepeatFlag Integer

Instances

Eq RepeatFlag

Methods

(==) :: RepeatFlag -> RepeatFlag -> Bool

(/=) :: RepeatFlag -> RepeatFlag -> Bool

Ord RepeatFlag

Methods

compare :: RepeatFlag -> RepeatFlag -> Ordering

(<) :: RepeatFlag -> RepeatFlag -> Bool

(<=) :: RepeatFlag -> RepeatFlag -> Bool

(>) :: RepeatFlag -> RepeatFlag -> Bool

(>=) :: RepeatFlag -> RepeatFlag -> Bool

max :: RepeatFlag -> RepeatFlag -> RepeatFlag

min :: RepeatFlag -> RepeatFlag -> RepeatFlag

Show RepeatFlag

Methods

showsPrec :: Int -> RepeatFlag -> ShowS

show :: RepeatFlag -> String

showList :: [RepeatFlag] -> ShowS

data Gate

The low-level representation of gates.

Constructors

QGate String InverseFlag [Wire] [Wire] Controls NoControlFlag

A named reversible quantum gate: Qbit^(m+n) -> Qbit^(m+n). The second [Wire] argument should be "generalized controls", i.e. wires not modified by the gate. The gate type is uniquely determined by: the name, the number of inputs, and the number of generalized controls. Gates that differ in one of these respects should be regarded as different gates.

QRot String InverseFlag Timestep [Wire] [Wire] Controls NoControlFlag

A named reversible quantum gate that also depends on a real parameter. This is typically used for phase and rotation gates. The gate name can contain '%' as a place holder for the parameter, e.g., "exp(-i%Z)". The remaining arguments are as for QGate.

GPhase Timestep [Wire] Controls NoControlFlag

Global phase gate: '1' -> '1'. The list of wires is just a hint for graphical rendering.

CNot Wire Controls NoControlFlag

Classical not: Cbit -> Cbit.

CGate String Wire [Wire] NoControlFlag

Generic classical gate 1 -> Cbit.

CGateInv String Wire [Wire] NoControlFlag

Uncompute classical gate Cbit -> 1, asserting that the classical bit is in the state specified by the corresponding CGate.

CSwap Wire Wire Controls NoControlFlag

Classical swap gate: Cbit * Cbit -> Cbit * Cbit.

QPrep Wire NoControlFlag

Initialization: Cbit -> Qbit.

QUnprep Wire NoControlFlag

Measurement Qbit -> Cbit with an assertion that the qubit is already in a computational basis state. This kind of measurement loses no information, and is formally the inverse of QPrep.

QInit Bool Wire NoControlFlag

Initialization: Bool -> Qbit.

CInit Bool Wire NoControlFlag

Initialization: Bool -> Cbit.

QTerm Bool Wire NoControlFlag

Termination of a Qbit wire with assertion that the qubit is in the specified state: Qbit * Bool -> 1.

CTerm Bool Wire NoControlFlag

Termination of a Cbit wire with assertion that the bit is in the specified state: Cbit * Bool -> 1.

QMeas Wire

Measurement: Qbit -> Cbit.

QDiscard Wire

Termination of a Qbit wire without assertion: Qbit -> 1

CDiscard Wire

Termination of a Cbit wire without assertion: Cbit -> 1

DTerm Bool Wire

Termination of a Cbit wire, with a comment indicating what the observed state of that wire was. This is typically inserted in a circuit after a dynamic lifting is performed. Unlike CTerm, this is in no way an assertion, but simply a record of observed behavior during a particular run of the algorithm.

Subroutine BoxId InverseFlag [Wire] Arity [Wire] Arity Controls NoControlFlag ControllableFlag RepeatFlag

Reference to a subroutine, assumed to be bound to another circuit. Arbitrary input and output arities. The domain of a1 must include the range of ws1, and similarly for a2 and ws2.

Comment String InverseFlag [(Wire, String)]

A comment. Does nothing, but can be useful for marking a location or some wires in a circuit.

Instances

Show Gate

Methods

showsPrec :: Int -> Gate -> ShowS

show :: Gate -> String

showList :: [Gate] -> ShowS

Basic information about gates

gate_arity :: Gate -> ([(Wire, Wiretype)], [(Wire, Wiretype)])

Compute the incoming and outgoing wires of a given gate (excluding controls, comments, and anchors). This essentially encodes the type information of the basic gates. If a wire is used multiple times as an input or output, then gate_arity also returns it multiple times; this enables run-time type checking.

Note that gate_arity returns the logical wires, and therefore excludes things like labels, comments, and graphical anchors. This is in contrast to wires_of_gate, which returns the syntactic set of wires used by the gate.

gate_controls :: Gate -> Controls

Return the controls of a gate (or an empty list if the gate has no controls).

gate_ncflag :: Gate -> NoControlFlag

Return the NoControlFlag of a gate, or False if it doesn't have one.

gate_with_ncflag :: NoControlFlag -> Gate -> Gate

Apply the given NoControlFlag to the given Gate. This means, if the first parameter is True, set the gate's NoControlFlag, otherwise do nothing. Throw an error if attempting to set the NoControlFlag on a gate that can't support this flag.

gate_reverse :: Gate -> Gate

Reverse a gate. Throw an error if the gate is not reversible.

Auxiliary functions on gates and wires

wires_of_controls :: Controls -> IntSet

Return the set of wires used by a list of controls.

wires_of_gate :: Gate -> IntSet

Return the set of wires used by a gate (including controls, labels, and anchors).

Unlike gate_arity, the function wires_of_gate is used for printing, and therefore returns all wires that are syntactically used by the gate, irrespective of whether they have a logical meaning.

wirelist_of_gate :: Gate -> [Wire]

Like wires_of_gate, except return a list of wires.

Dynamic arities

type ExtArity = XIntMap Wiretype

Recall that an Arity is a set of typed wires, and it determines the external interfaces at which circuits and gates can be connected. The type ExtArity stores the same information as the type Arity, but in a format that is more optimized for efficient updating. Additionally, it also stores the set of wires ever used.

arity_append_safe :: Gate -> ExtArity -> ExtArity

Check whether the given gate is well-formed and can be legally applied in the context of the given arity. If successful, return the updated arity resulting from the gate application. If unsuccessful, raise an error. Properties checked are:

• that each gate has non-overlapping inputs, including controls;
• that each gate has non-overlapping outputs, including controls;
• that the inputs of the gate (including controls) are actually present in the current arity;
• that the types of the inputs (excluding controls) match those of the current arity;
• that the outputs of the gate (excluding controls) don't conflict with any wires already existing in the current arity.

arity_append_unsafe :: Gate -> ExtArity -> ExtArity

Like arity_append, but without type checking. This is potentially faster, but should only used in applications that have already been thoroughly tested or type-checked.

arity_append :: Gate -> ExtArity -> ExtArity

For now, we disable run-time type checking, because we have not yet implemented run-time types properly. Therefore, we define arity_append to be a synonym for arity_append_unsafe.

arity_empty :: ExtArity

Return an empty arity.

arity_unused_wire :: ExtArity -> Wire

Return a wire unused in the current arity.

arity_unused_wires :: Int -> ExtArity -> [Wire]

Return the next k wires unused in the current arity.

arity_alloc :: Wiretype -> ExtArity -> (Wire, ExtArity)

Add a new typed wire to the current arity. This returns a new wire and the updated arity.

arity_of_extarity :: ExtArity -> Arity

Convert an extended arity to an ordinary arity.

n_of_extarity :: ExtArity -> Int

Return the smallest wire id nowhere used in the circuit.

Circuit abstraction

type Circuit = (Arity, [Gate], Arity, Int)

A completed circuit (a1,gs,a2,n) has an input arity a1, a list of gates gs, and an output arity a2. We also record n, the total number of wires used by the circuit. Because wires are allocated consecutively, this means that the wire id's used are [0..n-1].

wirelist_of_circuit :: Circuit -> [Wire]

Return the set of all the wires in a circuit.

Reversing low-level circuits

reverse_gatelist :: [Gate] -> [Gate]

Reverse a gate list.

reverse_circuit :: Circuit -> Circuit

Reverse a circuit. Throw an error if the circuit is not reversible.

NoControlFlag on low-level circuits

circuit_to_nocontrol :: Circuit -> Circuit

Set the NoControlFlag on all gates of a circuit.

Ordered circuits

newtype OCircuit

An ordered circuit is a Circuit together with an ordering on (usually all, but potentially a subset of) the input and output endpoints.

This extra information is required when a circuit is used within a larger circuit (e.g. via a Subroutine gate), to identify which wires of the sub-circuit should be bound to which wires of the surrounding circuit.

Constructors

OCircuit ([Wire], Circuit, [Wire])

reverse_ocircuit :: OCircuit -> OCircuit

Reverse an OCircuit. Throw an error if the circuit is not reversible.

Annotated circuits

data CircuitTypeStructure a

One often wants to consider the inputs and outputs of a circuit as more structuredtyped than just lists of bitsqubits; for instance, a list of six qubits could be structured as a pair of triples, or a triple of pairs, or a six-bit QDInt.

While for the most part this typing information is not included in low-level circuits, we need to consider it in hierarchical circuits, so that the information stored in a subroutine is sufficient to call the subroutine in a typed context.

Specifically, the extra information needed consists of functions to destructure the input/output data as a list of typed wires, and restructure such a list of wires into a piece of data of the appropriate type.

Constructors

CircuitTypeStructure (a -> ([Wire], Arity)) (([Wire], Arity) -> a)

id_CircuitTypeStructure :: CircuitTypeStructure ([Wire], Arity)

The trivial CircuitTypeStructure on ([Wire],Arity).

destructure_with :: CircuitTypeStructure a -> a -> ([Wire], Arity)

Use a CircuitTypeStructure to destructure a piece of (suitably typed) data into a list of typed wires.

structure_with :: CircuitTypeStructure a -> ([Wire], Arity) -> a

Use a CircuitTypeStructure to structure a list of typed wires (of the appropriate length/arity) into a piece of structured data.

Boxed circuits

data TypedSubroutine

A typed subroutine consists of:

• a low-level circuit, ordered to allow binding of incoming and outgoing wires;
• functions for structuring/destructuring the inputs and outputs to and from lists of wires (these functions being dynamically typed, since the input/output type may vary between subroutines);
• a ControllableFlag, recording whether the circuit is controllable.

Constructors

(Typeable a, Typeable b) => TypedSubroutine OCircuit (CircuitTypeStructure a) (CircuitTypeStructure b) ControllableFlag

circuit_of_typedsubroutine :: TypedSubroutine -> Circuit

Extract just the Circuit from a TypedSubroutine.

type Namespace = Map BoxId TypedSubroutine

A name space is a map from names to subroutine bindings. These subroutines can reference each other; it is the programmer’s responsibility to ensure there is no circular dependency, and no clash of names.

namespace_empty :: Namespace

The empty namespace.

showNames :: Namespace -> String

A function to display the names of all the subroutines in a Namespace.

type BCircuit = (Circuit, Namespace)

A boxed circuit is a distinguished simple circuit (analogous to a “main” function) together with a namespace.

Ordered circuits

type OBCircuit = (OCircuit, Namespace)

An ordered boxed circuit is a BCircuit together with an ordering on the input and output endpoints, or equivalently, an OCircuit together with a namespace.

ob_circuit :: [Wire] -> BCircuit -> [Wire] -> OBCircuit

Construct an OBCircuit from a BCircuit and an ordering on the input and output endpoints.

Basic functions lifted to boxed circuits

reverse_bcircuit :: BCircuit -> BCircuit

Reverse a simple boxed circuit, or throw an error if not reversible.

The ReadWrite monad

The ReadWrite monad encapsulates the interaction with a (real or simulated) low-level quantum device.

data ReadWrite a

The ReadWrite monad describes a standard read-write computation, here specialized to the case where writes are Gates, prompts are Bits, and reads are Bools. Thus, a read-write computation can do three things:

• terminate with a result. This is the case RW_Return.
• write a single Gate and continue. This is the case RW_Write.
• issue a prompt, which is a Wire, then read a Bool, then continue. This is the case RW_Read.

Constructors

RW_Return a

RW_Write !Gate (ReadWrite a)

RW_Read !Wire (Bool -> ReadWrite a)

RW_Subroutine BoxId TypedSubroutine (ReadWrite a)

Instances

Monad ReadWrite

Methods

(>>=) :: ReadWrite a -> (a -> ReadWrite b) -> ReadWrite b

(>>) :: ReadWrite a -> ReadWrite b -> ReadWrite b

return :: a -> ReadWrite a

fail :: String -> ReadWrite a

Functor ReadWrite

Methods

fmap :: (a -> b) -> ReadWrite a -> ReadWrite b

(<$) :: a -> ReadWrite b -> ReadWrite a

Applicative ReadWrite

Methods

pure :: a -> ReadWrite a

(<*>) :: ReadWrite (a -> b) -> ReadWrite a -> ReadWrite b

(*>) :: ReadWrite a -> ReadWrite b -> ReadWrite b

(<*) :: ReadWrite a -> ReadWrite b -> ReadWrite a

readwrite_wrap :: ReadWrite a -> ReadWrite ([Gate], a)

Transforms a read-write computation into one that behaves identically, but also returns the list of gates generated.

This is used as a building block, for example to allow a read-write computation to be run in a simulator while simultaneously using a static backend to print the list of generated gates.

readwrite_unwind_static :: ErrMsg -> ReadWrite a -> a

Extract the contents of a static ReadWrite computation. A ReadWrite computation is said to be static if it contains no RW_Read instructions, or in other words, no dynamic lifting. If an RW_Read instruction is encountered, issue an error message using the given stub.

gatelist_of_readwrite :: ErrMsg -> ReadWrite a -> Namespace -> ([Gate], Namespace, a)

Turn a static read-write computation into a list of gates, while also updating a namespace. "Static" means that the computation may not contain any RW_Read operations. If it does, the message "dynamic lifting" is passed to the given error handler.

Important usage note: This function returns a triple (gates, ns, x). The list of gates is generated lazily, and can be consumed one gate at a time. However, the values ns and x are only computed at the end of the computation. Any function using them should not apply a strict pattern match to ns or x, or else the whole list of gates will be generated in memory. For example, the following will blow up the memory:

(gates, ns, (a, n, x)) = gatelist_of_readwrite errmsg comp

whereas the following will work as intended:

(gates, ns, ~(a, n, x)) = gatelist_of_readwrite errmsg comp

Dynamic boxed circuits

type DBCircuit a = (Arity, ReadWrite (Arity, Int, a))

The type of dynamic boxed circuits. The type DBCircuit a is the appropriate generalization of (BCircuit, a), in a setting that is dynamic rather than static (i.e., with dynamic lifting or "interactive measurement").

bcircuit_of_static_dbcircuit :: ErrMsg -> DBCircuit a -> (BCircuit, a)

Convert a dynamic boxed circuit to a static boxed circuit. The dynamic boxed circuit may not contain any dynamic liftings, since these cannot be performed in a static setting. In case any output liftings are encountered, try to issue a meaningful error via the given stub error message.

dbcircuit_of_bcircuit :: BCircuit -> a -> DBCircuit a

Convert a boxed circuit to a dynamic boxed circuit. The latter, of course, contains no RW_Read instructions.