Brevet FR3072483A1 — publié le 20190419

OPTIMISATION D'UN CIRCUIT QUANTIQUE PAR INSERTION DE PORTES SWAP DOMAINE DE L'INVENTION La présente invention est relative au domaine des calculateurs quantiques. Plus précisément elle concerne l'optimisation d'un circuit quantique afin que l'ensemble des portes qui le composent soient locales. CONTEXTE DE L'INVENTION Un circuit quantique utilise les propriétés quantiques de la matière, telles que la superposition des états, et l'intrication, afin d'effectuer des opérations sur les données. À la différence d'un ordinateur conventionnel qui manipule des données binaires, ou « bits », un circuit quantique travaille sur des qubits, ou bits quantiques, qui représentent un état quantique basé sur deux états de base : 0> et 1> Des premiers algorithmes pour calculateur quantique sont apparus dans les années 1990. Un algorithme particulièrement emblématique et à l'origine d'un certain engouement pour les calculateurs quantiques et l'algorithme de Shor pour le calcul polynomial, décrit dans P. W. Shor, « Algorithm for quantum computation: discrete logarithms and factoring », in Foundations of Computer Science, pages 124-134, 1994. D'une façon générale, un circuit quantique peut être considéré comme une succession d'opérateurs, ou portes, quantiques, appliqués à un ensemble de qubit. La figure 1 illustre un circuit quantique simple, comportant 4 qubits ql, q2, q3, q4 agencés en ligne. De façon classique, on représente les qubits 2 sous forme de lignes horizontales. La direction allant de la gauche vers la droite correspond au temps. Sur cette figure, 3 portes quantiques Gl, G2, G3 composent le circuit et s'appliquent donc dans le temps selon une succession G1 puis G2 puis G3. 5 Une porte quantique peut impliquer un ou plusieurs qubits. Elles peuvent être de différentes natures et chacune représente un coût lié au temps de propagation du signal ou temps de passage de la porte. Parmi les portes classiquement utilisées, on peut citer la porte « swap » G1 qui intervertie 10 deux valeurs de qubit. L'agencement d'un circuit quantique implémentant un algorithme résulte souvent de procédé de génération plus au moins automatiques. Le circuit quantique est alors décrit de façon indépendante de l'infrastructure 15 physique sur lequel il est amené à être déployé. Notamment, ces circuits quantiques « théoriques » supposent que l'ensemble des qubits peuvent interagir les uns avec les autres. Ainsi, dans l'exemple de la figure 1, la porte quantique G2 fait interagir les qubits ql et 20 q4, alors que ceux-ci ne sont pas voisins, mais séparés par les qubits q2 et q3. Or, les développements des infrastructures physiques pour calculateurs quantiques ont fait apparaître des contraintes physiques. Parmi ces contraintes, une contrainte dite « de localité » fait que la plupart des 25 infrastructures ne permettent des interactions qu'entre qubits proches. Cette contrainte de localité est imposée par les mécanismes physiques utilisés pour mettre en oeuvre les qubits et les portes quantiques. Ainsi, une porte quantique telle que la porte G2 de la figure 1 ne peut donc pas être déployée directement sur ce type d'architecture. 30 3 Certains algorithmes classiques, tels que l'algorithme de Shor, peuvent faire l'objet d'optimisation permettant de les mettre en oeuvre sur des architectures soumises à la contrainte de localité. Mais en général, comme dit précédemment, les agencements des circuits quantiques sont synthétisés, 5 automatiquement ou non, sans prendre en compte cette contrainte. Il apparaît donc un besoin de méthodes et d'outils pour optimiser un tel circuit, c'est-à-dire pour transformer un circuit indépendant de l'infrastructure physique en un circuit respectant la contrainte de localité et fonctionnellement équivalent. 10 Plusieurs méthodes ont été proposées qui, de façon générale, reposent sur l'insertion de portes « SWAP » locales. En insérant des portes Swap en amont d'une porte quantique donnée, on peut intervertir des qubits de sorte à rendre locale cette porte quantique. Les figures 2a, 2b, 2c illustrent cette problématique. La figure 2a 15 représente un circuit quantique non-optimisé pour un réseau de 4 qubits ql, q2, q3, q4 disposés en ligne. On remarque plusieurs portes quantiques non-locales, c'est-à-dire faisant interagir des qubits non-adjacents. Les portes quantiques sont représentées par des segments verticaux comportant des jonctions avec les lignes de qubits sur lesquels elles agissent. Le graphisme 20 de ces jonctions (points, cercles...) correspond au type de porte quantique et au type d'action sur ces qubits. Les figures 2b, 2c représentent deux optimisations possibles de ce circuit par insertions de portes quantiques Swap rendant chacune des portes locales. Les portes Swap sont représentés avec des croix comme bornes. 25 Toutefois, chaque insertion de porte quantique Swap représente un coût, lié notamment au temps de traversé de la porte. Il importe donc de déterminer, parmi les circuits équivalents optimisé, celui qui minimise ce coût supplémentaire. 4 Trouver une solution optimale à ce problème n'est pas envisageable car il y a un nombre exponentiel, en nombres de qubits et de portes, de circuits fonctionnellement équivalents et respectant la contrainte de localité. 5 Différentes solutions à ce problème ont été proposées. On peut par exemple citer la thèse de master « An efficient method to convert arbitrary quantum circuits to ones on a Linear Nearest Neighbor architecture » de Yuichi Hirata, soutenue auprès de l'Institut de Sciences et de Technologies de Nara, Japon. On peut également citer l'article « Depth-optimal Quantum 10 Circuit Placement for Arbitrary Topologies » de Bhattacharjee Debyoti et Chattopadhyay Anupam, du Nanyang Technological University, 2017. On peut aussi citer l'article « Optimal SWAP Gate Insertion for Nearest Neighbor Quantum Circuits » de Robert Wille, Aaron Lye et Rolf Drechsler, in Design Automation Conference (ASP-DAC), 2014 19th Asia and South 15 Pacific Mais ces solutions sont très lentes à trouver une solution puisqu'elles recherchent la solution globalement optimale. Dès que le nombre de qubits et de portes devient important, dans le cadre d'un algorithme complexe, la 20 détermination du circuit quantique représentant un optimum global n'est plus possible, et ces solutions de l'état de la technique deviennent inutilisables en pratique. D'autres familles de solutions visent à résoudre le problème dans un 25 cadre plus encadré d'une infrastructure physique donnée, afin d'en limiter la combinatoire. L'invention vise à proposer une solution qui puisse à la fois s'appliquer à des infrastructures physiques variées et atteindre des performances 5 permettant son exploitation effective pour des grands nombres de qubits et de portes quantiques. RE SUME DE L'INVENTION 5 Le but de la présente invention est de fournir une solution palliant au moins partiellement les inconvénients précités. A cette fin, la présente invention propose un procédé d'optimisation 10 d'un circuit quantique composée d'une série ordonnée de portes quantiques appliqué à un agencement initial de valeurs de qubits, consistant à insérer un ensemble de portes SWAP locales, afin que l'ensemble des portes dudit circuit quantique soient locales, comprenant des étapes consistant à : - pour chaque porte quantique de ladite série, si celle-ci n'est pas 15 locale, insérer un ensemble de porte SWAP locales avant ladite porte quantique, - déterminer l'ensemble de permutations, chacune constituée d'une succession d'interversions de valeurs de qubits le long des plus courts chemins entre les positions de qubits associés à ladite porte 20 quantique, et - choisir parmi ledit ensemble de permutations, une permutation minimisant un coût représentant le nombre d'interversions nécessaires pour rendre locales les portes quantiques d'une suite au sein de ladite série, la taille de ladite suite étant 25 substantiellement inférieure à celle de ladite série; - rétablir ledit agencement initial en établissant un arbre couvrant d'un graphe représentatif de l'agencement des qubits dudit circuit, et en effectuant des interversions de valeurs de qubits le long de chemins dudit arbre. 30 6 Suivant des modes de réalisation préférés, l'invention comprend une ou plusieurs des caractéristiques suivantes qui peuvent être utilisées séparément ou en combinaison partielle entre elles ou en combinaison totale entre elles : - la détermination de l'ensemble des permutations est effectuée en 5 déterminant une succession d'interversions de valeurs de qubits le long des premiers plus courts chemins entre une paire quelconque de positions de qubits associés à ladite porte quantique, puis, si ladite porte quantique possède plus de deux positions de qubits associés, en déterminant une succession d'interversions de valeurs 10 de qubits le long des deuxièmes plus courts chemins entre chacune des positions en dehors de celles déjà considérées, et les premiers plus courts chemins déjà déterminés ; - ledit coût prend en compte un coût attaché à chaque type de portes SWAP à insérer ; 15 - le rétablissement de l'agencement initial est effectuée selon un l'algorithme d'échanges de jetons ; - ledit algorithme d'échanges de jetons est celui décrit au paragraphe 3.2 de l'article « Swapping Labeled Tokens on Graph » de K. Yamanaka, E. D. Demaine, T. Ito, J. Kawahara, M. 20 Kiyomi, Y. Okamoto, T. Saitô, A. Suzuki, K. Uchizawa et T. Uno ; - la taille de ladite suite est entre 3 et 5 ; - lesdits plus courts chemins sont déterminés en précalculant les plus courts chemins entre l'ensemble des sommets dudit graphe ; 25 - lesdits plus courts chemins sont calculés par l'algorithme de Dij sktra. Un autre objet de l'invention est relatif à un programme d'ordinateur comportant des instructions qui, lorsqu'exécutés par un processeur d'un 7 système informatique, entraînent la mise en oeuvre d'un procédé tel que précédemment défini. Un autre objet de l'invention est une plateforme d'optimisation d'un 5 circuit quantique composée d'une série ordonnée de portes quantiques appliqué à un agencement initial de valeurs de qubits, consistant à insérer un ensemble de portes SWAP locales, afin que l'ensemble des portes dudit circuit quantique soient locales, comprenant des moyens de traitement numérique pour 10 - pour chaque porte quantique de ladite série, si celle-ci n'est pas locale, insérer un ensemble de porte SWAP locales avant ladite porte quantique, - déterminer l'ensemble de permutations, chacune constituée d'une succession d'interversions de valeurs de qubits le long des plus 15 courts chemins entre les positions de qubits associés à ladite porte quantique, et - choisir parmi ledit ensemble de permutations, une permutation minimisant un coût représentant le nombre d'interversions nécessaires pour rendre locales les portes quantiques d'une suite 20 au sein de ladite série, la taille de ladite suite étant substantiellement inférieure à celle de ladite série; - rétablir ledit agencement initial en établissant un arbre couvrant d'un graphe représentatif de l'agencement des qubits dudit circuit, et en effectuant des interversions de valeurs de qubits le long de 25 chemins dudit arbre. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit d'un mode de réalisation préféré de l'invention, donnée à titre d'exemple et en référence aux dessins annexés. 30 8 BREVE DESCRIPTION DES DESSINS La figure 1 représente schématiquement un exemple de circuit quantique simple selon l'état de la technique. 5 Les figure 2a, 2b, 2c représentent schématiquement un exemple de circuit quantique et deux optimisations possibles par insertions de portes Swap, selon l'état de la technique. La figure 3 représente un organigramme schématique selon un mode de réalisation de l'invention. 10 Les figures 4a et 4b représentent schématiquement un exemple de graphe relatif à l'agencement des qubits, selon un mode de réalisation de l'invention. Les figures 5a à 5d représentent également un graphe relatif à l'agencement des qubits, et différentes étapes pour passer d'un agencement à 15 un autre. La figure 6 schématise une arborescence pour le choix d'une permutation, selon un mode de réalisation de l'invention. Les figures 7a à 7d illustrent schématiquement des étapes d'interversions de qubits sur un exemple de graphe, pour passer d'un 20 agencement à un autre. DESCRIPTION DETAILLEE DE L'INVENTION L'invention vise donc à optimiser un circuit quantique constitué d'un 25 ensemble de qubit et d'un ensemble de portes quantiques afin de rendre ces dernières locales. Bien évidemment, le circuit quantique résultant de cette optimisation doit être fonctionnellement équivalent, c'est-à-dire fournir une sortie identique pour une entrée donnée. 9 Un tel circuit quantique optimisé peut ensuite être implémenté sur une infrastructure physique quelconque, notamment une infrastructure imposant une contrainte de localité pour les portes quantiques. 5 On rappelle qu'on appelle « local » le fait de faire interagir des qubits adjacents dans l'arrangement topologique sous-jacent. Une porte est donc dite « locale » si les qubits qu'elle doit faire interagir en entrée sont adjacents. Cette notion est bien connue de l'homme du métier, désormais, et on peut trouver des définitions plus ou moins formalisées dans la littérature. 10 La figure 3 représente un organigramme décrivant une mise en oeuvre possible du procédé selon l'invention. En entrée de cet algorithme, on dispose d'une structure de données décrivant un agencement de qubit. Cette structure peut être un simple vecteur 15 dans le cas d'un agencement sous forme de lignes, mais elle peut être plus complexe pour d'autres types d'agencement. La première étape Si est une étape préliminaire consistant à déterminer un graphe représentatif de l'agencement des qubits du circuit quantique à 20 optimiser. Les sommets du graphe représentent les positions des qubits et les arêtes représentent les relations de voisinage entre les positions des qubits. On fait ici la distinction entre les positions des qubits et les valeurs des qubits. On appelle position l'emplacement au sein de l'infrastructure physique où se situe un qubit. A la suite d'une porte SWAP par exemple, les valeurs de 25 deux positions de qubits sont échangées. Dans les figures 1, 2a, 2b, 2c, les positions de qubits représentent les lignes horizontales, tandis peuvent se déplacer entre les lignes et être modifiés du fait des portes quantiques. Le graphe déterminé à l'issue de cette étape Si est invariant dans sa topologie tout au long du procédé d'optimisation qui va être décrit. Toutefois, 10 ces noeuds, qui correspondent à des positions de qubits, portent la valeur des qubits qui vont évoluer au fil des interversions. La figure 4a représente un graphe d'un agencement simple, en ligne, de 5 4 positions q0, ql, q2, q3 ayant respectivement pour valeurs 0, 1, 2, 3, 4. Après application d'une porte Swap entre les positions q0 et ql, le graphe devient tel que représenté en figure 4b. La figure 5a représente un graphe d'un agencement topologique plus 10 complexe de 10 positions q0, ql, q2, q3, q4, q5, q6, q7, q8, q9. Les arrêtes figurent les adjacences et les étiquettes en regard de chaque position les valeurs respectives. L'étape S2 est une autre étape de prétraitement consistant à calculer le 15 plus court chemin entre chaque paire de sommets du graphe. Ces plus courts chemins seront utilisés dans la suite du procédé. Dans la mesure où ce procédé est itératif, il est intéressant de pré-calculer ces plus courts chemins en amont de la boucle itérative S3-S6, plutôt que de les calculer à plusieurs itérations durant la boucle. 20 Le calcul du plus court chemin entre deux sommets peut être effectué selon plusieurs modes de réalisation. Un algorithme possible est l'algorithme de Dijsktra. Cet algorithme fait partie des algorithmes bien connus de l'homme du métier. Il est par exemple décrit dans tout manuel général 25 d'algorithmique, et sur la page wikipedia : https://fr.wikipedia.org/wiki/Algorithme de Dijkstra Ce pré-calcul a une complexité quadratique en nombre de qubits, mais permet d'ensuite pouvoir estimer en temps constant le coût de « localisation » 30 d'une porte quantique, c'est-à-dire le nombre et le coût des portes Swap à 11 insérer pour rendre locale une porte. On peut ainsi obtenir les portes à insérer en un temps linéaire en nombre de qubits. La boucle S3-S6 est itérée pour chacune des portes quantiques i du 5 circuit, que l'on cherche à rendre locale. Dans une étape S3, on vérifie que la porte quantique i est locale ou non. Si elle est déjà locale (ce qui est bien évidemment le cas pour les portes impliquant un unique qubit), alors on considère la porte quantique suivante, 10 en itérant le compteur i. Dans une étape S4, on créé un ensemble de permutations possibles permettant de rendre locale la porte quantique i. On appelle permutation une bijection de l'ensemble des qubits sur lui-même. 15 Chaque permutation peut être vue comme le résultat obtenu par une série d'interversions des valeurs de qubits de deux positions adjacentes (c'est-à-dire d'une série de portes Swap). Selon une mise en oeuvre de l'invention, on ne recherche pas l'ensemble 20 des permutations possibles, mais on utilise une heuristique permettant d'en déterminer un sous-ensemble. L'heuristique choisie permet de garantir que ce sous-ensemble est pertinent et permet de déterminer, in fine, une solution satisfaisante d'optimisation du circuit dans son entièreté. Pour ce faire, il est proposé de considérer les plus courts chemins entre 25 les qubits concernés par la porte quantique. Si la porte quantique est une porte impliquant deux qubits, on ne considère que ces deux-là ; si la porte quantique est une porte impliquant davantage de qubits, on considère une paire quelconque de qubits parmi ceux-ci. Pour une paire de qubits concernés par ladite porte quantique, on créé 30 l'ensemble des permutations possibles, obtenues, chacune, par une 12 succession d'interversions de qubits situés sur le plus court chemin entre les qubits de la paire, et permettant de rendre locale la porte quantique i. Dans l'exemple de la figure 5a, on considère une porte quantique i 5 impliquant les positions de qubits ql, q6, q8. On choisit, de façon arbitraire (qui peut être aléatoire) la paire ql, q8. On détermine ensuite le plus court chemin sur le graphe entre ces deux positions de qubit. Ce plus court chemin est [ql, q2, q3, q8] On créé ensuite l'ensemble des successions d'interversions entre deux 10 qubits adjacents situés sur ce chemin qui rendent adjacents les qubits situés initialement aux positions ql et q8. On détermine ainsi 3 permutations, pouvant chacune être obtenues par 2 interversions de qubits : - [1, 3, 9, 2, 6, 0, 8, 5, 7, 4], 15 résultant des interversions (q8, q3) et (q3, q2) ; - [1, 2, 3, 9, 6, 0, 8, 5, 7, 4], résultant des interversions (q8, q3) et (ql, q2) ; - [1, 2, 7, 3, 6, 0, 8, 5, 9, 4], résultant des interversions (ql, q2) et (q2, q3). 20 Les figures 5b, 5c et 5d schématisent des graphes d'agencement topologique correspondants, respectivement, à ces 3 permutations appliqués à l'exemple de la figure 5a comme état initial. On voit que les valeurs 3 et 9, situées initialement sur les positions ql 25 et q8, respectivement, sont maintenant adjacentes, dans ces 3 permutations. Si la porte quantique i ne concernait que ces 2 qubits, le travail s'arrêterait là : les valeurs 3 et 9 étant adjacentes, elles peuvent interagir de façon locale. Comme la porte quantique i de l'exemple concerne un troisième qubit 30 de valeur 8 en position q6, le processus se poursuit en considérant le plus 13 court chemin entre cette position q6 et la paire de positions où se situent les valeurs 3 et 9, c'est-à-dire ql, q2 en figure 5b, q2, q3 en figure 5c et q3, q8 en figure 5c. La paire de positions où se situent les valeurs de qubit 3 et 9 sont indiquées avec un trait plus épais. Les plus courts chemins sont 5 représentés par un trait pointillé. De la même façon que précédemment, on créé ensuite, pour chacune des permutations obtenues précédemment, l'ensemble des séries d'interversions entre deux qubits adjacents situés sur ce chemin qui rendent 10 adjacents les qubits situés initialement aux positions q6 et ql ou q2, ou bien q2 ou q3, ou bien q3 ou q8, respectivement. Ces séries supplémentaires peuvent être ajoutées aux successions d'interversions précédemment obtenues. 15 Dans cet exemple, le plus court chemin vers le troisième qubit est de longueur 2 : une seule interversion est donc nécessaire sur ce plus court chemin. Par conséquent, le nombre de successions d'interversions possibles permettant d'obtenir les 3 permutations reste de 3 : - [8, 3, 9, 2, 6, 0, 1, 5, 7, 4], 20 résultant des inversions (q8, q3), (q3, q2) et (q6, q0); - [1, 2, 3, 9, 6, 0, 5, 8, 7, 4], résultant des inversions (q8, q3), (ql, q2) et (q6, q7); - [1, 2, 7, 3, 6, 0, 5, 8, 9, 4], résultant des inversions (ql, q2), (q2, q3) et (q6, q7). 25 Chaque interversion ainsi appliquée correspond à l'insertion d'une porte quantique de type Swap. Il est à noter que cette façon de faire n'est pas optimale, dans le cas 30 d'une porte quantique à 3 qubits. Toutefois, dans la pratique, la plupart des 14 portes sont à un ou deux qubits, de sorte que la contribution non-optimale d'une porte à 3 qubits a peu d'impacts sur l'optimisation d'un circuit quantique dans son entièreté. 5 Pour chaque permutation, on peut calculer un coût associé. Ce coût correspond au nombre de portes quantiques Swap inséré (donc au nombre d'interversions), éventuellement pondéré par un coût par arrête du graphe. Il peut en effet être plus ou moins coûteux d'intervertir des qubits sur certaines positions, du fait de l'infrastructure physique sous-jacente. 10 Ce coût C peut alors s'exprimer : C = p,,, - S,11 , dans lequel Si j représente un nombre de portes Swap insérés entre les positions i et j, et pi,j le poids ed l'arrête i, j du graphe, c'est-à-dire un coûte associé à une interversion des qubits entre les positions i et j. 15 Dans une étape S5, on effectue un choix parmi les permutations déterminées à l'étape précédente S3, en fonction d'un coût. Ce coût peut être le coût C calculé ci-dessus que l'on cherche à minimiser Préférentiellement, toutefois, le coût utilisé pour le choix est un coût 20 cumulant les coûts déterminés pour les différentes interversions pour un certain nombre w de portes quantiques suivantes dans la série ordonnée constituant le circuit quantique. Autrement dit, il s'agit dans cette étape S5 de choisir la permutation minimisant ce coût représentant le nombre d'interversions nécessaires (et 25 éventuellement leurs poids pij) pour rendre locales les portes quantiques d'une suite de longueur w au sein de la série de portes constituant le circuit quantique. La figure 6 schématise une telle arborescence. Le noeud NO représente 30 la position initiale de la figure 5a. Les noeuds Ni, N2, N3 représentent les 3 15 permutations correspondant respectivement aux figures 5b, 5c, 5d et à la porte i. Afin de choisir entre ces 3 permutations, des permutations N4, Ns, N6 NZ sont déterminées pour les portes i+1 et i+2 (avec w=2), et les coûts associés à chaque permutation sont calculés. Ces coûts sont cumulés sur 5 chaque chemin allant d'une des permutations Ni, N2 N3 à une feuille. In fine, le noeud associé au coût le plus faible peut être choisi. Il est à noter qu'il s'agit là d'une heuristique basée sur une optimisation locale, dont le nombre de portes w est un paramètre. 10 Il définit une profondeur d'une arborescence de calcul, et sa détermination fixe un compromis entre une recherche globale sur l'ensemble du circuit qui aboutirait à une explosion combinatoire, et une recherche très locale (w=1, w=2..) qui aboutirait à une solution finale satisfaisante mais non optimale. 15 Le choix du paramètre w peut être fait de façon expérimentale, par exemple en effectuant une batterie de tests sur un ensemble varié de circuits quantiques. Cet ensemble de circuits peut notamment être obtenu de façon aléatoire. Expérimentalement, il est apparu que w=4 formait un bon compromis 20 entre la qualité de l'approximation et le temps de calcul. Dans une étape S6, une fois ce choix effectué, on peut mettre à jour des structures de données, stockées dans une mémoire, qui fourniront la solution en fin de la dernière itération de la boucle S3-S5. 25 Notamment, on peut mettre à jour Le nombre de portes Swap à insérer : nbtotal(i) = nbtotal(i1)+nb(i). Il s'agit simplement d'ajouter le nombre d'interversions de qubits dans la permutation choisie pour la porte quantique i courante, à un nombre total calculé pour la porte précédente i-1. 16 - La liste des portes Swap à insérer : listetotale(i)=listetotale (i- 1)+1iste(i). Il s'agit d'ajouter de façon ensembliste les interversions qui viennent d'être déterminer pour la porte i courante à l'ensemble des interversions déjà déterminées jusqu'à 5 la porte précédente i-1. - Enfin, la structure de données représentant l'agencement du circuit est mise à jour en prenant pour nouvelle valeur celle de la permutation choisie. 10 Si par exemple la permutation correspondant au noeud N2 est choisie, ces structures de valeurs auront pour valeur, à l'issue de l'itération correspondant à cette porte i : Nb(i) = 3 Liste(i) = {(q8, q3), (ql, q2) et (q6, q7)} 15 Et comme agencement : [1, 2, 3, 9, 6, 0, 5, 8, 7, 4], Le procédé peut alors reboucler vers l'étape S3, dans une itération suivante consistant à considérer la porte quantique i+1. Les portes quantiques de la série ordonnée formant le circuit quantique sont ainsi traitées de façon 20 successive. Lorsque la dernière porte quantique du circuit est atteinte (i=N), le procédé se poursuit avec une étape S7 consistant à rétablir l'agencement initial des qubits. Cette étape S7 vise donc à déterminer un ensemble d'insertions de 25 portes quantiques de type Swap permettant à partir de la permutation obtenue pour la dernière itération de la boucle S3-S5 (soit, pour i=N) de revenir à l'agencement initial. On cherche en outre à minimiser le coût associé à ces insertions. Ce problème technique peut être résolu de différentes façons peut être 30 ramené à un problème général d'échange de jetons (« token swapping») qui 17 a été décrit dans l'article « Swapping Labeled Tokens on Graph » de K. Yamanaka, E. D. Demaine, T. Ito, J. Kawahara, M. Kiyomi, Y. Okamoto, T. Saitô, A. Suzuki, K. Uchizawa et T. Uno, in Fun with Algorithms (2014): 364-375. 5 Il a été démontré qu'il existe au moins une solution quelque soient les états d'arrivée et de départ, qui consiste en une succession d'échanges de jetons (ou d'interversions de qubits adjacents, dans le cas équivalent qui nous intéresse). En revanche, ce problème d'échange de jetons étant NP-complet, il 10 s'avère en pratique impossible de déterminer la solution optimale, c'est-à- dire qui donne un nombre minimal d'échanges. L'article mentionné ci-dessus propose donc une méthode heuristique pour déterminer une solution approchée. 15 Les figures 7a, 7b, 7c, 7d sont tirées de l'article mentionné ci-dessus et illustre une série d'échanges de jetons (ou interversions de qubits) entre un agencement fo et un agencement ft au sein d'un graphe comprenant 6 noeuds (ou positions de qubits), vo, vt, v2, v3, v4, vs, v6, sur lesquels sont placés 6 jetons (ou qubits). 20 La figure 7a représente un agencement initial fo pouvant correspondre à la dernière permutation de qubits obtenue en étape S5 pour la dernière porte quantique i=N. La figure 7d représente un agencement final ft qui correspond à l'agencement initial des qubits, tel que par exemple disponible aux étapes 51, S2. 25 Les figures 7b, 7c représentent deux étapes intermédiaires dans le plus court chemin parmi les successions d'échanges de jeton menant de fo à ft : un premier échange entre les jetons 5 et 1 (initialement en positions v4 et vs respectivement) donne l'agencement illustré en figure 7b ; un second échange entre les jetons 4 et 1 (respectivement en positions vi et v4 donne 18 l'agencement illustré en figure 7c ; puis le dernier échange entre les jetons 16 et 3 (respectivement en positions v3 et v6) donne l'agencement final ft. L'article « Swapping Labeled Tokens on Graph » mentionné ci-dessus 5 décrit plusieurs méthodes de résolution d'un tel problème technique dans le cas d'échanges de jetons. Dans le cadre d'interversions de qubits, on créé dans un premier temps un arbre couvrant (ou « recouvrant ») du graphe représentatif de l'agencement 10 des qubits. On appelle arbre couvrant un arbre inclus dans ce graphe et qui connecte tous les sommets du graphe. Parmi les différents arbres couvrant, on choisit celui dont la somme des poids sur les arrêtes est minimale. Si les interversions de qubit, c'est-à-dire 15 les portes quantiques Swap, ont toutes le même poids (ou coût), alors ce problème revient à un simple parcours en profondeur dans le graphe. Dans le cas contraire, on peut utiliser les algorithmes existants de détermination d'arbres couvrants. Parmi ces algorithmes, on peut notamment citer l'algorithme de Kruskal. Une présentation de cet algorithme peut être 20 trouvée sur le site web de l'encyclopédie participative Wikipedia : htt s s://fr.wiki - edia.or - /wiki/Alcorithme de Kruskal et https://en.wikipedia.org/wiki/Kruskal%27salgorithm en langue anglaise. 25 L'algorithme a été présenté pour la première fois dans l'article de Kruskal, J. B., "On the shortest spanning subtree of a graph and the traveling salesman problem" in Proceedings of the American Mathematical Society. 7: 48-50. JSTOR 2033241 (1956). 19 A partir de cet arbre couvrant, on peut utiliser l'algorithme décrit au paragraphe 3.2 de l'article susmentionné. On appelle graphe de conflits D d'un arrangement f, un graphe D=(VD, ED ) répondant aux caractéristiques (1) et (2) : 5 (1) VD= {1)1 E V (G); f (1), ) ft (v1)} (2) Il existe un arc (vt, vt) dans le graphe D, si et seulement si gviKt(vi) avec : - VD l'ensemble des noeuds, c'est-à-dire des positions des 10 qubits, pour le graphe de conflit D, et V celui pour le graphe G représentatif de l'agencement des qubits. - ED et E, les ensembles des arrêtes pour, respectivement le graphe des conflits D et le graphe G. - f(vt) est la valeur du qubit en position vt dans l'arrangement f. 15 Ainsi, sur les figures 7a et 7d, on remarque notamment que f(vt)=4 et ft(v4)=4, que f(v5) = et ft(vt)=1 ; et que f(v4)=5 et ft(v5)=5. On obtient dans un sous-graphe du graphe de conflit D liant vt, v4 et vs. 20 Autrement dit, chaque valeur de qubit f(vt) située à la position vt à l'agencement f, dans le graphe de conflits D doit être déplacé jusqu'à la position vt dans ce même graphe de conflit, de sorte que (vt, vt)GED. On considère ensuite l'arbre couvrant T précédemment déterminé. 25 Pour deux noeuds u et v de cet arbre T, on note P(u, v) un chemin unique de l'arbre entre ces noeuds u et y. On considère le graphe de conflits D pour un agencement initial fo que qubits de l'arbre T. On considère également C=(wi, W2, W3...wq) un cycle orienté arbitraire du graphe de conflit D avec wewi. 30 Si on pose ek=f0(wk) pour tout k tel que 1 k q-1, alors ft(wk+0= ek 20 L'algorithme itératif proposé vise alors à déplacer les qubits e1, e2, e3...eq_i du cycle C vers leurs positions cibles le long des chemins uniques. Plus précisément, on construit une sous-séquence d'interversions Sc comme suit : 5 - on initialise f1,0=f0 - à l'étape k, avec 1 k q - 2 , on considère la valeur de qubit ek=f0(wk) dont la position courante est fj(e k), et on la déplace sur la position sur le chemin P(fk-:°-(? k) k+i) ) qui est adjacent à la position f (k+1,) k,ol.e - 10 On note fk-pt, 0 l'agencement de l'arbre T qui résulte de ce mouvement. - A l'étape (q-1), on déplace la valeur de qubit eq_i=f0(wq_1) depuis sa position courante J1 (e 1) vers la position wq (=w1). 15 A l'issue du déroulement de l'algorithme, cette séquence Sc possède deux caractéristiques : (a) Len(S c) où MIrq-1SPT (Wk, Wk+1) Len(Sc) est la longueur de la séquence Sc et spT(wk, wk-pt) est le nombre d'arrête sur l'arbre T entre les positions wk et Wk+1 20 (b) L'agencement f de l'arbre T obtenu par la séquence Sc satisfait, pour chaque position vt de V(T) : Si y, G C sinon Jf (v,) 1 f(v) Ces caractéristiques nous permettent de pouvoir choisir les cycles 25 orientés du graphe de conflits D dans un ordre arbitraire. Et ainsi, en construisant des séquences d'interversion pour tous les cycles orientés de D, de façon itérative, on aboutit finalement à l'agencement ft de l'arbre T. 21 Autrement dit, on arrive à reconstituer, de façon itérative et en un temps polynomial, l'agencement initial des qubits. Bien entendu, la présente invention n'est pas limitée aux exemples et aux modes de réalisation décrits et représentés, mais elle est susceptible de 5 nombreuses variantes accessibles à l'homme de l'art. REVENDICATIONS1. Procédé d'optimisation d'un circuit quantique composée d'une série ordonnée de portes quantiques appliqué à un agencement initial de valeurs de qubits, consistant à insérer un ensemble de portes SWAP locales, afin que l'ensemble des portes dudit circuit quantique soient locales, comprenant des étapes consistant à: - pour chaque porte quantique de ladite série, si celle-ci n'est pas locale, insérer un ensemble de porte SWAP locales avant ladite porte quantique, - déterminer (S4) l'ensemble de permutations, chacune constituée d'une succession d'interversions de valeurs de qubits le long des plus courts chemins entre les positions de qubits associés à ladite porte quantique, et - choisir (S5) parmi ledit ensemble de permutations, une permutation minimisant un coût représentant le nombre d'interversions nécessaires pour rendre locales les portes quantiques d'une suite au sein de ladite série, la taille de ladite suite étant substantiellement inférieure à celle de ladite série; - rétablir (S7) ledit agencement initial en établissant un arbre couvrant d'un graphe représentatif de l'agencement des qubits dudit circuit, et en effectuant des interversions de valeurs de qubits le long de chemins dudit arbre. 2. Procédé selon la revendication précédente, sans lequel la détermination (S4) de l'ensemble des permutations est effectuée en - déterminant une succession d'interversions de valeurs de qubits le long des premiers plus courts chemins entre une paire quelconque de positions de qubits associés à ladite porte quantique, puis, 23 et si ladite porte quantique possède plus de deux positions de qubits associés, en déterminant une succession d'interversions de valeurs de qubits le long des deuxièmes plus courts chemins entre chacune des positions en dehors de celles déjà considérées, et les 5 premiers plus courts chemins déjà déterminés. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel ledit coût prend en compte un coût attaché à chaque type de portes SWAP à insérer. 10 4. Procédé selon l'une des revendications précédentes, dans lequel le rétablissement de l'agencement initial est effectuée selon un l'algorithme d'échanges de jetons. - " ----------------- -- K. T- ' ' ' ' - - - - --------------------- . , '.. - , Procédé selon l'une des revendications précédentes, dans lequel la taille de ladite suite est entre 3 et 5. -. Procédé selon l'une des revendications précédentes dans lequel lesdits plus courts chemins sont déterminés en précalculant les plus courts 25 chemins entre l'ensemble des sommets dudit graphe. 7. Procédé selon la revendication précédente, dans lequel lesdits plus courts chemins sont calculés par l'algorithme de Dijsktra. 15 20 30 24 8. Programme d'ordinateur comportant des instructions qui, lorsqu'exécutés par un processeur d'un système informatique, entraînent la mise en oeuvre d'un procédé selon l'une des revendications précédentes. 5 . Plateforme d'optimisation d'un circuit quantique composée d'une série ordonnée de portes quantiques appliqué à un agencement initial de valeurs de qubits, consistant à insérer un ensemble de portes SWAP locales, afin que l'ensemble des portes dudit circuit quantique soient locales, 10 comprenant des moyens de traitement numérique pour - pour chaque porte quantique de ladite série, si celle-ci n'est pas locale, insérer un ensemble de porte SWAP locales avant ladite porte quantique, - en déterminant l'ensemble de permutations, chacune constituée 15 d'une succession d'interversions de valeurs de qubits le long des plus courts chemins entre les positions de qubits associés à ladite porte quantique, et - en choisissant parmi ledit ensemble de permutations, une permutation minimisant un coût représentant le nombre 20 d'interversions nécessaires pour rendre locales les portes quantiques d'une suite au sein de ladite série, la taille de ladite suite étant substantiellement inférieure à celle de ladite série; - rétablir ledit agencement initial en établissant un arbre couvrant d'un graphe représentatif de l'agencement des qubits dudit circuit, 25 et en effectuant des interversions de valeurs de qubits le long de chemins dudit arbre. 22 REVENDICATIONS 1. Procédé d'optimisation d'un circuit quantique composée d'une série 5 ordonnée de portes quantiques appliqué à un agencement initial de valeurs de qubits, consistant à insérer un ensemble de portes SWAP locales, afin que l'ensemble des portes dudit circuit quantique soient locales, comprenant des étapes consistant à: - pour chaque porte quantique de ladite série, si celle-ci n'est pas 10 locale, insérer un ensemble de porte SWAP locales avant ladite porte quantique, - déterminer (S4) l'ensemble de permutations, chacune constituée d'une succession d'interversions de valeurs de qubits le long des plus courts chemins entre les positions de qubits associés à ladite 15 porte quantique, et - choisir (S5) parmi ledit ensemble de permutations, une permutation minimisant un coût représentant le nombre d'interversions nécessaires pour rendre locales les portes quantiques d'une suite au sein de ladite série, la taille de ladite 20 suite étant substantiellement inférieure à celle de ladite série; - rétablir (S7) ledit agencement initial en établissant un arbre couvrant d'un graphe représentatif de l'agencement des qubits dudit circuit, et en effectuant des interversions de valeurs de qubits le long de chemins dudit arbre. 25 2. Procédé selon la revendication précédente, sans lequel la détermination (S4) de l'ensemble des permutations est effectuée en - déterminant une succession d'interversions de valeurs de qubits le long des premiers plus courts chemins entre une paire quelconque 30 de positions de qubits associés à ladite porte quantique, puis, 23 - et si ladite porte quantique possède plus de deux positions de qubits associés, en déterminant une succession d'interversions de valeurs de qubits le long des deuxièmes plus courts chemins entre chacune des positions en dehors de celles déjà considérées, et les 5 premiers plus courts chemins déjà déterminés. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel ledit coût prend en compte un coût attaché à chaque type de portes SWAP à insérer. 10 4. Procédé selon l'une des revendications précédentes, dans lequel le rétablissement de l'agencement initial est effectuée selon un l'algorithme d'échanges de jetons. 5. Procédé selon l'une des revendications précédentes, dans lequel la 15 taille de ladite suite est entre 3 et 5. 6. Procédé selon l'une des revendications précédentes dans lequel lesdits plus courts chemins sont déterminés en précalculant les plus courts chemins entre l'ensemble des sommets dudit graphe. 20 7. Procédé selon la revendication précédente, dans lequel lesdits plus courts chemins sont calculés par l'algorithme de Dijsktra. 8. Programme d'ordinateur comportant des instructions qui, 25 lorsqu'exécutés par un processeur d'un système informatique, entraînent la mise en oeuvre d'un procédé selon l'une des revendications précédentes. 9. Plateforme d'optimisation d'un circuit quantique composée d'une série ordonnée de portes quantiques appliqué à un agencement initial de 30 valeurs de qubits, consistant à insérer un ensemble de portes SWAP locales, 24 afin que l'ensemble des portes dudit circuit quantique soient locales, comprenant des moyens de traitement numérique pour - pour chaque porte quantique de ladite série, si celle-ci n'est pas locale, insérer un ensemble de porte SWAP locales avant ladite 5 porte quantique, - en déterminant l'ensemble de permutations, chacune constituée d'une succession d'interversions de valeurs de qubits le long des plus courts chemins entre les positions de qubits associés à ladite porte quantique, et 10 - en choisissant parmi ledit ensemble de permutations, une permutation minimisant un coût représentant le nombre d'interversions nécessaires pour rendre locales les portes quantiques d'une suite au sein de ladite série, la taille de ladite suite étant substantiellement inférieure à celle de ladite série; 15 - rétablir ledit agencement initial en établissant un arbre couvrant d'un graphe représentatif de l'agencement des qubits dudit circuit, et en effectuant des interversions de valeurs de qubits le long de chemins dudit arbre. 20