bloqade-circuit 0.6.0__py3-none-any.whl → 0.6.2__py3-none-any.whl

bloqade/cirq_utils/noise/conflict_graph.py

@@ -0,0 +1,166 @@
+import cirq
+import numpy as np
+
+
+class OneZoneConflictGraph:
+    """
+    Representation of the AOD conflict graph for qubits to more to their entangling partners in a single zone setup.
+
+    Assumes the qubits are specified as cirq.GridQubits with a chosen geometry.
+    """
+
+    def __init__(self, moment: cirq.Moment):
+        """
+        Initializes the conflict graph for a given moment of a cirq circuit.
+
+        :param moment: A cirq.Moment object containing operations (gates) to be analyzed.
+        """
+
+        self.moment = moment
+        self.gates_in_moment = [op for op in moment.operations if len(op.qubits) == 2]
+
+    def _get_nodes(self):
+        """Each qubit participating in a two-qubit gate is a node.
+
+        Sets the self.nodes attribute.
+        """
+        nodes = set()
+        for gate in self.gates_in_moment:
+            nodes.add(gate.qubits[0])
+            nodes.add(gate.qubits[1])
+        self.nodes = nodes
+
+    def _get_edges(self):
+        """
+        Generate the edges of the conflict graph for a given moment.
+
+        Defines self.edges as a set of tuples, where each tuple represents an edge between two qubits.
+        """
+
+        edges = set()
+        for idx1, gate1 in enumerate(self.gates_in_moment):
+            edges.add(gate1.qubits)
+            for gate2 in self.gates_in_moment[idx1 + 1 :]:
+                # X one-to-many, ie. we can't split/merge AOD tones
+                if (gate1.qubits[0].row == gate2.qubits[0].row) ^ (
+                    gate1.qubits[1].row == gate2.qubits[1].row
+                ):
+                    edges.add((gate1.qubits[0], gate2.qubits[0]))
+                    edges.add((gate1.qubits[0], gate2.qubits[1]))
+                    edges.add((gate1.qubits[1], gate2.qubits[0]))
+                    edges.add((gate1.qubits[1], gate2.qubits[1]))
+                # Y one-to-many
+                if (gate1.qubits[0].col == gate2.qubits[0].col) ^ (
+                    gate1.qubits[1].col == gate2.qubits[1].col
+                ):
+                    edges.add((gate1.qubits[0], gate2.qubits[0]))
+                    edges.add((gate1.qubits[0], gate2.qubits[1]))
+                    edges.add((gate1.qubits[1], gate2.qubits[0]))
+                    edges.add((gate1.qubits[1], gate2.qubits[1]))
+                # X ordering, ie. the ordering of AOD tones must be preserved.
+                if (gate1.qubits[0].row < gate2.qubits[0].row) ^ (
+                    gate1.qubits[1].row < gate2.qubits[1].row
+                ):
+                    edges.add((gate1.qubits[0], gate2.qubits[0]))
+                    edges.add((gate1.qubits[1], gate2.qubits[1]))
+                if (gate1.qubits[1].row < gate2.qubits[0].row) ^ (
+                    gate1.qubits[0].row < gate2.qubits[1].row
+                ):
+                    edges.add((gate1.qubits[0], gate2.qubits[1]))
+                    edges.add((gate1.qubits[1], gate2.qubits[0]))
+                # Y ordering
+                if (gate1.qubits[0].col < gate2.qubits[0].col) ^ (
+                    gate1.qubits[1].col < gate2.qubits[1].col
+                ):
+                    edges.add((gate1.qubits[0], gate2.qubits[0]))
+                    edges.add((gate1.qubits[1], gate2.qubits[1]))
+                if (gate1.qubits[1].col < gate2.qubits[0].col) ^ (
+                    gate1.qubits[0].col < gate2.qubits[1].col
+                ):
+                    edges.add((gate1.qubits[0], gate2.qubits[1]))
+                    edges.add((gate1.qubits[1], gate2.qubits[0]))
+
+        self.edges = edges
+
+    def _get_node_degrees(self):
+        """Sets the self.degrees attribute."""
+
+        deg_dict = {}
+        for node in self.nodes:
+            deg_dict[node] = np.sum([node in edge for edge in self.edges])
+        self.degrees = deg_dict
+
+    def get_move_schedule(self, mover_limit: int = 10000):
+        """Generates a move schedule by coloring the conflict graph greedily, first coloring nodes of highest degree.
+
+        Qubits that are the arguments of a single CZ gate are 'partners'. Only one partner need be moved to arrange the
+        atoms for the 2Q gate. Thus, in coloring the conflict graph, as soon as one partner is colored, the other can be
+        disregarded for the purpose of coloring the rest of the graph.
+
+        This sets the self.move_schedule attribute, which is a dictionary where the keys are the indices of the move moments.
+
+        :param mover_limit: The maximum number of qubits that can be moved in a single moment. Added as a constraint
+            when coloring the conflict graph.
+        :returns a dictionary of idx:[cirq.Qid] where idx indexes the move moment where the list of qubits move.
+        """
+
+        self._get_nodes()
+        self._get_edges()
+        self._get_node_degrees()
+
+        self.ordered_nodes = sorted(self.degrees, key=self.degrees.get, reverse=True)
+
+        move_schedule = {}
+        colored_nodes = set()
+        partner_node = None
+        for node in self.ordered_nodes:
+            colored = False
+            for gate in self.gates_in_moment:
+                if node in gate.qubits:
+                    partners = set(gate.qubits)
+                    partners.remove(node)
+                    partner_node = list(partners)[0]
+            if node in colored_nodes:
+                # NOTE: if a node is colored, both it and its partner are added to colored_nodes
+                continue
+            else:
+                connected_nodes = set()
+                for edge in self.edges:
+                    if node in edge:
+                        connected_nodes.add(edge[0])
+                        connected_nodes.add(edge[1])
+                connected_nodes.remove(node)
+                for color in move_schedule.keys():
+                    has_colored_neighbor = False
+                    for connected_node in connected_nodes:
+                        # NOTE: loop through to make sure none of the connected nodes are already assigned to color.
+                        if connected_node in move_schedule[color]:
+                            has_colored_neighbor = True
+                            break
+                        else:
+                            continue
+                    mover_limit_reached = len(move_schedule[color]) >= mover_limit
+                    if not (has_colored_neighbor or mover_limit_reached):
+                        # NOTE: node needs color
+                        move_schedule[color].add(node)
+                        colored = True
+
+                        # NOTE: add this node and it's partner to the solved nodes.
+                        colored_nodes.add(node)
+
+                        if partner_node is not None:
+                            colored_nodes.add(partner_node)
+                        break
+                if not colored:
+                    move_schedule[len(move_schedule)] = {node}
+                    colored = True
+                    colored_nodes.add(node)
+
+                    if partner_node is not None:
+                        colored_nodes.add(partner_node)
+
+        assert set(colored_nodes) == set(self.ordered_nodes)
+
+        self.move_schedule = move_schedule
+
+        return self.move_schedule