logo

Expression of type Prob

from the theory of proveit.physics.quantum.QPE

In [1]:
import proveit
# Automation is not needed when building an expression:
proveit.defaults.automation = False # This will speed things up.
proveit.defaults.inline_pngs = False # Makes files smaller.
%load_expr # Load the stored expression as 'stored_expr'
# import Expression classes needed to build the expression
from proveit import ExprRange, Variable, VertExprArray, t
from proveit.linear_algebra import ScalarMult, VecAdd
from proveit.numbers import Add, Exp, Interval, Mult, Neg, e, frac, i, one, pi, sqrt, two, zero
from proveit.physics.quantum import ket0, ket1, ket_plus
from proveit.physics.quantum.QPE import QPE1, _U, _ket_u, _phase, _s
from proveit.physics.quantum.circuits import Gate, Input, MultiQubitElem, Output, Qcircuit
from proveit.statistics import Prob
In [2]:
# build up the expression from sub-expressions
sub_expr1 = Variable("_a", latex_format = r"{_{-}a}")
sub_expr2 = Add(t, one)
sub_expr3 = Add(t, _s)
sub_expr4 = Interval(sub_expr2, sub_expr3)
sub_expr5 = MultiQubitElem(element = Gate(operation = QPE1(_U, t), part = sub_expr1), targets = Interval(one, sub_expr3))
expr = Prob(Qcircuit(vert_expr_array = VertExprArray([ExprRange(sub_expr1, Input(state = ket_plus), one, t), ExprRange(sub_expr1, MultiQubitElem(element = Input(state = _ket_u, part = sub_expr1), targets = sub_expr4), one, _s)], [ExprRange(sub_expr1, sub_expr5, one, t), ExprRange(sub_expr1, sub_expr5, sub_expr2, sub_expr3)], [ExprRange(sub_expr1, Output(state = ScalarMult(frac(one, sqrt(two)), VecAdd(ket0, ScalarMult(Exp(e, Mult(two, pi, i, Exp(two, Neg(sub_expr1)), _phase)), ket1)))), Add(Neg(t), one), zero).with_decreasing_order(), ExprRange(sub_expr1, MultiQubitElem(element = Output(state = _ket_u, part = sub_expr1), targets = sub_expr4), one, _s)])))
expr:
In [3]:
# check that the built expression is the same as the stored expression
assert expr == stored_expr
assert expr._style_id == stored_expr._style_id
print("Passed sanity check: expr matches stored_expr")
Passed sanity check: expr matches stored_expr
In [4]:
# Show the LaTeX representation of the expression for convenience if you need it.
print(stored_expr.latex())
\textrm{Pr}\left(\begin{array}{c} \Qcircuit@C=1em @R=.7em{
\qin{\lvert + \rangle} & \multigate{4}{\textrm{QPE}_1\left(U, t\right)} & \qout{\frac{1}{\sqrt{2}} \cdot \left(\lvert 0 \rangle + \left(\mathsf{e}^{2 \cdot \pi \cdot \mathsf{i} \cdot 2^{t - 1} \cdot \varphi} \cdot \lvert 1 \rangle\right)\right)} \\
\qin{\lvert + \rangle} & \ghost{\textrm{QPE}_1\left(U, t\right)} & \qout{\frac{1}{\sqrt{2}} \cdot \left(\lvert 0 \rangle + \left(\mathsf{e}^{2 \cdot \pi \cdot \mathsf{i} \cdot 2^{t - 2} \cdot \varphi} \cdot \lvert 1 \rangle\right)\right)} \\
\qin{\begin{array}{c}:\\ \left(t - 3\right) \times \\:\end{array}} & \ghost{\textrm{QPE}_1\left(U, t\right)} & \qout{\vdots} \\
\qin{\lvert + \rangle} & \ghost{\textrm{QPE}_1\left(U, t\right)} & \qout{\frac{1}{\sqrt{2}} \cdot \left(\lvert 0 \rangle + \left(\mathsf{e}^{2 \cdot \pi \cdot \mathsf{i} \cdot 2^{0} \cdot \varphi} \cdot \lvert 1 \rangle\right)\right)} \\
\qin{\lvert u \rangle} & \ghost{\textrm{QPE}_1\left(U, t\right)} & \qout{\lvert u \rangle}
} \end{array}\right)
In [5]:
stored_expr.style_options()
no style options
In [6]:
# display the expression information
stored_expr.expr_info()
 core typesub-expressionsexpression
0Operationoperator: 1
operand: 3
1Literal
2ExprTuple3
3Operationoperator: 4
operands: 5
4Literal
5ExprTuple6, 7, 8
6ExprTuple9, 10
7ExprTuple11, 12
8ExprTuple13, 14
9ExprRangelambda_map: 15
start_index: 88
end_index: 71
10ExprRangelambda_map: 16
start_index: 88
end_index: 72
11ExprRangelambda_map: 17
start_index: 88
end_index: 71
12ExprRangelambda_map: 17
start_index: 57
end_index: 58
13ExprRangelambda_map: 18
start_index: 19
end_index: 78
14ExprRangelambda_map: 20
start_index: 88
end_index: 72
15Lambdaparameter: 101
body: 21
16Lambdaparameter: 101
body: 22
17Lambdaparameter: 101
body: 23
18Lambdaparameter: 101
body: 24
19Operationoperator: 65
operands: 25
20Lambdaparameter: 101
body: 26
21Operationoperator: 43
operands: 27
22Operationoperator: 32
operands: 28
23Operationoperator: 32
operands: 29
24Operationoperator: 48
operands: 30
25ExprTuple31, 88
26Operationoperator: 32
operands: 33
27NamedExprsstate: 34
28NamedExprselement: 35
targets: 41
29NamedExprselement: 36
targets: 37
30NamedExprsstate: 38
31Operationoperator: 99
operand: 71
32Literal
33NamedExprselement: 40
targets: 41
34Operationoperator: 84
operand: 52
35Operationoperator: 43
operands: 49
36Operationoperator: 44
operands: 45
37Operationoperator: 50
operands: 46
38Operationoperator: 75
operands: 47
39ExprTuple71
40Operationoperator: 48
operands: 49
41Operationoperator: 50
operands: 51
42ExprTuple52
43Literal
44Literal
45NamedExprsoperation: 53
part: 101
46ExprTuple88, 58
47ExprTuple54, 55
48Literal
49NamedExprsstate: 56
part: 101
50Literal
51ExprTuple57, 58
52Literal
53Operationoperator: 59
operands: 60
54Operationoperator: 81
operands: 61
55Operationoperator: 62
operands: 63
56Literal
57Operationoperator: 65
operands: 64
58Operationoperator: 65
operands: 66
59Literal
60ExprTuple67, 71
61ExprTuple88, 68
62Literal
63ExprTuple69, 70
64ExprTuple71, 88
65Literal
66ExprTuple71, 72
67Literal
68Operationoperator: 95
operands: 73
69Operationoperator: 84
operand: 78
70Operationoperator: 75
operands: 76
71Variable
72Literal
73ExprTuple97, 77
74ExprTuple78
75Literal
76ExprTuple79, 80
77Operationoperator: 81
operands: 82
78Literal
79Operationoperator: 95
operands: 83
80Operationoperator: 84
operand: 88
81Literal
82ExprTuple88, 97
83ExprTuple86, 87
84Literal
85ExprTuple88
86Literal
87Operationoperator: 89
operands: 90
88Literal
89Literal
90ExprTuple97, 91, 92, 93, 94
91Literal
92Literal
93Operationoperator: 95
operands: 96
94Literal
95Literal
96ExprTuple97, 98
97Literal
98Operationoperator: 99
operand: 101
99Literal
100ExprTuple101
101Variable