Classical bounds on two-outcome bipartite Bell expressions and linear prepare-and-measure witnesses: Efficient computation in parallel environments such as graphics processing units

Abstract

The presented program aims at speeding up the brute force computation of the so-called $L_d$ norm of a matrix $M$ using graphics processing units (GPUs). Alternatives for CPUs have also been implemented, and the algorithm is applicable to any parallel environment. The $n\times m$ matrix $M$ has real elements which may represent coefficients of a bipartite Bell expression or those of a linear prepare-and-measure (PM) witness. In this interpretation, the $L_1$ norm is the local bound of the given correlation-type Bell expression, and the $L_d$ norm for $d\ge 2$ is the classical $d$-dimensional bound of the given PM witness, which is associated with the communication of $d$-level classical messages in the PM scenario. The program is also capable of calculating the local bound of Bell expressions including marginals. In all scenarios, the output is assumed to be binary. The code for GPUs is written in CUDA C and can utilize one NVIDIA GPU in a computer. To illustrate the performance of our implementation, we refer to Brierley et al. [arXiv:1609.05011] who needed approximately three weeks to compute the local bound on a Bell expression defined by a $42\times 42$ matrix on a standard desktop using a single CPU core. In contrast, our efficient implementation of the brute force algorithm allows us to reduce this to three minutes using a single NVIDIA RTX 6000 Ada graphics card on a workstation. For CPUs, the algorithm was implemented with OpenMP and MPI according to the shared and distributed memory models, respectively, and achieves a comparable speedup at a number of CPU cores around 100.

34 pages, 3 figures, 3 tables

References in corpus (8)

• Testing the Hilbert space dimension
• Certified randomness in quantum physics
• Device-independent tests of classical and quantum dimensions
• Grothendieck's constant and local models for noisy entangled quantum states
• A lower bound on the dimension of a quantum system given measured data
• Efficiency of higher dimensional Hilbert spaces for the violation of Bell inequalities
• Non-classicality thresholds for multiqubit states - numerical analysis
• Qutrit witness from the Grothendieck constant of order four

Cited by in corpus (1)

• Certifying asymmetry in the configuration of three qubits