arxiv: v2 [quant-ph] 2 Aug 2012
|
|
- Magdalene Young
- 5 years ago
- Views:
Transcription
1 Qcomiler: quantum comilation with CSD method Y. G. Chen a, J. B. Wang a, a School of Physics, The University of Western Australia, Crawley WA 6009 arxiv: v2 [quant-h] 2 Aug 202 Abstract In this aer, we resent a general quantum comutation comiler, which mas any given quantum algorithm to a quantum circuit consisting a sequential set of elementary quantum logic gates based on recursive cosinesine decomosition. The resulting quantum circuit diagram is rovided by directly linking the ackage outut written in LaTex to Qcircuit.tex <htt:// We illustrate the use of the Qcomiler ackage through various examles with full details of the derived quantum circuits. Besides its generality and simlicity, Qcomiler roduces quantum circuits which reflect the symmetry of the systems under study. Corresonding author address: jingbo.wang@uwa.edu.au (J. B. Wang) Prerint submitted to August 3, 202
2 Program summary Program Title: Qcomiler Journal Reference: Catalogue identifier: Licensing rovisions: none Programming language: Fortran Comuter: any comuter with a Fortran comiler Oerating system: Linux, Mac OS X 0.5 (and later) RAM: deend on the size of the unitary matrix to be decomosed Keywords: Quantum comiler, quantum circuit, unitary matrix, quantum gate, quantum algorithm Classification: External routines/libraries: Laack Nature of roblem: Decomose any given unitary oeration into a quantum circuit with only elementary quantum logic gates. Solution method: This ackage decomoses an arbitrary unitary matrix, by alying the CSD algorithm recursively, into a series of block-diagonal matrices, which can then be readily associated with elementary quantum gates to form a quantum circuit. Restrictions: The only limitation is imosed by the available memory on the user s comuter. Comments: This ackage is alicable for any arbitrary unitary matrices, both real and comlex. If the unitary matrix is real, its corresonding quantum circuit is much simler with only half number of quantum gates in comarison with comlex matrices of the same size. Running time: Memory and CPU time requirements deend critically on the size of the unitary matrix to be decomosed. Most examles resented in this aer require a few minutes of CPU time on Intel Pentium Dual Core 2 Duo 2.2GHz. 2
3 . Introduction Quantum comuting exloits the nature of the quantum world in a way that romises to solve roblems which are intractable using conventional comuters [ 3]. At the heart of a quantum comuter lies a set of qubits whose states are maniulated by a series of elementary quantum logic gates, namely a quantum circuit, to rovide the ultimate comutational results. In our earlier work [4], we develoed a highly efficient quantum comutation simulator to assist on the analysis of comlicated quantum circuits comrised of qubit and qudit quantum gates. In this aer, we resent a quantum comutation comiler which mas any given quantum algorithm to a quantum circuit consisting a set of elementary quantum logic gates. In the seminal aers by Barenco et al and Deutsch et al [5, 6], it was roven that any arbitrarily comlex unitary oeration can be imlemented by a quantum circuit involving only one- or two-qubit elementary quantum logic gates. Earlier studies alied the standard triangularization or QR-factorization scheme with Givens rotations and Gray codes to ma a quantum algorithm to a series of elementary gate oerations[3, 5 7]. Subsequently, a more efficient and versatile scheme based on the cosine-sine decomosition was roosed and utilized [8 3]. More recently, De Vos et al [4, 5] examined another decomosition scheme, namely the Birkhoff decomosition, which was found to rovide simler quantum circuits for certain tyes of unitary matrices than the cosine-sine decomosition. However, the Birkhoff decomosition does not work for general unitary matrices. In this work we develo a general quantum comiler, named the Qcomiler, based on the cosine-sine decomosition scheme, because it works for arbitrary unitary matrices, it is extremely adatable, and the core CSD code is now available as art of the LAPACK ackage. Furthermore, we have managed to half the number of quantum gates in the circuit, if the quantum algorithm involves only real unitary matrices. For the Qcomiler, the inut file contains a user secified unitary matrix U. Qcomiler alies the cosine-sine decomosition (CSD) recursively and comiles U into a quantum circuit with a sequence of controlled/uncontrolled rotation and hase gates. Its outut contains the comlete information on these elementary quantum gates, including a searate LaTex document which can be directly linked to Qcircuit.tex <htt:// to roduce the resulting quantum circuit diagram. This aer is structured as follows. Section 2 describes the recursive cosine-sine decomosition scheme, which mas an arbitrary unitary matrix U into a quantum circuit with only one- and two-qubit logic gates. In Section 3, we rovide a more cost-efficient method for real unitary matrices, which significantly reduces the number of quantum gates in the final circuit. In Section 4, we discuss the general structure as well as the usage of the Qcomiler ackage. In Section 5, we resent various examles with full details and discussions 3
4 of the resulting quantum circuits. Section 6 concludes the aer. 2. The recursive CSD scheme: General Unitary Matrices The Cosine Sine Decomosition (CSD) algorithm as described in [6, 7] decomoses an arbitrary 2 n 2 n unitary matrix U as the following u C S x U = v S C y, () where cos θ cos θ 2 C = diag (cos θ l ) = l=,...,2 n... cos θ 2 n and sin θ sin θ 2 S = diag (sin θ l ) = l=,...,2 n... sin θ 2 n are 2 n 2 n diagonal matrices, while u, v, x and y are 2 n 2 n unitary matrices which can be further decomosed by the CSD algorithm recursively into a string of block diagonal unitary matrices [9]. To show the general structure of the decomosition, we denote all sub-matrices as u i dim ignoring their exlicitly different numerical values, where i indicates the level of recursion and dim the matrix dimension. At the first level, we have U = u 2 n 2 n u 2 n 2 n u 2 n 2 n u 2 n 2 n =U A U. 4 C 2 n 2 n S 2 n 2 n S 2 n 2 n C 2 n 2 n (2)
5 At the second level, u = 2 n 2 u 2 2 n 2 2 n 2 n u 2 2 n 2 2 n 2 C 2 2 n 2 2 n 2 S 2 2 n 2 2 n 2 S 2 2 n 2 2 n 2 C 2 2 n 2 2 n 2 u 2 2 n 2 2 n 2 u 2, 2 n 2 2 n 2 U = u 2 n 2 n u 2 n 2 n = u 2 2 n 2 2 n 2 u 2 2 n 2 2 n 2 u 2 C 2 2 n 2 2 n 2 S 2 2 n 2 2 n 2 S 2 2 n 2 2 n 2 C 2 2 n 2 2 n 2 C 2 2 n 2 2 n 2 u 2 2 n 2 2 n 2 (3) 2 n 2 2 n 2 S 2 2 n 2 2 n 2 S 2 2 n 2 2 n 2 C 2 2 n 2 2 n 2 u 2 2 n 2 2 n 2 u 2 2 n 2 2 n 2 u 2 =U 2 A 2 U 2, 2 n 2 2 n 2 u 2 2 n 2 2 n 2 and U = U A U = U 2 A 2 U 2 A U 2 A 2 U 2. (4) At the ith level of recursion, the matrix U i is decomosed as the following: U i = U i A i U i (5) 5
6 where u i 2 n i 2 n i U i = u i 2 n i 2 n i... (6) u i 2 n i 2 n i and C i 2 n i 2 n i S i 2 n i 2 n i A i = S i 2 n i 2 n i C i 2 n i 2 n i.... (7) C i 2 n i 2 n i S i 2 n i 2 n i S i 2 n i 2 n i C i 2 n i 2 n i At the end of the recursive rocess, i.e. the nth level, we obtain U =U A U =U 2 A 2 U 2 A U 2 A 2 U 2 =U 3 A 3 U 3 A 2 U 3 A 3 U 3 A U 3 A 3 U 3 A 2 U 3 A 3 U 3 = 2 n = U n A i() Un 2 n, = (8) where marks the osition of the matrix sequence, i is imlicitly determined by, A is given by Eq. (7), and U n = diag k=,...,2 n (u n,k ) = diag k=,...,2 n (ex(iϕ,k )). (9) The above described recursive CSD scheme works for 2 n 2 n unitary matrices. For an arbitrary N N unitary matrix U, where 2 n < N 2 n, we add an identity matrix to U to form a new unitary matrix W 2 n 2 n = U N N I (2n N) (2n N), (0) and then aly the recursive CSD decomosition to W as described above. The decomosed matrices given by Eq. (8) can be directly related to elementary quantum gates, in articular, the hase gate Φ and the controlled rotation gate [3, 9], where the rotation oeration is defined as 6
7 If the rotation axis is y or z, we have R a (ρ) = ex(ia σ ρ 2 ) = I cos ρ 2 + ia σ sin ρ 2. () (ρ) = (2θ) = ex(iσ y θ) = I cos θ + iσ y sin θ cos θ sin θ = sin θ cos θ (2) and (ρ) = (2φ) = ex(iσ z φ) = I cos φ + iσ z sin φ ex(iφ) 0 (3) = 0 ex( iφ), resectively. P i() To establish a maing to elementary quantum gates, Möttönen et al [9] inserted identity matrices I = (P i() ) after each A i() in Eq. (8) with which P i() commutes, and thus U n A i() = U n A i() P i() (P i() ) = U n P i() A i() (P i() ). (4) The final decomosition becomes 2 n U = = B i() A i() Ũn 2 n, (5) where Ũ n = (P i( ) ) U n, excet Ũ n = Un, and Bi() = Ũ n P i() is a 2 n 2 n diagonal unitary matrix. With a set of P i() gate secified by solving a set of linear equations, B i() will have the required symmetry to be equivalent to C n (i;,..., i, i +,..., n), (6) where i denotes the target qubit with n control qubits. Similarly, A i can be maed to gate C n (i;,..., i, i +,..., n), (7) 7
8 and Ũ n 2 n is equivalent to a series of gates Φ(; ) I 2 n 2 n, (; ) I 2 n 2 n, C (2; ) I 2 n 2 2 n 2,... (8) C n 2 (n ;,..., n 2) I 2 2, C n (n;,..., n ). As an examle, we decomose a unitary matrix by alying the CSD scheme recursively, i.e. U = B 3 A3 B2 2 A2 2 B3 3 A3 3 B 4 A 4 B3 5 A3 5 B2 6 A2 6 B3 7 A3 7Ũ3 8, (9) where A ( = 4) C 2 (; 2, 3), A 2 ( = 2, 6) C 2 (2;, 3), A 3 ( =, 3, 5, 7) C 2 (), B ( = 4) C 2 (; 2, 3), (20) B 2 ( = 2, 6) C 2 (2;, 3), B 3 ( =, 3, 5, 7) C 2 (), and U 3 8 Φ(; ), (; ), C (2; ), C 2 (). For simlicity, the identity matrices are omitted in the above exressions. The equivalent quantum circuit for U given by Eq. 9 is shown in Fig.. For a unitary matrix, we have U =B 4 A4 B3 2 A3 2 B4 3 A4 3 B2 4 A2 4 B4 5 A4 5 B3 6 A3 6 B4 7 A4 7 B 8 A 8 B 4 9 A4 9 B3 0 A3 0 B4 A4 B2 2 A2 2 B4 3 A4 3 B3 4 A3 4 B4 5 A4 5Ũ4 6. (2) The corresonding quantum circuit is shown in Fig. 2. 8
9 U = Ũ 3 8 Ũ 3 8 A 3 7 B 3 7 A 2 6 B 2 6 A 3 5 B 3 5 A 4 B 4 A 3 3 B 3 3 A 2 2 B 2 2 A 3 B 3 Φ Ũ 3 8 = Figure : Quantum circuit for an 8-by-8 comlex unitary matrix. = U Ũ 4 6 Ũ 4 6 A 4 5 B 4 5 A 3 4 B 3 4 A 4 3 B 4 3 A 2 2 B 2 2 A 4 B 4 A 3 0 B3 0 A 4 9 B 4 9 A 8 B 8 A 4 7 B 4 7 A 3 6 B 3 6 A 4 5 B 4 5 A 2 4 B 2 4 A 4 3 B 4 3 A 3 2 B 3 2 A 4 B 4 Φ = Ũ 4 6 Figure 2: Quantum circuit for a 6-by-6 comlex unitary matrix. 3. The recursive CSD scheme: Real Unitary Matrices If the 2 n 2 n unitary matrix U is real, we have u n,k = ex(iϕ,k) = or from Eq. 9 and U n becomes a diagonal matrix consisting of only or. In this case, we can insert the identity matrix I = Ũ n (Ũ n ) = Ũ n Ũn after each A i() in Eq. 8, i.e. U n A i() = U n A i() Ũ n Ũn (22) 9
10 where Ũ n = Uq. n The decomosition for a real unitary U then becomes where à i() q= 2 n U = à i() = Ũn 2 n (23) = Ũ n A i() Ũ n. Note for general unitary matrices, it requires the solution of a set of linear equations to determine each inserted P i() in Eq. 4 [9]. Here, for real unitary matrices, we only need to calculate the roduct of U n, Un 2,..., Un to obtain the inserted Ũ n. Furthermore, we end u with only half number of decomosed matrices in comarison with Eq. 5. Since Ũ n is a diagonal matrix with only two ossible values, or, à i() A i() can be readily maed to the gate C n (i;,..., i, i +,..., n). The gates for à i() are almost the same excet for the signs of some rotation angles. and The last matrix Ũ n 2 n can be written as Ũ n 2 n = n diag (D m,k ) I 2 m 2m (24) k=,...,2 n m m= where D m,k are either ( ) ( 0 0 or 0 ) 0, with the former being an identity gate and the later a Π gate. Therefore, U n 2 n is equivalent to a subset of the following quantum gates: Π(; ) I 2 n 2 n C Π(2; ) I 2 n 2 2 n 2... C n 2 Π(n ;,..., n 2) I 2 2 C n Π(n;,..., n ). As an examle we decomose a 6-by-6 real unitary matrix U and get U = à 4 Ã3 2Ã4 3Ã2 4Ã4 5Ã3 6Ã4 7à 8Ã4 9Ã3 0Ã4 Ã2 2Ã4 3Ã3 4Ã4 5Ũ4 6 (25) Its corresonding circuit is, as shown in Fig. 3, much simler than that for a 6-by-6 comlex unitary matrix. The number of gates is reduced to half. 0
11 = U Ũ 4 6 Ũ 4 6 Ã 4 5 Ã 3 4 Ã 4 3 Ã 2 2 Ã 4 Ã 3 0 Ã 4 9 Ã 8 Ã 4 7 Ã 3 6 Ã 4 5 Ã 2 4 Ã 4 3 Ã 3 2 Ã 4 Ũ 4 6 = the subset of Figure 3: The circuit of a 6-by-6 real U. 4. Overview of the Software The Qcomiler ackage, written in Fortran, consists of main rogram and 2 subroutines. ackage utilizes recursively Sutton s CSD subroutines [7], which is now art of Laack available at htt:// under a ermissive free software license. As described in Section 3, the Qcomiler ackage decomoses U as U = ( 2 n ( = 2 n = ) Ã i() Ũ n 2 n B i() A i() ) Ũ n 2 n if U is real, if U is comlex. It then generates an outut file containing a comlete descrition of each quantum gate in the circuit. The rincial flow chart of Qcomiler is shown in Fig. 4 with all subroutines calls summarized in Fig. 5. This (26)
12 Start Main Program CYG_CSD Real Case Main Program CYGR_CSD Read in the real unitary matrix CYGR_READF Create a useful table called index CYG_INDEXTABLE Choose the matrix style Comlex Case Main Program CYGC_CSD Read in the comlex unitary matrix CYGC_READF Create a useful table called index CYG_INDEXTABLE CSD for one ste LAPACK: DORCSD CSD recursively CYGR_BLKCSD CSD recursively CYGC_BLKCSD CSD for one ste LAPACK: ZUNCSD Process the CSD result in the real case to match each decomosed matrix with a quantum gate CYGR_CUTGATE Create the outut file (details of each gate in the circuit) CYGR_BLKCSD Process the CSD result in the comlex case to match each decomosed matrix with a quantum gate CYGC_CSDPHASE Create the outut file (details of each gate in the circuit) CYGC_BLKCSD End Figure 4: Flow chart for the Qcomiler ackage. CYGR_READF CYGR_CSD CYGR_BLKCSD CYGR_CUTGATE CYGR_WRITEF LAPACK: DORCSD CYG_CSD CYG_INDEXTABLE CYGC_CSD CYGC_READF CYGC_BLKCSD CYGC_CSDGATEU CYGC_CSDPHASE CYGC_WRITEF LAPACK: ZUNCSD CYGC_COEFF Figure 5: Flow chart for all subroutine calls. 2
13 5. Examle 5.. Comlex unitary matrix An 8-by-8 comlex unitary matrix U is generated randomly as the following: U = i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i Recursive CSD decomosition gives U = B 3 A3 B2 2 A2 2 B3 3 A3 3 B 4 A 4 B3 5 A3 5 B2 6 A2 6 B3 7 A3 7Ũ3 8. Their corresonding quantum gates are detailed in the outut file roduced by Qcomiler, as exlained below.. (27) Outut of the ackage Matrix Gate GATEPHASE ; Φ(; ) Φ = Φ GATEZ ; Ũ 3 8 (; ) 2φ = 2 ( ) GATEZ 2; C (2; ) 2φ = 2 ( ) 2φ 2 = = R z R 2 z GATEZ C 2 () 2φ = φ 2 = 2 ( ) 2φ 3 = 2 ( ) 2φ 4 = = R z R 2 z R 3 z R 4 z 3
14 Outut of the ackage Matrix Gate A 3 7 C 2 () 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y GATEZ B 3 7 C 2 () 2φ = 2 ( ) 2φ 2 = φ 3 = 2 ( ) 2φ 4 = 2 ( 0.209) = R z R 2 z R 3 z R 4 z 2;, A 2 6 C 2 (2;, 3) 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y GATEZ 2;, B 2 6 C 2 (2;, 3) 2φ = 2 ( ) 2φ 2 = 2 ( 0.227) 2φ 3 = φ 4 = = R z R 2 z R 3 z R 4 z A 3 5 C 2 () 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y 4
15 Outut of the ackage Matrix Gate GATEZ B 3 5 C 2 () 2φ = 2 ( 0.346) 2φ 2 = φ 3 = 2 ( 0.530) 2φ 4 = 2 ( 0.047) = R z R 2 z R 3 z R 4 z ; 2, A 4 C 2 (; 2, 3) 2θ = θ 2 = θ 3 = θ 4 = R y R 2 y R 3 y R 4 y = GATEZ ; 2, B 4 C 2 (; 2, 3) 2φ = φ 2 = φ 3 = φ 4 = 2 ( ) R z R 2 z R 3 z R 4 z = A 3 3 C 2 () 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y GATEZ B 3 3 C 2 () 2φ = 2 ( ) 2φ 2 = 2 ( 0.300) 2φ 3 = 2 ( ) 2φ 4 = = R z R 2 z R 3 z R 4 z 5
16 Outut of the ackage Matrix Gate 2;, A 2 2 C 2 (2;, 3) 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y GATEZ 2;, B 2 2 C 2 (2;, 3) 2φ = 2 ( 0.607) 2φ 2 = 2 ( ) 2φ 3 = φ 4 = = R z R 2 z R 3 z R 4 z A 3 C 2 () 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y GATEZ B 3 C 2 () 2φ = 2 ( ) 2φ 2 = 2 ( ) 2φ 3 = 2 ( 0.305) 2φ 4 = 2 ( 0.447) = R z R 2 z R 3 z R 4 z The final quantum circuit of the comlex unitary matrix U given by Eq. 27 is shown in Fig. 6. Φ U = Ũ 3 8 A 3 7 B 3 7 A 2 6 B 2 6 A 3 5 B 3 5 A 4 B 4 A 3 3 B 3 3 A 2 2 B 2 2 A 3 B 3 Figure 6: The circuit of U, the random 8-by-8 comlex unitary matrix. 6
17 It requires 8 gates for Ũ 3 8, 4 7 gates for A i, and 4 7 gates for B i, thus 64 gates in total to build the equivalent quantum circuit Real unitary matrix An 8-by-8 real unitary matrix U is generated randomly as the following: U = (28) Recursive CSD decomosition gives U = Ã 3 Ã2Ã3Ã4Ã3 5Ã26Ã3 7Ũ3 8. Their corresonding quantum gates are detailed in the outut file roduced by Qcomiler, as exlained below. Outut of the ackage Matrix Gate GATEPI ; Y Π(; ) GATEPI 2; Y N Ũ 3 8 C Π(2; 0 ) GATEPI Y Y N Y C 2 Π(3; (, 2) 00 ) C 2 Π(3; (, 2) 0 ) C 2 Π(3; (, 2) ) Note: This art secifies the subset of gates: Ã 3 7 C 2 () 2θ = 2 ( 0.388) 2θ 2 = θ 3 = 2 ( ) 2θ 4 = = R y R 2 y R 3 y R 4 y 7
18 Outut of the ackage Matrix Gate 2;, à 2 6 C 2 (2;, 3) 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y à 3 5 C 2 () 2θ = 2 ( 0.026) 2θ 2 = θ 3 = 2 ( ) 2θ 4 = 2 ( ) = R y R 2 y R 3 y R 4 y ; 2, à 4 C 2 (; 2, 3) 2θ = 2 ( 0.052) 2θ 2 = θ 3 = 2 ( ) 2θ 4 = 2 ( ) R y R 2 y R 3 y R 4 y = à 3 3 C 2 () 2θ = θ 2 = 2 ( 0.357) 2θ 3 = 2 ( 0.277) 2θ 4 = = R y R 2 y R 3 y R 4 y 2;, à 2 2 C 2 (2;, 3) 2θ = θ 2 = θ 3 = θ 4 = = R y R 2 y R 3 y R 4 y 8
19 Outut of the ackage Matrix Gate à 3 C 2 () 2θ = θ 2 = 2 ( 0.267) 2θ 3 = 2 ( ) 2θ 4 = 2 ( 0.28) = R y R 2 y R 3 y R 4 y The final quantum circuit of the real unitary matrix U given by Eq. 28 is shown in Fig. 7, which is significantly simler than that shown in Fig. 6 for a comlex unitary matrix of the same size. U = Ũ 3 8 à 3 7 à 2 6 à 3 5 à 4 à 3 3 à 2 2 à 3 Figure 7: The circuit of U, the random 8-by-8 real unitary matrix. It requires 5 Π gates for Ũ 3 8, and 4 7 gates for à i, thus 33 gates in total to build the equivalent quantum circuit, significantly simler than its comlex counterart Random walk on a square grah Quantum walks have recently been exlored for their non-intuitive dynamics, which may hold the key to radically new quantum algorithms [8 20]. Douglas and Wang [2] recently resented efficient quantum circuits for several families of highly symmetrical grahs. Here we use Qcomiler to build quantum circuits for quantum walk on general grahs. The walker s quantum state is given : ψ = N A i, j i, j, (29) i= j S where N is the number of nodes of the given grah, i reresents the ith node state, j the coin state which connects the ith and jth node, and S is determined by the grah adjacent matrix. 9
20 To realize one ste quantum walk, we first aly a coin oerator Ĉ = diag (Ĉ i ) on the walker, where i=,...,n Ĉ i = j S k S c i j,k i, k i, j. (30) In this work, we choose the Grover Coin [22] so that the quantum walk can be reresented by a real matrix. We then aly the translation oerator ˆT = N j, i i, j. (3) i= j S The unitary matrix for the comlete ste is Û = ˆTĈ. For a simle square grah, shown in Fig. 8, the quantum walk oerator is Û = ˆTĈ = = Figure 8: The square grah The outut file of Qcomiler is given as the following with its quantum circuit shown in Fig. 9. Outut of the ackage Matrix Gate GATEPI ; Y Π(; ) GATEPI 2; Y N Ũ 3 8 C Π(2; 0 ) 20
21 Outut of the ackage Matrix Gate GATEPI N Y Y N C 2 Π(3; (, 2) 0 ) C 2 Π(3; (, 2) 0 ) Note: This art secifies the subset of gates: à 3 7 C 2 () 2θ = = 2θ 4 = 2 ( 0.5) = R y R 4 y 2;, à 2 6 C 2 (2;, 3) 2θ 2 = = 2θ 4 = = R 2 y R 4 y à 3 5 C 2 () 2θ = 2 ( 0.5) = 2θ 2 = 2 ( 0.5) = R y R 2 y ; 2, à 4 C 2 (; 2, 3) 2θ 3 = 2 ( 0.5) = 2θ 4 = = R 3 y R 4 y à 3 3 C 2 () 2θ 2 = = 2θ 3 = = R 2 y R 3 y 2
22 Outut of the ackage Matrix Gate 2;, à 2 2 C 2 (2;, 3) 2θ 2 = 2 ( 0.5) = 2θ 4 = = R 2 y R 4 y à 3 C 2 () 2θ 2 = 2 ( 0.5) = 2θ 3 = 2 ( 0.5) = R 2 y R 3 y R 3 y R 4 y U = R 2 y R 4 y R 2 y R 4 y R y R 4 y R y R 2 y R 2 y R 3 y R 2 y R 3 y Ũ 3 8 à 3 7 à 2 6 à 3 5 à 4 à 3 3 à 2 2 à 3 Figure 9: The circuit of U, the random walk evolution on the square grah (Fig. 8). Comare to the random 8-by-8 real matrix given by Eq. 28, which requires 33 gates, the circuit of the square grah consists of only 8 gates, reflecting the simle symmetry of this grah Random walk on the 8-star grah Fig. 0 shows another simle grah, namely the 8-star grah 22
23 Figure 0: The 8-star grah For this grah, the quantum walk oerator is Û = ˆTĈ = , which is a 6 6 unitary matrix. The comlete result is dislayed below. Note that since all of the rotation angles of à 4, Ã3 2, Ã4 3, Ã2 4, Ã4 5, Ã3 6, Ã4 7 are zeros, these gates vanish and are omitted here. Outut of the ackage Matrix Gate GATEPI ; N GATEPI 2; N N GATEPI N N N N GATEPI 4;, 2, 3 N N N N N N N Y Ũ 4 6 C 3 Π(4; (, 2, 3) ) 23
24 Outut of the ackage Matrix Gate Note: This art secifies the subset of gates: 4;, 2, Ã 4 5 C 3 (4;, 2, 3) 2θ 5 = 2 ( 0.53) 2θ 6 = θ 7 = 2 ( ) 2θ 8 = 2 ( ) R 5 y R 6 y R 7 y R 8 y, Ã 3 4 C 3 (, 4) 2θ 5 = 2 ( 0.289) 2θ 6 = 2 ( ) 2θ 7 = θ 8 = R 5 y R 6 y R 7 y R 8 y 4;, 2, Ã 4 3 C 3 (4;, 2, 3) 2θ 5 = θ 6 = 2 ( ) 2θ 7 = 2 ( ) 2θ 8 = R 5 y R 6 y R 7 y R 8 y 2;, 3, Ã 2 2 C 3 (2;, 3, 4) 2θ 8 = = R 8 y 24
25 Outut of the ackage Matrix Gate 4;, 2, Ã 4 C 3 (4;, 2, 3) 2θ 5 = 2 ( ) 2θ 6 = θ 7 = θ 8 = R 5 y R 6 y R 7 y R 8 y, Ã 3 0 C 3 (, 4) 2θ 5 = θ 6 = θ 7 = 2 ( 0.074) 2θ 8 = R 5 y R 6 y R 7 y R 8 y 4;, 2, Ã 4 9 C 3 (4;, 2, 3) 2θ 5 = θ 6 = 2 ( 0.573) 2θ 7 = θ 8 = 2 ( ) R 5 y R 6 y R 7 y R 8 y ; 2, 3, Ã 8 C 3 (; 2, 3, 4) 2θ = 2θ 2 = 2θ 3 = 2θ 4 = 2θ 5 = 2θ 6 = 2θ 7 = 2θ 8 = 2 ( 0.5) = R y R 2 y R 8 y =... Putting these gates together, we obtain its quantum circuit as shown in Fig.. It requires only 34 quantum gates to imlement the 6 6 unitary matrix. As a comarison, for a random 6 6 real matrix, it needs around = 35 gates, while for a random 6 6 comlex matrix, it 25
26 U = Ũ 4 6 Ã 4 5 Ã 3 4 Ã 4 3 Ã 2 2 Ã 4 Ã 3 0 Ã 4 9 Ã 8 Figure : The circuit of U, the random walk evolution on the star grah (Fig. 0). needs = 256 gates. The last 7 gate combinations (56 subgates) comletely disaear. This is another examle demonstrating the efficiency of Qcomiler Random walk on a comlicated grah Qcomiler also works for very comlicated grahs. Fig. 2 is an adjacent matrix of a 00-node grah with white dots standing for and black dots for 0. Figure 2: The adjacent matrix of a comlicated grah. The quantum walk oerator for this grah is a unitary matrix U. We first exand the size of 26
27 the unitary matrix with an identity matrix, i.e. U W = I (32) where U is the unitary matrix, I is a identity matrix, and W is a unitary matrix. The resulting quantum circuit contains 2 qubits. 6. Conclusions We have develoed an efficient and versatile ackage, Qcomiler, which mas any unitary matrix U of arbitrary size into a quantum circuit with only one- and two-qubit logic gates by alying cosine-sine decomosition recursively. For real unitary matrices, Qcomiler rovides a much simler quantum circuit with only half number of elementary gates in comarison with earlier work by Mottonen et al [9] which deals with general unitary matrices. The quantum circuits roduced by Qcomiler also reflects the symmetry of the systems under study. In articular, we examined the resulting quantum circuits corresonding to quantum walks on grahs with certain degree of symmetry [9]. In this case, many gates in such a circuit turned out to be an identity gate or a simle gate, which can be easily eliminated or combined with other gates to further reduce the comlexity of the final quantum circuit. For unitary matrices of certain structure or symmetry, it may be ossible to include qutris and qudits to roduce simler and more efficient quantum circuits. This will be an interesting subject for further study. Acknowledgements The authors would like to thank Sergey Koulikov, Blake Segler, and Gary Allwood for their earlier involvement in this roject. Thanks are also due to Kia Manouchehri, Brendan Douglas and Thomas Loke for several valuable discussions. References [] Deutsch, D., Proceedings of the Royal Society of London Series A - Mathematical Physical and Engineering Sciences 425 (989) 73. [2] Shor, P. W., SIAM Journal on Comuting 26 (997) 484. [3] Nielsen, M. and Chuang, I., Quantum comutation and quantum information, Cambridge, Cambridge,
28 [4] Loke, T. and Wang, J. B., Comuter Physics Communications 82 (20) [5] Barenco, A. et al., Physical Review A 52 (995) [6] Deutsch, D., Barenco, A., and Ekert, A., Proceedings of the Royal Society of London Series a-mathematical Physical and Engineering Sciences 449 (995) 669. [7] Cybenko, G., Comuting in Science and Engineering 3 (200) 27. [8] Tucci, R. R., quant-h/ (999), quant-h/04027 (2004), quant-h/ (2004). [9] Mottonen, M., Vartiainen, J. J., Bergholm, V., and Salomaa, M. M., Physical Review Letters 93 (2004) [0] Bergholm, V., Vartiainen, J. J., Mottonen, M., and Salomaa, M. M., Physical Review A 7 (2005) [] Shende, V. V. and Markov, I. L., Quantum Information & Comutation 5 (2005) 49. [2] Khan, F. S. and Perkowski, M., Theoretical Comuter Science 367 (2006) 336. [3] Manouchehri, K. and Wang, J. B., Physical Review A 80 (2009) [4] De Vos, A. and Van Rentergem, Y., Journal of Multile-Valued Logic and Soft Comuting 5 (2009) 489. [5] De Vos, A., Boes, M., and De Baerdemacker, S., Journal of Multile-Valued Logic and Soft Comuting 8 (20) 67. [6] Paige, C. C. and Wei, M., Linear Algebra and Its Alications 209 (994) 303. [7] Sutton, B., Numerical Algorithms 50 (2009) 33. [8] Keme, J., Contemorary Physics 44 (2003) 307. [9] Douglas, B. L. and Wang, J. B., J. Phys. A: Math. Theor 4 (2008) [20] Berry, S. D. and Wang, J. B., Physical Review A 83 (20) [2] Douglas, B. L. and Wang, J. B., Physical Review A 79 (2009) [22] Shenvi, N., Keme, J., and Whaley, K. B., Physical Review A 67 (2003)
arxiv: v1 [quant-ph] 3 Feb 2015
From reversible comutation to quantum comutation by Lagrange interolation Alexis De Vos and Stin De Baerdemacker 2 arxiv:502.0089v [quant-h] 3 Feb 205 Cmst, Imec v.z.w., vakgroe elektronica en informatiesystemen,
More informationTransformation of quantum states using uniformly controlled rotations
Transformation of quantum states using uniformly controlled rotations Mikko Möttönen, Juha J. Vartiainen, Ville Bergholm, and Martti M. Salomaa Materials Physics Laboratory, P.O. Box 2200, FIN-02015 Helsinki
More informationComputer arithmetic. Intensive Computation. Annalisa Massini 2017/2018
Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationFault Tolerant Quantum Computing Robert Rogers, Thomas Sylwester, Abe Pauls
CIS 410/510, Introduction to Quantum Information Theory Due: June 8th, 2016 Sring 2016, University of Oregon Date: June 7, 2016 Fault Tolerant Quantum Comuting Robert Rogers, Thomas Sylwester, Abe Pauls
More information1-way quantum finite automata: strengths, weaknesses and generalizations
1-way quantum finite automata: strengths, weaknesses and generalizations arxiv:quant-h/9802062v3 30 Se 1998 Andris Ambainis UC Berkeley Abstract Rūsiņš Freivalds University of Latvia We study 1-way quantum
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationarxiv:quant-ph/ v1 14 Apr 2005
Decompositions of general quantum gates Mikko Möttönen and Juha J. Vartiainen arxiv:quant-ph/0504100v1 14 Apr 2005 Abstract Materials Physics Laboratory, POB 2200 (Technical Physics) FIN-02015 HUT, Helsinki
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationarxiv: v2 [quant-ph] 5 Dec 2013
Decomposition of quantum gates Chi Kwong Li and Diane Christine Pelejo Department of Mathematics, College of William and Mary, Williamsburg, VA 23187, USA E-mail: ckli@math.wm.edu, dcpelejo@gmail.com Abstract
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationProof Nets and Boolean Circuits
Proof Nets and Boolean Circuits Kazushige Terui terui@nii.ac.j National Institute of Informatics, Tokyo 14/07/04, Turku.1/44 Motivation (1) Proofs-as-Programs (Curry-Howard) corresondence: Proofs = Programs
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationA MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION
O P E R A T I O N S R E S E A R C H A N D D E C I S I O N S No. 27 DOI:.5277/ord73 Nasrullah KHAN Muhammad ASLAM 2 Kyung-Jun KIM 3 Chi-Hyuck JUN 4 A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationFourier Series Tutorial
Fourier Series Tutorial INTRODUCTION This document is designed to overview the theory behind the Fourier series and its alications. It introduces the Fourier series and then demonstrates its use with a
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
More informationFinding Shortest Hamiltonian Path is in P. Abstract
Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs
More informationCSE 599d - Quantum Computing When Quantum Computers Fall Apart
CSE 599d - Quantum Comuting When Quantum Comuters Fall Aart Dave Bacon Deartment of Comuter Science & Engineering, University of Washington In this lecture we are going to begin discussing what haens to
More informationALTERNATIVE SOLUTION TO THE QUARTIC EQUATION by Farid A. Chouery 1, P.E. 2006, All rights reserved
ALTERNATIVE SOLUTION TO THE QUARTIC EQUATION b Farid A. Chouer, P.E. 006, All rights reserved Abstract A new method to obtain a closed form solution of the fourth order olnomial equation is roosed in this
More informationA Recursive Block Incomplete Factorization. Preconditioner for Adaptive Filtering Problem
Alied Mathematical Sciences, Vol. 7, 03, no. 63, 3-3 HIKARI Ltd, www.m-hiari.com A Recursive Bloc Incomlete Factorization Preconditioner for Adative Filtering Problem Shazia Javed School of Mathematical
More informationEfficient algorithms for the smallest enclosing ball problem
Efficient algorithms for the smallest enclosing ball roblem Guanglu Zhou, Kim-Chuan Toh, Jie Sun November 27, 2002; Revised August 4, 2003 Abstract. Consider the roblem of comuting the smallest enclosing
More informationDeriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.
Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More information1. INTRODUCTION. Fn 2 = F j F j+1 (1.1)
CERTAIN CLASSES OF FINITE SUMS THAT INVOLVE GENERALIZED FIBONACCI AND LUCAS NUMBERS The beautiful identity R.S. Melham Deartment of Mathematical Sciences, University of Technology, Sydney PO Box 23, Broadway,
More informationPHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationPretest (Optional) Use as an additional pacing tool to guide instruction. August 21
Trimester 1 Pretest (Otional) Use as an additional acing tool to guide instruction. August 21 Beyond the Basic Facts In Trimester 1, Grade 8 focus on multilication. Daily Unit 1: Rational vs. Irrational
More informationPell's Equation and Fundamental Units Pell's equation was first introduced to me in the number theory class at Caltech that I never comleted. It was r
Pell's Equation and Fundamental Units Kaisa Taiale University of Minnesota Summer 000 1 Pell's Equation and Fundamental Units Pell's equation was first introduced to me in the number theory class at Caltech
More informationRealization of Two-Qutrit Quantum Gates with Control Pulses
Commun. Theor. Phys. Beijing, China 51 pp. 65 65 c Chinese Physical Society and IOP Publishing Ltd Vol. 51, No., April 15, Realization of Two-Qutrit Quantum Gates with Control Pulses ZHANG Jie, DI Yao-Min,
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationModel checking, verification of CTL. One must verify or expel... doubts, and convert them into the certainty of YES [Thomas Carlyle]
Chater 5 Model checking, verification of CTL One must verify or exel... doubts, and convert them into the certainty of YES or NO. [Thomas Carlyle] 5. The verification setting Page 66 We introduce linear
More informationJosé Alberto Mauricio. Instituto Complutense de Análisis Económico. Universidad Complutense de Madrid. Campus de Somosaguas Madrid - SPAIN
A CORRECTED ALGORITHM FOR COMPUTING THE THEORETICAL AUTOCOVARIANCE MATRICES OF A VECTOR ARMA MODEL José Alberto Mauricio Instituto Comlutense de Análisis Económico Universidad Comlutense de Madrid Camus
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationA New Method of DDB Logical Structure Synthesis Using Distributed Tabu Search
A New Method of DDB Logical Structure Synthesis Using Distributed Tabu Search Eduard Babkin and Margarita Karunina 2, National Research University Higher School of Economics Det of nformation Systems and
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationLIMITATIONS OF RECEPTRON. XOR Problem The failure of the perceptron to successfully simple problem such as XOR (Minsky and Papert).
LIMITATIONS OF RECEPTRON XOR Problem The failure of the ercetron to successfully simle roblem such as XOR (Minsky and Paert). x y z x y z 0 0 0 0 0 0 Fig. 4. The exclusive-or logic symbol and function
More informationBeginnings of the Helicity Basis in the (S, 0) (0, S) Representations of the Lorentz Group
Beginnings of the Helicity Basis in the (S, 0 (0, S Reresentations of the Lorentz Grou Valeriy V. Dvoeglazov UAF, Universidad Autónoma de Zacatecas, México E-mail: valeri@fisica.uaz.edu.mx Abstract We
More informationLilian Markenzon 1, Nair Maria Maia de Abreu 2* and Luciana Lee 3
Pesquisa Oeracional (2013) 33(1): 123-132 2013 Brazilian Oerations Research Society Printed version ISSN 0101-7438 / Online version ISSN 1678-5142 www.scielo.br/oe SOME RESULTS ABOUT THE CONNECTIVITY OF
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationSolving Cyclotomic Polynomials by Radical Expressions Andreas Weber and Michael Keckeisen
Solving Cyclotomic Polynomials by Radical Exressions Andreas Weber and Michael Keckeisen Abstract: We describe a Male ackage that allows the solution of cyclotomic olynomials by radical exressions. We
More informationScaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering
More informationMetrics Performance Evaluation: Application to Face Recognition
Metrics Performance Evaluation: Alication to Face Recognition Naser Zaeri, Abeer AlSadeq, and Abdallah Cherri Electrical Engineering Det., Kuwait University, P.O. Box 5969, Safat 6, Kuwait {zaery, abeer,
More informationAn Overview of Witt Vectors
An Overview of Witt Vectors Daniel Finkel December 7, 2007 Abstract This aer offers a brief overview of the basics of Witt vectors. As an alication, we summarize work of Bartolo and Falcone to rove that
More informationTitle. Author(s)Okamoto, Ryo; Hofmann, Holger F.; Takeuchi, Shigeki; CitationPhysical Review Letters, 95: Issue Date
Title Demonstration of an Otical Quantum Controlled-NOT G Author(s)Okamoto, Ryo; Hofmann, Holger F.; Takeuchi, Shigeki; CitationPhysical Review Letters, 95: 210506 Issue Date 2005-11-18 Doc URL htt://hdl.handle.net/2115/5550
More informationRelaxed p-adic Hensel lifting for algebraic systems
Relaxed -adic Hensel lifting for algebraic systems Jérémy Berthomieu Laboratoire de Mathématiques Université de Versailles Versailles, France jeremy.berthomieu@uvsq.fr Romain Lebreton Laboratoire d Informatique
More informationUncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition
TNN-2007-P-0332.R1 1 Uncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition Haiing Lu, K.N. Plataniotis and A.N. Venetsanooulos The Edward S. Rogers
More informationImpossibility of a Quantum Speed-up with a Faulty Oracle
Imossibility of a Quantum Seed-u with a Faulty Oracle Oded Regev Liron Schiff Abstract We consider Grover s unstructured search roblem in the setting where each oracle call has some small robability of
More informationPROFIT MAXIMIZATION. π = p y Σ n i=1 w i x i (2)
PROFIT MAXIMIZATION DEFINITION OF A NEOCLASSICAL FIRM A neoclassical firm is an organization that controls the transformation of inuts (resources it owns or urchases into oututs or roducts (valued roducts
More informationUniversal Finite Memory Coding of Binary Sequences
Deartment of Electrical Engineering Systems Universal Finite Memory Coding of Binary Sequences Thesis submitted towards the degree of Master of Science in Electrical and Electronic Engineering in Tel-Aviv
More information1 More General von Neumann Measurements [10 Points]
Ph5c Sring 07 Prof. Sean Carroll seancarroll@gmail.com Homework - 3 Solutions Assigned TA: Thom Bohdanowicz thom@caltech.edu Due: 5:00m, 5/03/07 More General von Neumann Measurements [0 Points] Recall
More informationSAT based Abstraction-Refinement using ILP and Machine Learning Techniques
SAT based Abstraction-Refinement using ILP and Machine Learning Techniques 1 SAT based Abstraction-Refinement using ILP and Machine Learning Techniques Edmund Clarke James Kukula Anubhav Guta Ofer Strichman
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationAn Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices
Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationAny AND-OR formula of size N can be evaluated in time N 1/2+o(1) on a quantum computer
Any AND-OR formula of size N can be evaluated in time N /2+o() on a quantum comuter Andris Ambainis, ambainis@math.uwaterloo.ca Robert Šalek salek@eecs.berkeley.edu Andrew M. Childs, amchilds@uwaterloo.ca
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
More informationA Parallel Algorithm for Minimization of Finite Automata
A Parallel Algorithm for Minimization of Finite Automata B. Ravikumar X. Xiong Deartment of Comuter Science University of Rhode Island Kingston, RI 02881 E-mail: fravi,xiongg@cs.uri.edu Abstract In this
More informationq-ary Symmetric Channel for Large q
List-Message Passing Achieves Caacity on the q-ary Symmetric Channel for Large q Fan Zhang and Henry D Pfister Deartment of Electrical and Comuter Engineering, Texas A&M University {fanzhang,hfister}@tamuedu
More informationarxiv: v2 [quant-ph] 17 Nov 2011
Three-layer quantum Kolkata restaurant roblem under decoherence M. Ramzan Deartment of Physics Quaid-i-Azam University Islamabad 4, Pakistan (Dated: Setember 11, 17) arxiv:1111.91v [quant-h] 17 Nov 11
More informationAutomatic Generation and Integration of Equations of Motion for Linked Mechanical Systems
Automatic Generation and Integration of Equations of Motion for Linked Mechanical Systems D. Todd Griffith a, John L. Junkins a, and James D. Turner b a Deartment of Aerosace Engineering, Texas A&M University,
More informationThe Noise Power Ratio - Theory and ADC Testing
The Noise Power Ratio - Theory and ADC Testing FH Irons, KJ Riley, and DM Hummels Abstract This aer develos theory behind the noise ower ratio (NPR) testing of ADCs. A mid-riser formulation is used for
More informationTRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES
Khayyam J. Math. DOI:10.22034/kjm.2019.84207 TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES ISMAEL GARCÍA-BAYONA Communicated by A.M. Peralta Abstract. In this aer, we define two new Schur and
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6.2 Introduction We extend the concet of a finite series, met in Section 6., to the situation in which the number of terms increase without bound. We define what is meant by an infinite
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationThe Arm Prime Factors Decomposition
The Arm Prime Factors Decomosition Arm Boris Nima arm.boris@gmail.com Abstract We introduce the Arm rime factors decomosition which is the equivalent of the Taylor formula for decomosition of integers
More information1 Properties of Spherical Harmonics
Proerties of Sherical Harmonics. Reetition In the lecture the sherical harmonics Y m were introduced as the eigenfunctions of angular momentum oerators lˆz and lˆ2 in sherical coordinates. We found that
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationBlame, coercion, and threesomes: Together again for the first time
Blame, coercion, and threesomes: Together again for the first time Draft, 19 October 2014 Jeremy Siek Indiana University jsiek@indiana.edu Peter Thiemann Universität Freiburg thiemann@informatik.uni-freiburg.de
More informationLecture 1.2 Pose in 2D and 3D. Thomas Opsahl
Lecture 1.2 Pose in 2D and 3D Thomas Osahl Motivation For the inhole camera, the corresondence between observed 3D oints in the world and 2D oints in the catured image is given by straight lines through
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationOn the Toppling of a Sand Pile
Discrete Mathematics and Theoretical Comuter Science Proceedings AA (DM-CCG), 2001, 275 286 On the Toling of a Sand Pile Jean-Christohe Novelli 1 and Dominique Rossin 2 1 CNRS, LIFL, Bâtiment M3, Université
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
More informationD.1 Deutsch-Jozsa algorithm
4 CHAPTER III. QUANTUM COMPUTATION Figure III.: Quantum circuit for Deutsch algorithm. [fig. from Nielsen & Chuang (00)] D Quantum algorithms D. Deutsch-Jozsa algorithm D..a Deutsch s algorithm In this
More informationConstructive quantum scaling of unitary matrices
Quantum Inf Process (016) 15:5145 5154 DOI 10.1007/s1118-016-1448-z Constructive quantum scaling of unitary matrices Adam Glos 1, Przemysław Sadowski 1 Received: 4 March 016 / Accepted: 1 September 016
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More informationHotelling s Two- Sample T 2
Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test
More informationThe Graph Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule
The Grah Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule STEFAN D. BRUDA Deartment of Comuter Science Bisho s University Lennoxville, Quebec J1M 1Z7 CANADA bruda@cs.ubishos.ca
More informationDesign of NARMA L-2 Control of Nonlinear Inverted Pendulum
International Research Journal of Alied and Basic Sciences 016 Available online at www.irjabs.com ISSN 51-838X / Vol, 10 (6): 679-684 Science Exlorer Publications Design of NARMA L- Control of Nonlinear
More informationTopic: Lower Bounds on Randomized Algorithms Date: September 22, 2004 Scribe: Srinath Sridhar
15-859(M): Randomized Algorithms Lecturer: Anuam Guta Toic: Lower Bounds on Randomized Algorithms Date: Setember 22, 2004 Scribe: Srinath Sridhar 4.1 Introduction In this lecture, we will first consider
More informationPrincipal Components Analysis and Unsupervised Hebbian Learning
Princial Comonents Analysis and Unsuervised Hebbian Learning Robert Jacobs Deartment of Brain & Cognitive Sciences University of Rochester Rochester, NY 1467, USA August 8, 008 Reference: Much of the material
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationarxiv: v4 [quant-ph] 8 Mar 2013
Decomposition of unitary matrices and quantum gates Chi-Kwong Li, Rebecca Roberts Department of Mathematics, College of William and Mary, Williamsburg, VA 2387, USA. (E-mail: ckli@math.wm.edu, rlroberts@email.wm.edu)
More informationarxiv: v1 [quant-ph] 20 Jun 2017
A Direct Couling Coherent Quantum Observer for an Oscillatory Quantum Plant Ian R Petersen arxiv:76648v quant-h Jun 7 Abstract A direct couling coherent observer is constructed for a linear quantum lant
More informationQuantum Game Beats Classical Odds Thermodynamics Implications
Entroy 05, 7, 7645-7657; doi:0.3390/e77645 Article OPEN ACCESS entroy ISSN 099-4300 www.mdi.com/journal/entroy Quantum Game Beats Classical Odds Thermodynamics Imlications George Levy Entroic Power Cororation,
More informationarxiv: v1 [quant-ph] 22 Apr 2017
Quaternionic Quantum Particles SERGIO GIARDINO Institute of Science and Technology, Federal University of São Paulo (Unifes) Avenida Cesare G. M. Lattes 101, 147-014 São José dos Camos, SP, Brazil arxiv:1704.06848v1
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
More informationAlgorithms for Air Traffic Flow Management under Stochastic Environments
Algorithms for Air Traffic Flow Management under Stochastic Environments Arnab Nilim and Laurent El Ghaoui Abstract A major ortion of the delay in the Air Traffic Management Systems (ATMS) in US arises
More informationSpecialized and hybrid Newton schemes for the matrix pth root
Secialized and hybrid Newton schemes for the matrix th root Braulio De Abreu Marlliny Monsalve Marcos Raydan June 15, 2006 Abstract We discuss different variants of Newton s method for comuting the th
More informationComputations in Quantum Tensor Networks
Comutations in Quantum Tensor Networks T Huckle a,, K Waldherr a, T Schulte-Herbrüggen b a Technische Universität München, Boltzmannstr 3, 85748 Garching, Germany b Technische Universität München, Lichtenbergstr
More informationObserver/Kalman Filter Time Varying System Identification
Observer/Kalman Filter Time Varying System Identification Manoranjan Majji Texas A&M University, College Station, Texas, USA Jer-Nan Juang 2 National Cheng Kung University, Tainan, Taiwan and John L. Junins
More informationThe University of the State of New York REGENTS HIGH SCHOOL EXAMINATION COURSE III. Wednesday, August 16, :30 to 11:30 a.m.
The University of the State of New York REGENTS HIGH SCHOOL EXAMINATION THREE-YEAR SEQUENCE FOR HIGH SCHOOL MATHEMATICS COURSE III Wednesday, August 6, 000 8:0 to :0 a.m., only Notice... Scientific calculators
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More information