Next Article in Journal
Graph Partitions in Chemistry
Next Article in Special Issue
Advances in Quantum Computing
Previous Article in Journal
Improving Existing Segmentators Performance with Zero-Shot Segmentators
Previous Article in Special Issue
Principal Component Analysis and t-Distributed Stochastic Neighbor Embedding Analysis in the Study of Quantum Approximate Optimization Algorithm Entangled and Non-Entangled Mixing Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Quantum Convolution for Multidimensional Data

Department of Electrical Engineering and Computer Science, University of Kansas, Lawrence, KS 66045, USA
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(11), 1503; https://doi.org/10.3390/e25111503
Submission received: 18 September 2023 / Revised: 20 October 2023 / Accepted: 27 October 2023 / Published: 31 October 2023
(This article belongs to the Special Issue Advances in Quantum Computing)

Abstract

:
The convolution operation plays a vital role in a wide range of critical algorithms across various domains, such as digital image processing, convolutional neural networks, and quantum machine learning. In existing implementations, particularly in quantum neural networks, convolution operations are usually approximated by the application of filters with data strides that are equal to the filter window sizes. One challenge with these implementations is preserving the spatial and temporal localities of the input features, specifically for data with higher dimensions. In addition, the deep circuits required to perform quantum convolution with a unity stride, especially for multidimensional data, increase the risk of violating decoherence constraints. In this work, we propose depth-optimized circuits for performing generalized multidimensional quantum convolution operations with unity stride targeting applications that process data with high dimensions, such as hyperspectral imagery and remote sensing. We experimentally evaluate and demonstrate the applicability of the proposed techniques by using real-world, high-resolution, multidimensional image data on a state-of-the-art quantum simulator from IBM Quantum.

1. Introduction

Convolution is a common operation that is leveraged in a wide variety of practical applications, such as signal processing [1], image processing [2], and most recently, convolutional neural networks [3]. However, leveraging the widespread utility of convolution operations in quantum algorithms is limited by the lack of a systematic, generalized implementation of quantum convolution. Specifically, contemporary quantum circuits for performing quantum convolution with a given filter are designed on a case-by-case basis [4,5,6,7,8]. In other words, implementing a novel convolution filter on a quantum computer is arduous and time consuming, requiring substantial human effort. Such a workflow is impractical for applications, such as quantum convolutional neural networks, which require a generalized, parameterized quantum circuit to iteratively test thousands of unique filters per training cycle.
In this work, we propose a generalizable algorithm for quantum convolution compatible with amplitude-encoded multidimensional data that is able to implement arbitrary multidimensional filters. Furthermore, our proposed technique implements unity stride, which is essential for capturing the totality of local features in input data. We experimentally verify our technique by applying multiple filters on high-resolution, multidimensional images and report the fidelity of the quantum results against the classically computed expectations. The quantum circuits are implemented on a state-of-the-art quantum simulator from IBM Quantum [9] in both noise-free (as a statevector) and noisy environments. Compared to classical CPU- and GPU-based implementations of convolution, we achieve an exponential improvement in time complexity with respect to data size. Additionally, when compared to existing quantum implementations, we achieve improved circuit depth complexity when factoring in the data encoding.
The work is structured as follows. In Section 2, we cover important background information and review the related work. In Section 3, we introduce the proposed quantum convolution circuits and provide analyses of the corresponding circuit depth. In Section 4, we present the experimental setup and results, while in Section 5, we provide discussions of the results and comparisons to related work. Finally, in Section 6, we present our conclusions and potential avenues for future explorations and extensions.

2. Background

In this section, we discuss related work pertinent to quantum convolution. Quantum operations that are relevant to convolution, such as quantum data encoding, and quantum shift operation, are also presented.

2.1. Related Work

Classically, the convolution operation is implemented either directly or by leveraging fast Fourier transform (FFT). On CPUs, the direct implementation has a time complexity of O ( N 2 ) [10], where N is the data size, while FFT-based implementation has a time complexity of O ( N log N ) [10]. On GPUs, the FFT-based implementation has a similar O ( N log N ) complexity [11]. It is also common to take advantage of the parallelism offered by GPUs to implement convolution using general matrix multiplications (GEMMs) with O ( N F N ) FLOPS [12,13], where N F is the filter size.
Techniques for performing quantum convolution have previously been reported [4,5,6,7,8]. However, these techniques only use fixed sizes of filter windows for specific filters, e.g., edge detection [4,5,6,7,8]. We will denote such methods as fixed-filter quantum convolution. Reportedly, these methods possess a quadratic circuit depth complexity of O ( n 2 ) in terms of the number of qubits n = log 2 N , where N is the size of the input data [4,5,6,7,8]. Because the shortest execution time of classical convolution is in the order of O ( N ) or O ( 2 n ) [12,13] with respect to data size N, authors of the quantum counterparts often claim a quantum advantage [4]. However, the reported depth complexity of fixed-filter quantum convolution does not include the unavoidable overhead of data encoding. Furthermore, there does not exist, to the best of our knowledge, a method for performing generalized, multidimensional quantum convolution.
In reported related work [4,5,6,7,8], data encoding is performed with either the flexible representation of quantum images (FRQI) [14] or novel enhanced quantum representation (NEQR) [15] methods. In these encoding techniques, positional information is stored in the basis quantum states of n qubits, while color information is stored via angle encoding and basis encoding for FRQI and NEQR, respectively. FRQI and NEQR require a total of n + 1 and n + q qubits, respectively, where q is the number of qubits used to represent color values, e.g.,  q = 8 for standard grayscale pixel representation. The reported circuit depth complexities of FRQI and NEQR are O ( 4 n ) and O ( q n 2 n ) , respectively. When factoring in the depth complexities of either data-encoding technique, it is evident that the referenced fixed-filter quantum convolution techniques should be expected to perform worse than classical implementations.
In [16], the authors propose a method of edge detection based on amplitude encoding and the quantum wavelet transform (QWT), which they denote as quantum Hadamard edge detection (QHED). Although the work utilizes grayscale two-dimensional images, the QHED technique is highly customized for those data and does not easily scale or generalize to data of higher dimensions, such as colored and/or multispectral images. For example, the quantum discriminator operation in their technique is applied over all qubits in the circuit, without distinguishing between qubits representing each dimension, i.e., image rows or columns. Such a procedure not only forgoes parallelism and increases circuit depth but inhibits the algorithm’s ability to be generalized beyond capturing one-dimensional features.
In our proposed work, we achieve an exponential improvement in time complexity compared to classical implementations of convolution with respect to data size. Additionally, when compared to existing quantum convolution implementations, we achieve improved circuit depth complexity when factoring in the data encoding. The contribution of our work is analyzed, experimentally verified, and discussed in detail in Section 5.

2.2. Classical to Quantum (C2Q)

Our method of quantum convolution leverages amplitude encoding, which encodes N data values directly in the complex probability amplitudes c i C of the positional basis state | i for an n-qubit state | ψ , where n = log 2 N and 0 i < n , see (1):
| ψ = i = 0 2 n 1 c i | i : c i C .
We use the classical-to-quantum (C2Q) [17] data-encoding technique to encode the amplitude encoded state | ψ 0 from the ground state | 0 n , see Figure 1 and (2). The C2Q operation U C 2 Q has a circuit depth complexity of O ( 2 n ) , a quadratic and linear improvement over FRQI and NEQR, respectively:
U C 2 Q | 0 n = | ψ 0 U C 2 Q = | ψ 0 | × | × , where | × = don t care

2.3. Quantum Shift Operation

A fundamental operation for quantum convolution is the quantum shift operation, denoted in this work as U shift k , which shifts the basis states of the state vector by k positions when applied to an m-qubit state | ψ , see (3). The quantum shift operation is critical for performing the cyclic rotations needed to prepare strided windows when performing convolution. It is also common for the operation to be described as a quantum incrementer when k > 0 , see Figure 2a, and a quantum decrementer when k < 0 [16,18], see Figure 2b:
U shift k | ψ = i = 0 2 m 1 c i | j , where j = ( i k ) mod 2 m

3. Materials and Methods

In general, a convolution operation can be performed using a sequence of shift and multiply-and-accumulate operations. In our proposed methods, we implement the generalized convolution operations as follows:
  • Shift: Auxiliary filter qubits and controlled quantum decrementers are used to create shifted (unity-strided) replicas of input data.
  • Multiply-and-accumulate: Arbitrary state synthesis and classical-to-quantum (C2Q) encoding are applied to create generic multidimensional filters.
  • Data rearrangement: Quantum permutation operations are applied to restructure the fragmented data into one contiguous output datum.
In Section 3.1, we present our quantum convolution technique in detail for one-dimensional data. In the following sections, we illustrate optimizations to improve circuit depth and generalize our method for multidimensional data. For evaluating our proposed methods, we used real-world, high-resolution, black-and-white (B/W) and RGB images, ranging in a resolution from ( 8 × 8 ) pixels to ( 512 × 512 ) pixels and ( 8 × 8 × 3 ) pixels to ( 512 × 512 × 3 ) pixels, respectively. We also performed experiments on 1-D real-world audio data and 3-D real-world hyperspectral data to demonstrate our method’s applicability to data and filters of any dimensionality. Further details about our experimental setup and dataset can be found in Section 4.

3.1. Quantum Convolution for One-Dimensional Data

The proposed structure of quantum convolution for one-dimensional (1-D) data is shown in Figure 3. The following sections show the details of the five steps of the convolution operation procedure to transform the initial encoded data | ψ 0 to the final state | ψ 5 , see Figure 3.

3.1.1. Shift Operation

To perform convolution with unity stride with a filter of size N f terms, N f replicas of the input data must be made, strided for 0 k < N f . To store these replicas, we add n f = log 2 N f auxiliary qubits, which we denote as “filter qubits”, to the most significant positions of the initial quantum state | ψ 0 , see (4) and Figure 3:
| ψ 1 = | 0 n f | ψ 0 = | ψ 0 0 0 | ψ 0 2 n 0 0 2 n + n f
Placing the filter qubits in superposition using Hadamard (H) gates creates 2 n f identical replicas of the initial data | ψ 0 , as shown in (5):
| ψ 2 = H n f I n | ψ 1 = 1 2 n f | ψ 0 | ψ 0 | ψ 0 2 n | ψ 0 2 n 2 n + n f
Finally, multiplexed quantum shift operations can be used to generate the strided replicas, see (6):
| ψ 3 = U mux | ψ 2 = 1 2 n f U shift 0 | ψ 0 U shift ( 2 n f 1 ) | ψ 0 U shift 0 | ψ 0 2 n U shift ( 2 n f 1 ) | ψ 0 2 n 2 n + n f , where U mux = U shift 0 U shift ( 2 n f 1 )

3.1.2. Multiply-and-Accumulate Operation

For the traditional convolution operation, applying a filter F R N f to an array of data W R N f produces a scalar output x R , which can be expressed as x = F T W . In the quantum domain, we can represent an array of data as the partial state | ϕ and the normalized filter as | F . Accordingly, the output state can be expressed as shown in (7):
| ψ out = i = 0 2 n 1 F | ϕ i · | i , where | ϕ i = j = 0 2 n f 1 k | ψ 3 · | k , and k = ( 2 n f · i ) + j
To calculate | ψ 4 from | ψ 3 , it is necessary to embed 〈F| into a unitary operation U F as shown in (8). Since 〈F| is a normalized row vector, we can define U F as a matrix such that its first row is 〈F| and the remaining rows are arbitrarily determined to preserve the unitariness of U F such that U F U F = U F U F = I n f . From (2), we can realize U F using an inverse C2Q operation, see (8):
| ψ 4 = I n U F | ψ 3 = U F | ϕ 0 U F | ϕ 2 n 1 U F | ϕ 0 2 n f U F | ϕ 2 n 1 2 n f 2 n + n f , where U F = F | × | × | 2 n f = U C 2 Q

3.1.3. Data Rearrangement

As of | ψ 4 , the desired values of | ψ out are fragmented among undesired/extraneous values, which we denote using the symbol “×”. We apply SWAP permutations to rearrange and coalesce our desired values to be contiguous in the final statevector | ψ 5 , see (9) and Figure 3:
| ψ 5 = U SWAP 1 D | ψ 4 = F | ϕ 0 F | ϕ 2 n 1 × × F ϕ 0 F ϕ 2 n 1 2 n × × 2 n + n f = | ψ o u t × × | ψ o u t 2 n × × 2 n + n f , where U SWAP 1 D = j = n f 1 0 I ( n f 1 j ) SWAP j , j + n I j

3.1.4. Circuit Depth Analysis of 1-D Quantum Convolution

When considering the circuit depth complexity of the proposed 1-D quantum convolution technique, it is evident from Figure 3 that the operations described by (5) and (9) are performed using parallel Hadamard and SWAP operations, respectively, and thus are of constant depth complexity, i.e.,  O ( 1 ) . In contrast, the  U mux and U F operations incur the largest circuit depth, as they are both serial operations that scale with the data size and/or filter size, see Figure 3.
For the implementation of U mux in Figure 3, there are a total of 2 n f controlled quantum shift operations, where the i-th shift operation is a quantum decrementer U shift i = U shift 1 i . Since all the U shift 1 operations are performed in series, the circuit depth of U mux depends on the total number of unity quantum decrementers, N U shift 1 , see (10):
N U shift 1 ( n f ) = i = 0 2 n f 1 i = 2 n f ( 2 n f 1 ) 2 = 4 n f 2 2 n f 1
As shown in Figure 2b, each quantum decrementer U shift 1 acting on m qubits can be realized using m multi-controlled CNOT (MCX) gates, where the i-th MCX gate is controlled by i qubits and 0 i < m . Accordingly, for each quantum decrementer U shift 1 that is controlled by c qubits, its i-th MCX gate is controlled by a total of i + c qubits. Therefore, the depth of the quantum decrementer circuits can be expressed in terms of the MCX gate count as shown in (11):
D U shift 1 ( m , c ) = i = c m + c 1 D MCX ( i )
The depth of an MCX gate with a total of m qubits can be expressed with a linear function in terms of fundamental single-qubit rotation gates and CNOT gates [19] as shown in (12), where α represents the first-order coefficient and β represents the constant bias term. Thus, the depth complexity of U shift 1 can be expressed as shown in (13):
D MCX ( m ) = α m + β : α , β R
D U shift 1 ( m , c ) = i = 0 m 1 α ( i + c ) + β = α 2 m 2 + α c 1 2 + β m = O m 2
To derive the circuit depth complexity of U mux , we leverage the definitions of N U shift 1 ( n f ) and D U shift 1 ( m , c ) as shown in (14), where m = n and c = n f :
D U mux ( n , n f ) = N U shift 1 ( n f ) · D U shift 1 ( n , n f ) = 4 n f 2 2 n f 1 · α 2 n 2 + α n f 1 2 + β n = 4 n f 1 2 n f 2 · α n 2 + 2 α n f n ( α 2 β ) n = O 4 n f n 2 + 4 n f n f n
As discussed in Section 3.1.2, we implement the filter operation U F by leveraging the C2Q arbitrary synthesis operation [17]. Although C2Q incurs a circuit depth of exponential complexity in terms of fundamental quantum gates, as shown in (15), U F is only applied to n f qubits, a small subset of the total number of qubits, which somewhat mitigates the circuit depth. Furthermore, in most practical scenarios, the dimensions of the filter are typically much smaller than the dimensions of the input data, i.e.,  n f n . As a result, U F should not incur overly large circuit depth relative to other circuit components, e.g.,  U mux . Altogether, the overall circuit depth complexity of the 1-D quantum convolution operation can be expressed according to (16):
D U F ( n f ) = O ( 2 n f )
D 1 - D conv ( n , n f ) = O 4 n f n 2 + 4 n f n f n + 2 n f , where n n f

3.2. Depth-Optimized 1-D Quantum Convolution

In Figure 4, we present an optimized implementation of U mux that greatly reduces the circuit depth.
In Section 3.1, we implemented U mux with 2 n f controlled quantum decrementers U shift k , where 0 k < 2 n f . We can represent each k in binary notation, as shown in (17), to express U shift k as a sequence of controlled shift operations by powers of 2. As shown in (18), we can denote such operations with the notation U shift b j 2 j ( n ) , where 0 j < n f , and ( n ) reflects that the shift operation is applied to an n-qubit state.
k = b n f 1 b n f 2 b j b 1 b 0 2 = j = 0 n f 1 b j 2 j : b j { 0 , 1 }
U shift k ( n ) = U shift j = 0 n f 1 b j 2 j ( n ) = j = 0 n f 1 U shift b j 2 j ( n )
The binary decomposition of the uniformly controlled U shift k operations is conducive to several circuit depth optimizations. As shown in (19), the value of b j is dependent only on the state of the j-th filter qubit q n + j . In other words, each U shift 2 j ( n ) can only be controlled by one qubit q n + j , independently from the other control qubits. Accordingly, it is possible to coalesce the multiplexed U shift 2 j ( n ) operations across the k control indices. The resultant implementation of U mux , therefore, becomes a sequence of 2 n f single-controlled U shift 2 j ( n ) operations, where 0 j < n f , which comparatively has a smaller circuit depth by a factor of 2 n f . Furthermore, each U shift 2 j ( n ) operation can be implemented using a single U shift 1 ( n j ) operation in lieu of sequential U shift 1 ( n ) operations, see (20) and Figure 4, further reducing the depth by a factor of 2 j per operation:
b j = 1 , | q n + j = | 1 0 , otherwise , k [ 0 , 2 n f ]
U shift 2 j ( n ) 2 j U shift 1 ( n ) = U shift 1 ( n j ) I j

Circuit Depth Analysis of Optimized 1-D Quantum Convolution

With the aforementioned optimizations, the depth of the updated U mux operation can be expressed in terms of D U shift 1 ( m , c ) as described by (21), where m = n j and c = 1 for all 0 j < n f . In comparison with the depth complexity of the unoptimized U mux , see (14), the dominant term remains quadratic, i.e.,  n 2 , in terms of the data qubits n. However, its coefficient is improved exponentially, from  4 n f to n f , see (14) and (21). Note that a cubic term n f 3 in terms of the number of filter qubits is introduced in the optimized U mux implementation, see (21). However, when considering the total depth for the optimized 1-D quantum convolution circuit D 1 - D conv opt , the  n f 3 term becomes negligible because of U F , whose complexity O ( 2 n f ) is exponential in terms of the number of filter qubits, see (15), (21), and (22):
D U mux opt ( n , n f ) = j = 0 n f 1 D U shift 1 ( n j , 1 ) = j = 0 n f 1 α 2 ( n j ) 2 + α 2 + β ( n j ) = α 2 n f n 2 α 2 n f 2 n + ( α + β ) n f n + α 6 n f 3 α + β 2 n f 2 + 2 α + 3 β 6 n f = O n f n 2 n f 2 n + n f 3 , where n n f
D 1 - D conv opt ( n , n f ) = O n f n 2 n f 2 n + n f 3 + 2 n f = O n f n 2 n f 2 n + 2 n f , where n n f

3.3. Generalized Quantum Convolution for Multidimensional Data and Filters

In this section, we present the quantum circuit of our proposed quantum convolution technique generalized for multidimensional data and filters. Although quantum statevectors are one-dimensional, it is possible to map multidimensional data to a 1-D vector in either row- or column-major order. In this work, we represent multidimensional input data and convolutional filters in a quantum circuit in column-major order. In other words, for d-dimensional data of size ( N 0 × × N i × N d 1 ) data values, the positional information of the i-th dimension is encoded in the j = 0 i 1 n j to ( j = 0 i n j ) 1 qubits, where n i = log 2 N i . Using this representation, the optimized 1-D quantum convolution circuit shown in Figure 4 can be generalized for d dimensions by “stacking” d 1-D circuits as shown in Figure 5.
The “stacked” quantum circuit in Figure 5 is based on the assumption that the overall (lumped) d-dimensional filter operator U F is separable and decomposable into d one-dimensional filters U F i for 0 i < d . However, it would be more practically useful to generalize our multidimensional quantum convolution technique independently from the separability/decomposability of U F . For this purpose, the identity in (23), which could be easily derived from either Figure 3 or Figure 4 for 1-D convolution, can be leveraged and generalized for multidimensional convolution circuits, see (24). The identity in (24) allows us to reverse the order of multiply-and-accumulate and data rearrangement steps and, therefore, generate one generic lumped U F that acts on the contiguous n f = i = 0 d 1 ( n f i ) filter qubits, where n f i is the number of qubits representing the filter dimension i for 0 i < d , see Figure 6. U F can be derived based on the given arbitrary multidimensional filter F using the method discussed in Section 3.1.2 when F is represented as a normalized 1-D vector | F in a column major ordering:
U SWAP 1 D · I n U F = U F I n · U SWAP 1 D
U SWAP d D · I n f i = d 1 0 I ( n i n f i ) U F i = i = d 1 0 U F i I n · U SWAP d D = U F I n · U SWAP d D , where U F = i = d 1 0 U F i U F d 1 U F d 2 U F 1 U F 0 , U SWAP d D = i = d 1 0 j = n f i 1 0 I n f 1 j q f i SWAP j + q i , j + n + q f i I ( j + q i ) , q f i = k = 0 i 1 n f k , q i = k = 0 i 1 n k , n f = k = 0 d 1 n f k , and n = k = 0 d 1 n k

Circuit Depth Analysis of Generalized Multidimensional Quantum Convolution

As a result of the “stacked" structure, the data of all d dimensions could be concurrently processed in parallel. Consequently, the circuit depth of the multidimensional quantum circuit becomes dependent on the largest data dimension N max = 2 n max , where n max = max i = 0 d 1 ( n i ) , in lieu of the total data size N. The circuit component of the optimized 1-D circuit with the greatest depth contribution U mux is performed in parallel on each dimension in the generic d-D circuit. Specifically, U mux scales with the number of qubits used to represent the largest data dimension n max = log 2 N max . Note that the parallelization across dimensions applies to the Hadamard and SWAP operations from (5) and (9), see Figure 6, and therefore these operations are of constant depth complexity, i.e.,  O ( 1 ) . The depth complexity of the multidimensional U F operation is determined by the total number of elements in the filter N F , and therefore the C2Q-based implementation of U F does not benefit from multidimensional stacking. Accordingly, the circuit depth of the generalized multidimensional quantum convolution operation could be derived from (22) and expressed in (25), where n f max = max i = 0 d 1 ( n f i ) is the number of qubits representing the maximum filter dimension N f max = 2 n f max . It is worth mentioning that the generic multidimensional formula in (25) reduces to the 1-D formula in (22) when n = n max and n f = n f max :
D d - D conv opt ( n , n f ) = O n f max n max 2 n f max 2 n max + 2 n f , where n max = max i = 0 d 1 ( n i ) , n f max = max i = 0 d 1 ( n f i ) , and n max n f max

4. Experimental Setup and Results

We experimentally demonstrate our proposed technique for generalized, multidimensional quantum convolution with unity stride on real-world, high-resolution, multidimensional image data, see Figure 7. By leveraging the Qiskit SDK (v0.39.4) from IBM Quantum [9], we simulate our quantum circuits in the following formats: (1) classically (to present the ideal/theoretical expectation), (2) noise-free (using statevector simulation), and (3) noisy (using 1,000,000 “shots” or circuit samples). Moreover, we present a quantitative comparison of the obtained results using fidelity [20] as a similarity metric between compared results ρ and σ , see (26). Experiments were performed on a 16-core AMD EPYC 7302P CPU with frequencies up to 3.3 GHz, 128 MB of L3 cache, and access to 256 GB of DDR4 RAM. In our analysis, we evaluated the correctness of the proposed techniques by comparing classical results with noise-free results. We also evaluated the scalability of the proposed techniques for higher-resolution images by comparing the classical results with both the noise-free and noisy results. Finally, we plotted the circuit depth of our techniques with respect to the data size and filter size as shown in Figure 8.
Fidelity ( ρ , σ ) = ρ · σ · ρ
In our experiments, we evaluated our techniques using well-known ( 3 × 3 ) and ( 5 × 5 ) filters, i.e., Averaging F avg , Sobel edge-detection F Sx / F Sy , Gaussian blur F blur , and Laplacian of Gaussian blur (Laplacian) F L , see (27)–(30). We applied zero padding to maintain the size of the filter dimensions at a power of two for quantum implementation. In addition, we used wrapping to resolve the boundary conditions, and we restricted the magnitude of the output between [ 0 , 255 ] to mitigate quantization errors in the classical domain:
F avg 3 × 3 = 1 9 1 1 1 1 1 1 1 1 1 , F avg 5 × 5 = 1 25 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
F blur 3 × 3 = 1 16 1 2 1 2 4 2 1 2 1 , F blur 5 × 5 = 1 273 1 4 7 4 1 4 16 26 16 4 7 26 41 26 7 4 16 26 16 4 1 4 7 4 1
F Sx = 1 4 1 0 1 2 0 2 1 0 1 , F Sy = 1 4 1 2 1 0 0 0 1 2 1
F L 3 × 3 = 1 6 1 1 1 1 8 1 1 1 1 , F L 5 × 5 = 1 20 1 1 1 1 1 1 1 1 1 1 1 1 24 1 1 1 1 1 1 1 1 1 1 1 1
We applied 2-D convolution filters to black-and-white (B/W) and RGB images, see Figure 7, ranging in resolution from ( 8 × 8 ) pixels to ( 512 × 512 ) pixels and ( 8 × 8 × 3 ) pixels to ( 512 × 512 × 3 ) pixels, respectively. The number of filter qubits can be obtained by the size of filter dimensions, i.e.,  n f = log 2 3 + log 2 3 = 4 qubits for ( 3 × 3 ) filters and n f = log 2 5 + log 2 5 = 6 qubits for ( 5 × 5 ) filters. Therefore, our simulated quantum circuits ranged in size ( n + n f ) from 10 qubits to 26 qubits. Figure 9 and Figure 10 present the reconstructed output images for classical, noise-free, and noisy experiments using ( 128 × 128 ) and ( 128 × 128 × 3 ) -pixel input images, respectively.
To demonstrate our method’s applicability to data and filters of any dimensionality, we also performed experiments applying the 1-D and 3-D averaging filter to 1-D real-world audio data and 3-D real-world hyperspectral data, respectively. The audio files were sourced from the publicly available sound quality assessment material published by the European Broadcasting Union and modified to be single channel with data sizes ranging from 2 8 values to 2 20 values when sampled at 44.1 kHz [21]. Figure 11 and Figure 12 present the reconstructed output images and fidelity, respectively, from applying (3) and (5) averaging filters. The hyperspectral images were sourced from the Kennedy Space Center (KSC) dataset [22] and resized to range from ( 8 × 8 × 8 ) pixels to ( 128 × 128 × 128 ) pixels. Figure 13 and Figure 14 present the reconstructed output images and fidelity, respectively, from applying ( 3 × 3 × 3 ) and ( 5 × 5 × 5 ) averaging filters.
Comparison of the noise-free quantum results against the ideal classical results demonstrates a 100 % fidelity across all trials. Thus, in a noise-free (statevector) environment, our proposed quantum convolution technique correctly performs an identical operation to classical convolution given the same input parameters and boundary conditions.
When considering the behavior of noisy (sampled) environments, Figure 12, Figure 14 and Figure 15 plot the fidelity of the noisy quantum results against the ideal classical results. We observe a monotonic decrease in fidelity as the data size (image resolution) increases, consistent with previously reported behavior [23]. Such behavior derives from how the number of shots required to properly characterize a quantum state increases with the corresponding number of qubits in order to reduce the effects of statistical noise. Notably, the fidelity varies dramatically depending on the filter category and size, where the largest discrepancy occurs between the ( 5 × 5 ) Averaging and ( 5 × 5 ) Laplacian filters for a data size of 65 , 536 values. Specifically, for a B/W image of ( 256 × 256 ) pixels, the Averaging filter had a fidelity of 94.84%, while the Laplacian filter had a fidelity of 8.82%—a difference of 86.02%, see Figure 15a. In general, we observed that the average and blur filters perform better than the edge detection methods (Sobel/Laplacian). Since the output data are reconstructed from only a portion of the final state | ψ 5 , it is likely that sparse filters, represented in Figure 9 and Figure 10 as being mostly black, are significantly less likely to be recorded during sampled measurement, resulting in reduced fidelity. For practical applications, dimension reduction techniques, such as pooling, can be used to mitigate information loss [23].

5. Discussion

In the following section, we compare our proposed method of quantum convolution to the related work discussed in Section 2.1 in terms of filter generalization and circuit depth.

5.1. Arbitrary Multidimensional Filtering

Our generalizable and parameterized technique of quantum convolution with unity stride offers distinct workflow advantages over existing fixed-filter quantum convolution techniques in variational applications, such as quantum machine learning. Our technique does not require extensive development for each new filter design. For instance, current quantum convolutional filters are primarily two dimensional, only focusing on image processing. However, the development of even similar filters targeting audio and video processing, for example, would require extensive development and redesign. Our method offers a systematic and straightforward approach for generating practical quantum circuits given fundamental input variables.

5.2. Circuit Depth

Figure 8a,b show the circuit depth for our proposed technique of generalized quantum convolution with respect to the total number of data qubits n = i = 0 d 1 ( n i ) and the total number of filter qubits n f = i = 0 d 1 ( n f i ) , respectively. The results were gathered using the QuantumCircuit.depth() method built into Qiskit for a QuantumCircuit transpiled to fundamental single-qubit and CNOT quantum gates. Figure 8a illustrates quadratic circuit depth complexity with respect to the data qubits n for a fixed filter size N F = 2 n f , aligning with our theoretical expectation in (25). Note that n = i = 0 d 1 n i and n max = max i = 0 d 1 ( n i ) for d-dimensional data. Similarly, Figure 8b (plotted on a log-scale) illustrates exponential circuit depth complexity with respect to n f for a fixed data size N = 2 n , which also aligns with our theoretical expectation in (25).
The time complexity comparison of our proposed quantum convolution technique against related work is shown in Table 1. Compared to classical direct implementations on CPUs, our proposed technique for generalized quantum convolution demonstrates an exponential improvement with respect to data size N = 2 n , i.e.,  O ( n 2 ) vs. O ( N 2 ) , see (25) and Table 1c. As discussed in Section 2.1, the fastest classical GEMM implementation of convolution on GPUs (excluding data I/O overhead) [12,13] has a complexity of O ( N F N ) , see Table 1c. Even when including quantum data encoding, which is equivalent to data I/O overhead, our method remains to demonstrate a linear improvement with respect to data size N by a factor of the filter size N F , see (31), over classical GEMM GPUs:
D proposed ( n ) = D C 2 Q ( n ) + D d - D conv opt ( n ) = O 2 n + n max 2 = O 2 n = O N , for fixed n f
Compared to fixed-filter quantum convolution techniques [4,5,6,7,8], our proposed arbitrary filter quantum convolution technique (for unity stride) demonstrates improved circuit depth complexity with respect to data size when factoring in the circuit depth contribution from data encoding. For a fixed filter, i.e.,  n f is constant, the depth of the proposed method scales quadratically with the largest data dimension n max , see (25) and (31). As described in Section 2.1, fixed-filter quantum convolution techniques similarly show quadratic depth scaling with respect to the number of qubits n, see Table 1b. For data encoding, the fixed-filter techniques use either the FRQI [14] or NEQR [15] algorithms, which have circuit depth complexities of O ( 4 n ) or O ( q n 2 n ) , respectively. In contrast, our proposed technique uses C2Q data encoding [17], which has a depth complexity of O ( 2 n ) —a quadratic and linear improvement over FRQI and NEQR, respectively, see Table 1a.

6. Conclusions

In this work, we proposed and evaluated a method for generalizing the convolution operation with arbitrary filtering and unity stride for multidimensional data in the quantum domain. We presented the corresponding circuits and their performance analyses along with experimental results that were collected using the IBM Qiskit development environment. In our experimental work, we validated the correctness of our method by comparing classical results to noise-free quantum results. We also demonstrated the practicality of our method for various convolution filters by evaluating the noisy quantum results. Furthermore, we presented experimentally verified analyses that highlight our technique’s advantages in terms of time complexity and/or circuit depth complexity compared to existing classical and quantum methods, respectively. Future work will focus on adapting our proposed technique for arbitrary strides. In addition, we will investigate multidimensional quantum machine learning as a potential application of our proposed technique.

Author Contributions

Conceptualization: M.J., V.J., D.K. and  E.E.-A.; Methodology: M.J., V.J., D.K. and  E.E.-A.; Software: M.J., D.L., D.K. and  E.E.-A.; Validation: M.J., V.J., D.L., D.K. and E.E.-A.; Formal analysis: M.J., V.J., D.L., D.K. and E.E.-A.; Investigation: M.J., A.N., V.J., D.L., D.K., M.C., I.I., M.M.R. and E.E.-A.; Resources: M.J., A.N., V.J., D.L., D.K., M.C., I.I., M.M.R. and  E.E.-A.; Data curation: M.J., A.N., V.J., D.L., D.K., M.C., I.I., M.M.R. and  E.E.-A.; Writing—original draft preparation: M.J., A.N., V.J., D.L., D.K., M.C., I.I. and  E.E.-A.; Writing—review and editing: M.J., A.N., V.J., D.L., D.K., M.C., I.I., M.M.R. and  E.E.-A.; Visualization: M.J., A.N., V.J., D.L., D.K., M.C., I.I., M.M.R. and E.E.-A.; Supervision: E.E.-A.; Project administration: E.E.-A.; Funding acquisition: E.E.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The audio samples used in this work are publicly available from the European Broadcasting Union at https://tech.ebu.ch/publications/sqamcd (accessed on 19 October 2023) as file 64.flac [21]. The hyperspectral data used in this work are publicly available from the Grupo de Inteligencia Computacional (GIC) at https://www.ehu.eus/ccwintco/index.php/Hyperspectral_Remote_Sensing_Scenes#Kennedy_Space_Center_(KSC) (accessed on 19 October 2023) under the heading Kennedy Space Center (KSC) [22].

Acknowledgments

This research used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-00OR22725.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
1-DOne Dimensional
C2QClassical to Quantum
FFTFast Fourier Transform
FRQIFlexible Representation of Quantum Images
GEMMGeneral Matrix Multiplication
NEQRNovel Enhanced Quantum Representation
QHEDQuantum Hadamard Edge Detection
QWTQuantum Wavelet Transform

References

  1. Rhodes, W. Acousto-optic signal processing: Convolution and correlation. Proc. IEEE 1981, 69, 65–79. [Google Scholar] [CrossRef]
  2. Keys, R. Cubic convolution interpolation for digital image processing. IEEE Trans. Acoust. Speech Signal Process. 1981, 29, 1153–1160. [Google Scholar] [CrossRef]
  3. LeCun, Y.; Kavukcuoglu, K.; Farabet, C. Convolutional networks and applications in vision. In Proceedings of the 2010 IEEE International Symposium on Circuits and Systems, Paris, France, 30 May–2 June 2010; pp. 253–256. [Google Scholar] [CrossRef]
  4. Fan, P.; Zhou, R.G.; Hu, W.; Jing, N. Quantum image edge extraction based on classical Sobel operator for NEQR. Quantum Inf. Process. 2018, 18, 24. [Google Scholar] [CrossRef]
  5. Ma, Y.; Ma, H.; Chu, P. Demonstration of Quantum Image Edge Extration Enhancement Through Improved Sobel Operator. IEEE Access 2020, 8, 210277–210285. [Google Scholar] [CrossRef]
  6. Zhang, Y.; Lu, K.; Gao, Y. QSobel: A novel quantum image edge extraction algorithm. Sci. China Inf. Sci. 2015, 58, 1–13. [Google Scholar] [CrossRef]
  7. Zhou, R.G.; Yu, H.; Cheng, Y.; Li, F.X. Quantum image edge extraction based on improved Prewitt operator. Quantum Inf. Process. 2019, 18, 261. [Google Scholar] [CrossRef]
  8. Li, P. Quantum implementation of the classical Canny edge detector. Multimed. Tools Appl. 2022, 81, 11665–11694. [Google Scholar] [CrossRef]
  9. IBM Quantum. Qiskit: An Open-source Framework for Quantum Computing. 2021. Available online: https://zenodo.org/records/7416349 (accessed on 19 October 2023).
  10. Burrus, C.S.; Parks, T.W. DFT/FFT and Convolution Algorithms: Theory and Implementation, 1st ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1991. [Google Scholar]
  11. Podlozhnyuk, V. FFT-based 2D convolution. NVIDIA 2007, 32. Available online: https://developer.download.nvidia.com/compute/cuda/1.1-Beta/x86_64_website/projects/convolutionFFT2D/doc/convolutionFFT2D.pdf (accessed on 19 October 2023).
  12. NVIDIA. Convolution Algorithms. 2023. Available online: https://docs.nvidia.com/deeplearning/performance/dl-performance-convolutional/index.html#conv-algo (accessed on 19 October 2023).
  13. NVIDIA. CUTLASS Convolution. 2023. Available online: https://github.com/NVIDIA/cutlass/blob/main/media/docs/implicit_gemm_convolution.md (accessed on 19 October 2023).
  14. Le, P.Q.; Dong, F.; Hirota, K. A flexible representation of quantum images for polynomial preparation, image compression, and processing operations. Quantum Inf. Process. 2011, 10, 63–84. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Lu, K.; Gao, Y.; Wang, M. NEQR: A novel enhanced quantum representation of digital images. Quantum Inf. Process. 2013, 12, 2833–2860. [Google Scholar] [CrossRef]
  16. Yao, X.W.; Wang, H.; Liao, Z.; Chen, M.C.; Pan, J.; Li, J.; Zhang, K.; Lin, X.; Wang, Z.; Luo, Z.; et al. Quantum Image Processing and Its Application to Edge Detection: Theory and Experiment. Phys. Rev. X 2017, 7, 031041. [Google Scholar] [CrossRef]
  17. El-Araby, E.; Mahmud, N.; Jeng, M.J.; MacGillivray, A.; Chaudhary, M.; Nobel, M.A.I.; Islam, S.I.U.; Levy, D.; Kneidel, D.; Watson, M.R.; et al. Towards Complete and Scalable Emulation of Quantum Algorithms on High-Performance Reconfigurable Computers. IEEE Trans. Comput. 2023, 72, 2350–2364. [Google Scholar] [CrossRef]
  18. Li, X.; Yang, G.; Torres, C.M.; Zheng, D.; Wang, K.L. A Class of Efficient Quantum Incrementer Gates for Quantum Circuit Synthesis. Int. J. Mod. Phys. B 2014, 28, 1350191. [Google Scholar] [CrossRef]
  19. Balauca, S.; Arusoaie, A. Efficient Constructions for Simulating Multi Controlled Quantum Gates. In Proceedings of the Computational Science—ICCS 2022, London, UK, 21–23 June 2022; Groen, D., de Mulatier, C., Paszynski, M., Krzhizhanovskaya, V.V., Dongarra, J.J., Sloot, P.M.A., Eds.; Springer: Cham, Switzerland, 2022; pp. 179–194. [Google Scholar]
  20. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information: 10th Anniversary Edition; Cambridge University Press: Cambridge, UK, 2010; p. 409. [Google Scholar] [CrossRef]
  21. Geneva, S. Sound Quality Assessment Material: Recordings for Subjective Tests. 1988. Available online: https://tech.ebu.ch/publications/sqamcd (accessed on 19 October 2023).
  22. Graña, M.; Veganzons, M.A.; Ayerdi, B. Hyperspectral Remote Sensing Scenes. Available online: https://www.ehu.eus/ccwintco/index.php/Hyperspectral_Remote_Sensing_Scenes#Kennedy_Space_Center_(KSC) (accessed on 19 October 2023).
  23. Jeng, M.; Islam, S.I.U.; Levy, D.; Riachi, A.; Chaudhary, M.; Nobel, M.A.I.; Kneidel, D.; Jha, V.; Bauer, J.; Maurya, A.; et al. Improving quantum-to-classical data decoding using optimized quantum wavelet transform. J. Supercomput. 2023, 79, 20532–20561. [Google Scholar] [CrossRef]
Figure 1. Quantum circuit for classical-to-quantum (C2Q) arbitrary state synthesis [17].
Figure 1. Quantum circuit for classical-to-quantum (C2Q) arbitrary state synthesis [17].
Entropy 25 01503 g001
Figure 2. Quantum shift operation using quantum incrementers/decrementers.
Figure 2. Quantum shift operation using quantum incrementers/decrementers.
Entropy 25 01503 g002
Figure 3. The 1-D quantum convolution circuit.
Figure 3. The 1-D quantum convolution circuit.
Entropy 25 01503 g003
Figure 4. Depth-optimized 1-D quantum convolution circuit.
Figure 4. Depth-optimized 1-D quantum convolution circuit.
Entropy 25 01503 g004
Figure 5. Multidimensional quantum convolution circuit with distributed/stacked 1-D filters.
Figure 5. Multidimensional quantum convolution circuit with distributed/stacked 1-D filters.
Entropy 25 01503 g005
Figure 6. Generalized multidimensional quantum convolution circuit.
Figure 6. Generalized multidimensional quantum convolution circuit.
Entropy 25 01503 g006
Figure 7. Real-world, high-resolution, multidimensional images used in experimental trials.
Figure 7. Real-world, high-resolution, multidimensional images used in experimental trials.
Entropy 25 01503 g007
Figure 8. Circuit depth of quantum convolution with respect to data and filter qubits.
Figure 8. Circuit depth of quantum convolution with respect to data and filter qubits.
Entropy 25 01503 g008
Figure 9. The 2-D convolution filters applied to a ( 128 × 128 ) B/W image.
Figure 9. The 2-D convolution filters applied to a ( 128 × 128 ) B/W image.
Entropy 25 01503 g009
Figure 10. The 2-D convolution filters applied to a ( 128 × 128 × 3 ) RGB image.
Figure 10. The 2-D convolution filters applied to a ( 128 × 128 × 3 ) RGB image.
Entropy 25 01503 g010
Figure 11. The 1-D convolution (averaging) filters applied to 1-D audio samples.
Figure 11. The 1-D convolution (averaging) filters applied to 1-D audio samples.
Entropy 25 01503 g011
Figure 12. Fidelity of 1-D convolution (averaging) filters with unity stride on 1-D audio data (sampled with 1,000,000 shots).
Figure 12. Fidelity of 1-D convolution (averaging) filters with unity stride on 1-D audio data (sampled with 1,000,000 shots).
Entropy 25 01503 g012
Figure 13. The 3-D convolution (averaging) filters applied to 3-D hyperspectral images (bands 0, 1, and 2).
Figure 13. The 3-D convolution (averaging) filters applied to 3-D hyperspectral images (bands 0, 1, and 2).
Entropy 25 01503 g013
Figure 14. Fidelity of 3-D convolution (averaging) filters with unity stride on 3-D hyperspectral data (sampled with 1,000,000 shots).
Figure 14. Fidelity of 3-D convolution (averaging) filters with unity stride on 3-D hyperspectral data (sampled with 1,000,000 shots).
Entropy 25 01503 g014
Figure 15. Fidelity of 2-D convolution with unity stride (sampled with 1,000,000 shots).
Figure 15. Fidelity of 2-D convolution with unity stride (sampled with 1,000,000 shots).
Entropy 25 01503 g015
Table 1. Comparison of depth/time complexity of proposed generalized quantum convolution technique against related work.
Table 1. Comparison of depth/time complexity of proposed generalized quantum convolution technique against related work.
a Depth complexity of quantum data encoding (I/O) techniques
FRQI [14]NEQR [15]C2Q [17]
O I / O ( 4 n ) O I / O ( q n 2 n ) O I / O ( 2 n )
b Depth complexity of quantum convolution algorithms for a fixed filter
ProposedRelated Work [4,5,6,7,8]
O alg ( n max 2 ) O alg ( n 2 )
c Complexity of proposed technique compared to classical convolution
ProposedDirect (CPU) [10]FFT (CPU/GPU) [10,11]GEMM (GPU) [12,13]
O alg n f max n max 2 n f max 2 n max + 2 n f O alg + I / O n f max n max 2 n f max 2 n max + 2 n + 2 n f O alg N 2 O alg 4 n O alg N log N O alg n 2 n O alg N F N O alg 2 ( n + n f )
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jeng, M.; Nobel, A.; Jha, V.; Levy, D.; Kneidel, D.; Chaudhary, M.; Islam, I.; Rahman, M.M.; El-Araby, E. Generalized Quantum Convolution for Multidimensional Data. Entropy 2023, 25, 1503. https://doi.org/10.3390/e25111503

AMA Style

Jeng M, Nobel A, Jha V, Levy D, Kneidel D, Chaudhary M, Islam I, Rahman MM, El-Araby E. Generalized Quantum Convolution for Multidimensional Data. Entropy. 2023; 25(11):1503. https://doi.org/10.3390/e25111503

Chicago/Turabian Style

Jeng, Mingyoung, Alvir Nobel, Vinayak Jha, David Levy, Dylan Kneidel, Manu Chaudhary, Ishraq Islam, Muhammad Momin Rahman, and Esam El-Araby. 2023. "Generalized Quantum Convolution for Multidimensional Data" Entropy 25, no. 11: 1503. https://doi.org/10.3390/e25111503

APA Style

Jeng, M., Nobel, A., Jha, V., Levy, D., Kneidel, D., Chaudhary, M., Islam, I., Rahman, M. M., & El-Araby, E. (2023). Generalized Quantum Convolution for Multidimensional Data. Entropy, 25(11), 1503. https://doi.org/10.3390/e25111503

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop