examples/2_Expressing_classical_control_flow.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "135014b8",
   "metadata": {},
   "source": [
    "# Expressing classical control flow in quantum programs\n",
    "In the getting-started notebook, you learn how to compose a basic quantum program with only quantum gates. A quantum program should also be able to express classical control flow, such as if-else branches, for-loops and subroutines. With AutoQASM, you write control flow the same way you would do in a Python program. The Pythonic experience in expressing classical control flow makes the quantum programs more readable, and makes it easier to build more complex programs.\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "15cdfaee",
   "metadata": {},
   "outputs": [],
   "source": [
    "# general imports\n",
    "from collections import Counter\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# AutoQASM imports\n",
    "import autoqasm as aq\n",
    "from autoqasm.instructions import cnot, h, measure\n",
    "\n",
    "# AWS imports: Import Braket SDK modules\n",
    "from braket.devices import LocalSimulator"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2cccd8c2",
   "metadata": {},
   "source": [
    "## Subroutine\n",
    "When your quantum program uses the same logic multiple times, subroutines can be a helpful tool. With subroutines, you can generate a more efficient program representation. Your programs become easier to read and reason about, and less work to maintain.\n",
    "\n",
    "As an example, we want to prepare two Bell states, each on different qubit pair. Because the gates to prepare Bell states are the same, we can implement them as a subroutine. Similar to a quantum program, a subroutine is marked by the `@aq.subroutine` decorator. Any number of functions in your program can be marked as a quantum program with this decorator, and they can call one another. Here, we define the Bell state preparation as a subroutine."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6e012372",
   "metadata": {},
   "outputs": [],
   "source": [
    "@aq.subroutine\n",
    "def bell(q0: int, q1: int):\n",
    "    h(q0)\n",
    "    cnot(q0, q1)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e2a884e6",
   "metadata": {},
   "source": [
    "In `two_bell`, the subroutine is called twice to prepare Bell states on different qubit pairs."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "c268be19",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "OPENQASM 3.0;\n",
      "def bell(int[32] q0, int[32] q1) {\n",
      "    h __qubits__[q0];\n",
      "    cnot __qubits__[q0], __qubits__[q1];\n",
      "}\n",
      "qubit[4] __qubits__;\n",
      "bell(0, 1);\n",
      "bell(2, 3);\n"
     ]
    }
   ],
   "source": [
    "@aq.main(num_qubits=4)\n",
    "def two_bell():\n",
    "    bell(0, 1)\n",
    "    bell(2, 3)\n",
    "\n",
    "\n",
    "print(two_bell.build().to_ir())"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c6c59f55",
   "metadata": {},
   "source": [
    "With classical control flow, sometimes the number of qubits used could be undetermined before the execution time (i.e., when program is executed on a QPU or a simulator). Therefore, AutoQASM does not know how many qubits the program may need. In cases like subroutine and using variables as qubit indices, we must define the qubit count using the keyword argument `num_qubits` in `@aq.main`."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1e86b3d5",
   "metadata": {},
   "source": [
    "## If-else logic branch"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bf201c13",
   "metadata": {},
   "source": [
    "In this example, we demonstrate if-else logic. We use a Bell state again, but this time we prepare the state on a dynamically determined qubit pair. We measure qubit `0` and apply the Bell state preparation onto qubits `1` and `2` if the result is `1` or else `3` and `4` if the result is `0`. Finally, we measure all five qubits."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "850bb524",
   "metadata": {},
   "outputs": [],
   "source": [
    "@aq.main(num_qubits=5)\n",
    "def conditioned_bell():\n",
    "    h(0)\n",
    "    if measure(0):\n",
    "        bell(1, 2)\n",
    "    else:\n",
    "        bell(3, 4)\n",
    "\n",
    "    return measure()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "db483f0e",
   "metadata": {},
   "source": [
    "Let's run the quantum program on the Braket simulator."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "c6603968",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "measurement counts:  Counter({'11100': 131, '00000': 131, '00011': 124, '10000': 114})\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGwCAYAAABPSaTdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAozklEQVR4nO3df1TVdZ7H8ddXflyRAMUfXElK2qVswnQW05GmJEvM8cfOeiYr3bLRSse0iBqVsRI7Cepsyg5sOraNsrVme3ayabdZF3QUM3RExErb1JkhpeIOZQioCAjf/aPj93RD0uDCvXx6Ps6553i/38+9vC/fsKff+wUs27ZtAQAAGKqHvwcAAADoTMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIwW7O8BAkFLS4s+/fRTRUREyLIsf48DAAAug23bqqurU2xsrHr0aPv8DbEj6dNPP1VcXJy/xwAAAO1QUVGhQYMGtbmf2JEUEREh6ctPVmRkpJ+nAQAAl6O2tlZxcXHO/8fbQuxIzltXkZGRxA4AAN3MpS5B4QJlAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGC/b3AKYbvPgtf4/wnfXRion+HgHdFF+3/sPXLToDZ3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYLRgfw8AdFeDF7/l7xG+sz5aMdHfIwDoRjizAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACj+TV2du3apcmTJys2NlaWZemNN95w9jU1NWnRokUaOnSowsPDFRsbq/vvv1+ffvqp13M0NDRowYIF6tevn8LDwzVlyhR9/PHHXfxKAABAoPJr7Jw5c0bDhg1TXl5eq31nz57VgQMH9PTTT+vAgQN6/fXXdfToUU2ZMsVrXVpamrZs2aLNmzdr9+7dOn36tCZNmqTm5uauehkAACCA+fWHCk6YMEETJky46L6oqCgVFhZ6bcvNzdXIkSN14sQJXXXVVaqpqdFLL72kl19+WXfccYck6ZVXXlFcXJy2bdum8ePHX/S5Gxoa1NDQ4Nyvra310SsCAACBpltds1NTUyPLstS7d29JUmlpqZqampSamuqsiY2NVWJiooqLi9t8nuzsbEVFRTm3uLi4zh4dAAD4SbeJnXPnzmnx4sWaPn26IiMjJUkej0ehoaHq06eP19qYmBh5PJ42nysjI0M1NTXOraKiolNnBwAA/tMtfjdWU1OT7rnnHrW0tOiFF1645HrbtmVZVpv7XS6XXC6XL0cEAAABKuDP7DQ1NWnatGkqLy9XYWGhc1ZHktxutxobG1VdXe31mKqqKsXExHT1qAAAIAAFdOxcCJ1jx45p27Zt6tu3r9f+pKQkhYSEeF3IXFlZqUOHDik5ObmrxwUAAAHIr29jnT59Wn/605+c++Xl5Tp48KCio6MVGxurn/zkJzpw4ID++7//W83Nzc51ONHR0QoNDVVUVJRmz56tJ554Qn379lV0dLSefPJJDR061PnuLAAA8N3m19jZv3+/brvtNud+enq6JGnmzJnKzMzUm2++KUkaPny41+N27NihlJQUSdKaNWsUHBysadOmqb6+Xrfffrs2btyooKCgLnkNAAAgsPk1dlJSUmTbdpv7v2nfBT179lRubq5yc3N9ORoAADBEQF+zAwAA0FHEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACjdYvfeg4AQEcNXvyWv0f4zvpoxUS/fnzO7AAAAKMROwAAwGjEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIxG7AAAAKMROwAAwGjEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIxG7AAAAKMROwAAwGjEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACj+TV2du3apcmTJys2NlaWZemNN97w2m/btjIzMxUbG6uwsDClpKTo8OHDXmsaGhq0YMEC9evXT+Hh4ZoyZYo+/vjjLnwVAAAgkPk1ds6cOaNhw4YpLy/vovtXrVql1atXKy8vTyUlJXK73Ro3bpzq6uqcNWlpadqyZYs2b96s3bt36/Tp05o0aZKam5u76mUAAIAAFuzPDz5hwgRNmDDhovts21ZOTo6WLFmiqVOnSpLy8/MVExOjTZs2ac6cOaqpqdFLL72kl19+WXfccYck6ZVXXlFcXJy2bdum8ePHd9lrAQAAgSlgr9kpLy+Xx+NRamqqs83lcmnMmDEqLi6WJJWWlqqpqclrTWxsrBITE501F9PQ0KDa2lqvGwAAMFPAxo7H45EkxcTEeG2PiYlx9nk8HoWGhqpPnz5trrmY7OxsRUVFObe4uDgfTw8AAAJFwMbOBZZled23bbvVtq+71JqMjAzV1NQ4t4qKCp/MCgAAAk/Axo7b7ZakVmdoqqqqnLM9brdbjY2Nqq6ubnPNxbhcLkVGRnrdAACAmQI2duLj4+V2u1VYWOhsa2xsVFFRkZKTkyVJSUlJCgkJ8VpTWVmpQ4cOOWsAAMB3m1+/G+v06dP605/+5NwvLy/XwYMHFR0drauuukppaWnKyspSQkKCEhISlJWVpV69emn69OmSpKioKM2ePVtPPPGE+vbtq+joaD355JMaOnSo891ZAADgu82vsbN//37ddtttzv309HRJ0syZM7Vx40YtXLhQ9fX1mjdvnqqrqzVq1CgVFBQoIiLCecyaNWsUHBysadOmqb6+Xrfffrs2btyooKCgLn89AAAg8Pg1dlJSUmTbdpv7LctSZmamMjMz21zTs2dP5ebmKjc3txMmBAAA3V3AXrMDAADgC8QOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIxG7AAAAKMROwAAwGjEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIxG7AAAAKMROwAAwGjEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIxG7AAAAKMROwAAwGjEDgAAMBqxAwAAjEbsAAAAoxE7AADAaMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIwW0LFz/vx5PfXUU4qPj1dYWJiuueYaPfvss2ppaXHW2LatzMxMxcbGKiwsTCkpKTp8+LAfpwYAAIEkoGNn5cqVWrdunfLy8vR///d/WrVqlX75y18qNzfXWbNq1SqtXr1aeXl5Kikpkdvt1rhx41RXV+fHyQEAQKAI6NjZs2eP/v7v/14TJ07U4MGD9ZOf/ESpqanav3+/pC/P6uTk5GjJkiWaOnWqEhMTlZ+fr7Nnz2rTpk1+nh4AAASCgI6dH/7wh9q+fbuOHj0qSXr33Xe1e/du/ehHP5IklZeXy+PxKDU11XmMy+XSmDFjVFxc3ObzNjQ0qLa21usGAADMFOzvAb7JokWLVFNToyFDhigoKEjNzc1avny57r33XkmSx+ORJMXExHg9LiYmRsePH2/zebOzs7Vs2bLOGxwAAASMgD6z89prr+mVV17Rpk2bdODAAeXn5+uf/umflJ+f77XOsiyv+7Ztt9r2VRkZGaqpqXFuFRUVnTI/AADwv4A+s/Pzn/9cixcv1j333CNJGjp0qI4fP67s7GzNnDlTbrdb0pdneAYOHOg8rqqqqtXZnq9yuVxyuVydOzwAAAgIAX1m5+zZs+rRw3vEoKAg51vP4+Pj5Xa7VVhY6OxvbGxUUVGRkpOTu3RWAAAQmAL6zM7kyZO1fPlyXXXVVbrhhhtUVlam1atXa9asWZK+fPsqLS1NWVlZSkhIUEJCgrKystSrVy9Nnz7dz9MDAIBAENCxk5ubq6efflrz5s1TVVWVYmNjNWfOHD3zzDPOmoULF6q+vl7z5s1TdXW1Ro0apYKCAkVERPhxcgAAECgCOnYiIiKUk5OjnJycNtdYlqXMzExlZmZ22VwAAKD7COhrdgAAADqK2AEAAEYjdgAAgNHaFTsVFRX6+OOPnfv79u1TWlqa1q9f77PBAAAAfKFdsTN9+nTt2LFD0pc/0G/cuHHat2+ffvGLX+jZZ5/16YAAAAAd0a7YOXTokEaOHClJ+o//+A8lJiaquLhYmzZt0saNG305HwAAQIe0K3aampqcX7ewbds2TZkyRZI0ZMgQVVZW+m46AACADmpX7Nxwww1at26d3n77bRUWFurOO++UJH366afq27evTwcEAADoiHbFzsqVK/XrX/9aKSkpuvfeezVs2DBJ0ptvvum8vQUAABAI2vUTlFNSUvT555+rtrZWffr0cbY//PDDCg8P99lwAAAAHdWuMztjx45VXV2dV+hIUnR0tO6++26fDAYAAOAL7YqdnTt3qrGxsdX2c+fO6e233+7wUAAAAL7yrd7Geu+995w/f/DBB/J4PM795uZmbd26VVdeeaXvpgMAAOigbxU7w4cPl2VZsixLY8eObbU/LCxMubm5PhsOAACgo75V7JSXl8u2bV1zzTXat2+f+vfv7+wLDQ3VgAEDFBQU5PMhAQAA2utbxc7VV18tSWppaemUYQAAAHytXd96LklHjx7Vzp07VVVV1Sp+nnnmmQ4PBgAA4Avtip0XX3xRP/vZz9SvXz+53W5ZluXssyyL2AEAAAGjXbHz3HPPafny5Vq0aJGv5wEAAPCpdv2cnerqat11112+ngUAAMDn2hU7d911lwoKCnw9CwAAgM+1622sv/3bv9XTTz+tvXv3aujQoQoJCfHa/+ijj/pkOAAAgI5qV+ysX79eV1xxhYqKilRUVOS1z7IsYgcAAASMdsVOeXm5r+cAAADoFO26ZgcAAKC7aNeZnVmzZn3j/t/85jftGgYAAMDX2hU71dXVXvebmpp06NAhnTp16qK/IBQAAMBf2hU7W7ZsabWtpaVF8+bN0zXXXNPhoQAAAHzFZ9fs9OjRQ48//rjWrFnjq6cEAADoMJ9eoPznP/9Z58+f9+VTAgAAdEi73sZKT0/3um/btiorK/XWW29p5syZPhkMAADAF9oVO2VlZV73e/Toof79++v555+/5HdqAQAAdKV2xc6OHTt8PQcAAECnaFfsXPDZZ5/pyJEjsixL1157rfr37++ruQAAAHyiXRconzlzRrNmzdLAgQN166236pZbblFsbKxmz56ts2fP+npGAACAdmtX7KSnp6uoqEj/9V//pVOnTunUqVP63e9+p6KiIj3xxBO+nhEAAKDd2vU21m9/+1v953/+p1JSUpxtP/rRjxQWFqZp06Zp7dq1vpoPAACgQ9p1Zufs2bOKiYlptX3AgAG8jQUAAAJKu2Jn9OjRWrp0qc6dO+dsq6+v17JlyzR69GifDQcAANBR7XobKycnRxMmTNCgQYM0bNgwWZalgwcPyuVyqaCgwNczAgAAtFu7Ymfo0KE6duyYXnnlFX344YeybVv33HOPZsyYobCwMF/PCAAA0G7tip3s7GzFxMTooYce8tr+m9/8Rp999pkWLVrkk+EAAAA6ql3X7Pz617/WkCFDWm2/4YYbtG7dug4PBQAA4Cvtih2Px6OBAwe22t6/f39VVlZ2eKiv+uSTT/SP//iP6tu3r3r16qXhw4ertLTU2W/btjIzMxUbG6uwsDClpKTo8OHDPp0BAAB0X+2Knbi4OL3zzjuttr/zzjuKjY3t8FAXVFdX6+abb1ZISIj+53/+Rx988IGef/559e7d21mzatUqrV69Wnl5eSopKZHb7da4ceNUV1fnszkAAED31a5rdh588EGlpaWpqalJY8eOlSRt375dCxcu9OlPUF65cqXi4uK0YcMGZ9vgwYOdP9u2rZycHC1ZskRTp06VJOXn5ysmJkabNm3SnDlzfDYLAADontoVOwsXLtQXX3yhefPmqbGxUZLUs2dPLVq0SBkZGT4b7s0339T48eN11113qaioSFdeeaXmzZvnXBhdXl4uj8ej1NRU5zEul0tjxoxRcXFxm7HT0NCghoYG535tba3PZgYAAIGlXW9jWZallStX6rPPPtPevXv17rvv6osvvtAzzzzj0+H+8pe/aO3atUpISND//u//au7cuXr00Uf1b//2b5K+vHZIUquf5hwTE+Psu5js7GxFRUU5t7i4OJ/ODQAAAke7zuxccMUVV+imm27y1SyttLS0aMSIEcrKypIkff/739fhw4e1du1a3X///c46y7K8HmfbdqttX5WRkaH09HTnfm1tLcEDAICh2nVmp6sMHDhQ3/ve97y2XX/99Tpx4oQkye12S1KrszhVVVUX/d1dF7hcLkVGRnrdAACAmQI6dm6++WYdOXLEa9vRo0d19dVXS5Li4+PldrtVWFjo7G9sbFRRUZGSk5O7dFYAABCYOvQ2Vmd7/PHHlZycrKysLE2bNk379u3T+vXrtX79eklfvn2VlpamrKwsJSQkKCEhQVlZWerVq5emT5/u5+kBAEAgCOjYuemmm7RlyxZlZGTo2WefVXx8vHJycjRjxgxnzcKFC1VfX6958+apurpao0aNUkFBgSIiIvw4OQAACBQBHTuSNGnSJE2aNKnN/ZZlKTMzU5mZmV03FAAA6DYC+podAACAjiJ2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGC0bhU72dnZsixLaWlpzjbbtpWZmanY2FiFhYUpJSVFhw8f9t+QAAAgoHSb2CkpKdH69et14403em1ftWqVVq9erby8PJWUlMjtdmvcuHGqq6vz06QAACCQdIvYOX36tGbMmKEXX3xRffr0cbbbtq2cnBwtWbJEU6dOVWJiovLz83X27Flt2rSpzedraGhQbW2t1w0AAJipW8TOI488ookTJ+qOO+7w2l5eXi6Px6PU1FRnm8vl0pgxY1RcXNzm82VnZysqKsq5xcXFddrsAADAvwI+djZv3qwDBw4oOzu71T6PxyNJiomJ8doeExPj7LuYjIwM1dTUOLeKigrfDg0AAAJGsL8H+CYVFRV67LHHVFBQoJ49e7a5zrIsr/u2bbfa9lUul0sul8tncwIAgMAV0Gd2SktLVVVVpaSkJAUHBys4OFhFRUX61a9+peDgYOeMztfP4lRVVbU62wMAAL6bAjp2br/9dr3//vs6ePCgcxsxYoRmzJihgwcP6pprrpHb7VZhYaHzmMbGRhUVFSk5OdmPkwMAgEAR0G9jRUREKDEx0WtbeHi4+vbt62xPS0tTVlaWEhISlJCQoKysLPXq1UvTp0/3x8gAACDABHTsXI6FCxeqvr5e8+bNU3V1tUaNGqWCggJFRET4ezQAABAAul3s7Ny50+u+ZVnKzMxUZmamX+YBAACBLaCv2QEAAOgoYgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGC2gYyc7O1s33XSTIiIiNGDAAP34xz/WkSNHvNbYtq3MzEzFxsYqLCxMKSkpOnz4sJ8mBgAAgSagY6eoqEiPPPKI9u7dq8LCQp0/f16pqak6c+aMs2bVqlVavXq18vLyVFJSIrfbrXHjxqmurs6PkwMAgEAR7O8BvsnWrVu97m/YsEEDBgxQaWmpbr31Vtm2rZycHC1ZskRTp06VJOXn5ysmJkabNm3SnDlz/DE2AAAIIAF9ZuframpqJEnR0dGSpPLycnk8HqWmpjprXC6XxowZo+Li4jafp6GhQbW1tV43AABgpm4TO7ZtKz09XT/84Q+VmJgoSfJ4PJKkmJgYr7UxMTHOvovJzs5WVFSUc4uLi+u8wQEAgF91m9iZP3++3nvvPb366qut9lmW5XXftu1W274qIyNDNTU1zq2iosLn8wIAgMAQ0NfsXLBgwQK9+eab2rVrlwYNGuRsd7vdkr48wzNw4EBne1VVVauzPV/lcrnkcrk6b2AAABAwAvrMjm3bmj9/vl5//XX94Q9/UHx8vNf++Ph4ud1uFRYWOtsaGxtVVFSk5OTkrh4XAAAEoIA+s/PII49o06ZN+t3vfqeIiAjnOpyoqCiFhYXJsiylpaUpKytLCQkJSkhIUFZWlnr16qXp06f7eXoAABAIAjp21q5dK0lKSUnx2r5hwwY98MADkqSFCxeqvr5e8+bNU3V1tUaNGqWCggJFRER08bQAACAQBXTs2LZ9yTWWZSkzM1OZmZmdPxAAAOh2AvqaHQAAgI4idgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGIHAAAYjdgBAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGI3YAQAARiN2AACA0YgdAABgNGNi54UXXlB8fLx69uyppKQkvf322/4eCQAABAAjYue1115TWlqalixZorKyMt1yyy2aMGGCTpw44e/RAACAnxkRO6tXr9bs2bP14IMP6vrrr1dOTo7i4uK0du1af48GAAD8LNjfA3RUY2OjSktLtXjxYq/tqampKi4uvuhjGhoa1NDQ4NyvqamRJNXW1vp8vpaGsz5/TlyezjieX8Wx9R+Orbk689hyXP2ns47rhee1bfsb13X72Pn888/V3NysmJgYr+0xMTHyeDwXfUx2draWLVvWantcXFynzAj/iMrx9wToLBxbc3FszdTZx7Wurk5RUVFt7u/2sXOBZVle923bbrXtgoyMDKWnpzv3W1pa9MUXX6hv375tPua7qLa2VnFxcaqoqFBkZKS/x4EPcWzNxHE1F8f24mzbVl1dnWJjY79xXbePnX79+ikoKKjVWZyqqqpWZ3sucLlccrlcXtt69+7dWSN2e5GRkXxxGYpjayaOq7k4tq190xmdC7r9BcqhoaFKSkpSYWGh1/bCwkIlJyf7aSoAABAouv2ZHUlKT0/XfffdpxEjRmj06NFav369Tpw4oblz5/p7NAAA4GdGxM7dd9+tkydP6tlnn1VlZaUSExP1+9//XldffbW/R+vWXC6Xli5d2uotP3R/HFszcVzNxbHtGMu+1PdrAQAAdGPd/podAACAb0LsAAAAoxE7AADAaMQOAAAwGrFjkF27dmny5MmKjY2VZVl64403vPa//vrrGj9+vPr16yfLsnTw4MFWz7F+/XqlpKQoMjJSlmXp1KlTrdZUV1frvvvuU1RUlKKionTfffe1WnfixAlNnjxZ4eHh6tevnx599FE1Njb67sV+B73wwguKj49Xz549lZSUpLffftvZZ9u2MjMzFRsbq7CwMKWkpOjw4cNej29oaNCCBQvUr18/hYeHa8qUKfr444+91ixfvlzJycnq1atXmz9o87HHHlNSUpJcLpeGDx/u65f5ndQVx5av2651qb+POa5di9gxyJkzZzRs2DDl5eW1uf/mm2/WihUr2nyOs2fP6s4779QvfvGLNtdMnz5dBw8e1NatW7V161YdPHhQ9913n7O/ublZEydO1JkzZ7R7925t3rxZv/3tb/XEE0+0/8V9x7322mtKS0vTkiVLVFZWpltuuUUTJkzQiRMnJEmrVq3S6tWrlZeXp5KSErndbo0bN051dXXOc6SlpWnLli3avHmzdu/erdOnT2vSpElqbm521jQ2Nuquu+7Sz372szZnsW1bs2bN0t133915L/g7pKuOLV+3XetSfx9zXLuYDSNJsrds2XLRfeXl5bYku6ysrM3H79ixw5ZkV1dXe23/4IMPbEn23r17nW179uyxJdkffvihbdu2/fvf/97u0aOH/cknnzhrXn31Vdvlctk1NTXtfk3fZSNHjrTnzp3rtW3IkCH24sWL7ZaWFtvtdtsrVqxw9p07d86Oioqy161bZ9u2bZ86dcoOCQmxN2/e7Kz55JNP7B49ethbt25t9fE2bNhgR0VFfeNMS5cutYcNG9b+FwXbtrvm2PJ1619f//uY49r1OLODb2XPnj2KiorSqFGjnG0/+MEPFBUVpeLiYmdNYmKi1y9mGz9+vBoaGlRaWtrlM3d3jY2NKi0tVWpqqtf21NRUFRcXq7y8XB6Px2u/y+XSmDFjnGNSWlqqpqYmrzWxsbFKTEx01qDrddWx5es2sHBcux6xg2/F4/FowIABrbYPGDDA+WWsHo+n1S9h7dOnj0JDQ1v9wlZc2ueff67m5uZWn9OYmBh5PB7nc9rWfunLYxIaGqo+ffq0uQZdr6uOLV+3gYXj2vWIHXxrlmW12mbbttf2y1mDb+frn7tLfc4v5/PNMQkMXXFs+boNPBzXrkPs4Ftxu93661//2mr7Z5995vzrwe12t/oXQ3V1tZqamlr9CwOX1q9fPwUFBbX6nFZVVSkmJkZut1uS2twvfXlMGhsbVV1d3eYadL2uOrZ83QYWjmvXI3bwrYwePVo1NTXat2+fs+2Pf/yjampqlJyc7Kw5dOiQKisrnTUFBQVyuVxKSkrq8pm7u9DQUCUlJamwsNBre2FhoZKTkxUfHy+32+21v7GxUUVFRc4xSUpKUkhIiNeayspKHTp0yFmDrtdVx5av28DCcfUDf1wVjc5RV1dnl5WV2WVlZbYke/Xq1XZZWZl9/Phx27Zt++TJk3ZZWZn91ltv2ZLszZs322VlZXZlZaXzHJWVlXZZWZn94osv2pLsXbt22WVlZfbJkyedNXfeead944032nv27LH37NljDx061J40aZKz//z583ZiYqJ9++232wcOHLC3bdtmDxo0yJ4/f37XfTIMs3nzZjskJMR+6aWX7A8++MBOS0uzw8PD7Y8++si2bdtesWKFHRUVZb/++uv2+++/b9977732wIED7draWuc55s6daw8aNMjetm2bfeDAAXvs2LH2sGHD7PPnzztrjh8/bpeVldnLli2zr7jiCue/p7q6OmfNsWPH7LKyMnvOnDn2tdde66xpaGjouk+IQbrq2PJ127Uu9fcxx7VrETsGufDt4l+/zZw507btL7+d+GL7ly5d6jzH0qVLL7pmw4YNzpqTJ0/aM2bMsCMiIuyIiAh7xowZrb5F/fjx4/bEiRPtsLAwOzo62p4/f7597ty5zv8kGOxf/uVf7KuvvtoODQ21/+7v/s4uKipy9rW0tNhLly613W637XK57FtvvdV+//33vR5fX19vz58/346OjrbDwsLsSZMm2SdOnPBaM3PmzIse/x07djhrxowZc9E15eXlnfnyjdYVx5av2651qb+POa5dy7Jt2+7ss0cAAAD+wjU7AADAaMQOAAAwGrEDAACMRuwAAACjETsAAMBoxA4AADAasQMAAIxG7AAAAKMROwC6REpKitLS0trcP3jwYOXk5HTJLB999JEsy9LBgwe75OMB8C9iB0BAKCkp0cMPP+zctyxLb7zxxmU/fuPGjerdu/dlrY2Li1NlZaUSExO/5ZQAuqNgfw8AAJLUv3//Lvk4jY2NCg0Nldvt7pKPB8D/OLMDoMucP39e8+fPV+/evdW3b1899dRTuvDr+b76NtbgwYMlSf/wD/8gy7Kc+++++65uu+02RUREKDIyUklJSdq/f7927typn/70p6qpqZFlWbIsS5mZmc5zPffcc3rggQcUFRWlhx56qNXbWDt37pRlWdq+fbtGjBihXr16KTk5WUeOHPGa/7nnntOAAQMUERGhBx98UIsXL9bw4cOd/Tt37tTIkSMVHh6u3r176+abb9bx48c769MJ4DIROwC6TH5+voKDg/XHP/5Rv/rVr7RmzRr967/+a6t1JSUlkqQNGzaosrLSuT9jxgwNGjRIJSUlKi0t1eLFixUSEqLk5GTl5OQoMjJSlZWVqqys1JNPPuk83y9/+UslJiaqtLRUTz/9dJvzLVmyRM8//7z279+v4OBgzZo1y9n37//+71q+fLlWrlyp0tJSXXXVVVq7dq2z//z58/rxj3+sMWPG6L333tOePXv08MMPy7KsDn/eAHQMb2MB6DJxcXFas2aNLMvSddddp/fff19r1qzRQw895LXuwltavXv39nq76cSJE/r5z3+uIUOGSJISEhKcfVFRUbIs66JvT40dO9Yrfj766KOLzrd8+XKNGTNGkrR48WJNnDhR586dU8+ePZWbm6vZs2frpz/9qSTpmWeeUUFBgU6fPi1Jqq2tVU1NjSZNmqS/+Zu/kSRdf/313+rzA6BzcGYHQJf5wQ9+4HWmY/To0Tp27Jiam5sv6/Hp6el68MEHdccdd2jFihX685//fFmPGzFixGWtu/HGG50/Dxw4UJJUVVUlSTpy5IhGjhzptf6r96Ojo/XAAw9o/Pjxmjx5sv75n/9ZlZWVl/VxAXQuYgdAt5GZmanDhw9r4sSJ+sMf/qDvfe972rJlyyUfFx4eflnPHxIS4vz5QpS1tLS02nbBheuNLtiwYYP27Nmj5ORkvfbaa7r22mu1d+/ey/rYADoPsQOgy3z9f/x79+5VQkKCgoKCWq0NCQm56Bmfa6+9Vo8//rgKCgo0depUbdiwQZIUGhp62WeI2uO6667Tvn37vLbt37+/1brvf//7ysjIUHFxsRITE7Vp06ZOmwnA5SF2AHSZiooKpaen68iRI3r11VeVm5urxx577KJrBw8erO3bt8vj8ai6ulr19fWaP3++du7cqePHj+udd95RSUmJc13M4MGDdfr0aW3fvl2ff/65zp4969PZFyxYoJdeekn5+fk6duyYnnvuOb333nvO2Z7y8nJlZGRoz549On78uAoKCnT06FGu2wECABcoA+gy999/v+rr6zVy5EgFBQVpwYIFXj9I8Kuef/55paen68UXX9SVV16po0eP6uTJk7r//vv117/+Vf369dPUqVO1bNkySVJycrLmzp2ru+++WydPntTSpUudbz/3hRkzZugvf/mLnnzySZ07d07Tpk3TAw884Jzt6dWrlz788EPl5+fr5MmTGjhwoObPn685c+b4bAYA7WPZX3/TGQBwWcaNGye3262XX37Z36MA+Aac2QGAy3D27FmtW7dO48ePV1BQkF599VVt27ZNhYWF/h4NwCVwZgcALkN9fb0mT56sAwcOqKGhQdddd52eeuopTZ061d+jAbgEYgcAABiN78YCAABGI3YAAIDRiB0AAGA0YgcAABiN2AEAAEYjdgAAgNGIHQAAYDRiBwAAGO3/ASAmhWlZTD/rAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "device = LocalSimulator(\"autoqasm\")\n",
    "result = device.run(conditioned_bell, shots=500).result()\n",
    "counts = Counter(result.measurements[\"return_value\"])\n",
    "print(\"measurement counts: \", counts)\n",
    "\n",
    "plt.bar(counts.keys(), counts.values())\n",
    "plt.xlabel(\"bitstrings\")\n",
    "plt.ylabel(\"counts\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c1ded7b8",
   "metadata": {},
   "source": [
    "The simulator result shows that the Bell state is prepared on the qubit pair (1, 2) when the 0th (leftmost) qubit is measured as 1, resulting in measured bitstrings '10000' and '11100'. Otherwise, the Bell state is prepared on the qubit pair (3, 4), resulting in measured bitstrings '00000' and '00011'."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "edee6e9c",
   "metadata": {},
   "source": [
    "## For-loop"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ce1da10e",
   "metadata": {},
   "source": [
    "In this example, we demonstrate for-loop by creating many Bell states in a quantum program. Instead of writing the preparation steps explicitly, we can use a for-loop to make the program simpler and more readable. In the generated QASM script, you can see the code structure of for-loop is preserved."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "88403f29",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "OPENQASM 3.0;\n",
      "def bell(int[32] q0, int[32] q1) {\n",
      "    h __qubits__[q0];\n",
      "    cnot __qubits__[q0], __qubits__[q1];\n",
      "}\n",
      "qubit[6] __qubits__;\n",
      "for int i in [0:2:6 - 1] {\n",
      "    bell(i, i + 1);\n",
      "}\n"
     ]
    }
   ],
   "source": [
    "n_bell = 3\n",
    "\n",
    "\n",
    "@aq.main(num_qubits=n_bell * 2)\n",
    "def multiple_bell():\n",
    "    for i in aq.range(0, len(aq.qubits), 2):\n",
    "        bell(i, i + 1)\n",
    "\n",
    "\n",
    "print(multiple_bell.build().to_ir())"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9303a270",
   "metadata": {},
   "source": [
    "## Summary"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2619f38d",
   "metadata": {},
   "source": [
    "In this notebook, you have learned how to use AutoQASM to express classical control flow in your quantum program. AutoQASM not only enables writing program with real-time control flow based on measurement results, but also increases the readability by taking advantage of high level programming features such as for-loops and subroutines."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.18"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}