{ "cells": [ { "cell_type": "markdown", "id": "ca6d6102", "metadata": {}, "source": [ "(sec-tutorial-ops-qubits)=\n", "\n", "# Tutorial: Operations and Qubits\n", "\n", "```{seealso}\n", "The complete source code of this tutorial can be found in\n", "\n", "{nb-download}`Operations and Qubits.ipynb`\n", "```\n", "\n", "## Gates, measurements and qubits\n", "\n", "In the previous tutorials, experiments were created on the {ref}`quantum-device level`. On this level,\n", "operations are defined in terms of explicit signals and locations on the chip, rather than the qubit and the intended operation.\n", "To work at a greater level of abstraction, `quantify-scheduler` allows creating operations on the\n", "{ref}`quantum-circuit level`.\n", "Instead of signals, clocks, and ports, operations are defined by the effect they have on specific qubits. This representation of the schedules can be compiled to the quantum-device level to create the pulse schemes.\n", "\n", "In this tutorial we show how to define operations on the {ref}`quantum-circuit level`, combine them into schedules, and show their circuit-level visualization.\n", "We go through the configuration file needed to compile the schedule to the quantum-device level and show how these configuration files can be created automatically and dynamically.\n", "Finally, we showcase the hybrid nature of `quantify-scheduler`, allowing the scheduling of circuit-level and device-level operations side by side in the same schedule.\n", "\n", "Many of the gates used in the circuit layer description are defined in\n", "{class}`~quantify_scheduler.operations.gate_library` such as {class}`~quantify_scheduler.operations.gate_library.Reset`, {class}`~quantify_scheduler.operations.gate_library.X90` and\n", "{class}`~quantify_scheduler.operations.gate_library.Measure`.\n", "Operations are instantiated by providing them with the name of the qubit(s) on which\n", "they operate:" ] }, { "cell_type": "code", "execution_count": 1, "id": "699ef160", "metadata": {}, "outputs": [], "source": [ "from quantify_scheduler.operations.gate_library import CZ, Measure, Reset, X90\n", "\n", "q0, q1 = (\"q0\", \"q1\")\n", "X90(q0)\n", "Measure(q1)\n", "CZ(q0, q1)\n", "Reset(q0);\n" ] }, { "cell_type": "markdown", "id": "181aa232", "metadata": {}, "source": [ "Let's investigate the different components present in the circuit-level description of\n", "the operation. As an example, we create a 45 degree rotation operation over the\n", "x-axis." ] }, { "cell_type": "code", "execution_count": 2, "id": "e71605e4", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'acquisition_info': [],\n", " 'gate_info': {'operation_type': 'Rxy',\n", " 'phi': 0.0,\n", " 'plot_func': 'quantify_scheduler.schedules._visualization.circuit_diagram.gate_box',\n", " 'qubits': ['q0'],\n", " 'tex': '$R_{xy}^{45, 0}$',\n", " 'theta': 45.0,\n", " 'unitary': array([[0.92387953+0.j , 0. -0.38268343j],\n", " [0. -0.38268343j, 0.92387953+0.j ]])},\n", " 'logic_info': {},\n", " 'name': \"Rxy(45, 0, 'q0')\",\n", " 'pulse_info': []}\n" ] } ], "source": [ "from pprint import pprint\n", "from quantify_scheduler.operations.gate_library import Rxy\n", "\n", "rxy45 = Rxy(theta=45.0, phi=0.0, qubit=q0)\n", "pprint(rxy45.data)\n" ] }, { "cell_type": "markdown", "id": "9d07d091", "metadata": {}, "source": [ "As we can see, the structure of a circuit-level operation is similar to a pulse-level\n", "operation. However, the information is contained inside the {code}`gate_info` entry rather\n", "than the {code}`pulse_info` entry of the data dictionary.\n", "Importantly, there is no device-specific information coupled to the operation such that\n", "it represents the abstract notion of this qubit rotation, rather than how to perform it\n", "on any physical qubit implementation.\n", "\n", "The entries present above are documented in the `operation` schema.\n", "Generally, these schemas are only important when defining custom operations, which is\n", "not part of this tutorial. This schema can be inspected via:" ] }, { "cell_type": "code", "execution_count": 3, "id": "19c5580b", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'operation_type': {'description': 'Defines what class of operations this gate '\n", " 'refers to (e.g. Rxy, CZ etc.).',\n", " 'type': 'string'},\n", " 'plot_func': {'description': 'reference to a function for plotting this '\n", " 'operation. If not specified, defaults to using '\n", " ':func:`~quantify_scheduler.schedules._visualization.circuit_diagram.gate_box`.',\n", " 'type': ['string', 'null']},\n", " 'qubits': {'description': 'A list of strings indicating the qubits the gate '\n", " 'acts on. Valid qubits are strings that appear in '\n", " 'the device_config.json file.',\n", " 'type': 'array'},\n", " 'symmetric': {'description': 'A boolean to indicate whether a two qubit gate '\n", " 'is symmetric. This is used in the device config '\n", " 'compilation stage. By default, it is set as '\n", " 'False',\n", " 'type': 'boolean'},\n", " 'tex': {'description': 'latex snippet for plotting', 'type': 'string'},\n", " 'unitary': {'description': 'A unitary matrix describing the operation.'}}\n" ] } ], "source": [ "import importlib.resources\n", "import json\n", "from quantify_scheduler import schemas\n", "\n", "operation_schema = json.loads(importlib.resources.read_text(schemas, \"operation.json\"))\n", "pprint(operation_schema[\"properties\"][\"gate_info\"][\"properties\"])\n" ] }, { "cell_type": "markdown", "id": "70b91da6", "metadata": {}, "source": [ "## Schedule creation from the circuit layer (Bell)\n", "\n", "The circuit-level operations can be used to create a {class}`~quantify_scheduler.schedules.schedule.Schedule` within\n", "`quantify-scheduler` using the same method as for the pulse-level operations.\n", "This enables creating schedules on a more abstract level.\n", "Here, we demonstrate this extra layer of abstraction by creating a {class}`~quantify_scheduler.schedules.schedule.Schedule` for measuring\n", "Bell violations.\n", "\n", "```{note}\n", "Within a single {class}`~quantify_scheduler.schedules.schedule.Schedule`, high-level circuit layer operations can be mixed with quantum-device level operations. This mixed representation is useful for experiments where some pulses cannot easily be represented as qubit gates. An example of this is given by the `Chevron` experiment given in {ref}`Mixing pulse and circuit layer operations (Chevron)`.\n", "```\n", "\n", "As the first example, we want to create a schedule for performing the\n", "[Bell experiment](https://en.wikipedia.org/wiki/Bell%27s_theorem).\n", "The goal of the Bell experiment is to create a Bell state\n", "{math}`|\\Phi ^+\\rangle=\\frac{1}{2}(|00\\rangle+|11\\rangle)` which is a perfectly entangled state, followed by a measurement.\n", "By rotating the measurement basis, or equivalently one of the qubits, it is possible\n", "to observe violations of the CSHS inequality.\n", "\n", "We create this experiment using the\n", "{ref}`quantum-circuit level` description.\n", "This allows defining the Bell schedule as:" ] }, { "cell_type": "code", "execution_count": 4, "id": "d286adc5", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Schedule \"Bell experiment\" containing (66) 147 (unique) operations." ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import numpy as np\n", "from quantify_scheduler import Schedule\n", "from quantify_scheduler.operations.gate_library import CZ, Measure, Reset, Rxy, X90\n", "\n", "sched = Schedule(\"Bell experiment\")\n", "\n", "for acq_idx, theta in enumerate(np.linspace(0, 360, 21)):\n", " sched.add(Reset(q0, q1))\n", " sched.add(X90(q0))\n", " sched.add(X90(q1), ref_pt=\"start\") # Start at the same time as the other X90\n", " sched.add(CZ(q0, q1))\n", " sched.add(Rxy(theta=theta, phi=0, qubit=q0))\n", "\n", " sched.add(Measure(q0, acq_index=acq_idx), label=\"M q0 {:.2f} deg\".format(theta))\n", " sched.add( \n", " Measure(q1, acq_index=acq_idx),\n", " label=\"M q1 {:.2f} deg\".format(theta),\n", " ref_pt=\"start\", # Start at the same time as the other measure\n", " )\n", "\n", "sched\n" ] }, { "cell_type": "markdown", "id": "52c601fc", "metadata": {}, "source": [ "By scheduling 7 operations for 21 different values for {code}`theta` we indeed get a schedule containing 7\\*21=147 operations. To minimize the size of the schedule, identical operations are stored only once. For example, the {class}`~quantify_scheduler.operations.gate_library.CZ` operation is stored only once but used 21 times, which leaves only 66 unique operations in the schedule.\n", "\n", "```{note}\n", "The acquisitions are different for every iteration due to their different {code}`acq_index`. The {class}`~quantify_scheduler.operations.gate_library.Rxy`-gate rotates over a different angle every iteration and must therefore also be different for every iteration (except for the last since {math}`R^{360}=R^0`). Hence the number of unique operations is 3\\*21-1+4=66.\n", "```\n", "\n", "## Visualizing the quantum circuit\n", "\n", "We can directly visualize the created schedule on the\n", "{ref}`quantum-circuit level`.\n", "This visualization shows every operation on a line representing the different qubits." ] }, { "cell_type": "code", "execution_count": 5, "id": "0d51c398", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzEAAADECAYAAACvIa+3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAtKElEQVR4nO3deXxU1f3/8deZyQIIhCWETTAssgURK4gKAgqaQCi7pVZcIl8VLVRLfyoColVoRamCyxcoilasllJQR7CoLAFExW/RGsEgiGHVBAISNiFk5vz+mBCTkECWSWZu8n4+HvMIc++Zcz8ZZu4nn3vPPddYaxEREREREXEKV7ADEBERERERKQ0VMSIiIiIi4igqYkRERERExFFUxIiIiIiIiKOoiBEREREREUdRESMiIiIiIo6iIkZERERERBxFRYyIiIiIiDiKihgREREREXEUFTEiItWEMSbZGJOc73msMcYaY24PXlQVr/DvLSIizqciRkQkxBhjbs8tLvI/9htj1hhjBgQ7Pgk8Y0wzY8xjxpiuwY5FRMQJwoIdgIiIFGsqkAYYoDFwO/CeMeaX1tplwQzMYW4IdgAl0Ax4FNgJ/DeokYiIOICKGBGR0PVva+1/zjwxxrwMZAA3ASpizsMYU8tae8Jamx3sWEREJLA0nExExDkOAz8BOfkXGmNcxpj7jTFbjDEnjTEZxph5xpj6gdqwMaaeMWaWMWaPMeaUMeZbY8xDxhhX7nqTO9ztgDEmJt/rIowxXxljdhhjLshd9ljuELkOxph/GmOOGGMOGmNmG2NqFLHt0caYTcaYn4wxh4wx/zDGtCjUJtkYs9kYc7kxZp0x5gTwp3zrkvO17Zu7/V8ZYx41xuwzxhw1xvzLGBNljInM/V33G2OOGWNeMcZEljOuTrnvz4nc7T2YPx7g/3KfvpJvCOHtpfxvEhGpNnQmRkQkdEUZY6LxDyeLAcYDtYHXC7Wbh3+o2SvAc0ArYBxwmTGmp7X2dHmCMMbUAtYCzXO3tRu4Gvgz0BS431prjTF3ACnAXGB47sv/CMQBfa21xwt1/U/8w6ceBq4EfgfUB27Nt+3JwBO5bV8CGuF/H9YZYy6z1h7O119D4N/AP/C/Rxnn+dUexl8UPgm0ze33NODLjeOx3Lhuxz+s7/EyxlUfWAEszW0/EphhjPnKWvtvIBX/0MHHgb8C63Nf9/F54hcRqb6stXrooYceeoTQA/8fzbaIx0ngtkJte+Wu+02h5fGFlwPJQHK+57G5bW4/TzxTgGPAxYWW/xn/WaEW+ZbdldvnzUCP3PXPFnrdY7lt3im0/MXc5V1yn1+U+/pJhdp1xl9sTCr0u1ng7iLiL/x7981t+xUQnm/5G/gLmPcKvf5jYGe+52WJ65Z8yyKAH4B/5VvWrST/F3rooYceevgfGk4mIhK6fgtcn/sYDawBXjLGDM/X5kYgC/jQGBN95gFswl94XBuAOG7Ef3bgx0LbWAm4gd5nGlpr/wq8DzwPLAR2AJOK6ffFQs+fz/05MPfncPzDnv9ZaLvpwHbO/t1O4T8bVVKv2YJnqTbiP+u1oFC7jUALY8yZ0QuljesY+c6eWf81Op8BrUsRq4iI5KPhZCIioeszW/DC/jeBL4AXjDHLcv8YvhiIAvYX00dMMctL42KgC3CghNsYg794uRi42lr7UzGv217o+Q78Z0Ji823XFNHujMLD5PbZ0l3Ev7vQ86zcn3uKWO7C/z4fLENce621ttCyH/G/pyIiUgYqYkREHMJa6zPGrAHuw/+H9Bb8f1zvxz98qyjFFR6l4QI+BJ4qZv22Qs/7AmcuhL8E+KSE2yn8h74rd9kAwFtE+2OFnhdXLBWnqD7PtdyUMa7z9SciIqWkIkZExFnO7Ldr5/7cAfQHNpzjjEd57QBqW2tXnq+hMaYp/mFhHwDZwExjzPvW2l1FNL8Y/wXzZ7TFXyDszLddA6RZawsXSsFUEXEVLuBEROQcdE2MiIhDGGPC8d+4MRv/jFbgn+3KDTxSRPswY0y9AGz6n8BVxpj4IrZRL9+1IgDz8eeWMfgv8s8BXjbGFHXW4beFno/P/fnv3J9L8Z/FeLTw63OndG5Y6t8kMCoirjMzt9UrZ2wiItWCzsSIiISuAcaYDrn/jgF+g//sxZPW2iMA1tq1xph5wMPGmK74z4Cczm13I/6hZ/8qZxxPA4OBZcaYV/FPGnAB/qFiI/Ffw5JpjEkCEvHPsLUXwBgzHv9F7fcA/1uo31bGGA/+6Yevwj95wRvW2i9zf7cdxpgp+GdBizXGvA0cxT+F9DD80xHPLOfvVmoVFNcO/PcBGmuMOYq/qNlorU0756tERKopFTEiIqHr8Xz/PglsxV8MzMvfyFo71hizCbgb/w0ec/APyXod2FDeIKy1J4wxffDPMnYj/vu4HMF/LcyjQJYx5kLgWeBda+3f8r3278aYEcBTxph/F/qjfFTu7/hkbswvAA8U2vaTxphtwO9ztwX+C+8/ADzl/d3KKtBxWWtPG2Nuw18YzcWfn5MoONxORERymbMnTBEREak4xpjH8P/h38hamxnkcERExIF0TYyIiIiIiDiKihgREREREXEUFTEiIiIiIuIouiZGREREREQcRWdiRERERETEUVTEiIiIiIiIo6iIERERERERRwnVm13qQh0RERERkerHlKSRzsSIiIiIiIijqIgRERERERFHUREjIiIiIiKOoiJGREREREQcRUWMiIiIiIg4iooYERERERFxFBUxIiIiIiLiKCpiRERERETEUVTEiIiIiIiIo6iIERERERERR1ERIyIiIiIijqIiRkREREREHEVFjIiIiIiIOIqKGBERERERcRQVMSIiIiIi4igqYkRERERExFFUxIiIiIiIiKOoiBEREREREUcJC3YAUjmeX70z2CEExfjrYoMdgsg5bd26NdghBEWHDh2CHUK1ohwgEpqUA8pOZ2JERERERMRRVMSIiIiIiIijqIgRERERERFHUREjIiIiIiKOoiJGRESKdPz4cXbt2hXsMEREJAhCPQeoiBERkbP4fD4efPBBEhISSE1NDXY4IiJSiZyQA1TEiIjIWebPn8/q1asBcLvdQY5GREQqkxNygIoYEREpYN26dcyePTvveXh4eBCjERGRyuSUHKCbXYqEAGstaZnZfLv/FHt+zCY7xxIRZmhRP4K2MZG0io7AGBPsMKUa2L17NxMmTCiwLCIiIkjRiFQPygESKpyUA1TEiARRjteyMvUIy1Ky+CErBwC3C7CAAa/P365pVBiDukTRv2NdwtxKZFIxTpw4wb333svJkyex1uYtD9UEJuJ0ygESSpyWA1TEiATJ7oPZzFq1n10HswssP5O0+Hn/wQ9ZOcxff5APvj7K/f1iaNkwNHco4lzWWiZPnkxaWho+n6/AulBNYCJOphwgocSJOUDXxIgEwZbvf+KhJfvYcyj7/I3z2XMom4eW7GPL9z9VUGRSXS1cuJAVK1aclbwgdBOYiFMpB0iocWIOUBFTSHJyMrGxsaV6zYIFC9i6dWvFBBSifD4f027vx7KXny6wfMunq/l9fDu+XL8iSJGFvt0Hs5m2LJ3TPovPnr99fj4Lp32WacvS2X2wdMlP5FwOHToEFD0LTagmsIqgHFAyygFlpxwgociJOUBFTAm8+OKLxMbGUqNGDXr06MFnn31WYP3GjRuZN29ekKILDpfLxfU33cv6dxby07EjAOzZtplXnhjP4Dsf4tJrEoIcYWjK8VpmrdpPjs9iS5m8zrAWcnyW2av2k+MtYycihdx3331MmDABr9cL+L/jxhhcLlfITq9ZWZQDzqYcUDbKARKqnJgDAn5NjDFmODAWuBxoAFxmrf1voLdTWRYtWsSECROYO3cuPXr0YNasWcTHx/PNN98QExMDwJAhQxg/fjzPPvtskKOtXN36D2HFwtmse+c1uvcfxrwpY7gyYSTXjhwT7NBC1srUI2eNfy4Ln4WdB7NZmXqEhM5RAYhMqiJrLenp6Wzbti3vsXfvXo4ePcqxY8cwxlC7dm3q1KlDixYtSE1NxeVy4fP5aNCgAZ07d6Zp06bB/jWCSjmgeMoBpaccIJWpqueAiriw/wLgI+CfwPwK6L9SPfPMM9x5550kJSUBMHfuXJYvX86CBQuYOHEiAP369SMjI4PNmzfTuXPnYIZbqdzuMPqPGsuyV/7C56vfpWX7Lgz/7aPBDitkWWtZlpIVsP4MsDzlCPFxdTX1puQ5ceIEK1euxOPxsGHDhhK/7osvvijwPDMzE4Du3btz8uRJatSoEdA4nUI5oHjKAaWjHCCVoTrlgFIXMcaYC4A5wHDgKDAT+CXwX2vt/dbahbntYgMYZ1BkZ2ezadMmHn744bxlLpeL/v3788knn+Qti4yM5IYbbsDj8VSrBAb+I3FL//dxaNCI2yc/h8tVshGKr077HdcMHs3yV5/h6I8HcbvdxN8ynsv6JFZwxMGTlpmdN4VmIFjg+6zTpGVm07pRZMD6Fefxer18+umneDwePB5PwPpNTk4mOTkZgBEjRjB48GC6detW4u+50ykHnJ9yQMkpB0hFqa45oCxnYp4G+gBDgP3An4BfAP8NXFihITMzE6/XS+PGjQssb9y48VkXcQ4ZMoQ5c+YwadKkygwx6BY/5z/qdjzrEKaEH2qfz0fm97uIbh7L8HuncmHbThw5dICnx/6STldcS2TNWhUZctB8u/9UxfR74JQSWDVlreXDDz9k8uTJHDt2rMg2HTt2JC4ujnbt2tG+fXtiY2OpVatW3oWa2dnZHD9+nDFjxvDdd99hrSUyMpJTpwp+XpcsWcKSJUuIjo5m2rRp9OnTp8J/v2BTDjg/5YCSUw6QQKvuOaBURYwxpjYwBhhtrV2Vu+w2YG9ZAzDGRAIFvn3ffvtthc6E0KJFi4D3OXDgQJKSksjMzCQ6OrrINnv27An4doNp+St/YcvGNUx4YSkvPHALn7y3iN5Db81bn7F7B8/+biQ1a9elRq0LOLBvFxe27cSwe6bQ4uLORDWMIaqhf0x53QaNuCCqPieOHg5oAgul9/zr3V7cBgJ5HabbQOquQ3SsE7ghCuIMW7du5YknnuDzzz8/a12/fv0YOnQovXv3Pu++NCIigsOHD7Njxw4AjDFMnDiRoUOHsmbNGt555x3Wrl2b1z4zM5OxY8fSs2dPJk2aROvWrcsUfzC/m8oBgaEcUDrKARJIVTkHlHQfXdozMW2ACGDjmQXW2kPGmG9K2U9+DwMFBtH27t27HN2d3759+0rULjo6GrfbTUZGRoHlGRkZNGnSpMCytLQ06tWrR7169Yrt78orryx1rIEycWHJx0WWxMfL32T14pcYN/MNmrfpRN8RSaxaNI+eg27CHRYOQOOWbYjt2JXBdz5Es9YdmHbbddzz5N9IXrqADt2uKdDf7m1fYX0+6sc0C2icwXzPC2uV+CCNugzA5Q7cpWinc3JY/Na/eep/ngpYn1K5Vq1aVar2P/74I7Nnz2bRokUFlrds2ZKkpCQSEhLOuR8qisfjybuY0+12k5CQQI0aNRgwYAADBgzg0KFDvPfee8yfP5/9+/cDsGHDBhITE7n11lsZN24cderUKdU2g/ndVA4oP+WA0lMOkKIoB5ytpPvoiriwv7T+DDyTf8G6deuyQmFO6oiICC6//HJWrVrF0KFDAf9p8FWrVjFu3LgCbT0eDwMHDiQsrPi39NNPP63IcM/p7e3egPW1ZeMaFj//GLdPeY5WnS4DoPeQW1m96K989uFbXDXgV3ltM/Z8R0yL1mSf/Amfz0dkzVps+/xjeg/5+Wjd8SOHef3JP/DrP/w5YDGeEcz3vLAlm718vMsX0KNw4WFh3DhsCMMfGR64TqVSHT9+vMRtU1JSuPnmm8nJKTiu/vHHH2f48OFlmgbTWsvSpUvzklefPn3OSoANGjRg9OjRjBo1ikWLFjF9+vS8da+99hpvvfUWCxcupH379iXebih9N4ujHFA05YCyUQ6QoigHlF1pi5gdwGmgB7AbwBhTH2gHrD3H64plrT0FVMxA0QCYMGECt912G926deOKK65g1qxZHD9+PG+mmjM8Hg+PPPLIOfuqiCEMJbZ9Z0C62b3tK155fBxD7pzIpb3i85bXrF2X3sNuZ+Wbc+hxwwhcbjcnTxwjLDyCsPAIdu34ksYt23DyxDGs9VGzdl0ATmef4qWpd9H/prG0jrs8IDHmF9T3vJBOR46wfmdmQPv0Wuh4UQNatKgb0H6l8pT0JomLFy9m6tSpBZbdcccdjB07ttRHwPJLSUnJO+rl9Xrz/lgvSnh4OKNHjyYxMZEXXniBN954A4CjR48ydOhQZs6cSWJiyS7MDqXv5rkoBxSkHFB2ygFSFOWAsitVEWOtPWaMeRl42hhzEP+F/dMB35k2xpgGQEvgzDnh9rlT/6Vba9PLHXElGzVqFAcOHGDq1Kmkp6fTtWtXVqxYUeBCz7S0NL755hsSEqr+zb1atruEmcu3FLkuMWkCiUkT8p6n79pOk4vaApCx+1tiLmzFti8+5uKuVwH+6v/vTz3AxZddzRXXV/2jSG1jKubCy7a6oLPK279/PzNmzMh73rRpU+bPn0+bNm3K3bfH48HtduP1eqldu3aJhvPWr1+fRx55hJEjR5KUlERWln88/vTp0+nVqxdRUVXnvhXKAQUpB5SdcoCUlXJA0coyR9oDwHrgXWAl/nvCbMq3fjDwBbA89/k/cp+PLXuYwTVu3Dh27drFqVOn2LhxIz169Ciw3uPx0Ldv33JVwlXRDzu30zS2HQARNWqxddNHpP7fOjpc7h8L/d3m//BF8jK+2vABM+4ayIy7BvL9dyU7IuFEraIjaBoVuBGcBmgWFU6r6OAPvZSKFRMTwzPP+EfdxsfH8+677wYkeWVnZ+PxePB6vbjdbgYNGlSqSVU6duzI+++/T69evQB47rnnqlQBc4ZyQNkoBxSkHCBlpRxQtFJ/m6y1x4Bbch8AGGMS861/FXg1ALE5hsfjYfjwqn8UqbTyj4vu2nsAXXsP4Ol7BnPj+D8C0OaS7sxe+V2wwqt0xhgGdYli/vqDAenPAolddJOz6qJ3794sXryYuLi4gP2fr1+/Pm9aTq/Xy+DBg0vdR1RUFPPmzSM1NZW4uLiAxOU0ygFFUw4oSDlAykM54GyhcGG/45V2Zonq7IE5gbsJkxP171iXD74+yp5D2fjKcXGny0DLBhFc30njoKuTQN9I8e23384bRtC8eXO6du1apn5cLle1LWBAOaA0lAOUA6TslAMKbbfStxjiYmNjuf/++4MdhlRRYW7D/f1iCHMZynogxRgIcxnu6xeD26UjcFI2WVlZJCcn4/V6cblcDBs2TEd0UQ6QiqUcIKGiKuSAgBQx1tq+1tr7A9FXsCmBSUVr2TCCKYOaEO4ylDb/uAyEuwxTBjWhZUONg5ayW7FiRd40nT6fr0zDCKoi5QCpaMoBEgqqQg7QmRiRIIhrVpMZI5rTskHJktCZPNeyQQQzRjQnrlnNigtOqoWlS5dijMEYQ5cuXUJqKlqRqk45QIKtKuQAXRMjEiQtG/qT0crUIyxLyeKHrII3r3K7wJs7eXnTqHASu9Slf8e6hLmddbpXQs+ePXtISUnJez5ixIggRiNSPSkHSLBUlRygIkYkiMLchoTOUcTH1SUtM5ut6Sd5+SP/zDX9O9YhNjqSto0iaRUd4bixqhK6VqxYkffvsLAw4uPjz9FaRCqKcoAEQ1XJASpiREKAMYbWjSJpVi88L4HdelVDaoRrxKcEXvfu3WnYsCEHDx7kuuuuq5L3dhFxEuUAqUxVJQeoiBERqWa6du3K2rVr2bJlCx06dAh2OCIiUomqSg5QESMiUg253W66dOkS7DBERCQIqkIO0HlKERERERFxFBUxIiIiIiLiKCpiRERERETEUVTEiIiIiIiIo+jC/mpi/HWxwQ5BRIrg5JlhxDmUA0RCk3JA2elMjIiIiIiIOIqKGBERERERcRQVMSIiIiIi4igqYkRERERExFFUxIiIiIiIiKOoiBEREREREUdRESMiIiIiIo6iIkZERERERBxFRYyIiIiIiDiKihgREREREXEUFTEiIiIiIuIoKmJERERERMRRVMSIiIiIiIijqIgRERERERFHUREjIiIiIiKOoiJGREREREQcRUWMiIiIiIg4iooYERERERFxFBUxIiIiIiLiKCpiRERERETEUVTEiIiIiIiIo6iIERERERERR1ERIyIiIiIijqIiRkREREREHEVFjIiIiIiIOIqKGBERERERcRQVMSIiIiIi4igqYkRERERExFFUxIiIiIiIiKOoiBEREREREUdRESMiIiIiIo6iIkZERERERBxFRYyIiIiIiDiKihgREREREXEUFTEiIiIiIuIoKmJERERERMRRVMSIiIiIiIijhAU7AKkcz6/eGewQgmL8dbHBDqFa0ees9LZu3Rq4QBykQ4cOwQ6hWtF3UyqDPmelpxxQdjoTIyIiIiIijqIiRkREREREHEVFjIiIiIiIOIqKGBERERERcRQVMSIiUqTjx4+za9euYIchIiJBEOo5QEWMiIicxefz8eCDD5KQkEBqamqwwxERkUrkhBygIkZERM4yf/58Vq9eDYDb7Q5yNCIiUpmckANUxIiISAHr1q1j9uzZec/Dw8ODGI2IiFQmp+QA3exSJARYa0nLzCb1h5N5y/728UFaRUfSNiaSVtERGGOCGGFo2vzJKt6eOx1rLf1G3c3Vib8OSNvqbPfu3UyYMKHAsoiIiCBFI1I9KAeU3EtT72b7l5/S7rKrGfPYnLzla/71Mp+8twhrLe0v78mI3z5a7HumfFA8J+UAFTEiQZTjtaxMPcKylCx+yMopsG7V1qN4fUcBaBoVxqAuUfTvWJcwtxIZgNebw1tzpjH+L29Ss3Ydnh77Sy7tFc8FUfXL1bY6O3HiBPfeey8nT57EWpu3PFQTmIjTKQeUXp/hSfRIuJHPPliSt+zo4YOse/s1Ji34ALc7jNm/H8XOr7+gVdwvznq98kHxnJYDNJxMJEh2H8zmwSX7mL/+4FnJC8Dr+/nfP2TlMH/9QR5cso/dB7MrMcrQtWvrlzSNbUe9Rk2IrHkBHa/oy9ZN68vdtrqy1jJ58mTS0tLwer0F1oVqAhNxMuWAsrm465XUqFX7rOU+bw6ns0/h9ebgzTlN7foNi3y98kHRnJgDdCZGJAi2fP8T05alk+Oz52+cz55D2Ty0ZB9TBjUhrlnNCoouNMz+/Sh2pHwGgDssnOhmLYkfPZ5u/YYAkJWZQVR047z29aIbczgzvci+StO2ulq4cCErVqwocl2oJjARp1IOKNr59vvFqVOvIdf96k4evaknbrebnoNuplGzi4psq3xQNCfmAJ2JKSQ5OZnY2NhSvWbBggVs3bq1YgIKUT6fj2m392PZy08XWL7l09X8Pr4dX64v+osg/qNv05alc9pnKWX+wmfhtM8ybVl6lT4aZ61l77dfM/TuSUxb/BlTXl1Fx+59eP3JP3Dwhz3BDq9KOnToEFD0LDShmsAqgnJAySgHlJ1yQNHKs98/cTSLLZ+u5rG/r+fxRZ+S9vUmvk3ZWEmRVw1OzAEqYkrgxRdfJDY2lho1atCjRw8+++yzAus3btzIvHnzghRdcLhcLq6/6V7Wv7OQn44dAWDPts288sR4Bt/5EJdekxDkCENTjtcya9V+cnwWW8rkdYa1kOOzzF61nxxvGTsJcQf2pnHqxDE6XtGHug0a0bBpC64c8Ct8Pi8Ze74DICq6MVmZGXmvOZyZQVTDxkX2V5q21dV9993HhAkT8oYRuFwujDG4XK6QnV6zsigHnE05oGyUA4pXkv1+cb75/COim13EBXXrERFZg7ge17Lz6y+KbKt8UDQn5oCADycz/qkg/gjcCdQDNgD3WGu3B3pblWHRokVMmDCBuXPn0qNHD2bNmkV8fDzffPMNMTExAAwZMoTx48fz7LPPBjnaytWt/xBWLJzNundeo3v/YcybMoYrE0Zy7cgxwQ4tZK1MPcKuABw981nYeTCblalHSOgcFYDIQsue7ZupVSeKJhddDMCPB35g2YKZhIVH0Lx1BwAu6nApP+zcxuED6dSsXYfUz5JJGD0egBf+382Mfugv1GvU5LxtqzJrLenp6Wzbti3vsXfvXo4ePcqxY8cwxlC7dm3q1KlDixYtSE1NxeVy4fP5aNCgAZ07d6Zp06bB/jWCSjmgeMoBpaccULyS7PeLU69RM9K+/pzT2adwu8PY/uXGvBnHqnM+qOo5oCKuiXkQ+B1wG5AGPAG8b4zpZK09ec5XhqBnnnmGO++8k6SkJADmzp3L8uXLWbBgARMnTgSgX79+ZGRksHnzZjp37hzMcCuV2x1G/1FjWfbKX/h89bu0bN+F4b99NNhhhSxrLctSsgLWnwGWpxwhPq5ulZt6c8/2zfx0/CgPDOqM9Xk5nX2K8MgajLp/et5YZrc7jKFjJ/H8H27KnSbzLi6Iqo+1lgP7dlKrbr28/oprWxWdOHGClStX4vF42LBhQ4lf98UXBY9aZmZmAtC9e3dOnjxJjRo1AhqnUygHFE85oHSUA86tJPt9gBceGM2+HalknzzBI6Ou4o6pL9Iq7hfEXdGXp+5OxBgX7X5xNZdcfX21zAfVKQeUuogxxlwAzAGGA0eBmcAvgf8CvwfuB6ZZa9/JbX8rkAEMBf4RgJgrTXZ2Nps2beLhhx/OW+Zyuejfvz+ffPJJ3rLIyEhuuOEGPB5PtUpg4D8St/R/H4cGjbh98nO4XCUbofjqtN9xzeDRLH/1GY7+eBC32038LeO5rE9iBUccPGmZ2UXOQFNWFvg+6zRpmdm0bhQZsH5Dwd7tW7hmyC30GXY7Px07wtvz/kTruMvpkTCyQLtLrr6eS66+vsCyjN07uPSaAURE1jhv26rC6/Xy6aef4vF48Hg8Aes3OTmZ5ORkAEaMGMHgwYPp1q1bib/nTqcccH7KASWnHHBuJd3vj3v69SJfP2jMAwwa80CBZem7vq0W+aC65oCynIl5GugDDAH2A38CfoG/iGkFNAFWnmlsrc0yxmwErqKERUzhqd0CraRj+zIzM/F6vTRuXHCsZOPGjc+6iHPIkCHMmTOHSZMmFdtfRf9ewbD4Of9Rt+NZhzAl/FD7fD4yv99FdPNYht87lQvbduLIoQM8PfaXdLriWiJr1gpYfKH0nm9L/6lC+t2e8RMXNahaEw3u2b6Zqwb+mkbNYwG48XdPMOPOAVyVeBMNm1x4ztc2uagtw++dUglR/ixYnzNrLR9++CGTJ0/m2LFjRbbp2LEjcXFxtGvXjvbt2xMbG0utWrXyLtTMzs7m+PHjjBkzhu+++w5rLZGRkZw6dapAP0uWLGHJkiVER0czbdo0+vTpU67Yg/ndVA4IHOWAklMOOLfy7PeLU1n5QDmg9M71npV0H12qT70xpjYwBhhtrV2Vu+w2YG9ukya5PzMKvTQj37rCfUYCBQ4hzJkzp6imATN27NiA9zlw4ECSkpLIzMwkOjq6yDbz588P+HZLrF3gL7Jc/spf2LJxDRNeWMoLD9zCJ+8tovfQW/PWZ+zewbO/G0nN2nWpUesCDuzbxYVtOzHsnim0uLgzUQ1jiGroH1Net0EjLoiqz4mjhwOawIL6nhey1bbH0AIbwPk0DD7+vf4L0j76JmB9lksAPmeZ3+/mp2NHaNqqXd6yprEXE92sJZtWvcMNN/+23NsItPJ8zvr27Vum123dupUnnniCzz///Kx1/fr1Y+jQofTu3fu8s8pERERw+PBhduzYAYAxhokTJzJ06FDWrFnDO++8w9q1a/PaZ2ZmMnbsWHr27MmkSZNo3bp1meIP5ndTOSAwlANKRzmgeE7c7+enHFB653rPSrqPLm3p3gaIAPLmrbPWHjLGlOfb8zBQYBDtE088UY7uzq+kb050dDRut5uMjII1WUZGBk2aFKzJ0tLSqFevHvXq1Su2v4r+vc5l4sLAJrCPl7/J6sUvMW7mGzRv04m+I5JYtWgePQfdhDssHIDGLdsQ27Erg+98iGatOzDttuu458m/kbx0AR26XVOgv93bvsL6fNSPaRbQOIP5nhfWKvFBGnVpjssduATm9fr4T0oKi5c/FbA+yyMQn7M927/CHRZOzIWtCixv94uepGz4IC+ZPfu7kQy7ZzKxHS/jjacfommrdrS99ErefWkG9854DYCUDR/w9cY1/HrCn8sd17mU53NW2gT2448/Mnv2bBYtWlRgecuWLUlKSiIhIeGc+6GieDyevIs53W43CQkJ1KhRgwEDBjBgwAAOHTrEe++9x/z589m/fz8AGzZsIDExkVtvvZVx48ZRp06dUm0zmN9N5YDyUw4oPeWA4jlxv5+fckBgc0BFFTHnc+ZuQY2BH/Itb4x/uFlR/gw8k3/Bm2++mRUeHh7g0EovIiKCyy+/nFWrVjF06FDAfxp81apVjBs3rkBbj8fDwIEDCQsr/i3917/+VZHhntN/AngWe8vGNSx+/jFun/IcrTpdBkDvIbeyetFf+ezDt7hqwK/y2mbs+Y6YFq3JPvkTPp+PyJq12Pb5x/Qe8vPRuuNHDvP6k3/g138I/A4nmO95Yav21uHLTDe+8zctsTC3mwHX9+O6pCsC2GvZBeJztmf7Fho1jyUsvODRo/a/6MmGd//Ojwd+oH6jpsSPHs/KN+fS+pLuGJfh2pFj8HpzSN/1LQBebw4rFj7P3dNfLn9Q51FZn7OUlBRuvvlmcnIKjqt//PHHGT58eJmmwbTWsnTp0rzk1adPn7MSYIMGDRg9ejSjRo1i0aJFTJ8+PW/da6+9xltvvcXChQtp3759ibcbSt/N4igHFE05oGyUA4rnxP1+fsoBwckBpS1idgCngR7AbgBjTH2gHbAW/2xk6UA/cosWY0zd3PZFjhGz1p4CThW1LhRMmDCB2267jW7dunHFFVcwa9Ysjh8/njdTzRkej4dHHnnknH1dddVVFRnqOf1n9c6A9LN721e88vg4htw5kUt7xectr1m7Lr2H3c7KN+fQ44YRuNxuTp44Rlh4BGHhEeza8SWNW7bh5IljWOujZu26AJzOPsVLU++i/01jaR13eUBizC+Y73lhR7cc4Yt1mQHt04ehe1xLrupUN6D9llUgPmeD/+dBBv/Pg2ctv/SaBGav/PleAZ2u6MOyBTM5uXEN9zz5KuCfdaZeo6YcTN/L1xvXENfj2rzhKhWpPJ+zkt4kcfHixUydOrXAsjvuuIOxY8eW+ghYfikpKezbtw/wj1E+88d6UcLDwxk9ejSJiYm88MILvPHGGwAcPXqUoUOHMnPmTBITS3Zhdih9N89FOaAg5YCyUw4onhP3+/kpBwQnB5SqiLHWHjPGvAw8bYw5iP/C/ungP7BgrbXGmFnAFGPMdn6eYvl74O1yRxsEo0aN4sCBA0ydOpX09HS6du3KihUrClzomZaWxjfffENCQtW/uVfLdpcwc/mWItclJk0gMWlC3vP0XdtpclFbADJ2f0vMha3Y9sXHXNzV/8G11vL3px7g4suu5orrh1d88EHWNqZiZo9pWwVmpSmLXVu/5MTRwzRo3Dxv+Ar47wGw/b+f8NG7f+f3z4XOUdjy2L9/PzNmzMh73rRpU+bPn0+bNm3K3bfH48HtduP1eqlduza9e/c+72vq16/PI488wsiRI0lKSiIryz9t7PTp0+nVqxdRUVXjvhWgHFCYckDZKQeUX3Xa7+enHFC0sgzMfABYD7yLfxayj4BN+dY/BTwP/BX4P6A2kODEe8ScMW7cOHbt2sWpU6fYuHEjPXr0KLDe4/HQt2/fclXCVdEPO7fTNNZ/kV5EjVps3fQRqf+3jg6X+8dCf7f5P3yRvIyvNnzAjLsGMuOugXz/XcmOSDhRq+gImkYFbgSnAZpFhdMq+twX7VVFhw+k8+ZfJjJ+5hscSt/L92k/X5Z3UYdLeXvun+g56DfUqFU7iFEGTkxMDM884x91Gx8fz7vvvhuQ5JWdnY3H48Hr9eJ2uxk0aNB5LwLNr2PHjrz//vv06tULgOeee65KFTBnKAeUjXJAQcoB5VPd9vv5KQcUrdTfJmvtMeCW3AcAxpjEfOstMDX3US14PB6GD6/6R5FKK/+46K69B9C19wCevmcwN47/IwBtLule4DRxVWeMYVCXKOavPxiQ/iyQ2KVq3OSsNLJPnWTB479l5PjHaNi0Bdf/5h7ef/15kh55AYCYC1txQd169Pzlb4IcaWD17t2bxYsXExcXF7D/8/Xr1+dNy+n1ehk8eHCp+4iKimLevHmkpqYSFxcXkLicRjmgaMoBBSkHlF113e/npxxwNudPLB4CVq1aFewQHOOBOYG7CZMT9e9Ylw++PsqeQ9n4bNn7cRlo2SCC60NkHHRlioiswYTnl+Q9v6xPYoEb5K33vM6Qux/G7a56u7dA30jx7bffzhtG0Lx5c7p27VqmflwuV7UtYEA5oDSUA5QDyqI67/fzUw4otN1K32KIi42N5f777w92GFJFhbkN9/eLIcxlKOuBFGMgzGW4r18MblfVPwJXUge+38W0264jskZNuvS8IdjhhLysrCySk5Pxer24XC6GDRtWLY7ono9ygFQk5YDA0n6/7KpCDghIyWqt7RuIfkKBEphUtJYNI5gyqAnTlqWT47OlOhrnyk1eUwY1oWXD6jEOuqQaNbuIKX9bHewwHGPFihV503T6fL4yDSOoipQDpKIpBwSO9vtlVxVygM7EiARBXLOazBjRnJYNSpaEzhwbadkgghkjmhPXrGbFBSfVwtKlSzHGYIyhS5cutGjRItghiVQbygESbFUhB1TtwYMiIaxlQ38yWpl6hGUpWfyQ5T8i4nbhv2LTgDf3rmhNo8JJ7FKX/h3rEuZ21uleCT179uwhJSUl7/mIESOCGI1I9aQcIMFSVXKAihiRIApzGxI6RxEfV5e0zGy+PXCKPYeyyc6xRIQZWjSIoG2jSFpFRzhurKqErhUrVuT9OywsjPj4+HO0FpGKohwgwVBVcoCKGJEQYIyhdaNIWlejm5ZJ8HTv3p2GDRty8OBBrrvuuip5bxcRJ1EOkMpUVXKAihgRkWqma9eurF27li1bttChQ4dghyMiIpWoquQAFTEiItWQ2+2mS5cuwQ5DRESCoCrkAM1OJiIiIiIijqIiRkREREREHEVFjIiIiIiIOIqKGBERERERcRRjrQ12DEUJyaBERERERKRCleimSDoTIyIiIiIijqIiRkREREREHEVFjIiIiIiIOIqKGBERERERcRQVMSIiIiIi4igqYkRERERExFFUxIiIiIiIiKOoiBEREREREUdRESMiIiIiIo6iIkZERERERBxFRYyIiIiIiDiKihgREREREXEUFTEiIiIiIuIoKmJERERERMRRVMSIiIiIiIijqIgRERERERFHUREjIiIiIiKOoiJGREREREQcJSzYARTDBDsAEREREREJTToTIyIiIiIijqIiRkREREREHEVFjIiIiIiIOIqKGBERERERcRQVMSIiIiIi4igqYkRERERExFFUxIiIiIiIiKOoiBEREREREUdRESMiIiIiIo7y/wGU2lZ37LD6FgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "\n", "_, ax = sched.plot_circuit_diagram()\n", "# all gates are plotted, but it doesn't all fit in a matplotlib figure.\n", "# Therefore we use :code:`set_xlim` to limit the number of gates shown.\n", "ax.set_xlim(-0.5, 9.5)\n", "plt.show()\n" ] }, { "cell_type": "markdown", "id": "981e15a7", "metadata": {}, "source": [ "In previous tutorials, we visualized the `schedules` on the pulse level using {meth}`~quantify_scheduler.schedules.schedule.ScheduleBase.plot_pulse_diagram` .\n", "Up until now, however, all gates have been defined on the\n", "{ref}`quantum-circuit level` without defining the\n", "corresponding pulse shapes.\n", "Therefore, trying to run {meth}`~quantify_scheduler.schedules.schedule.ScheduleBase.plot_pulse_diagram` will raise an error which\n", "signifies no {code}`pulse_info` is present in the schedule:" ] }, { "cell_type": "code", "execution_count": 6, "id": "b3c23bea", "metadata": { "tags": [ "raises-exception" ] }, "outputs": [ { "ename": "RuntimeError", "evalue": "Attempting to sample schedule Bell experiment, but the schedule does not contain any `pulse_info`. Please verify that the schedule has been populated and device compilation has been performed.", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mRuntimeError\u001b[0m Traceback (most recent call last)", "Cell \u001b[0;32mIn[6], line 1\u001b[0m\n\u001b[0;32m----> 1\u001b[0m \u001b[43msched\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mplot_pulse_diagram\u001b[49m\u001b[43m(\u001b[49m\u001b[43m)\u001b[49m\n", "File \u001b[0;32m~/.anaconda3/envs/tmp/lib/python3.9/site-packages/quantify_scheduler/schedules/schedule.py:357\u001b[0m, in \u001b[0;36mScheduleBase.plot_pulse_diagram\u001b[0;34m(self, port_list, sampling_rate, modulation, modulation_if, plot_backend, plot_kwargs)\u001b[0m\n\u001b[1;32m 350\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m plot_backend \u001b[38;5;241m==\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mmpl\u001b[39m\u001b[38;5;124m\"\u001b[39m:\n\u001b[1;32m 351\u001b[0m \u001b[38;5;66;03m# NB imported here to avoid circular import\u001b[39;00m\n\u001b[1;32m 352\u001b[0m \u001b[38;5;66;03m# pylint: disable=import-outside-toplevel\u001b[39;00m\n\u001b[1;32m 353\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01mquantify_scheduler\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mschedules\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01m_visualization\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mpulse_diagram\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m (\n\u001b[1;32m 354\u001b[0m pulse_diagram_matplotlib,\n\u001b[1;32m 355\u001b[0m )\n\u001b[0;32m--> 357\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43mpulse_diagram_matplotlib\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m 358\u001b[0m \u001b[43m \u001b[49m\u001b[43mschedule\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m,\u001b[49m\n\u001b[1;32m 359\u001b[0m \u001b[43m \u001b[49m\u001b[43msampling_rate\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msampling_rate\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 360\u001b[0m \u001b[43m \u001b[49m\u001b[43mport_list\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mport_list\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 361\u001b[0m \u001b[43m \u001b[49m\u001b[43mmodulation\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmodulation\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 362\u001b[0m \u001b[43m \u001b[49m\u001b[43mmodulation_if\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmodulation_if\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 363\u001b[0m \u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mplot_kwargs\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 364\u001b[0m \u001b[43m \u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 365\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m plot_backend \u001b[38;5;241m==\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mplotly\u001b[39m\u001b[38;5;124m\"\u001b[39m:\n\u001b[1;32m 366\u001b[0m \u001b[38;5;66;03m# NB imported here to avoid circular import\u001b[39;00m\n\u001b[1;32m 367\u001b[0m \u001b[38;5;66;03m# pylint: disable=import-outside-toplevel\u001b[39;00m\n\u001b[1;32m 368\u001b[0m \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01mquantify_scheduler\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mschedules\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01m_visualization\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mpulse_diagram\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m (\n\u001b[1;32m 369\u001b[0m pulse_diagram_plotly,\n\u001b[1;32m 370\u001b[0m )\n", "File \u001b[0;32m~/.anaconda3/envs/tmp/lib/python3.9/site-packages/quantify_scheduler/schedules/_visualization/pulse_diagram.py:517\u001b[0m, in \u001b[0;36mpulse_diagram_matplotlib\u001b[0;34m(schedule, port_list, sampling_rate, modulation, modulation_if, x_range, ax)\u001b[0m\n\u001b[1;32m 478\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mpulse_diagram_matplotlib\u001b[39m(\n\u001b[1;32m 479\u001b[0m schedule: Union[Schedule, CompiledSchedule],\n\u001b[1;32m 480\u001b[0m port_list: Optional[List[\u001b[38;5;28mstr\u001b[39m]] \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;01mNone\u001b[39;00m,\n\u001b[0;32m (...)\u001b[0m\n\u001b[1;32m 485\u001b[0m ax: Optional[matplotlib\u001b[38;5;241m.\u001b[39maxes\u001b[38;5;241m.\u001b[39mAxes] \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;01mNone\u001b[39;00m,\n\u001b[1;32m 486\u001b[0m ) \u001b[38;5;241m-\u001b[39m\u001b[38;5;241m>\u001b[39m Tuple[matplotlib\u001b[38;5;241m.\u001b[39mfigure\u001b[38;5;241m.\u001b[39mFigure, matplotlib\u001b[38;5;241m.\u001b[39maxes\u001b[38;5;241m.\u001b[39mAxes]:\n\u001b[1;32m 487\u001b[0m \u001b[38;5;250m \u001b[39m\u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m 488\u001b[0m \u001b[38;5;124;03m Plots a schedule using matplotlib.\u001b[39;00m\n\u001b[1;32m 489\u001b[0m \n\u001b[0;32m (...)\u001b[0m\n\u001b[1;32m 515\u001b[0m \u001b[38;5;124;03m The matplotlib ax.\u001b[39;00m\n\u001b[1;32m 516\u001b[0m \u001b[38;5;124;03m \"\"\"\u001b[39;00m\n\u001b[0;32m--> 517\u001b[0m times, pulses \u001b[38;5;241m=\u001b[39m \u001b[43msample_schedule\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m 518\u001b[0m \u001b[43m \u001b[49m\u001b[43mschedule\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 519\u001b[0m \u001b[43m \u001b[49m\u001b[43msampling_rate\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43msampling_rate\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 520\u001b[0m \u001b[43m \u001b[49m\u001b[43mport_list\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mport_list\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 521\u001b[0m \u001b[43m \u001b[49m\u001b[43mmodulation\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmodulation\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 522\u001b[0m \u001b[43m \u001b[49m\u001b[43mmodulation_if\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmodulation_if\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 523\u001b[0m \u001b[43m \u001b[49m\u001b[43mx_range\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mx_range\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m 524\u001b[0m \u001b[43m \u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 525\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m ax \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n\u001b[1;32m 526\u001b[0m _, ax \u001b[38;5;241m=\u001b[39m plt\u001b[38;5;241m.\u001b[39msubplots()\n", "File \u001b[0;32m~/.anaconda3/envs/tmp/lib/python3.9/site-packages/quantify_scheduler/schedules/_visualization/pulse_diagram.py:402\u001b[0m, in \u001b[0;36msample_schedule\u001b[0;34m(schedule, port_list, modulation, modulation_if, sampling_rate, x_range)\u001b[0m\n\u001b[1;32m 399\u001b[0m max_x \u001b[38;5;241m=\u001b[39m \u001b[38;5;28mmin\u001b[39m(\u001b[38;5;28mmax\u001b[39m(t0 \u001b[38;5;241m+\u001b[39m pulse_info[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mduration\u001b[39m\u001b[38;5;124m\"\u001b[39m], max_x), x_range[\u001b[38;5;241m1\u001b[39m])\n\u001b[1;32m 401\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m np\u001b[38;5;241m.\u001b[39misinf(min_x) \u001b[38;5;129;01mor\u001b[39;00m np\u001b[38;5;241m.\u001b[39misinf(max_x):\n\u001b[0;32m--> 402\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mRuntimeError\u001b[39;00m(\n\u001b[1;32m 403\u001b[0m \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mAttempting to sample schedule \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mschedule\u001b[38;5;241m.\u001b[39mname\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m, \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 404\u001b[0m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mbut the schedule does not contain any `pulse_info`. \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 405\u001b[0m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mPlease verify that the schedule has been populated and \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 406\u001b[0m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mdevice compilation has been performed.\u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 407\u001b[0m )\n\u001b[1;32m 409\u001b[0m time_window \u001b[38;5;241m=\u001b[39m [min_x, max_x]\n\u001b[1;32m 411\u001b[0m logger\u001b[38;5;241m.\u001b[39mdebug(\u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mtime_window \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mtime_window\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m, port_map \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mport_map\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m\"\u001b[39m)\n", "\u001b[0;31mRuntimeError\u001b[0m: Attempting to sample schedule Bell experiment, but the schedule does not contain any `pulse_info`. Please verify that the schedule has been populated and device compilation has been performed." ] } ], "source": [ "sched.plot_pulse_diagram()" ] }, { "cell_type": "markdown", "id": "3ee6df2c", "metadata": {}, "source": [ "And similarly for the {code}`timing_table`:" ] }, { "cell_type": "code", "execution_count": 7, "id": "cded5484", "metadata": { "tags": [ "raises-exception" ] }, "outputs": [ { "ename": "ValueError", "evalue": "Absolute time has not been determined yet.", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)", "Cell \u001b[0;32mIn[7], line 1\u001b[0m\n\u001b[0;32m----> 1\u001b[0m \u001b[43msched\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mtiming_table\u001b[49m\n", "File \u001b[0;32m~/.anaconda3/envs/tmp/lib/python3.9/site-packages/quantify_scheduler/schedules/schedule.py:432\u001b[0m, in \u001b[0;36mScheduleBase.timing_table\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m 429\u001b[0m \u001b[38;5;28;01mfor\u001b[39;00m schedulable \u001b[38;5;129;01min\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mschedulables\u001b[38;5;241m.\u001b[39mvalues():\n\u001b[1;32m 430\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mabs_time\u001b[39m\u001b[38;5;124m\"\u001b[39m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;129;01min\u001b[39;00m schedulable:\n\u001b[1;32m 431\u001b[0m \u001b[38;5;66;03m# when this exception is encountered\u001b[39;00m\n\u001b[0;32m--> 432\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mValueError\u001b[39;00m(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mAbsolute time has not been determined yet.\u001b[39m\u001b[38;5;124m\"\u001b[39m)\n\u001b[1;32m 433\u001b[0m operation \u001b[38;5;241m=\u001b[39m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39moperations[schedulable[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124moperation_repr\u001b[39m\u001b[38;5;124m\"\u001b[39m]]\n\u001b[1;32m 435\u001b[0m \u001b[38;5;28;01mfor\u001b[39;00m i, op_info \u001b[38;5;129;01min\u001b[39;00m chain(\n\u001b[1;32m 436\u001b[0m \u001b[38;5;28menumerate\u001b[39m(operation[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mpulse_info\u001b[39m\u001b[38;5;124m\"\u001b[39m]),\n\u001b[1;32m 437\u001b[0m \u001b[38;5;28menumerate\u001b[39m(operation[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124macquisition_info\u001b[39m\u001b[38;5;124m\"\u001b[39m]),\n\u001b[1;32m 438\u001b[0m ):\n", "\u001b[0;31mValueError\u001b[0m: Absolute time has not been determined yet." ] } ], "source": [ "sched.timing_table" ] }, { "cell_type": "markdown", "id": "f600c83f", "metadata": {}, "source": [ "## Device configuration\n", "\n", "Up until now, the schedule is not specific to any qubit implementation.\n", "The aim of this section is to add device-specific information to the schedule.\n", "This knowledge is contained in the {ref}`device configuration`, which we introduce in this section.\n", "By compiling the schedule to the quantum-device layer, we incorporate the device configuration into the schedule (for example by adding pulse information to every gate) and thereby enable it to run on a specific qubit implementation.\n", "\n", "To start this section, we will unpack the structure of the device configuration.\n", "Here we will use an example device configuration for a transmon-based system that is used in the\n", "`quantify-scheduler` test suite." ] }, { "cell_type": "code", "execution_count": 8, "id": "6a81d4fe", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "['backend', 'clocks', 'elements', 'edges']" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from quantify_scheduler.backends.circuit_to_device import DeviceCompilationConfig\n", "from quantify_scheduler.schemas.examples.device_example_cfgs import (\n", " example_transmon_cfg,\n", ")\n", "\n", "device_cfg = DeviceCompilationConfig.parse_obj(example_transmon_cfg)\n", "\n", "list(device_cfg.dict())" ] }, { "cell_type": "markdown", "id": "06935383", "metadata": {}, "source": [ "Before explaining how this can be used to compile schedules, let us first investigate\n", "the contents of the device configuration." ] }, { "cell_type": "code", "execution_count": 9, "id": "29d2bdab", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " 'Schedule'>" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "device_cfg.backend\n" ] }, { "cell_type": "markdown", "id": "10afeda1", "metadata": {}, "source": [ "The backend of the device configuration specifies what function will be used to add\n", "pulse information to the gates. In other words, it specifies how to interpret the\n", "qubit parameters present in the device configuration and achieve the required gates.\n", "Let us briefly investigate the backend function:" ] }, { "cell_type": "code", "execution_count": 10, "id": "6f5d4fc4", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Help on function compile_circuit_to_device in module quantify_scheduler.backends.circuit_to_device:\n", "\n", "compile_circuit_to_device(schedule: 'Schedule', config: 'CompilationConfig | DeviceCompilationConfig | Dict | None' = None, *, device_cfg: 'DeviceCompilationConfig | Dict | None' = None) -> 'Schedule'\n", " Add pulse information to all gates in the schedule.\n", " \n", " Before calling this function, the schedule can contain abstract operations (gates or\n", " measurements). This function adds pulse and acquisition information with respect to\n", " `config` as they are expected to arrive to device (latency or distortion corrections\n", " are not taken into account).\n", " \n", " From a point of view of :ref:`sec-compilation`, this function converts a schedule\n", " defined on a quantum-circuit layer to a schedule defined on a quantum-device layer.\n", " \n", " Parameters\n", " ----------\n", " schedule\n", " The schedule to be compiled.\n", " config\n", " Compilation config for\n", " :class:`~quantify_scheduler.backends.graph_compilation.QuantifyCompiler`, of\n", " which only the :attr:`.CompilationConfig.device_compilation_config`\n", " is used in this compilation step.\n", " device_cfg\n", " (deprecated) Device compilation config. Pass a full compilation config instead\n", " using `config` argument. Note, if a dictionary is passed, it will be parsed to a\n", " :class:`~.DeviceCompilationConfig`.\n", " \n", " Returns\n", " -------\n", " :\n", " A copy of `schedule` with pulse information added to all gates.\n", " \n", " Raises\n", " ------\n", " ValueError\n", " When both `config` and `device_cfg` are supplied.\n", "\n" ] } ], "source": [ "help(device_cfg.backend)\n" ] }, { "cell_type": "markdown", "id": "727beb1a", "metadata": {}, "source": [ "The {ref}`device configuration ` also contains the\n", "parameters required by the backend for all qubits and edges." ] }, { "cell_type": "code", "execution_count": 11, "id": "2ab99cc8", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "['q0', 'q1']\n", "['q0_q1']\n", "['q0.01', 'q0.ro', 'q1.01', 'q1.ro']\n" ] } ], "source": [ "print(list(device_cfg.elements))\n", "print(list(device_cfg.edges))\n", "print(list(device_cfg.clocks))\n" ] }, { "cell_type": "markdown", "id": "cc4838a6", "metadata": {}, "source": [ "For every qubit and edge, we can investigate the contained parameters." ] }, { "cell_type": "code", "execution_count": 12, "id": "3a0599d4", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'reset': OperationCompilationConfig(factory_func=, factory_kwargs={'duration': 0.0002}, gate_info_factory_kwargs=None), 'Rxy': OperationCompilationConfig(factory_func=, factory_kwargs={'amp180': 0.32, 'motzoi': 0.45, 'port': 'q0:mw', 'clock': 'q0.01', 'duration': 2e-08}, gate_info_factory_kwargs=['theta', 'phi']), 'Z': OperationCompilationConfig(factory_func=, factory_kwargs={'amp': 0.23, 'duration': 4e-09, 'port': 'q0:fl', 'clock': 'cl0.baseband'}, gate_info_factory_kwargs=None), 'measure': OperationCompilationConfig(factory_func=, factory_kwargs={'port': 'q0:res', 'clock': 'q0.ro', 'pulse_type': 'SquarePulse', 'pulse_amp': 0.25, 'pulse_duration': 1.6e-07, 'acq_delay': 1.2e-07, 'acq_duration': 3e-07, 'acq_channel': 0}, gate_info_factory_kwargs=['acq_index', 'bin_mode', 'acq_protocol'])}\n", "{'amp180': 0.32, 'motzoi': 0.45, 'port': 'q0:mw', 'clock': 'q0.01', 'duration': 2e-08}\n" ] } ], "source": [ "print(device_cfg.elements[\"q0\"])\n", "print(device_cfg.elements[\"q0\"][\"Rxy\"].factory_kwargs)\n" ] }, { "cell_type": "code", "execution_count": 13, "id": "00729431", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'q0_q1': {'CZ': OperationCompilationConfig(factory_func=, factory_kwargs={'square_port': 'q0:fl', 'square_clock': 'cl0.baseband', 'square_amp': 0.5, 'square_duration': 2e-08, 'virt_z_parent_qubit_phase': 44, 'virt_z_parent_qubit_clock': 'q0.01', 'virt_z_child_qubit_phase': 63, 'virt_z_child_qubit_clock': 'q1.01'}, gate_info_factory_kwargs=None)}}\n" ] } ], "source": [ "print(device_cfg.edges)\n" ] }, { "cell_type": "code", "execution_count": 14, "id": "9ac78f5c", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'q0.01': 6020000000.0, 'q0.ro': 7040000000.0, 'q1.01': 5020000000.0, 'q1.ro': 6900000000.0}\n" ] } ], "source": [ "print(device_cfg.clocks)\n", "\n" ] }, { "cell_type": "markdown", "id": "ecda2bc0", "metadata": {}, "source": [ "Lastly, the complete example device configuration (also see {class}`~quantify_scheduler.backends.graph_compilation.DeviceCompilationConfig`):" ] }, { "cell_type": "code", "execution_count": 15, "id": "662ce733", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'backend': 'quantify_scheduler.backends.circuit_to_device.compile_circuit_to_device',\n", " 'clocks': {'q0.01': 6020000000.0,\n", " 'q0.ro': 7040000000.0,\n", " 'q1.01': 5020000000.0,\n", " 'q1.ro': 6900000000.0},\n", " 'edges': {'q0_q1': {'CZ': {'factory_func': 'quantify_scheduler.operations.pulse_factories.composite_square_pulse',\n", " 'factory_kwargs': {'square_amp': 0.5,\n", " 'square_clock': 'cl0.baseband',\n", " 'square_duration': 2e-08,\n", " 'square_port': 'q0:fl',\n", " 'virt_z_child_qubit_clock': 'q1.01',\n", " 'virt_z_child_qubit_phase': 63,\n", " 'virt_z_parent_qubit_clock': 'q0.01',\n", " 'virt_z_parent_qubit_phase': 44}}}},\n", " 'elements': {'q0': {'Rxy': {'factory_func': 'quantify_scheduler.operations.pulse_factories.rxy_drag_pulse',\n", " 'factory_kwargs': {'amp180': 0.32,\n", " 'clock': 'q0.01',\n", " 'duration': 2e-08,\n", " 'motzoi': 0.45,\n", " 'port': 'q0:mw'},\n", " 'gate_info_factory_kwargs': ['theta', 'phi']},\n", " 'Z': {'factory_func': 'quantify_scheduler.operations.pulse_library.SoftSquarePulse',\n", " 'factory_kwargs': {'amp': 0.23,\n", " 'clock': 'cl0.baseband',\n", " 'duration': 4e-09,\n", " 'port': 'q0:fl'}},\n", " 'measure': {'factory_func': 'quantify_scheduler.operations.measurement_factories.dispersive_measurement',\n", " 'factory_kwargs': {'acq_channel': 0,\n", " 'acq_delay': 1.2e-07,\n", " 'acq_duration': 3e-07,\n", " 'clock': 'q0.ro',\n", " 'port': 'q0:res',\n", " 'pulse_amp': 0.25,\n", " 'pulse_duration': 1.6e-07,\n", " 'pulse_type': 'SquarePulse'},\n", " 'gate_info_factory_kwargs': ['acq_index',\n", " 'bin_mode',\n", " 'acq_protocol']},\n", " 'reset': {'factory_func': 'quantify_scheduler.operations.pulse_library.IdlePulse',\n", " 'factory_kwargs': {'duration': 0.0002}}},\n", " 'q1': {'Rxy': {'factory_func': 'quantify_scheduler.operations.pulse_factories.rxy_drag_pulse',\n", " 'factory_kwargs': {'amp180': 0.4,\n", " 'clock': 'q1.01',\n", " 'duration': 2e-08,\n", " 'motzoi': 0.25,\n", " 'port': 'q1:mw'},\n", " 'gate_info_factory_kwargs': ['theta', 'phi']},\n", " 'measure': {'factory_func': 'quantify_scheduler.operations.measurement_factories.dispersive_measurement',\n", " 'factory_kwargs': {'acq_channel': 1,\n", " 'acq_delay': 1.2e-07,\n", " 'acq_duration': 3e-07,\n", " 'clock': 'q1.ro',\n", " 'port': 'q1:res',\n", " 'pulse_amp': 0.21,\n", " 'pulse_duration': 1.6e-07,\n", " 'pulse_type': 'SquarePulse'},\n", " 'gate_info_factory_kwargs': ['acq_index',\n", " 'bin_mode',\n", " 'acq_protocol']},\n", " 'reset': {'factory_func': 'quantify_scheduler.operations.pulse_library.IdlePulse',\n", " 'factory_kwargs': {'duration': 0.0002}}}}}\n" ] } ], "source": [ "pprint(example_transmon_cfg)\n" ] }, { "cell_type": "markdown", "id": "0fb10cde", "metadata": {}, "source": [ "## Quantum Devices and Elements\n", "\n", "The {ref}`device configuration` contains all knowledge\n", "of the physical device under test (DUT).\n", "To generate these device configurations on the fly, `quantify-scheduler` provides the\n", "{class}`~quantify_scheduler.device_under_test.quantum_device.QuantumDevice` and\n", "{class}`~quantify_scheduler.device_under_test.device_element.DeviceElement` classes.\n", "\n", "These classes contain the information necessary to generate the device configs and allow\n", "changing their parameters on-the-fly.\n", "The {class}`~quantify_scheduler.device_under_test.quantum_device.QuantumDevice` class\n", "represents the DUT containing different {class}`~quantify_scheduler.device_under_test.device_element.DeviceElement` s.\n", "Currently, `quantify-scheduler` contains the\n", "{class}`~quantify_scheduler.device_under_test.transmon_element.BasicTransmonElement` class\n", "to represent a fixed-frequency transmon qubit connected to a feedline. We show their interaction below:" ] }, { "cell_type": "code", "execution_count": 16, "id": "7746d06f", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(, ['qubit'])" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from quantify_scheduler.device_under_test.quantum_device import QuantumDevice\n", "from quantify_scheduler.device_under_test.transmon_element import BasicTransmonElement\n", "\n", "# First create a device under test\n", "dut = QuantumDevice(\"DUT\")\n", "\n", "# Then create a transmon element\n", "qubit = BasicTransmonElement(\"qubit\")\n", "\n", "# Finally, add the transmon element to the QuantumDevice\n", "dut.add_element(qubit)\n", "dut, dut.elements()\n" ] }, { "cell_type": "markdown", "id": "e3d8dc8e", "metadata": {}, "source": [ "The different transmon properties can be set through attributes of the {class}`~quantify_scheduler.device_under_test.transmon_element.BasicTransmonElement` class instance, e.g.:" ] }, { "cell_type": "code", "execution_count": 17, "id": "20486027", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "['reset', 'rxy', 'measure', 'ports', 'clock_freqs']\n", "\n", "qubit.reset: ['duration']\n", "qubit.rxy: ['amp180', 'motzoi', 'duration']\n", "qubit.measure: ['pulse_type', 'pulse_amp', 'pulse_duration', 'acq_channel', 'acq_delay', 'integration_time', 'reset_clock_phase', 'acq_weights_a', 'acq_weights_b', 'acq_weights_sampling_rate', 'acq_weight_type']\n", "qubit.ports: ['microwave', 'flux', 'readout']\n", "qubit.clock_freqs: ['f01', 'f12', 'readout']\n" ] } ], "source": [ "qubit.clock_freqs.f01(6e9)\n", "\n", "print(list(qubit.submodules.keys()))\n", "print()\n", "for submodule_name, submodule in qubit.submodules.items():\n", " print(f\"{qubit.name}.{submodule_name}: {list(submodule.parameters.keys())}\")\n" ] }, { "cell_type": "markdown", "id": "56498c2d", "metadata": {}, "source": [ "The device configuration is now simply obtained using {code}`dut.generate_device_config()`.\n", "In order for this command to provide a correct device configuration, the different properties need to be set to applicable values in the {class}`~quantify_scheduler.device_under_test.transmon_element.BasicTransmonElement` and {class}`~quantify_scheduler.device_under_test.quantum_device.QuantumDevice` objects." ] }, { "cell_type": "code", "execution_count": 18, "id": "88d6de8e", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "DeviceCompilationConfig(backend=, clocks={'qubit.01': 6000000000.0, 'qubit.12': nan, 'qubit.ro': nan}, elements={'qubit': {'reset': OperationCompilationConfig(factory_func=, factory_kwargs={'duration': 0.0002}, gate_info_factory_kwargs=None), 'Rxy': OperationCompilationConfig(factory_func=, factory_kwargs={'amp180': nan, 'motzoi': 0, 'port': 'qubit:mw', 'clock': 'qubit.01', 'duration': 2e-08, 'reference_magnitude': None}, gate_info_factory_kwargs=['theta', 'phi']), 'measure': OperationCompilationConfig(factory_func=, factory_kwargs={'port': 'qubit:res', 'clock': 'qubit.ro', 'pulse_type': 'SquarePulse', 'pulse_amp': 0.25, 'pulse_duration': 3e-07, 'acq_delay': 0, 'acq_duration': 1e-06, 'acq_channel': 0, 'acq_protocol_default': 'SSBIntegrationComplex', 'reset_clock_phase': True, 'reference_magnitude': None, 'acq_weights_a': None, 'acq_weights_b': None, 'acq_weights_sampling_rate': None}, gate_info_factory_kwargs=['acq_index', 'bin_mode', 'acq_protocol'])}}, edges={})\n" ] } ], "source": [ "pprint(dut.generate_device_config())\n" ] }, { "cell_type": "markdown", "id": "b5f3bc36", "metadata": {}, "source": [ "## Mixing pulse and circuit layer operations (Chevron)\n", "\n", "As well as defining our schedules in terms of gates, we can also mix the circuit layer\n", "representation with pulse-level operations.\n", "This can be useful for experiments involving pulses not easily represented by Gates,\n", "such as the Chevron experiment.\n", "In this experiment, we want to vary the length and amplitude of a square pulse between\n", "X gates on a pair of qubits." ] }, { "cell_type": "code", "execution_count": 19, "id": "7a302e1f", "metadata": {}, "outputs": [], "source": [ "from quantify_scheduler import Schedule\n", "from quantify_scheduler.operations.gate_library import Measure, Reset, X, X90\n", "from quantify_scheduler.operations.pulse_library import SquarePulse\n", "from quantify_scheduler.resources import ClockResource\n", "\n", "sched = Schedule(\"Chevron Experiment\")\n", "acq_idx = 0\n", "\n", "for duration in np.linspace(start=20e-9, stop=60e-9, num=6):\n", " for amp in np.linspace(start=0.1, stop=1.0, num=10):\n", " reset = sched.add(Reset(\"q0\", \"q1\"))\n", " sched.add(X(\"q0\"), ref_op=reset, ref_pt=\"end\") # Start at the end of the reset\n", " # We specify a clock for tutorial purposes, Chevron experiments do not necessarily use modulated square pulses\n", " square = sched.add(SquarePulse(amp=amp, duration=duration, port=\"q0:mw\", clock=\"q0.01\"))\n", " sched.add(X90(\"q0\"), ref_op=square) # Start at the end of the square pulse\n", " sched.add(X90(\"q1\"), ref_op=square)\n", " sched.add(Measure(q0, acq_index=acq_idx), label=f\"M q0 {acq_idx}\")\n", " sched.add(\n", " Measure(q1, acq_index=acq_idx),\n", " label=f\"M q1 {acq_idx}\",\n", " ref_pt=\"start\", # Start at the same time as the other measure\n", " )\n", "\n", " acq_idx += 1\n", "\n", "\n", "# Specify the frequencies for the clocks; this can also be done via the DeviceElement (BasicTransmonElement) instead\n", "sched.add_resources([ClockResource(\"q0.01\", 6.02e9), ClockResource(\"q1.01\", 6.02e9), ClockResource(\"q0.ro\", 5.02e9), ClockResource(\"q1.ro\", 5.02e9)]) " ] }, { "cell_type": "code", "execution_count": 20, "id": "01ea0b6f", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzEAAADECAYAAACvIa+3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA4VklEQVR4nO3dd3xUVf7/8ddJSKGHFgzNCEiXolRpCiggAhFYWX+yFFG/yAK67IrYsMG6NgTLIuICig0RxFE0uwIGEQWVIoIUpSsSQGqCJDA5vz/uJGRCAikzmUzyfj4e87hzzy3zmWEyHz4z55xrrLWIiIiIiIgEi5BAByAiIiIiIpIXKmJERERERCSoqIgREREREZGgoiJGRERERESCiooYEREREREJKipiREREREQkqKiIERERERGRoKIiRkREREREgoqKGBERERERCSoqYkRECpkxxhpjXgp0HJJ3nn+7RwMdh4hISaciRkTER4wx9YwxM40xO40xp40xJ4wxq4wxdxtjSgc6vqLAGLPbUwhkd4sPdHzFgTHmamPMo8aYqEDHIiLiL6UCHYCISHFgjOkDLABSgDeATUA40Al4BmgK3BmwAIuWDcBz2bTvL+Q48qM0cDbQQVzE1cAjwFzgWEAjERHxExUxIiIFZIy5DHgX2AN0s9b+lmnzy8aY+kCfgAR3EcaYstba5EJ+2F+ttW8W8mPmmzEmBAi31p621p4OdDwiIqLuZCIivjABKAeMzFLAAGCt/dlaOz1ruzEmzhizyRiTYozZbIzplc0+NY0xs40xiZn2uy3T9urGmLPGmEeyObahp5vWGM/6cM96V2PMv40xB4FfMu0/2nP+FGPMfmPMy1m7JBljEjwxNzHGfG6MOWWM+dUYMyFvL1nOjDHRxphDnscymdrrG2OSjTHzs4nnKmPMV8aYP4wxu4wxo7I5b4Qx5jFjzM+e57jPGPO0MSYiy37WGPOSMeZWY8xmnF/XemXa9mimfR/1tDUwxrxpjDnuif0J46htjPnQ07XwgDHm7z6IK8f3jSe2ZzyruzJ11YvN0z+CiEgRp19iREQKri+w01r7VR6O6QQMAP4NnATGAQuNMXWstb+DU6AAqwELvAQcAnoD/zHGVLDWTrPWJhpjVgA3A49leYzBgBunm1tm//ac63GgrOexHsXpgrQUmAE0BO4C2hhjOlprz2Q6vhIQDywC3gMGAU8ZY36w1n6ai+ceZoypmk17srX2D2vtQWPMXZ64xwIveH4Nmet5rUZnOa4S8IknlndwXosZxphUa+1sz/MLAVw4r/urwBbgCuBvQAMgLss5u3nO8xJwGNh9kec033POiTi/uj0EHAH+D1gO3AfcCjxrjPnWWvtFPuO62Ptmkee4WzznOOw57tBF4hcRCS7WWt1000033fJ5AyrgFBmL83CMxfl2v16mtuae9jGZ2l7DGSdSJcvx7+CMdSjtWb/Tc2yzLPttBpZlWh/u2W8lEJqpvZonnv8CIZna/+rZf0SmtgRP218ytYUDvwHv5+K57/Ycn91tYpZ93waSgcuBf3j26Z9ln/R4xmeJZz2QCIR52obgFHSdshz/f57jr87y7+MGmuTwb/dopvVHPW0zM7WFAvuANOC+TO1RwClgbqa2vMaVm/dN+msVG+i/D9100003f93UnUxEpGAqeJYn83jcUmvtjvQVa+1G4ARQF8DTjWog8JFntWr6DafYqAhc6Tl8Ec5g88Hp5zPGNAOa4PxCkNUsa60703oPnP/4T7PWpmXezxNT1vE8SUDGmBZrbSrwTXrsubAGuC6b2ztZ9hsDHAfeB54A5llrP8zmfGeBmVnimQlEA1d5mv+E8yvH1iyv5XLP9muznHOFtfbHXD4fcArO9Md3A98BBvhPpvZjwDa8X6e8xnXB942ISEmh7mQiIgVzwrMsn8fj9mbTdhSnaxQ4v45E4fzKktOsZtEA1trDxphlON2fHvZsG4zzn/tF2Ry3K8v6pZ7ltsyN1tpUY8zOTNvT/WKttdnE3jyHOLM6bK1derGdrLVHjDHjcLqVJeJ0ncrOfnv+5ATbPctYnC55lwONyblbVXSW9ayv0cVk/fc8Dpy21h7Opr1KpvW8xnWx942ISImgIkZEpACstSeMMfuBZnk81J1De/pA9vRfyt8EXs9h342Z7r8LzDHGtLTWbsApaJZl859ogD/yGGtWF4vdl3p6lpWAWuR/yuAQ4AdgfA7b92VZz+trlN1rkpvXKa9xFeZrLyJSZKmIEREpuI+BO40xHay1X/vonIdwuqiF5uZXC2AxTheqwZ4JvRoAT+bysfZ4lg2BnemNxphw4DKcwf6FzjPr1u3A0ziD4l83xrSz1ma9TksNc/5U0Q08y92e5Q6gBU5hl/VXpEDyR1xF6fmJiPiFxsSIiBTc0zgD0F/zzCjmxRhTzxhzd15O6BlXsRAY6BnfkvWc1bLsfwxnrMzNwJ+BVJzCJjeWevYfl3lKY2AkztibJXmJ3Rc8Uzu/hjPW5gGcYuZKz/2sSuEMhE8/NtyzfghY62l+D6gJ3JHNY5U2xpT1Yfh54Y+40ou5qALEJSJSpOmXGBGRArLW7jDG/D880+waY94ANuEMlr8aZ/D23HyceiLOwO41xphZwI9AZZz/zPfw3M9sPk73s9HAfz2FTW7iP2SMeRJniuV4Y4wL51eZ0cC3ZBrE7yM1jTFDsmlPstYu9tyfjjN2pIenoIs3xrwGPGSM+dBa+32m4/YD93muhbIdZzxQS+BOe25q6Hk4Bd4rxphrgVU4s4g18rT3xBmMX9j8EVd64TbFGPMucAb4KJtxQyIiQUtFjIiID1hrXcaY5sC9QH+ca6yk4Ixb+TvOTF95PWeiMaYtMAnn2iCjgd9xpk6+L5tDXDhjOcqT/axkF3qsR40xh3BmBHse5xonrwIPWO9rxPhCS5z/vGe1B1hsjOkHDAX+bq3dmmn7eJxZzF43xrTJFNdRYBjwIs4vGok4Uw5nvObW2jRjTBzOtVOGAjfhTHe8E6dgSp8IoFD5Iy5r7bfGmIeBUTgX6QzB6RaoIkZEig1TtLoGi4iI5J4xJgGoaq3N68QKIiISxDQmRkREREREgoqKGBERERERCSoqYkREREREJKhoTIyIiIiIiAQV/RIjIiIiIiJBRUWMiIiIiIgEFRUxIiIiIiISVIrqxS41UEdEREREpOQxudlJv8SIiIiIiEhQUREjIiIiIiJBRUWMiIiIiIgEFRUxIiIiIiISVFTEiIiIiIhIUFERIyIiIiIiQUVFjIiIiIiIBBUVMSIiIiIiElRUxIiIiIiISFBRESMiIiIiIkFFRYyIiIiIiAQVFTEiIiIiIhJUVMSIiIiIiEhQUREjIiIiIiJBRUWMiIiIiIgEFRUxIiIiIiISVFTEiIiIiIhIUFERIyIiIiIiQaVUoAMQkeLjxeW7Ax1CQIztFpvvY7du3eq7QIJIo0aNAh2CiPiYckDeKQfkn36JERERERGRoKIiRkREREREgoqKGBERERERCSoqYkREREREJKioiBERkWwlJyezZ8+eQIchIiIBUNRzgIoYERE5T1paGhMmTKBXr15s2bIl0OGIiEghCoYcoCJGRETOM2vWLJYvXw5AaGhogKMREZHCFAw5QEWMiIh4+eKLL5g+fXrGelhYWACjERGRwhQsOUBFjIiIZNi7dy/jx4/3agsPDw9QNCIiUpiCKQeoiBEREQBOnTrF6NGjOX36NNbajPaimsBERMR3gi0HlAp0ACIiEnjWWh588EF27dpFWlqa17aimsBERMQ3gjEHqIgRERHmzZtHfHx8ttuKagITERHfCMYcoO5kWSQkJBAbG5unY2bPns3WrVv9E5BICZaWlsbk4d35+D/PeLVvXr2cv/VswPcrs//Albw7cuQIkP0sNEU1gfmDcoBI0aEcUHiCMQeoiMmFl19+mdjYWCIjI2nXrh3ffPON1/Y1a9Ywc+bMAEUnUnyFhIRw3S2jWfnhPP5IOgHAvu2bmPPEWPrdcR8tOvcKcITFx91338348eNxu92A89obYwgJCSmy02sWFuUAkcBQDig8wZgDfN6dzBgzABgFXAVUBlpZazf4+nEKy/z58xk/fjyvvPIK7dq1Y9q0afTs2ZNt27YRHR0NQP/+/Rk7dizPP/98gKMVKX5a9+hP/LzpfPHhG7TpcRMzHxpJ+16DuHbQyECHVuRZazlw4ADbt2/PuP3yyy+cPHmSpKQkjDGUK1eO8uXLU7t2bbZs2UJISAhpaWlUrlyZZs2aERMTE+inEVDKASKBpRyQf8U9B/hjTExZ4EvgPWCWH85fqKZOncodd9zBiBEjAHjllVdYsmQJs2fPZuLEiQB0796dxMRENm3aRLNmzQIZrkixExpaih6DR/HxnOdYt/wj6jRszoC/PhLosIqsU6dOsXTpUlwuF6tWrcr1cevXr/daP3z4MABt2rTh9OnTREZG+jTOYKEcIBJYygF5U5JyQJ6LGGNMWWAGMAA4CTwL9AU2WGvvsdbO8+wX68M4AyI1NZW1a9dy//33Z7SFhITQo0cPvv7664y2iIgIrr/+elwulxKYiB+07tGfRf9+HCpXY/iDLxASkruesHMnj6NzvyEsmTuVk0d/JzQ0lJ5/GUurrn38HHHhcrvdrF69GpfLhcvl8tl5ExISSEhIAGDgwIH069eP1q1b5/r1D3bKASJFg3LAhZXUHJCfX2KeAboC/YGDwD+BK4ENvguraDh8+DBut5vq1at7tVevXv28QZz9+/dnxowZPPDAA4UZokiJsOAF51u35ONHMLn88ExLS+Pw/j1UrRnLgNGTqFW/CSeOHOKZUX1p0vZaIkqX8WfIhcJay2effcaDDz5IUlJStvs0btyYpk2b0qBBAxo2bEhsbCxlypTJGKiZmppKcnIyI0eOZOfOnVhriYiIICUlxes8CxcuZOHChVStWpXJkyfTtWtXvz+/QFMOECkalAOyV9JzQJ6KGGNMOWAkMMRau8zTNgz4Jb8BGGMigIjMbT///LNfZ0KoXbu2z895ww03MGLECA4fPkzVqlWz3Wffvn0+f1yR4m7JnOfYvOZzxr+0iJfu/QtffzKfLnFDM7Yn7t3B8+MGUbpcBSLLlOXQr3uoVb8JN931ELUvb0bFKtFUrOKMXahQuRplK1bi1MljPk1ggfjb3rp1K0888QTr1q07b1v37t2Ji4ujS5cuF/0sDQ8P59ixY+zYsQMAYwwTJ04kLi6Ozz//nA8//JAVK1Zk7H/48GFGjRpFx44deeCBB6hbt26+4g/k56FygEjwUA7IXnHOAbn9jM7rLzH1gHBgTXqDtfaIMWZbHs+T2f2AV+fGLl26FOB0F/frr7/mar+qVasSGhpKYmKiV3tiYiKXXHKJV9uuXbuIiooiKioqx/O1b98+z7GKBJOJ83Lf/zY3vlryDssXvMaYZ9+mZr0mXDNwBMvmz6TjjbcQWioMgOp16hHbuCX97riPGnUbMXlYN+761+skLJpNo9advc63d/sP2LQ0KkXX8GmcBfnbXrZsWZ72P3r0KNOnT2f+/Ple7XXq1GHEiBH06tXrgp9D2XG5XBmDOUNDQ+nVqxeRkZH07t2b3r17c+TIET755BNmzZrFwYMHAVi1ahV9+vRh6NChjBkzhvLly+fpMQP5eagcIOIfygF5pxxwvtx+RheFi10+CUzN3PDFF18cLwpzUoeHh3PVVVexbNky4uLiAOfnyWXLljFmzBivfV0uFzfccAOlSuX8kq5evdqf4YoE3OKf3D471+Y1n7PgxUcZ/tALXNakFQBd+g9l+fxX+eazD+jQ++aMfRP37SS6dl1ST/9BWloaEaXLsH3dV3Tpf+7buuQTx3jzX3/nz39/0mcxpivI33ZycnKu9924cSO33norZ8+e9Wp//PHHGTBgQL6mwbTWsmjRoozk1bVr1/MSYOXKlRkyZAiDBw9m/vz5TJkyJWPbG2+8wQcffMC8efNo2LBhrh83GD4PlQNE8kY5IO+UA/Ivr0XMDuAM0A7YC2CMqQQ0AFZc4LgcWWtTgJSL7hgg48ePZ9iwYbRu3Zq2bdsybdo0kpOTM2aqSedyuXj44YcveC5/dGEQKVJ+2u2T0+zd/gNzHh9D/zsm0qJTz4z20uUq0OWm4Sx9Zwbtrh9ISGgop08lUSosnFJh4ezZ8T3V69Tj9KkkrE2jdLkKAJxJTeG1SXfS45ZR1G16lU9izKwgf9u5vUjiggULmDRpklfbbbfdxqhRo/L8DVhmGzduzPjWy+12Z/xnPTthYWEMGTKEPn368NJLL/H2228DcPLkSeLi4nj22Wfp0yd3A2aD5fNQOUAkD5QD8kw5IP/yVMRYa5OMMf8BnjHG/I4zsH8KkJa+jzGmMlAHSP+trqExBuCAtfZAgSMuZIMHD+bQoUNMmjSJAwcO0LJlS+Lj470Geu7atYtt27bRq5cuuiTiC3UaXMGzSzZnu63PiPH0GTE+Y/3Anp+45NL6ACTu/ZnoWpexff1XXN6yA+B8y/TW0/dyeauraXvdAP8H7wcHDx7kqaeeyliPiYlh1qxZ1KtXr8DndrlchIaG4na7KVeuXK6681aqVImHH36YQYMGMWLECI4fPw7AlClT6NSpExUrVixwXEWFcoBI4VMO8KYckL38zJF2L7AS+AhYinNNmLWZtvcD1gNLPOvvetZH5T/MwBozZgx79uwhJSWFNWvW0K5dO6/tLpeLa665pkCVsIjkz2+7fyImtgEA4ZFl2Lr2S7Z8+wWNrnL6Qu/c9B3rEz7mh1X/46k7b+CpO29g/87cffNVVERHRzN1qtPrtmfPnnz00Uc+SV6pqam4XC7cbjehoaHceOONeZpUpXHjxvz3v/+lU6dOALzwwgvFqoBJpxwgUnQpB+RfsOeAPI+JsdYmAX/x3AAwxvTJtH0uMNcHsQUNl8vFgAHBWd2LBLvM/aJbdulNyy69eeaufvxp7GMA1LuiDdOX7gxUeD7TpUsXFixYQNOmTfH8ul1gK1euzJiW0+12069fvzyfo2LFisycOZMtW7bQtGlTn8QVbJQDRAJHOSD/gj0HFIWB/UEvrzNLiIh/3TvDdxf7Kkp8fSHFxYsXZ3QjqFmzJi1btszXeUJCQkpsAQPKASJFjXJA7gR7Digal9wsQmJjY7nnnnsCHYaIiF8dP36chIQE3G43ISEh3HTTTT77di+YKQeISElQHHKAT36JsdZe44vzFAVKYCJSEsTHx2dM05mWlpavbgTFkXKAiJQExSEHqDuZiEgJtGjRooxv3a644gpN/ysiUoIUhxygIkZEpITZt28fGzduzFgfOHBgAKMREZHCVFxygMbEiIiUMPHx8Rn3S5UqRc+ePS+wt4iIFCfFJQeoiBERKWHatGlDlSpVAOjWrVuxvLaLiIhkr7jkAHUnExEpYVq2bMmKFSvYvHkzjRo1CnQ4IiJSiIpLDlARIyJSAoWGhtK8efNAhyEiIgFQHHKAupOJiIiIiEhQUREjIiIiIiJBRUWMiIiIiIgEFRUxIiIiIiISVIy1NtAxZKdIBiUiIiIiIn5lcrOTfokREREREZGgoiJGRERERESCiooYEREREREJKipiREREREQkqKiIERERERGRoKIiRkREREREgoqKGBERERERCSoqYkREREREJKioiBERERERkaCiIkZERERERIKKihgREREREQkqKmJERERERCSoqIgREREREZGgoiJGRERERESCiooYEREREREJKipiRLKRkJCAMYZjx44FOpQi5dFHH6Vly5aBDkNExK+UA7KnHCBFiYoYCUrDhw/HGMOoUaPO2/bXv/4VYwzDhw/36WPOmjWLFi1aUK5cOaKiomjVqhVPPvmkTx/Dn9KTcvqtevXqDBw4kJ07dwY6NBGRPFEOyDvlACluVMRI0Kpduzbvvvsuf/zxR0bb6dOnefvtt6lTp45PH2v27Nncc889jBs3jg0bNrBq1SomTJhAUlKSTx8nP9xuN2lpabnef9u2bezfv58FCxawefNm+vbti9vt9mOEIiK+pxzgUA6QkkpFjAStK6+8ktq1a7No0aKMtkWLFlGnTh1atWqV0ZaSksK4ceOIjo4mMjKSTp068e2333qd65NPPqFBgwaULl2aa6+9lt27d3ttd7lc3HzzzYwcOZL69evTtGlTbrnlFqZMmZKxj9vtZvz48URFRVGlShUmTJjAsGHDiIuLy9gnNjaWadOmeZ27ZcuWPProoxnrU6dO5YorrqBs2bLUrl2b0aNHeyXKuXPnEhUVhcvlokmTJkRERLB3715SUlL4xz/+Qc2aNSlbtizt2rUjISHhvNctOjqamJgYunTpwqRJk/jxxx/5+eefM86b2eLFizHG5PAv4Hyz17ZtW8qWLUtUVBQdO3Zkz549Gds//PBDrrzySiIjI6lbty6PPfYYZ8+ezfF8IiK5pRygHCAlm4oYCWq33XYbc+bMyVifPXs2I0aM8NpnwoQJLFy4kNdff51169ZRv359evbsyZEjRwDYt28fAwYMoG/fvmzYsIHbb7+diRMnep3jkksuYfXq1V4fzlk999xzzJ07l9mzZ/Pll19y5MgRPvjggzw/p5CQEF544QU2b97M66+/zvLly5kwYYLXPqdOneKpp57itddeY/PmzURHRzNmzBi+/vpr3n33XTZu3Mif/vQnevXqxU8//ZTjY5UuXRqA1NTUPMd59uxZ4uLi6Nq1Kxs3buTrr7/mzjvvzEh4K1euZOjQodx99938+OOPzJw5k7lz53olfRGRglAOUA6QEsxaWxRvIhc0bNgw279/f3vw4EEbERFhd+/ebXfv3m0jIyPtoUOHbP/+/e2wYcNsUlKSDQsLs2+99VbGsampqbZGjRr26aefttZae//999smTZp4nf++++6zgD169Ki11tr9+/fb9u3bW8A2aNDADhs2zM6fP9+63e6MY2JiYjLOaa21Z86csbVq1bL9+/fPaLv00kvt888/7/VYLVq0sI888kiOz3XBggW2SpUqGetz5syxgN2wYUNG2549e2xoaKj99ddfvY7t3r27vf/++6211n7++efnPaerr77a1qxZ06akpNg5c+bYihUreh3/wQcfWOdjwvHII4/YFi1aWGut/f333y1gExISso27e/fu9p///KdX27x582xMTEyOz1VEJDeUA5QDpFjLVb1QKjClk4hvVKtWjT59+jB37lystfTp04eqVatmbN+xYwdnzpyhY8eOGW1hYWG0bduWLVu2ALBlyxbatWvndd4OHTp4rcfExPD111+zadMmvvjiC7766iuGDRvGa6+9Rnx8PCdPnuS3337zOk+pUqVo3bo11to8PaelS5fy5JNPsnXrVk6cOMHZs2c5ffo0p06dokyZMgCEh4fTvHnzjGN++OEH3G43DRo08DpXSkoKVapU8WqrVasW1lpOnTpFixYtWLhwIeHh4XmKEaBy5coMHz6cnj17ct1119GjRw9uvvlmYmJiAPj+++9ZtWrVed0tsj4XEZH8Ug5wKAdISaQiRoLebbfdxpgxYwB4+eWX/fpYzZo1o1mzZowePZpRo0bRuXNnVqxYwZVXXpmr40NCQs5LaGfOnMm4v3v3bm688UbuuusupkyZQuXKlfnyyy8ZOXIkqampGR/6pUuX9uqnnJSURGhoKGvXriU0NNTr/OXKlfNaX7lyJRUqVCA6Opry5cvnOrbszJkzh3HjxhEfH8/8+fN56KGH+Oyzz2jfvj1JSUk89thjDBgw4LzjIiMjL3heEZHcUg5QDpCSSUWMBL1evXqRmpqKMYaePXt6batXrx7h4eGsWrWKSy+9FHA+lL/99lvuueceABo3bozL5fI6bvXq1Rd93CZNmgCQnJxMxYoViYmJYc2aNXTp0gVw+guvXbvWK7lVq1aN3377LWP9xIkT7Nq1K2N97dq1pKWl8dxzzxES4gxZe++99y4aS6tWrXC73Rw8eJDOnTtfcN/LLrvsvMGb6bGdPHmS5ORkypYtC8CGDRty9ditWrXi/vvvp0OHDrz99tu0b9+eK6+8km3btlG/fv2LnkNEJL+UA5QDpGRSESNBLzQ0NKNbQNZvoMqWLctdd93FvffeS+XKlalTpw5PP/00p06dYuTIkQCMGjWK5557jnvvvZfbb7+dtWvXMnfuXK/z3HXXXdSoUYNu3bpRq1YtfvvtNyZPnky1atUyuh3cfffd/Otf/+Lyyy+nUaNGTJ069bwLpXXr1o25c+fSt29foqKimDRpklfM9evX58yZM7z44ov07duXVatW8corr1z0NWjQoAG33norQ4cO5bnnnqNVq1YcOnSIZcuW0bx5c/r06XPRc7Rr144yZcrwwAMPMG7cONasWXPe65DZrl27ePXVV+nXrx81atRg27Zt/PTTTwwdOhSASZMmceONN1KnTh0GDRpESEgI33//PZs2bWLy5MkXjUdEJDeUA5QDpITK7eCZQr6JXFD6oM6cpA/qtNbaP/74w44dO9ZWrVrVRkRE2I4dO9pvvvnGa/+PPvrI1q9f30ZERNjOnTvb2bNnew2AfP/99+0NN9xgY2JibHh4uK1Ro4YdOHCg3bhxY8Y5zpw5Y++++25boUIFGxUVZcePH2+HDh3qFefx48ft4MGDbYUKFWzt2rXt3LlzzxvUOXXqVBsTE2NLly5te/bsad944w2vWLIbfGmtM1h10qRJNjY21oaFhdmYmBh70003ZcSYdVBndj744ANbv359W7p0aXvjjTfaV199NcdBnQcOHLBxcXEZr8mll15qJ02a5DXQNT4+3l599dW2dOnStkKFCrZt27b21VdfzfHxRURyQzmg4nnPWTlAipFc1QvG5nHAWSEpkkGJ5NXw4cM5duwYixcvDnQoIiJSyJQDRPIl54sTZaLrxIiIiIiISFBRESMiIiIiIkFF3clERERERKSoyFV3Ms1OJlJUHNwFrmfgxwQwIXBFd+h7L1SpFejIRETE35QDRPJEv8SIFBVLnocP/undVqYijJkHDTpkf4yIiBQPygEi6XL1S4yKGJGiwlr4dDpc3h7cblj4OOxaB6XC4Z750KhToCMUERF/UQ4QSaciRqTIS0mG8DJgsvl7TTkFf28Kp5OcfV7aBSGai0NEpNhQDhDJjqZYFiny5oyD6bfA8YPnb/tls5O8AFJPwbYvCzc2ERHxL+UAkXxTESMSKBv/B9+5nEGcx347f/sPS7OsLyuUsEREpBAoB4gUiIoYkUA4kwLvPOjcv24UXNri/H12f+8s67Vxlj+tgR9XQJq7cGIUERH/UA4QKTBNsVxCvLh8d6BDCIix3WIDHUL2lv8HDu2GitHQ9x/Z77N3o7PsNhJ2fAu718PUQTB0KnT5S6GFmhd6n+Xd1q1bfRdIEGnUqFGgQyhR9LdZxCgHFCvKAXnnixygX2JEClvyMWcqTYC4+yGy3Pn7JB2FE4ec+y16QkQZsGnO+odPOYNBRUQk+CgHiPiEihiRwvbpC3DqGNRsDB1vyX6fgzudZaUYJ8HVbOKsl68KxxNh2axCCVVERHxMOUDEJ1TEiBQm91nY8Klzf8BDEBKa/X6JngQWXddZ1mjgLC9v7yw/fdH5pk5ERIKHcoCIz6iIESlMoaXgkc/h9hnQ/Lqc9zu8x1lGX+Ysq9dzlqXCoVZT+OOE822eiIgED+UAEZ9RESNS2MIiof2g7C9ulu73fc6ySi1neUl9Z5m4AwY84Nxf/hoc2e+/OKXES05OZs+ePYEOQ6R4UQ6QIFHUc4CKGJHCsmm505UgN37/xVlW9iSw9C4FB3dBsx5Ol4LKNeFENhdIE/GBtLQ0JkyYQK9evdiyZUugwxEJfsoBEkSCIQdoimWRwvDTapg22Bmc+dD/ICziwvtn/Rau2qXO8o8TzoDQ/5sF5apAqTC/hSwl26xZs1i+fDkAoaE59NsXkdxRDpAgEww5QL/EiPhbmhte+6tz/5J6F09e1sJRz9WbK9d0luGlIeoS5/7BXc59JS/xky+++ILp06dnrIeF6b0mkm/KARJkgiUHqIgR8be598Dve537276CP05eeP/ko3DmtHM/PWkBVIt1lod2n2s7ewb+NwO+c/koWCnp9u7dy/jx473awsPDAxSNSDGgHCBBJJhygIoYEX86vBe+mn9uPel3WPnWhY9J/wauXBVnAGi6jASWaZDdF6/De5PgrYmQdMQnIUvJderUKUaPHs3p06ex1ma0F9UEJlLkKQdIEAm2HKAiRsRfrIUX/wJYCA2DPz3qtK/76MLHHfMksEox3u3pCexwpgTW+S9QoyGcPARv3ec8pkg+WGt58MEH2bVrF26322tbUU1gIkWacoAEkWDMASpiRPxl+Wvw64/O/a7DoHV/5/7OtXA6Kefj0r+Fy9yNAM4N7Dy4+1xbWASMeNG5YNq3i2Hlm76IXEqgefPmER8fT1pa2nnbimoCEynSlAMkiARjDlARk0VCQgKxsbF5Omb27Nls3brVPwEFibS0NCYP787H/3nGq33z6uX8rWcDvl8ZH6DIAii2peeOgV5jnFlmqtR2Bnnu+C7n4457psw8L4HFOsvDWeZsv6wV3OS5bsBb98H2rwoYeNGh91XhOXLE6YqS3Sw0RTWB+YNyQP7obzUbygEFpvdV4QnGHKAiJhdefvllYmNjiYyMpF27dnzzzTde29esWcPMmTMDFF3REBISwnW3jGblh/P4I+kEAPu2b2LOE2Ppd8d9tOjcK8ARBsC2Vc6yYcdzM8zUa+Ms92zI+bhjB5xlxere7enfwh3dD2dSvLf1HANX9QX3GXhxiPNNXzGg91Xhufvuuxk/fnxGN4KQkBCMMYSEhBTZ6TULi3LAxelvNRvKAQWm91XhCcYc4PPrxBhjDPAYcAcQBawC7rLW/uTrxyoM8+fPZ/z48bzyyiu0a9eOadOm0bNnT7Zt20Z0dDQA/fv3Z+zYsTz//PMBjjawWvfoT/y86Xzx4Ru06XETMx8aSfteg7h20MhAh1Z4Vr3rzOHf4//g6wVOW4dB57bXbgbfLIK9m3I+x/FEZ5n1W7jyVSGiLKQkO4NFYy4/ty0kBEa+DCcPw85152a2KQb0vso/ay0HDhxg+/btGbdffvmFkydPkpSUhDGGcuXKUb58eWrXrs2WLVsICQkhLS2NypUr06xZM2JiYi7+QMWYckDu6W8V5QA/0Psq/4p7DvDHxS4nAOOAYcAu4Angv8aYJtbaoPurmjp1KnfccQcjRowA4JVXXmHJkiXMnj2biRMnAtC9e3cSExPZtGkTzZo1C2S4ARUaWooeg0fx8ZznWLf8I+o0bM6Avz4S6LAKx4EdsPAJWL8EjIHIcvDbdmdu/yv7ntuvdlNnuS8XCSzrt3DGQHQs7NvsXCcgcwID57HGvQ17f4DSFeDTF5w46raGS5sX+CkGSol+X+XDqVOnWLp0KS6Xi1WrVuX6uPXr13utHz58GIA2bdpw+vRpIiMjszus2FMOyL0S/beqHOA3Jfp9lQ8lKQfkuYgxxpQFZgADgJPAs0BfYAPwN+AeYLK19kPP/kOBRCAOeNcHMRea1NRU1q5dy/3335/RFhISQo8ePfj6668z2iIiIrj++utxuVwlOoGB843Jon8/DpWrMfzBFwgJKaY9FtPS4NAu2P41rP0YNi93ZoUJLQU3/gN+WuPsd+WNUKbCueNqNHKWh3bD2VQolU0/05y6EgBE1/UksJ3ZxxVZDn75Ed6533uWmsq1oNtIuOpGqHqpkwyDSH7fV3Mnj6NzvyEsmTuVk0d/JzQ0lJ5/GUurrn38HHHhcrvdrF69GpfLhcvlu+tFJCQkkJCQAMDAgQPp168frVu3Lr5/11koB+SdcoBygD8oB1xYSc0B+fkl5hmgK9AfOAj8E7gSp4i5DLgEWJq+s7X2uDFmDdCBXBYxWad287Xc9u07fPgwbreb6tW9P0iqV69+3iDO/v37M2PGDB544IEcz+fv51UULHjB+XYk+fgRTBF4k/vlNXefwbw0lJAfP/dqtldcR9pND0DpCoQ81B4DuLsOh8wxVIgmJLIc5nQS7t9+dqbGzCwtjZATh5xjy1f1PhYw1WIJAdIO7MBm99x+XkPIO/djrMU26w6/bsEc3Q9HfoH3H4P3HyOt2+3YXmOzT5BFVH7eV2lpaRzev4eqNWMZMHoSteo34cSRQzwzqi9N2l5LROkyPosvUH/b1lo+++wzHnzwQZKSsp/tqHHjxjRt2pQGDRrQsGFDYmNjKVOmTMZAzdTUVJKTkxk5ciQ7d+7EWktERAQpKd597hcuXMjChQupWrUqkydPpmvXrgWKPZCfh8oB/qMcoBzgD8oB2SuuOSC3n9F5KmKMMeWAkcAQa+0yT9sw4BfPLukdOBOzHJqYaVvWc0YAEZnbZsyYkZew8mzUqFE+P+cNN9zAiBEjOHz4MFWrVs12n1mzZvn8cXOtgf8Hvy2Z8xyb13zO+JcW8dK9f+HrT+bTJW5oxvbEvTt4ftwgSperQGSZshz6dQ+16jfhnukL/BaTP17zqsm/MmjLueT1R2gZvqtxLZsj2sOSlVy/4x3qus/wa/nL+GjZOli2zuv4m0KjqE4SS9+cwa5KTb22RZ5JZrj7DAD/WeAiLcT7T7Th4X1cC+xft5KPk7M8N2sZuOUVqlnLtiqt+DyiK9S7hha/reTKA18Q4XZ6c4Ysf423D5bjRGQVH70imfjhfZbf99VNdz1E7cubUbFKNBWrOGMXKlSuRtmKlTh18phPE1hB3mfXXHNNvo7bunUrTzzxBOvWrTtvW/fu3YmLi6NLly4XnVUmPDycY8eOsWPHDgCMMUycOJG4uDg+//xzPvzwQ1asWJGx/+HDhxk1ahQdO3bkgQceoG7duvmKP5Cfh8oB/qEcoBygHJB3ygHny+1ndF5/iakHhANr0hustUeMMdvyeJ7M7ge8Ojc+8cQTBTjdxeX2xalatSqhoaEkJnrXZImJiVxyiXdNtmvXLqKiooiKisrxfP5+XhcycZ5/E9hXS95h+YLXGPPs29Ss14RrBo5g2fyZdLzxFkJLhQFQvU49Yhu3pN8d91GjbiMmD+vGXf963a9x+eM1rxFpOV0vhCF1nLnUS7tP0XnfEpruWMLOZEPdSpYzaTDi031sfu/8x6/W/CwDa8J3nyzg3zsXeW1rWC6N4Z3hSCo8NuXJ845tHZXGtR0g5OAOnshy7q5V0xjV5izJZ2HQe5s4krqZiQ3O0qGe95zvnx4w3Pvpy4DvuxP4+n1WkPdVwqLZNGrd2et8e7f/gE1Lo1J0DZ/GWZD3WV4T2NGjR5k+fTrz58/3aq9Tpw4jRoygV69eF/wcyo7L5coYzBkaGkqvXr2IjIykd+/e9O7dmyNHjvDJJ58wa9YsDh50pn9dtWoVffr0YejQoYwZM4by5cvn6TED+XmoHOB7ygHKAaAckB/KAefzVxFzMZ6OnFQHfsvUXh2nu1l2ngSmZm545513joeFhfk4tLwLDw/nqquuYtmyZcTFxQHOz5PLli1jzJgxXvu6XC5uuOEGSpXK+SV9//33/RnuBX33h//OvXnN5yx48VGGP/QClzVpBUCX/kNZPv9VvvnsAzr0vjlj38R9O4muXZfU03+Qlpbm029CsuO319ymsf7EAarsWkn1LZ8SfuoIlcOhcrjFYtjTdSyP3dEj20NrrXsH1r/D7X270a3zWK9tFX9dD/GPEBldh/fff+m8Y0NTkuDN/0eN0rDo7TdICz/3+jX8bDLs/YYTLfry6v/dAUC5xK3YJRM5Wrs1R2I7cKJ6U6LKV+d9P/WH9uX7rKDvq+3rvqJL/3Pf1iWfOMab//o7f/77+f8xKKjC+tveuHEjt956K2fPnvVqf/zxxxkwYEC+psG01rJo0aKM5NW1a9fzEmDlypUZMmQIgwcPZv78+UyZMiVj2xtvvMEHH3zAvHnzaNgwS9eYCwjk52FuKQfkjnKAckA65QD/Ug7wltciZgdwBmgH7AUwxlQCGgArcGYjOwB0x1O0GGMqePbPto+YtTYFSMluW1Ewfvx4hg0bRuvWrWnbti3Tpk0jOTk5Y6aadC6Xi4cffviC5+rQoYM/Q72g75bv9st5927/gTmPj6H/HRNp0alnRnvpchXoctNwlr4zg3bXDyQkNJTTp5IoFRZOqbBw9uz4nup16vklpsz8/5oPBKbBlpUw/2H4ZTMGS/1Ovaiffj2ArEL2w/p3qE4S1bPG99VeAMrEXJZz7B9Fw/GDtKsdBXWvctqOJ8Ic57oAMe17E5NxbAe4tg+VK0ZTuSBPM5d89T4r6Pvq9KkkrE2jdDlnQO2Z1BRem3QnPW4ZRd2mV/kkxswK8j7L7UUSFyxYwKRJk7zabrvtNkaNGpXnb8Ay27hxI7/++ivg9FFO/896dsLCwhgyZAh9+vThpZde4u233wbg5MmTxMXF8eyzz9KnT+4GzAby8zAvlAMuTDlAOSAz5YC8Uw7IvzwVMdbaJGPMf4BnjDG/4wzsnwKkebZbY8w04CFjzE+cm2J5P7C4wNEGwODBgzl06BCTJk3iwIEDtGzZkvj4eK+Bnrt27WLbtm306lXyLrpUp8EVPLtkc7bb+owYT58R4zPWD+z5iUsurQ9A4t6fia51WaHEWCgad4ZHE85dIyCn5AXnrrp8aPf529JnpYm6wIDLWk2dKzrv23wuga35wLkKdGQ5eP1vULmGc4E1gIrReXsuRUBB31fb13/F5S2dD0hrLW89fS+Xt7qattcN8H/wfnDw4EGeeuqpjPWYmBhmzZpFvXoF/0+gy+UiNDQUt9tNuXLl6NKly0WPqVSpEg8//DCDBg1ixIgRHD9+HIApU6bQqVMnKlasWOC4igrlgAtTDvBQDvAp5QBvygHZy8/UIfcCK4GPcGYh+xLIfGnYp4EXgVeBb4FyQK9gvEZMujFjxrBnzx5SUlJYs2YN7dq189rucrm45pprClQJlwS/7f6JmNgGAIRHlmHr2i9JPn40wFH5WMc/w3WZ+nKeSfGe5hLOJbCjv0Fqlt/ec7o+QGa1PVO47t14rm21Z2Ds6SQIC4ewojefu79k977a8u0XNLrK6Qu9c9N3rE/4mB9W/Y+n7ryBp+68gf07c/fNV1ERHR3N1KlOr9uePXvy0Ucf+SR5paam4nK5cLvdhIaGcuONN150EGhmjRs35r///S+dOnUC4IUXXihWBUw65QDfUA7wUA7wKeWA/Av2HJDnMTHW2iTgL54bAMaYPpm2W2CS51YiuFwuBgwIzuq+MGXuv9qyS29adukdwGgKwZFf4aWh0HkIXJup60m5ys63ZaeT4PA+qNHg3LaMb+GynczPcZnTN5gd3znL/dvOJbPQMBj75rlv50qA7N5Xz9zVjz+NfQyAele0YfrSHK6pEES6dOnCggULaNq0KcZH/dlXrlyZMS2n2+2mX79+eT5HxYoVmTlzJlu2bKFp06YXP6AYUg7IHeUAD+UAn1IOyL9gzwG+HthfIi1btizQIUhRtP4TJ7HMfxgadICanoucGeN8E7dvk9OdIHMCS/8W7kJdCeq3dZa//ghJR+GT6ee23foUNLjal88iKN07w3cX+ypKfH0hxcWLF2d0I6hZsyYtW7bM13lCQkJKbAEDygGSA+WAgFEOyJ1gzwGBvxJVERMbG8s999wT6DCkOOh2OzTrDmdTYM44p79yumqXOsvDe7yPyfgWLibn81asDjUbO10UVr8PaxY67fXbOd/4ieTC8ePHSUhIwO12ExISwk033eSzb/eCmXKA+IxygBRhxSEH+KSIsdZeY629xxfnCjQlMPEZY2D4NChdAXavh+X/ObctY2BnpgSWlub0kQaodIEEBnBVX2f57gNg08CEwOg5zmOK5EJ8fHzGNJ1paWn56kZQHCkHiM8oB0gRVhxygH6JEfGnqEtgkGd42If/ghOHnPvZzU6T9Dt4rtR8wUGd4PSvLl3h3Pr1d0GFar6IWEqIRYsWYYzBGEPz5s2pXbt2oEMSKX6UA6SIKg45QEWMiL91HgKXtoA/TsJiz0W20rsSZE5gxzzfwFWoBqUuMjtI+arwj0VOd4W7ZsOfHvV11FKM7du3j40bN2KtxVrLwIEDAx2SSPGlHCBFTHHJASpiRPwtJBT+PNm5v3OdM+VmtOf6CIf2OF0IAI7sd5aValz8nMcOQM0m8P+ePNetQCSX4uPjM+6XKlWKnj17XmBvESkQ5QApYopLDlARI1IYLm8Pf3sPHl4KYRFQuRaEloIzp88N5Px9n7OsUuvC57IWZt4Bj3SG3Rv8GrYUT23atKFKlSoAdOvWrVhe20WkSFEOkCKkuOQATbEsUliaXnvufmgpqFoHEnfCwZ3O1ZV//8XZVuUi/VJ/+Ax+Wu1czKxC8F2JWQKvZcuWrFixgs2bN9OoUaNAhyNSMigHSBFRXHKAfokRKWxnTjvTYkZ7rrabuMNZ5uZbuLQ0WPRP5363253EJ5IPoaGhNG/ePE9XZxYRH1AOkCKgOOQA/RIjUpjcZ+Gxa+HAz9DSc7XqAz87y4O7nGXVS3M+/ptF8MtmZ1aa3uP8G6uIiPiWcoCIz+iXGJHCFFrqXOLavd5Z7t/q9HE+uNNZr14v+2PPpMDifzn3e4+FcpX8G6uIiPiWcoCIz6iIESlsvcdBmahzgzl/2QLHEyHllHPBsmo5fAv3+Rzn6s4Vq0P3OwotXBER8SHlABGfUBEjUtjKRkGfv51bP54IW1Y696Mvy/n6ALvWOsv+90FEWb+GKCIifqIcIOITxlob6BiyUySDEvGZMykwqdO5C53VbQ07v4O2N8Gdr2Z/jLWwdSU07Ohcd0BERIKTcoDIhZjc7KRfYkQCISwCbplybn3nd87y0pY5H2MMNO6i5CUiEuyUA0QKTEWMSKA0vx7qtfFuu6KH93ryMXhzAiQdKbSwRESkECgHiBSIihiRQBr3FlSLde437Ag1GpzbduY0zLgNEubAy8MCEp6IiPiRcoBIvmlMjEigHT8IG/8HV/aFFXPg8vbgdsPCx2HXOogoA/d9DHWuCHSkIiLia8oBIlnlakyMihiRomLJ8/DBP73bylSEMfOgQYfAxCQiIoVDOUAkXa6KmFL+jkJEcqlNHPz2E/yY4Fwr4Iru0PdeqFIr0JGJiIi/KQeI5Il+iRERERERkaJCUyyLiIiIiEjxoyJGRERERESCiooYEREREREJKipiREREREQkqKiIERERERGRoKIiRkREREREgoqKGBERERERCSoqYkREREREJKioiBERERERkaCiIkZERERERIKKihgREREREQkqpQIdQA5MoAMQEREREZGiSb/EiIiIiIhIUFERIyIiIiIiQUVFjIiIiIiIBBUVMSIiIiIiElRUxIiIiIiISFBRESMiIiIiIkFFRYyIiIiIiAQVFTEiIiIiIhJUVMSIiIiIiEhQ+f+WTN3GZj/6IwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = sched.plot_circuit_diagram()\n", "ax.set_xlim(-0.5, 9.5)\n", "for t in ax.texts:\n", " if t.get_position()[0] > 9.5:\n", " t.set_visible(False)\n" ] }, { "cell_type": "markdown", "id": "08afde3d", "metadata": {}, "source": [ "This example shows that we add gates using the same interface as pulses. Gates are Operations, and\n", "as such support the same timing and reference operators as Pulses.\n", "\n", "\n", "## Device and Hardware compilation combined: Serial Compiler\n", "\n", "{class}`~quantify_scheduler.backends.graph_compilation.SerialCompiler` can be used to execute the device and hardware compilation separately, or execute both in one call. Here we will not set the hardware configuration thus only executing device compilation. The {ref}`Compiling to Hardware ` tutorial demonstrates how to set the hardware configuration.\n", "\n", "{class}`~quantify_scheduler.backends.graph_compilation.SerialCompiler` requires a {class}`~quantify_scheduler.backends.graph_compilation.CompilationConfig` and this holds both the device and hardware configurations (when set). In the example below, we generate a {class}`~quantify_scheduler.backends.graph_compilation.CompilationConfig` via {meth}`~quantify_scheduler.device_under_test.quantum_device.QuantumDevice.generate_compilation_config` of {class}`~quantify_scheduler.device_under_test.quantum_device.QuantumDevice`." ] }, { "cell_type": "code", "execution_count": 21, "id": "0792bdf2", "metadata": {}, "outputs": [], "source": [ "from quantify_scheduler.backends import SerialCompiler\n", "from quantify_scheduler.device_under_test.quantum_device import QuantumDevice\n", "from quantify_scheduler.device_under_test.transmon_element import BasicTransmonElement\n", "\n", "dut.close()\n", "dut = QuantumDevice(\"DUT\")\n", "q0 = BasicTransmonElement(\"q0\")\n", "q1 = BasicTransmonElement(\"q1\")\n", "dut.add_element(q0)\n", "dut.add_element(q1)\n", "dut.get_element(\"q0\").rxy.amp180(0.6)\n", "dut.get_element(\"q1\").rxy.amp180(0.6)\n", "\n", "compiler = SerialCompiler(name='compiler')\n", "compiled_sched = compiler.compile(schedule=sched, config=dut.generate_compilation_config())" ] }, { "cell_type": "markdown", "id": "a1d3d5e3", "metadata": {}, "source": [ "So, finally, we can show the timing table associated with the Chevron schedule and plot\n", "its pulse diagram:" ] }, { "cell_type": "code", "execution_count": 22, "id": "65340573", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/rsoko/.anaconda3/envs/tmp/lib/python3.9/site-packages/quantify_scheduler/schedules/schedule.py:451: FutureWarning: The behavior of DataFrame concatenation with empty or all-NA entries is deprecated. In a future version, this will no longer exclude empty or all-NA columns when determining the result dtypes. To retain the old behavior, exclude the relevant entries before the concat operation.\n", " timing_table = pd.concat(timing_table_list, ignore_index=True)\n" ] }, { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 portclockis_acquisitionabs_timedurationoperationwf_idx
0Nonecl0.basebandFalse0.0 ns200,000.0 nsReset('q0','q1')0
1Nonecl0.basebandFalse0.0 ns200,000.0 nsReset('q0','q1')1
2q0:mwq0.01False200,000.0 ns20.0 nsX(qubit='q0')0
3q0:mwq0.01False200,020.0 ns20.0 nsSquarePulse(amp=0.1,duration=2e-08,port='q0:mw',clock='q0.01',reference_magnitude=None,phase=0,t0=0)0
4q0:mwq0.01False200,040.0 ns20.0 nsX90(qubit='q0')0
5q1:mwq1.01False200,040.0 ns20.0 nsX90(qubit='q1')0
6Noneq0.roFalse200,060.0 ns0.0 nsMeasure('q0', acq_index=0, acq_protocol=\"None\", bin_mode=None)0
7q0:resq0.roFalse200,060.0 ns300.0 nsMeasure('q0', acq_index=0, acq_protocol=\"None\", bin_mode=None)1
8q0:resq0.roTrue200,060.0 ns1,000.0 nsMeasure('q0', acq_index=0, acq_protocol=\"None\", bin_mode=None)0
9Noneq1.roFalse200,060.0 ns0.0 nsMeasure('q1', acq_index=0, acq_protocol=\"None\", bin_mode=None)0
10q1:resq1.roFalse200,060.0 ns300.0 nsMeasure('q1', acq_index=0, acq_protocol=\"None\", bin_mode=None)1
\n" ], "text/plain": [ "" ] }, "execution_count": 22, "metadata": {}, "output_type": "execute_result" } ], "source": [ "compiled_sched.timing_table.hide(slice(11, None), axis=\"index\").hide(\n", " \"waveform_op_id\", axis=\"columns\"\n", ")\n" ] }, { "cell_type": "code", "execution_count": 23, "id": "37e1441e", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAGwCAYAAACtlb+kAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAABbC0lEQVR4nO3dd3wUZeI/8M/sZtMbIaRRQmhBjhKK5hCQYHKA51FORUR+AsqJIqAciByooJ4nRdGI3qGCNA/lDg+RrwKCSIIFQ42AKCUGg0ASWnqy9fn9sZlJlmzCtmQ3y+f9eu3L7Mzs7DPZxfnkqZIQQoCIiIjIi6ncXQAiIiKixsbAQ0RERF6PgYeIiIi8HgMPEREReT0GHiIiIvJ6DDxERETk9Rh4iIiIyOv5uLsAnsBkMuHChQsICQmBJEnuLg4RERHZQAiB0tJSxMXFQaVquA6HgQfAhQsX0LZtW3cXg4iIiBxw7tw5tGnTpsFjGHgAhISEADD/wkJDQ91cGiIiIrJFSUkJ2rZtq9zHG8LAAyjNWKGhoQw8REREzYwt3VHYaZmIiIi8HgMPEREReT0GHiIiIvJ6DDxERETk9Rh4iIiIyOsx8BAREZHXY+AhIiIir8fAQ0RERF6PgYeIiIi8HgMPEREReT0GHiIiIvJ6DDxERETk9Rh4yCUqy0vdXQQiIqJ6MfCQ0/atfAo+S+Nx6nCmu4tCRERkFQMPOS0y/xtoJCOuHv/S3UUhIiKyioGHnNbKmA8AkIrz3FwSIiIi6xh4yCnF1y4jHGUAgIDyc24uDRERkXUMPOSUS3knlZ9baC+6sSRERET1Y+Ahp5Tmn1F+jjYVwGgwuLE0RERE1jHwkFO0l35RfvaVjLh08az7CkNERFQPBh5yilRk2VH5yrlTbioJERFR/Rh4yCkBZZYdlcsLctxUEiIiovox8JBTwnTmIennpWgAgPEah6YTEZHnYeAhpwSYygEAV33jAACStsSdxSEiIrKKgYecEigqAQBV/uYaHpWuzJ3FISIisoqBhxwmTCYEogoAYAiOAQCoDeXuLBIREZFVDDzksIryEqgkAQBQhcYCAHwYeIiIyAMx8JDDKkuLAQBGIcEnxNykpTFUuLNIREREVjHwkMMqy4sAABXwhyYwFADgZ2LgISIiz8PAQw6rKjePyKqUAuDLwENERB6MgYccpq2u4alUBcIvyBx4AgUDDxEReR4fdxeAmi9DRSkAQKsKRHBQOICaYer1KSu+gq/f/BtM14oauXRERN5LCgxAv8efRVTbRHcXpdlg4CGHGSrNTVo6dSCCQsIBAP6SHga9Dj4aX6uvOfjft9D+w2+aqohERF7rgO4F3P3qR+4uRrPBwEMOM1aZa3j06iAEVgceACgvLUZYRCurr9GVmEd2FbbSoCz5lkYvIxGRt9H8+Ava5pZBlHMaEHsw8JDDhNYceAyaIGh8/aAVGvhJelSWFdUbeGA0AgDK2rfC3a/9p6mKSkTkNbYteBjI/R4wmtxdlGaFnZbJcVrzMhImTTAAoEIKAABUlRfX+xJRHXiEil89IiJHSCq1+QcTA489eNchx+ksA0+lLYHHZA48YOAhInKIpK7+/ydreOzCuw45TKWvXijULwgAUFUdePQV9a+YLtfwMPAQETmINTwO4V2HHKbWmzvMSX4hAACtOhBAzegtqwzVgUfNrx4RkSMkNQOPI3jXIYdpqhcKVfmbA49eba7pabCGR/4HysBDROQQycc83khik5ZdeNchh2mM5sDjUx14DD7mGh5T9XB1q9iHh4jIKey07BjedchhfkbzMhKawDAAgKG687I8essaIf9FIlfJEhGRXeROy5JJuLkkzQsDDznMv3oZCXmldHm0ltCxhoeIqLFI6uop9FjDYxfedchhAdWBxy/IXMMDX3PgUenqr+GRh1FKDDxERA6ROy1LDDx24V2HHCYvFOpfvVI6/KoDj77+6c7ZaZmIyDk1o7TYpGUP3nXIIUaDAX6SHgDgX92kJfmaR2mpDQ2smC4HHhX78BAROaKmhoeBxx4MPOSQivKaoecBQeZRWqrqCQh9jA0EnuqJByXW8BAROUSlqh6WziYtu/CuQw7Rlps7JpuEBD9/83B0tS2BhzU8RETOYZOWQzwy8Pzzn/9E+/bt4e/vj+TkZOzfv7/eY1euXIlBgwahRYsWaNGiBdLS0ho8nlxDW2XumFwFX6UDshx4NKaq+l/ITstERE5RyU1anHjQLh531/nPf/6DWbNmYeHChTh8+DB69eqFYcOGobCw0OrxGRkZGDduHPbs2YN9+/ahbdu2GDp0KM6fP9/EJb+5aCvMgadS8le2+fibOy37ihsHHviwhoeIyBHsw+MYjws8r7/+Oh599FE8/PDD6NatG9555x0EBgZi9erVVo/fsGEDnnjiCSQlJaFr165YtWoVTCYTdu/eXe97aLValJSUWDzIPrpKc5OWtlbg8Q0wBx6/hmp4qpu0JE48SETkELWPLwAGHnt5VODR6XQ4dOgQ0tLSlG0qlQppaWnYt2+fTeeoqKiAXq9HREREvccsWrQIYWFhyqNt27ZOl/1mo680Dz23Fnj8YUPgYR8eIiKHKDU8goHHHh4VeC5fvgyj0Yjo6GiL7dHR0cjPz7fpHHPnzkVcXJxFaLrevHnzUFxcrDzOnTvnVLlvRkatOfDoVDWBxy/QPForQGjrfZ3EGh4iIqewScsxPu4ugCstXrwYGzduREZGBvz9/es9zs/PD35+fk1YMu9jqF4vy1A78MhNWpIeRoMBah8rXy+upUVE5BRV9dISKgYeu3hUDU9kZCTUajUKCgosthcUFCAmJqbB17722mtYvHgxdu7ciZ49ezZmMQmAqbqGR68OULbJ8/EAQGWF9fW05L9I2KRFROQYeS0t1vDYx6MCj6+vL/r27WvR4VjugNy/f/96X7d06VL8/e9/x44dO9CvX7+mKOpNTw48Rp+awOPnHwiTkAAAVfUEHnZaJiJyjloJPG4uSDPjcU1as2bNwsSJE9GvXz/cdtttSE9PR3l5OR5++GEAwIQJE9C6dWssWrQIALBkyRIsWLAAH374Idq3b6/09QkODkZwcLDbrsPbCX0FAMBYq4ZHUqlQAT8EoQraegNPdQ0PAw8RkUPYadkxHhd4xo4di0uXLmHBggXIz89HUlISduzYoXRkzsvLg6rWpHUrVqyATqfDfffdZ3GehQsX4oUXXmjKot9cdObAI2rV8ABAlSQHHusrprPTMhGRc1Rqjfm/bNKyi8cFHgCYPn06pk+fbnVfRkaGxfOzZ882foGoDqm6hsekCbLYrpX8AVGszNNT53XVnZZVDDxERA5R+XCUliM8qg8PNR8qgznwwDfQYrs8L4+hqtz6C+V/oOy0TETkkJoaHjcXpJlh4CGHqAzmBUKl6wKPPC9PfYFHbtJSaTyycpGIyOOpfcyBh3147MPAQw5RK4HHsklLHqYuz9NzPbkKVqVi4CEickTNPDxuLkgzw8BDDvExmgOP6roaHnkiQnnY+vXkv0jYaZmIyDEqHwYeRzDwkEM01QuE+vhb1vDI8/LUG3iUYems4SEicoRazSYtRzDwkEM0JnMNj9rPcq4jeV4eeZ6e6ylNWqzhISJyCGt4HMPAQw7xNZkXCNUEWAYeoalu4tLdKPCwhoeIyBFyDY+KFTx2YeAhh/gJc5OW5romLZOPOfBIN6jhYR8eIiLHsIbHMQw85BB/mAOP73U1PPK8PMo8PdepqeHRNF7hiIi8mFxDrmYNj10YeMghAcLcpOUXGGKxXVICT6XV13GUFhGRc9S1/mA0Gg1uLEnzwsBDdtNpq6CRjAAA/8BQi33yvDzqegKPylhdw+PDPjxERI6o/f9Po07nxpI0Lww8ZLfKWguDBgRZ1vCo/cyBR56n53pyDY+aTVpERA7x0fgpPxuMDDy2YuAhu2krzAuDGoQKGo2vxT5V9TB1TX2Bp7qTHWt4iIgcU3uUq4lNWjZj4CG7VZUXAwAqpABIKsuvkDxM3ddUT5MWh6UTETnFx6fmD02jQe/GkjQvDDxkN215CQCgEv519mkCzE1c/qK+Jq3q/zLwEBE5xLKGh4HHVgw8ZDdddZNWlSqwzj6/oDAA9QceuYZHzcBDROQQNWt4HMLAQ3bTV5kDj1YVUGefX/WorYDqiQmvJ0+UJbEPDxGRQ1QqFeQ5B9mHx3YMPGQ3Q6W5SUtvJfAEBFUHHkkHg77u6AGO0iIicp6p+u7NGh7bMfCQ3UxV5mHpOp+gOvsCQ8KUnyvKS+vsl2t42GmZiMhxJqn6v6zhsRkDD9nNpDUHGaNP3T48vr7+0AvzLMqVZUV19suL3al9WMNDROQowRoeuzHwkN2EthwAYNTUreGRVCpUSubRW1Ws4SEiahRyDQ8Dj+0YeMhuks7cpCWs1PAAQAXMfXt0FSV19smBh6O0iIgcZ1KZE4+Rw9JtxsBDdpP05hoe4Rtsdb88ektbUWyx3WQyKV849XUzNBMRke3kTsvCwD48tmLgIbupqwOP5Ndw4DFUlllsNxpqRm2xDw8RkeNqangYeGzFwEN2UxsqANQfeHRqc1OXocqyD0/t0QQqDksnInKYqO7DIxh4bMbAQ3bTGM01PGr/EKv79dWBx1hpGXgMrOEhInIJuYaHw9Jtx8BDdpNXQvepJ/DIw9WFzrJJy7KGh52WiYgcJQ9LN7EPj80YeMhufiZzk5ZPgPUmLXm4utBe34enZjSBj5qdlomIHCUkuYbH6OaSNB/8M7sZKyypwn8PnkOV3nTjg+2UktgK/dpHWN3nbzKvk+UbGGZ1vzxcXdKVW2yv/ZeIimtpERE5TChNWhyWbivedZqxFZk5WPPt2UY598eHfsP381Ot7vOHuUnLL9B6k5Y8XF0evi6rPV8Em7SIiBwn9+ERrOGxGe86zdi1cnMn4NvaR6BbXKhLzlmhM+C/B3/D1Yq6C3/KAkUlIAH+wdZreOTRW2q99WHpJsm82i8RETlGruGpPd0HNYyBpxnTGsxNWSOS4vDQ7+Ndcs4rZVr89+Bv0BlMEEJAqm4nlum0VfCVzH9R+AfdIPBUD1+XyU1aJqnOS4iIyA5yHx7W8NiOf2Y3Y3Lg8VO77mP09ak5l3z+2irLamZPDgq2XqskD1fXGC0Dj9ykZeS3jojIKULNYen24q2nGdMazMneT+O6j9HPR13r/HUDT0X1CuhVQgOfepaHkIer+xrrqeHht46IyCk1nZZZw2Mr3nqaMZ1cw+Pjuo9Ro5Ygt2LprAQebYV5MsFKKaDec8jD1X1NlRbb5RoeIbFNi4jIGULptMwaHlsx8DRjSpNWrVoZZ0mSpAQouQaptqpyc5NWpeRf7znk4er+1wUe1vAQEbkG+/DYj7eeZkxbPf+OrwtreADAVy0Hnro1PLrSawCACpX1IekAEBBqnr8nWFw3D4+RgYeIyBWUJi0Ta3hsxVtPM6b04XFx4PHTmGuMtFYmNNSVXQEAVPnUH3hCwlsBAIKlSuh1WmW7HHjYpEVE5BxRPbWHsFITT9Yx8DRjjdGkZT6f+WuhM9YNPMbyq+Z9GutD0gEgJDxS+bnk2iXlZ9bwEBG5CDst2423nmZM6bTswlFaQE3g0err/kMyVZibtPR+4fW+Xu3jgxKY19MqK7pc81rW8BARuYQ8LB0mBh5bMfA0Y3INj68L5+EBAN/qGiNrfXikKnPgMTUQeACgVDKP1KosrlvDI7c9ExGRg+QmLdbw2IyBpxlrjHl4gFo1PFYCj0914JECWzR4jgq1eVLCqtIryjaT3jwsnU1aRETOUTotG9hp2Va89TRTJpOA3igANGIfHiuBR6M3D0tXB1pfSV1W6WMOPPpagcfIGh4iItdQmf+/zxoe2zHwNFO1OxQ32igtK73//fUlAACf4JYNnkNf3anZWHFV2SYYeIiIXEKZeJB9eGzGwNNM1R4y3pTz8AQazYHHL6ThwGOQ+/hUd3IGao/SYuAhInKK3HeTNTw2Y+BppuTaF5UE+Lg4QMh9gqyN0goWZQCAgLDIOvtqM/mHA6jp5AzUVL2yhoeIyDmCnZbtxsDTTNWeg0dy8TDv+ubhMRmNCBXmtbSCqycXrI8UYO7UrNbWrK6uBB4OSycico4ceNikZTMGnmZK20hz8AA1naCvn2m5rLQIasncUTqkRcOBRx1kbvLy1dcEHg5LJyJyEdbw2I2Bp5lqrGUlap/z+j48ZdWzJlcIP/gHBDV4Dk11Hx9/Q4myjU1aREQuIvfhMdXta0nWMfA0U8qkg40aeCz/ciivnkRQnlSwIf7VgSfIWLeGR/7LhIiIHMQaHrvxztNMyc1Nrp6Dx3xO6/PwVFUHnvIGVkqXBYSZm7xCqjs5A6zhISJyGfn/o6zhsRkDTzMldyhulCYtjfWlJbSl5nWxKhtYKV0WXL2AaCjKYdDrAAAmo3mmZQYeIiInqasnHuRq6TZj4Gmm5CHjjRF46puHx3DtNwBApX/0Dc8R3jIGOmH+B3n54q8AalW9skmLiMg58v9HOUrLZrzzNFON2odHY70Pj1RkDi76kLY3PIfaxweFqigAwJXfTgOo1aSlZg0PEZEzpOoaHhjZpGUrBp5mqvY8PK5WXx+egPLzAAB1RLxN57nmGwsAKC/IAVC7Dw+/dkRETlHL8/Aw8NiKd55mqjGHpfvWMyw9THcRABAY1cGm81QGtQYAGK/KTVocpUVE5BLVi4ey07LteOdppnTKxIONUcNTd+JBYTIh2lgAAGjRupNN5zGGtQMA+JTkmc9RXfXKTstERM6RuJaW3Rh4mimlD4+6aebhuVL4G/wlPYxCQlQb2wKPJjIBABBUeQFArSnQWcNDROQc1vDYjXeeZkqZh6cxl5ao1aR1+dwpAMAlKRIaXz+bzhMcbW76iqhuCqvptMyvHRGRMyT5D0d2WrYZ7zzNVFP04andabks39zx+KpvjM3niWzTBQAQJa5Ap62qGT7JJi0iIuf4sIbHXgw8zZSuCUZp1a7h0V86AwAoD2ht83laRrVGhfCDShK4mHuC8/AQEbmIUsPDwGMz3nmaqaaehyfs4ncAAFPrfjafR1KpkOuXCADI/2FnzYygDDxERE7hPDz2452nmWrc1dIt+/CUFl9FZ91PAIC2t46w61wlre8wn/PXDHZaJiJylepOyxJreGzmVXeef/7zn2jfvj38/f2RnJyM/fv3u7tIjaZm4sHGn4fnTNY2aCQjzklxiEvoate5WiX9EQDQufwIhL56Hp5GaIYjIrqZSD4+5h9Yw2Mzrwk8//nPfzBr1iwsXLgQhw8fRq9evTBs2DAUFha6u2iNonHn4anptGwymiAd+QAAcKFlf7vP1aH773EZ4QiSqlB2xTw8nTU8RETOkZQaHuHmkjQfPu4ugKu8/vrrePTRR/Hwww8DAN555x18/vnnWL16Nf72t7+5pUwXz/6E4wcy0KZFADQWQ7FrfUGFqLtd1PMFrrU96uyvGFRVDv+ffsaZ0jDr57a67cbnrjIYcWfVaUgQ+GzRWnQrPoKfRCB0vW7DmSN7rL++AT8FpeCWwu0IuPwbABW02grk/pgFlcprvn5ERE2qoqzY/ENpGb7butK9hbFRaFQbdP/9XW57f0mI+u6uzYdOp0NgYCA+/vhjjB49Wtk+ceJEFBUV4dNPP7U4XqvVQqvVKs9LSkrQtm1bFBcXIzQ01GXl+uTF/4euHx1y2fm8RU4vPf50yyV3F4OIqNn66nw4Yr8OdHcx7PJTogb3fHrUpecsKSlBWFiYTfdvr/gT+/LlyzAajYiOjrbYHh0djZ9//rnO8YsWLcKLL77Y6OUyqDUoDWj0t2kakgvmzhECWl8gJFaFKwiDqr7aJiIialB0Swn5LYGgCneXxHZGX/dGDq8IPPaaN28eZs2apTyXa3hcbcxzayCeFcg8dQlTPjgEncGEz2YMRPfWYTd+sRVHfyvCyLe/ha+PCu891BeDu7SC5IogQkREzUoLAL9zdyHsdJub398reo9GRkZCrVajoKDAYntBQQFiYurODOzn54fQ0FCLR2ORJAkpiVH4Qzdz7dPmw+cdPpf82mG/i0FKYhTDDhERkY28IvD4+vqib9++2L17t7LNZDJh9+7d6N/f/pFFjeGe3uYZirf+cB4GB4YR6o0m/N8PFyzORURERLbxmiatWbNmYeLEiejXrx9uu+02pKeno7y8XBm15W53dGmFED8fXC7T4cylMnSNsa9WKedSGa6U6xDi54NBnSMbqZRERETeyWsCz9ixY3Hp0iUsWLAA+fn5SEpKwo4dO+p0ZHYXjVqFhFZBOPpbMX69UmF34Dl72dwzrUNUMHy42jgREZFdvCbwAMD06dMxffp0dxejXu0iAqsDT7ndr5VfEx/RvIYhEhEReQJWFTSh9i2DAAC/XrF/HOHZ6te0b8nAQ0REZC8GnibUrjqsOBJ48q6WV58jyKVlIiIiuhkw8DQhpYbnqv1NWnIfHtbwEBER2Y+BpwnJYeX8tUpl8U9baA1GXCyuBADEs4aHiIjIbjZ1Wo6IiLDrpJIk4fDhw4iPj3eoUN6qVYgfAjRqVOqNOF9UiYRI28LLb9cqYRJAoK8akcG+jVxKIiIi72NT4CkqKkJ6ejrCwm68JIIQAk888QSMRqPThfM2kiShTYsAnC4sw/lrtgee89fMtTttWwRydmUiIiIH2Dws/YEHHkBUVJRNx86YMcPhAnm7FoHmGpriSr3Nr5GPDQ/UNEqZiIiIvJ1Ngcdksm8phNLSUocKczMIDTD/ykuqbA888rGhAQw8REREjrC50/Jnn31md/ChuuTQUmJHDU9JpcH8Wn8GHiIiIkfYHHhGjx6Ntm3b4tlnn8WZM2cas0xeTQ4t9jRpyTU8YazhISIicojNgSc3NxePPfYYNm7ciMTERAwePBgffPABKisrG7N8Xkep4bGjSUsOR3JzGBEREdnH5sDTtm1bLFiwADk5Ofjyyy/Rvn17TJ06FbGxsXj88cdx4MCBxiyn1whTmrQMNr9Gbv5ikxYREZFjHJp4cMiQIVi3bh0uXryIV199FceOHcPvf/979OrVy9Xl8zqh/uZaGvuatKr78LBJi4iIyCFOtZGEhIQgNTUVv/76K37++WecOHHCVeXyWs40abEPDxERkWMcquGprKzE+vXrkZKSgs6dO2Pjxo2YNWsWzp496+LieZ8wB0ZplSpNWuzDQ0RE5Ai77qDff/89Vq9ejf/+97/Q6XS455578OWXX2LIkCGNVT6vUzNKy44+PJyHh4iIyCk2B55u3brh5MmT6N27NxYtWoQHH3zQpqUmyJK9Ew8KIZQOzmzSIiIicozNgSctLQ0fffQROyY7SQ4tOoMJVXoj/DXqBo+v0pugM5onfGQNDxERkWNsDjzLly9vzHLcNIJ8faCSAJMw9+O5UeCRa4JUEhDk2/CxREREZJ1NnZb79OmDa9eu2XzSgQMH4vz58w4XypupVBJC/G0fqaXMwROg4UrpREREDrKphic7Oxs//PADIiIibDppdnY2tFqtUwXzZmEBGhRX6m3quMwh6URERM6zuUkrNTUVQgibjmVNRMOUjss2DE1XRmhxlmUiIiKH2RR4cnNz7T5xmzZt7H7NzSLUriYteZZlzsFDRETkKJvuovHx8Y1djptKsJ/5115adeMmrVKtweI1REREZD+HZlom58gjs6r0xhseq60+5kajuYiIiKh+DDxu4K8x/9q1BtMNj5VDkb8PAw8REZGjGHjcwJ4aniq9qfo1/KiIiIgcxbuoG9gXeNikRURE5CyHAk9RURFWrVqFefPm4erVqwCAw4cPc7JBG/n7mH/tcu1NQ6oM5sDjx8BDRETkMLuH/hw9ehRpaWkICwvD2bNn8eijjyIiIgKbN29GXl4e1q9f3xjl9Cp+bNIiIiJqUnbfRWfNmoVJkybh9OnT8Pf3V7b/8Y9/xN69e11aOG+lNGmx0zIREVGTsDvwHDhwAI899lid7a1bt0Z+fr5LCuXt5Noa+2p4GHiIiIgcZXfg8fPzQ0lJSZ3tp06dQqtWrVxSKG8n19bYNA+PQe60zCYtIiIiR9l9Fx05ciReeukl6PXmZREkSUJeXh7mzp2Le++91+UF9EZybY3Wlk7LHKVFRETkNLsDz7Jly1BWVoaoqChUVlZi8ODB6NSpE0JCQvCPf/yjMcrodZQmLQM7LRMRETUFu0dphYWFYdeuXfjmm29w9OhRlJWVoU+fPkhLS2uM8nklh+bhYadlIiIihzm8IuXAgQMxcOBAV5blplHTafnGTVqVes7DQ0RE5CybAs/y5cttPuGTTz7pcGFuFn52dFpmkxYREZHzbAo8b7zxhsXzS5cuoaKiAuHh4QDMMy8HBgYiKiqKgccGXC2diIioadlUbZCbm6s8/vGPfyApKQk//fQTrl69iqtXr+Knn35Cnz598Pe//72xy+sVAnztmHiwumNzAAMPERGRw+xuJ3n++efx1ltvITExUdmWmJiIN954A88995xLC+et5LW0dAYTTCZR73FGk4DeaN7PGh4iIiLH2R14Ll68CIPBUGe70WhEQUGBSwrl7WqHF20DtTy1m7zYh4eIiMhxdt9FU1NT8dhjj+Hw4cPKtkOHDmHq1Kkcmm6j2oGnoX48FoGHw9KJiIgcZnfgWb16NWJiYtCvXz/4+fnBz88Pt912G6Kjo7Fq1arGKKPXUaskaNQSgIYnH5T7+PiqVVCppCYpGxERkTeyex6eVq1aYdu2bTh16hR+/vlnAEDXrl3RpUsXlxfOm/n7qKE3Ghqci6dKmYOHzVlERETOcHjiwS5dujDkOMFPo0ap1mBTkxY7LBMRETnH7sDzyCOPNLh/9erVDhfmZlIz23JDgYeTDhIREbmC3YHn2rVrFs/1ej2OHz+OoqIi3HnnnS4rmLermXyw/iYtLdfRIiIicgm7A88nn3xSZ5vJZMLUqVPRsWNHlxTqZmBTDY+BTVpERESu4JK2EpVKhVmzZtVZgoLq52/Delps0iIiInINl91Jc3JyrE5ISNYpTVoNDUtnp2UiIiKXsLtJa9asWRbPhRC4ePEiPv/8c0ycONFlBfN2NU1aDQ1LN+/zYx8eIiIip9gdeI4cOWLxXKVSoVWrVli2bNkNR3BRDT8bVkyvVGp42KRFRETkDLsDz549exqjHDedmj48N554kE1aREREzrG76uDOO+9EUVFRne0lJSUclm4HW0ZpaVnDQ0RE5BJ230kzMjKg0+nqbK+qqsLXX3/tkkLdDOR+OQ2tlq41sg8PERGRK9jcpHX06FHl5xMnTiA/P195bjQasWPHDrRu3dq1pfNivj7mrGkw1h949AZhcSwRERE5xubAk5SUBEmSIEmS1aargIAAvPXWWy4tnDfzrV4tXddA4NEZzU1aGjUDDxERkTNsDjy5ubkQQqBDhw7Yv38/WrVqpezz9fVFVFQU1Go2vdhKDjF6W2p4qsMREREROcbmwBMfHw/AvIwEOU9T3Uylqw411shhiDU8REREzrEp8GzduhV33XUXNBoNtm7d2uCxI0eOdEnBvJ0tNTxycxf78BARETnHpsAzevRo5OfnIyoqCqNHj673OEmSYDTWP8yaasghpsEmLdbwEBERuYRNgad2MxabtFxD7pfTcOCR+/Aw8BARETmDd1I3kWttGpqHR1e9T+PDTstERETOsKmGZ/ny5Taf8Mknn3S4MEIILFy4ECtXrkRRUREGDBiAFStWoHPnzvW+ZtGiRdi8eTN+/vlnBAQE4Pbbb8eSJUuQmJjocDmagj19eNikRURE5BybAs8bb7xh08kkSXIq8CxduhTLly/HunXrkJCQgOeffx7Dhg3DiRMn4O/vb/U1mZmZmDZtGm699VYYDAbMnz8fQ4cOxYkTJxAUFORwWRpbTeC58SgtNmkRERE5x6bAk5ub29jlgBAC6enpeO655zBq1CgAwPr16xEdHY0tW7bggQcesPq6HTt2WDxfu3YtoqKicOjQIdxxxx2NXm5H+dnTaZmjtIiIiJzi1J1UCAEh6q+hsEdubi7y8/ORlpambAsLC0NycjL27dtn83mKi4sBABEREfUeo9VqUVJSYvFoanINj66BPjw1Ew8y8BARETnDoTvp+++/j+7du8Pf3x/+/v7o3r07Vq1a5VRB5LW5oqOjLbZHR0dbrNvVEJPJhJkzZ2LAgAHo3r17vcctWrQIYWFhyqNt27aOF9xBGptGabEPDxERkSvYfSddsGABnnrqKYwYMQKbNm3Cpk2bMGLECPz1r3/FggULbD7Phg0bEBwcrDz0er29Ralj2rRpOH78ODZu3NjgcfPmzUNxcbHyOHfunNPvbS9lpuUGAo88gkvDpSWIiIicYvPSErIVK1Zg5cqVGDdunLJt5MiR6NmzJ2bMmIGXXnrJpvOMHDkSycnJynOtVgsAKCgoQGxsrLK9oKAASUlJNzzf9OnT8dlnn2Hv3r1o06ZNg8f6+fnBz8/PpnI2FrmZSm/D0hKcaZmIiMg5dgcevV6Pfv361dnet29fGAwGm88TEhKCkJAQ5bkQAjExMdi9e7cScEpKSpCVlYWpU6fWex4hBGbMmIFPPvkEGRkZSEhIsP1i3MimxUM5SouIiMgl7L6TPvTQQ1ixYkWd7e+99x7Gjx/vcEEkScLMmTPx8ssvY+vWrTh27BgmTJiAuLg4i+UsUlNT8fbbbyvPp02bhn//+9/48MMPERISgvz8fOTn56OystLhsjQFXxuatOQh6+zDQ0RE5By7a3gAc6flnTt34ve//z0AICsrC3l5eZgwYQJmzZqlHPf666/bdd5nnnkG5eXlmDJlCoqKijBw4EDs2LHDYg6enJwcXL58WXkuh6+UlBSLc61ZswaTJk2y88qaji2dlnUclk5EROQSdgee48ePo0+fPgDM4QMAIiMjERkZiePHjyvHSZL9HW0lScJLL73UYD+gs2fPWjx31bD4puZ7g2HpQoiapSXYaZmIiMgpdgeePXv2NEY5bjpyM5VJAEaTgFplGWoMppog56dWN2nZiIiIvA3bStykdjOVtWat2tu4eCgREZFz7K7hqaqqwltvvYU9e/agsLAQJpPlzfrw4cMuK5w3qz3ySmc0wV9jWYtTe7g6Oy0TERE5x+7AM3nyZOzcuRP33XcfbrvtNof66pBlvxy9lX48tUdv+aj4OyYiInKG3YHns88+w7Zt2zBgwIDGKM9NQ5IkaNQS9EZhdcX02nPwMFQSERE5x+62ktatW1tMGEiOa2gBUXkbZ1kmIiJynt1302XLlmHu3Ln49ddfG6M8NxUl8DTQaZlD0omIiJxnd5NWv379UFVVhQ4dOiAwMBAajcZi/9WrV11WOG/X0PISOq6UTkRE5DJ2B55x48bh/PnzeOWVVxAdHc3+JU7w86k/8HBZCSIiItexO/B899132LdvH3r16tUY5bmpNLS8BFdKJyIich2776Zdu3b1+IU5m4uaTstWRmkZuFI6ERGRq9h9N128eDFmz56NjIwMXLlyBSUlJRYPsl1DnZa1ysKhbDIkIiJylt1NWsOHDwcApKamWmwXQkCSJBiNRteU7CYgLy9hbeJBvYGdlomIiFyFi4e6kW+DfXjYaZmIiMhV7A48gwcPrnff8ePHnSrMzUbukNzQPDzsw0NEROQ8p++mpaWleO+993Dbbbdx5JadaubhqdtpWcdRWkRERC7j8N107969mDhxImJjY/Haa6/hzjvvxPfff+/Ksnk9W5aW4EzLREREzrOrSSs/Px9r167F+++/j5KSEtx///3QarXYsmULunXr1lhl9Fq+Dcy0rOdMy0RERC5j8910xIgRSExMxNGjR5Geno4LFy7grbfeasyyeT2bJh5k4CEiInKazTU827dvx5NPPompU6eic+fOjVmmm0bDnZY5SouIiMhVbL6bfvPNNygtLUXfvn2RnJyMt99+G5cvX27Msnk9pdOylZmW5T487LRMRETkPJvvpr///e+xcuVKXLx4EY899hg2btyIuLg4mEwm7Nq1C6WlpY1ZTq/U0Grp7MNDRETkOnbfTYOCgvDII4/gm2++wbFjxzB79mwsXrwYUVFRGDlyZGOU0Ws11KSljNLi0hJEREROc6r6IDExEUuXLsVvv/2Gjz76yFVlumnInZatDUtnp2UiIiLXccndVK1WY/To0di6dasrTnfTaKhJS8dOy0RERC7Du6kbyU1aDQ5LZ6dlIiIip/Fu6ka+Dcy0LAceHxX78BARETmLgceN5DCjN9Udls4aHiIiItfh3dSNNHKTltW1tNiHh4iIyFV4N3Ujjcr86zdYqeExmNikRURE5CoMPG4kz7FjrdOygaO0iIiIXIZ3UzfyUTU0LJ0zLRMREbkK76ZuJE88KNfm1GaQR2mp2aRFRETkLAYeN2p4LS1zCOJMy0RERM7j3dSNfJTAU/+wdNbwEBEROY+Bx42UJi2TlU7L1SO35H4+RERE5DjeTd1IY0MNjy9XSyciInIaA48bKTMtNzAsnTU8REREzuPd1I0a7rTMYelERESuwrupG8lhxtqw9JrAwyYtIiIiZzHwuJEcZhps0mINDxERkdN4N3Wjhjot61jDQ0RE5DIMPG7kY8OwdPbhISIich7vpm5Uu4ZHiJpaHpNJwMjAQ0RE5DK8m7qRptaQc7lGBwD0tWp8ONMyERGR8xh43EhTa1LB2iO1av+s4Tw8RERETuPd1I1qTyqoqzVSq/aoLXZaJiIich4DjxvVDjMGi8BTU8OjVjHwEBEROYuBx40kSaq1vEStJi1TzZB0SWLgISIichYDj5v5WJl8UG/gCC0iIiJX4h3VzaytpyWP0vJhcxYREZFLMPC4mbKelqnuKC1fH348RERErsA7qpvV9OGpO0rLh0PSiYiIXIJ3VDeztp6WEng4JJ2IiMglGHjcTB6abm1Yui87LRMREbkE76huJtfw1J540MAaHiIiIpdi4HEzH7nTcu0mLS4cSkRE5FK8o7qZr9ykZao9D49cw8OPh4iIyBV4R3UzOdToDFZmWuY8PERERC7BwONm8rD02jU8OiObtIiIiFyJd1Q3kycX1LPTMhERUaNh4HEzq4uHsoaHiIjIpXhHdTONlVFa8hB1DWt4iIiIXMLH3QW42VlbPLSmSYt5lIjIUxiNRuj1encX46ai0WigVqtdci6PCjxCCCxcuBArV65EUVERBgwYgBUrVqBz5842vX7x4sWYN28ennrqKaSnpzduYV1E7qdjEXhMnGmZiMhTCCGQn5+PoqIidxflphQeHo6YmBhIknOtHh4VeJYuXYrly5dj3bp1SEhIwPPPP49hw4bhxIkT8Pf3b/C1Bw4cwLvvvouePXs2UWldw9paWjpl8VA2aRERuZscdqKiohAYGOj0jZdsI4RARUUFCgsLAQCxsbFOnc9jAo8QAunp6XjuuecwatQoAMD69esRHR2NLVu24IEHHqj3tWVlZRg/fjxWrlyJl19++YbvpdVqodVqleclJSXOX4CDrK2lJffnYZMWEZF7GY1GJey0bNnS3cW56QQEBAAACgsLERUV5VTzlsfcUXNzc5Gfn4+0tDRlW1hYGJKTk7Fv374GXztt2jTcfffdFq9tyKJFixAWFqY82rZt61TZnaHU8Jjqrpbuy07LRERuJffZCQwMdHNJbl7y797Z/lMeE3jy8/MBANHR0Rbbo6OjlX3WbNy4EYcPH8aiRYtsfq958+ahuLhYeZw7d86xQruAj6pup2U9a3iIiDwKm7Hcx1W/e7fdUTds2IDg4GDl4UhyO3fuHJ566ils2LDhhn18avPz80NoaKjFw12sN2nJw9IZeIiIiFzBbX14Ro4cieTkZOW53KemoKDAomNSQUEBkpKSrJ7j0KFDKCwsRJ8+fZRtRqMRe/fuxdtvvw2tVuuy4WyNxVqnZT3n4SEiInIpt1UhhISEoFOnTsqjW7duiImJwe7du5VjSkpKkJWVhf79+1s9R2pqKo4dO4bs7Gzl0a9fP4wfPx7Z2dkeH3YA68PS5f48cnMXERGRqwkhsGDBAsTGxiIgIABpaWk4ffq0u4vVaDzmjipJEmbOnImXX34ZW7duxbFjxzBhwgTExcVh9OjRynGpqal4++23AZhDU/fu3S0eQUFBaNmyJbp37+6mK7GPtZmW9YbqGh4f1vAQEVHjkKeCeeedd5CVlYWgoCAMGzYMVVVV7i5ao/CYwAMAzzzzDGbMmIEpU6bg1ltvRVlZGXbs2GHRPycnJweXL192YyldS9PAxIMa1vAQEZGDysvLMWHCBAQHByM2NhbLli1DSkoKZs6cWWcqmJ49e2L9+vW4cOECtmzZUu85X3jhBSQlJWH16tVo164dgoOD8cQTT8BoNGLp0qWIiYlBVFQU/vGPfyivefrpp/GnP/1JeZ6eng5JkrBjxw5lW6dOnbBq1apG+T3IPGYeHsBcy/PSSy/hpZdeqveYs2fPNniOjIwM1xaqkSmjtKwMS+dq6UREnkcIgUq9scnfN0CjtmvE0pw5c5CZmYlPP/0UUVFRmD9/Pg4fPoykpKQbTgUjz32XkpKC9u3bY+3atcpxOTk52L59O3bs2IGcnBzcd999+OWXX9ClSxdkZmbiu+++wyOPPIK0tDQkJydj8ODBWLVqFYxGI9RqNTIzMxEZGYmMjAwMHz4c58+fR05ODlJSUlz1q7LKowLPzUjjUx14DLWHpXOUFhGRp6rUG9FtwRdN/r4nXhqGQF/bbttlZWV4//338e9//xupqakAgHXr1qFNmzYAbJ8Kpl27dnVmODaZTFi9ejVCQkLQrVs3DBkyBCdPnsS2bdugUqmQmJiIJUuWYM+ePUhOTsagQYNQWlqKI0eOoG/fvti7dy/mzJmj1CRlZGSgdevW6NSpk0O/F1sx8LiZpnr5CIOp7kzLHKVFRESOyMnJgU6nsxgNHRERgcTERLvOs379+jrb2rdvj5CQEOV5dHQ01Go1VLW6YURHRytLQoSHh6NXr17IyMiAr68vfH19MWXKFCxcuBBlZWXIzMzE4MGD7b1EuzHwuFlDa2mxhoeIyPMEaNQ48dIwt7yvq8TExACwbyoYmUajsXguSZLVbaZaf8inpKQgIyMDfn5+GDx4MCIiInDLLbfgm2++QWZmJmbPnu3kFd0YA4+bWV0tnTMtExF5LEmSbG5acpeOHTtCo9EgKysL7dq1AwBcu3YNp06dwuDBg5GQkKBMBSMHHHkqmKlTp7q8PIMHD8bq1avh4+OD4cOHAzCHoI8++ginTp1q9P47gIeN0roZWRuWLjdvabhaOhEROSA4OBiTJ0/GnDlz8NVXX+H48eOYNGmS0uxk61QwEyZMwLx585wuzx133IHS0lJ89tlnSrhJSUnBhg0bEBsbiy5dujj9Hjfi2RH1JiAHHl2tGh4da3iIiMhJr776KsrKyjBixAiEhIRg9uzZKC4uVvY/88wzKC8vx5QpU1BUVISBAwfWmQomLy/Pom+Oo1q0aIEePXqgoKAAXbt2BWAOQSaTqUn67wCAJIQQNz7Mu5WUlCAsLAzFxcVNvq5WxslCTFpzAL+LC8XnTw4CAPzxza9x4mIJ1j1yGwZ3adWk5SEiohpVVVXIzc1FQkKCXWs2eqqUlBQkJSUhPT3d3UWxWUOfgT33b1YhuJlv9bB0naF2DY/5Z1/W8BAREbkE76hu5quM0qo7D48vl5YgIiJyCfbhcTOrNTwGuYbH8xc/JSKi5qO5rUbgSqzhcbOaTst1l5bg4qFERESuwcDjZjU1PDXrsmgN7MNDRETkSryjupmvlZmWa/rw8OMhIiJyBd5R3Uyp4TFa68PDj4eIiMgV2GnZzeRQYzQJGE0CQgiYqit7WMNDRETkGgw8bqapFWr0RhNqTwPJxUOJiIhcg4HHzWo3W2kNJqBW4GENDxERkWvwjupmGnXN0HO90WTRl8eHi4cSEVEj2bx5M4YOHYqWLVtCkiRkZ2e7u0iNioHHzSRJUmp5dIaawOPro4IkMfAQEVHjKC8vx8CBA7FkyRJ3F6VJMPB4ALmWR280Qc8RWkREnk0IQFfe9A871/ouLy/HhAkTEBwcjNjYWCxbtgwpKSmYOXMmAOChhx7CggULkJaWZvM5165di/DwcHz22WdITExEYGAg7rvvPlRUVGDdunVo3749WrRogSeffBJGo3l+ubfffhvdu3dXzrFlyxZIkoR33nlH2ZaWlobnnnvOruuzF/vweABfHxXKdUboDCalCw/77xAReSh9BfBKXNO/7/wLgG+QzYfPmTMHmZmZ+PTTTxEVFYX58+fj8OHDSEpKsvkckyZNwtmzZy2WpKioqMDy5cuxceNGlJaW4p577sGf//xnhIeHY9u2bfjll19w7733YsCAARg7diwGDx6MJ598EpcuXUKrVq2QmZmJyMhIZGRk4PHHH4der8e+ffvwt7/9zY5fhv0YeDxAzfISNaO0avftISIiskdZWRnef/99/Pvf/0ZqaioAYN26dWjTpo1d54mNjYXJZLLYptfrsWLFCnTs2BEAcN999+GDDz5AQUEBgoOD0a1bNwwZMgR79uzB2LFj0b17d0RERCAzMxP33XcfMjIyMHv2bLz55psAgP3790Ov1+P22293wZXXj4HHA9ReQJQ1PEREHk4TaK5tccf72ignJwc6nQ7JycnKtoiICCQmJtr1losWLaqzLTAwUAk7ABAdHY327dsjODjYYlthYSEAc1/VO+64AxkZGUhLS8OJEyfwxBNPYOnSpfj555+RmZmJW2+9FYGBtl+fIxh4PIDVwMM+PEREnkmS7Gpa8jYajcbiuSRJVrfVrhlKSUnBe++9h6+//hq9e/dGaGioEoIyMzMxePDgRi8376oeoPZ6WspK6Qw8RETkoI4dO0Kj0SArK0vZdu3aNZw6dcot5Rk8eDBOnDiBTZs2ISUlBYA5BH355Zf49ttvlW2NiTU8HqBmPS2j0ofHj01aRETkoODgYEyePBlz5sxBy5YtERUVhWeffRYqVc295erVq8jLy8OFC+bmuZMnTwIAYmJiEBMTAwCYN28ezp8/j/Xr1ztVnp49e6JFixb48MMP8dlnnwEwB56nn34akiRhwIABTp3fFryregCl07KBNTxEROQar776KgYNGoQRI0YgLS0NAwcORN++fZX9W7duRe/evXH33XcDAB544AH07t3bYrj4xYsXkZeX53RZJEnCoEGDIEkSBg4cCMAcgkJDQ9GvXz8EBTV+E6EkhJ0D+71QSUkJwsLCUFxcjNDQ0CZ//3HvfY99v1zB8nG9IYTAUxuzcXvHlvjw0d83eVmIiKhGVVUVcnNzkZCQAH9/f3cXx2kpKSlISkpCenq6u4tis4Y+A3vu32zS8gDyAqL6Wp2WWcNDRETkOgw8HsDXyjw8HJZORETkOgw8HsDXp2ZpCSXwsIaHiIhcrPaMyTcbBh4PUHvxUNbwEBERuR4Djwfg0hJERESNi4HHA1jMtMwaHiIiIpdj4PEA1peWULuvQERERF6GgccD1CwtUatJy4dNWkRERK7CwOMBrNXw+HGUFhERkcsw8HiAmk7LAqiOPJx4kIiIyHV4V/UAtWt4tAaTxTYiIqLGIITAggULEBsbi4CAAKSlpeH06dPuLlaj4V3VA2hq9eHRG1nDQ0REjW/p0qVYvnw53nnnHWRlZSEoKAjDhg1DVVWVzefQ6XSNWELX4l3VA9Su4dEZjBbbiIjIswghUKGvaPKHvWt9l5eXY8KECQgODkZsbCyWLVuGlJQUzJw5E0IIpKen47nnnsOoUaPQs2dPrF+/HhcuXMCWLVvqPecLL7yApKQkrFq1ymIxz6KiIvzlL39Bq1atEBoaijvvvBM//PCD8roffvgBQ4YMQUhICEJDQ9G3b18cPHjQod+/o9iHxwP4qmstLaFsY+AhIvJElYZKJH+Y3OTvm/VgFgI1gTYfP2fOHGRmZuLTTz9FVFQU5s+fj8OHDyMpKQm5ubnIz89HWlqacnxYWBiSk5Oxb98+PPDAAwDMq6u3b98ea9euVY47c+YM/ve//2Hz5s1QV0+hMmbMGAQEBGD79u0ICwvDu+++i9TUVJw6dQoREREYP348evfujRUrVkCtViM7OxsajcY1vxgbMfB4AKWGh4uHEhGRC5SVleH999/Hv//9b6SmpgIA1q1bhzZt2gAA8vPzAQDR0dEWr4uOjlb2AUC7du0QGxtrcYxOp8P69evRqlUrAMA333yD/fv3o7CwEH5+fgCA1157DVu2bMHHH3+MKVOmIC8vD3PmzEHXrl0BAJ07d26Eq24YA48HkPvraA0meZAW+/AQEXmoAJ8AZD2Y5Zb3tVVOTg50Oh2Sk2tqoiIiIpCYmGjXe65fv77Otvj4eCXsAObmqrKyMrRs2dLiuMrKSuTk5AAAZs2ahb/85S/44IMPkJaWhjFjxqBjx452lcVZDDwewNrEg6zhISLyTJIk2dW05IliYmIAAAUFBRY1OAUFBUhKSmrwtUFBQRbPy8rKEBsba3Ul9vDwcADmvj8PPvggPv/8c2zfvh0LFy7Exo0b8ec//9mp67AH76oeQGPRadk8LJ2LhxIRkaM6duwIjUaDrKyamqhr167h1KlTAICEhATExMRg9+7dyv6SkhJkZWWhf//+dr1Xnz59kJ+fDx8fH3Tq1MniERkZqRzXpUsX/PWvf8XOnTtxzz33YM2aNU5epX0YeDyAn8WwdM7DQ0REzgkODsbkyZMxZ84cfPXVVzh+/DgmTZoElcp8b5EkCTNnzsTLL7+MrVu34tixY5gwYQLi4uIwevRo5TwTJkzAvHnzGnyvtLQ09O/fH6NHj8bOnTtx9uxZfPfdd3j22Wdx8OBBVFZWYvr06cjIyMCvv/6Kb7/9FgcOHMAtt9zSmL+COtik5QGsLi3BwENERE549dVXUVZWhhEjRiAkJASzZ89GcXGxsv+ZZ55BeXk5pkyZgqKiIgwcOBA7duxQhpoDQF5enhKS6iNJErZt24Znn30WDz/8MC5duoSYmBjccccdiI6OhlqtxpUrVzBhwgQUFBQgMjIS99xzD1588cVGu3ar5RT2Duz3QiUlJQgLC0NxcTFCQ0Ob/P1/OFeEUf/8Fq3DAyCEwIXiKmydPgA924Q3eVmIiKhGVVUVcnNzLeacac5SUlKQlJSE9PR0dxfFZg19Bvbcv1nD4wHkGh5t9aSDtbcRERGR8xh4PECwn/ljKK0yKE1a8jYiIiJyHu+qHiAs0DzbpLxwKACEBTTtDJREROT9rA0dv1mw3cQDBPv6QFVrFLpKYg0PERGRKzHweACVSkKIf02NTmiABpLEeXiIiIhchYHHQ9RuwmJzFhERkWsx8HgIBh4iIqLGw8DjIRh4iIiIGg8Dj4eoHXJCGXiIiIhcioHHQ4QG1IzKCvVn4CEiInIlBh4PEcomLSIiakKbN2/G0KFD0bJlS0iShOzsbHcXqVEx8HgI9uEhIqKmVF5ejoEDB2LJkiUOn8NoNMJkMt34QA/AwOMhGHiIiJoHIQRMFRVN/rB3re/y8nJMmDABwcHBiI2NxbJly5CSkoKZM2cCAB566CEsWLAAaWlpNp9z7dq1CA8Px9atW9GtWzf4+fkhLy8PWq0WTz/9NFq3bo2goCAkJydbzOr866+/YsSIEWjRogWCgoLwu9/9Dtu2bbPrepzF6Xw9BAMPEVHzICorcbJP3yZ/38TDhyAFBtp8/Jw5c5CZmYlPP/0UUVFRmD9/Pg4fPoykpCSbzzFp0iScPXvWIrxUVFRgyZIlWLVqFVq2bImoqChMnz4dJ06cwMaNGxEXF4dPPvkEw4cPx7Fjx9C5c2dMmzYNOp0Oe/fuRVBQEE6cOIHg4GA7rt55HlXDI4TAggULEBsbi4CAAKSlpeH06dM3fN358+fx//7f/0PLli0REBCAHj164ODBg01QYtcJtZhpmTmUiIgcV1ZWhvfffx+vvfYaUlNT0aNHD6xbtw4Gg8Gu88TGxqJdu3YW2/R6Pf71r3/h9ttvR2JiIi5fvow1a9Zg06ZNGDRoEDp27Iinn34aAwcOxJo1awAAeXl5GDBgAHr06IEOHTrgT3/6E+644w6XXa8tPOrOunTpUixfvhzr1q1DQkICnn/+eQwbNgwnTpyAv7+/1ddcu3YNAwYMwJAhQ7B9+3a0atUKp0+fRosWLZq49M5hDQ8RUfMgBQQg8fAht7yvrXJycqDT6ZCcnKxsi4iIQGJiol3vuWjRojrbfH190bNnT+X5sWPHYDQa0aVLF4vjtFotWrZsCQB48sknMXXqVOzcuRNpaWm49957Lc7RFDwm8AghkJ6ejueeew6jRo0CAKxfvx7R0dHYsmULHnjgAauvW7JkCdq2baukSABISEhokjK7EgMPEVHzIEmSXU1L3iYgIMBivceysjKo1WocOnQIarXa4li52eovf/kLhg0bhs8//xw7d+7EokWLsGzZMsyYMaPJyu0xTVq5ubnIz8+36DwVFhaG5ORk7Nu3r97Xbd26Ff369cOYMWMQFRWF3r17Y+XKlQ2+l1arRUlJicXD3Rh4iIjIVTp27AiNRoOsrCxl27Vr13Dq1CmXv1fv3r1hNBpRWFiITp06WTxiYmKU49q2bYvHH38cmzdvxuzZs294r3Y1jwk8+fn5AIDo6GiL7dHR0co+a3755ResWLECnTt3xhdffIGpU6fiySefxLp16+p9zaJFixAWFqY82rZt65qLcEJ4oAapXaOQdksUAw8RETklODgYkydPxpw5c/DVV1/h+PHjmDRpElSqmtv+1atXkZ2djRMnTgAATp48iezsbIt77rx58zBhwoQG36tLly4YP348JkyYgM2bNyM3Nxf79+/HokWL8PnnnwMAZs6ciS+++AK5ubk4fPgw9uzZg1tuuaURrrx+bgs8GzZsQHBwsPLQ6/UOncdkMqFPnz545ZVX0Lt3b0yZMgWPPvoo3nnnnXpfM2/ePBQXFyuPc+fOOXoZLiNJEt6fdCtWTbzVoqqQiIjIEa+++ioGDRqEESNGIC0tDQMHDkTfvjWjy7Zu3YrevXvj7rvvBgA88MAD6N27t8X98+LFi8jLy7vhe61ZswYTJkzA7NmzkZiYiNGjR+PAgQNKh2ej0Yhp06bhlltuwfDhw9GlSxf861//cvEVN0wS9g7sd5HS0lIUFBQoz7VaLbp3744jR45YDJkbPHgwkpKS8Oabb1o9T3x8PP7whz9g1apVyrYVK1bg5Zdfxvnz520qS0lJCcLCwlBcXIzQ0FDHLoiIiLxOVVUVcnNzkZCQUO/gmeYkJSUFSUlJSE9Pd3dRbNbQZ2DP/dttnZZDQkIQEhKiPBdCICYmBrt371YCT0lJCbKysjB16tR6zzNgwACcPHnSYtupU6cQHx/fKOUmIiKi5sdj+vBIkoSZM2fi5ZdfxtatW3Hs2DFMmDABcXFxGD16tHJcamoq3n77beX5X//6V3z//fd45ZVXcObMGXz44Yd47733MG3aNDdcBREREXkijxmWDgDPPPMMysvLMWXKFBQVFWHgwIHYsWOHRRVWTk4OLl++rDy/9dZb8cknn2DevHl46aWXkJCQgPT0dIwfP94dl0BEROSxas+YfLNxWx8eT8I+PEREZI239eFpjlzVh8djmrSIiIg8FesG3MdVv3sGHiIionpoNOZ50SoqKtxckpuX/LuXPwtHeVQfHiIiIk+iVqsRHh6OwsJCAEBgYCDnSmsiQghUVFSgsLAQ4eHhdZatsBcDDxERUQPk5RHk0ENNKzw83GKJCkcx8BARETVAkiTExsYiKirK4VUByDEajcbpmh0ZAw8REZEN1Gq1y26+1PTYaZmIiIi8HgMPEREReT0GHiIiIvJ67MODmkmNSkpK3FwSIiIispV837ZlckIGHgClpaUAgLZt27q5JERERGSv0tJShIWFNXgM19ICYDKZcOHCBYSEhLh8QqmSkhK0bdsW586d88p1urz9+gDvv0Zvvz7A+6/R268P8P5r9PbrAxrnGoUQKC0tRVxcHFSqhnvpsIYHgEqlQps2bRr1PUJDQ732Swx4//UB3n+N3n59gPdfo7dfH+D91+jt1we4/hpvVLMjY6dlIiIi8noMPEREROT1GHgamZ+fHxYuXAg/Pz93F6VRePv1Ad5/jd5+fYD3X6O3Xx/g/dfo7dcHuP8a2WmZiIiIvB5reIiIiMjrMfAQERGR12PgISIiIq/HwENERERe76YLPHq9HnPnzkWPHj0QFBSEuLg4TJgwARcuXGjwdaWlpZg5cybi4+MREBCA22+/HQcOHLA4ZvPmzRg6dChatmwJSZKQnZ3tkjL/+OOPuPfee9G+fXtIkoT09PQ6x+zduxcjRoxAXFwcJEnCli1bXPLeshdeeAGSJFk8unbtquyvqqrCtGnT0LJlSwQHB+Pee+9FQUGBS8vgSjf6fQkhsGDBAsTGxiIgIABpaWk4ffq0xTFXr17F+PHjERoaivDwcEyePBllZWVNeBUNu9E1Tpo0qc5nOnz4cItjPPkaFy1ahFtvvRUhISGIiorC6NGjcfLkSYtjbPle5uXl4e6770ZgYCCioqIwZ84cGAyGprwUq2y5vpSUlDqf4eOPP25xjKdeHwCsWLECPXv2VCai69+/P7Zv367sb86fH3Dj62vun9/1Fi9eDEmSMHPmTGWbJ32GN13gqaiowOHDh/H888/j8OHD2Lx5M06ePImRI0c2+Lq//OUv2LVrFz744AMcO3YMQ4cORVpaGs6fP68cU15ejoEDB2LJkiUuL3OHDh2wePFixMTEWD2mvLwcvXr1wj//+U+Xvndtv/vd73Dx4kXl8c033yj7/vrXv+L//u//sGnTJmRmZuLChQu45557Gq0szrrR72vp0qVYvnw53nnnHWRlZSEoKAjDhg1DVVWVcsz48ePx448/YteuXfjss8+wd+9eTJkypaku4YZs+U4MHz7c4jP96KOPLPZ78jVmZmZi2rRp+P7777Fr1y7o9XoMHToU5eXlyjE3+l4ajUbcfffd0Ol0+O6777Bu3TqsXbsWCxYscMclWbDl+gDg0UcftfgMly5dquzz5OsDgDZt2mDx4sU4dOgQDh48iDvvvBOjRo3Cjz/+CKB5f37Aja8PaN6fX20HDhzAu+++i549e1ps96jPUJDYv3+/ACB+/fVXq/srKiqEWq0Wn332mcX2Pn36iGeffbbO8bm5uQKAOHLkSJ19165dE5MnTxaRkZEiJCREDBkyRGRnZ9tc1vj4ePHGG280eAwA8cknn9h8TlssXLhQ9OrVy+q+oqIiodFoxKZNm5RtP/30kwAg9u3b59JyNIbrf18mk0nExMSIV199VdlWVFQk/Pz8xEcffSSEEOLEiRMCgDhw4IByzPbt24UkSeL8+fNNVnZbWftOTJw4UYwaNare1zS3aywsLBQARGZmphDCtu/ltm3bhEqlEvn5+coxK1asEKGhoUKr1TbtBdzA9dcnhBCDBw8WTz31VL2vaU7XJ2vRooVYtWqV131+Mvn6hPCez6+0tFR07txZ7Nq1y+KaPO0zvOlqeKwpLi6GJEkIDw+3ut9gMMBoNMLf399ie0BAgEUthy3GjBmDwsJCbN++HYcOHUKfPn2QmpqKq1evOlr8JnP69GnExcWhQ4cOGD9+PPLy8gAAhw4dgl6vR1pamnJs165d0a5dO+zbt89dxXVYbm4u8vPzLa4nLCwMycnJyvXs27cP4eHh6Nevn3JMWloaVCoVsrKymrzMjsrIyEBUVBQSExMxdepUXLlyRdnX3K6xuLgYABAREQHAtu/lvn370KNHD0RHRyvHDBs2DCUlJRZ/hXuC669PtmHDBkRGRqJ79+6YN28eKioqlH3N6fqMRiM2btyI8vJy9O/f3+s+v+uvT+YNn9+0adNw9913W3xWgOf9G7zpFw+tqqrC3LlzMW7cuHoXMwsJCUH//v3x97//Hbfccguio6Px0UcfYd++fejUqZPN7/XNN99g//79KCwsVGaafO2117BlyxZ8/PHHHtNUYE1ycjLWrl2LxMREXLx4ES+++CIGDRqE48ePIz8/H76+vnUCY3R0NPLz891TYCfIZa79D1B+Lu/Lz89HVFSUxX4fHx9EREQ0m2sePnw47rnnHiQkJCAnJwfz58/HXXfdhX379kGtVjerazSZTJg5cyYGDBiA7t27A4BN38v8/Hyrn7O8z1NYuz4AePDBBxEfH4+4uDgcPXoUc+fOxcmTJ7F582YAzeP6jh07hv79+6OqqgrBwcH45JNP0K1bN2RnZ3vF51ff9QHe8flt3LgRhw8frtOnFfC8f4NeH3g2bNiAxx57THm+fft2DBo0CIC5A/P9998PIQRWrFjR4Hk++OADPPLII2jdujXUajX69OmDcePG4dChQzaX5YcffkBZWRlatmxpsb2yshI5OTnIy8tT/iEAwPz58zF//nybz9+Y7rrrLuXnnj17Ijk5GfHx8fjvf/+LgIAAN5aMHPXAAw8oP/fo0QM9e/ZEx44dkZGRgdTUVDeWzH7Tpk3D8ePH7a5xbS7qu77afyT16NEDsbGxSE1NRU5ODjp27NjUxXRIYmIisrOzUVxcjI8//hgTJ05EZmamu4vlMvVdX7du3Zr953fu3Dk89dRT2LVrV50WEE/k9YFn5MiRSE5OVp63bt0aQE3Y+fXXX/HVV1/dcKn6jh07IjMzE+Xl5SgpKUFsbCzGjh2LDh062FyWsrIyxMbGIiMjo86+8PBwhIeHW4zsur7q2pOEh4ejS5cuOHPmDP7whz9Ap9OhqKjIIskXFBTU28nak8llLigoQGxsrLK9oKAASUlJyjGFhYUWrzMYDLh69WqzvGYA6NChAyIjI3HmzBmkpqY2m2ucPn260qG6TZs2yvaYmJgbfi9jYmKwf/9+i/PJI0g85Rrruz5r5P/XnTlzBh07dmwW1+fr66vUlPft2xcHDhzAm2++ibFjx3rF51ff9b377rt1jm1un9+hQ4dQWFiIPn36KNuMRiP27t2Lt99+G1988YVHfYZe34cnJCQEnTp1Uh4BAQFK2Dl9+jS+/PLLOjUuDQkKCkJsbCyuXbuGL774AqNGjbL5tX369EF+fj58fHwsytSpUydERkbW2e7JgaesrAw5OTmIjY1F3759odFosHv3bmX/yZMnkZeXZ9FW3VwkJCQgJibG4npKSkqQlZWlXE///v1RVFRkUcP31VdfwWQyWQTs5uS3337DlStXlJDn6dcohMD06dPxySef4KuvvkJCQoLFflu+l/3798exY8csgt2uXbsQGhpqUdvqDje6PmvkP5hqf4aeen31MZlM0Gq1zf7zq498fdY0t88vNTUVx44dQ3Z2tvLo168fxo8fr/zsUZ+hS7tANwM6nU6MHDlStGnTRmRnZ4uLFy8qj9o9wu+8807x1ltvKc937Nghtm/fLn755Rexc+dO0atXL5GcnCx0Op1yzJUrV8SRI0fE559/LgCIjRs3iiNHjoiLFy8KIcyjfwYOHCh69eolvvjiC5Gbmyu+/fZbMX/+fIuRMNfTarXiyJEj4siRIyI2NlY8/fTT4siRI+L06dPKMaWlpcoxAMTrr78ujhw5Uu/IM3vNnj1bZGRkKGVOS0sTkZGRorCwUAghxOOPPy7atWsnvvrqK3Hw4EHRv39/0b9/f5e8d2O40e9r8eLFIjw8XHz66afi6NGjYtSoUSIhIUFUVlYq5xg+fLjo3bu3yMrKEt98843o3LmzGDdunLsuqY6GrrG0tFQ8/fTTYt++fSI3N1d8+eWXok+fPqJz586iqqpKOYcnX+PUqVNFWFiYyMjIsPh3XFFRoRxzo++lwWAQ3bt3F0OHDhXZ2dlix44dolWrVmLevHnuuCQLN7q+M2fOiJdeekkcPHhQ5Obmik8//VR06NBB3HHHHco5PPn6hBDib3/7m8jMzBS5ubni6NGj4m9/+5uQJEns3LlTCNG8Pz8hGr4+b/j8rLl+5JknfYY3XeCRh4xbe+zZs0c5Lj4+XixcuFB5/p///Ed06NBB+Pr6ipiYGDFt2jRRVFRkce41a9ZYPW/t85SUlIgZM2aIuLg4odFoRNu2bcX48eNFXl6e3WUePHiwcsyePXusHjNx4kQnf2NmY8eOFbGxscLX11e0bt1ajB07Vpw5c0bZX1lZKZ544gnRokULERgYKP785z8rQc8T3ej3ZTKZxPPPPy+io6OFn5+fSE1NFSdPnrQ4x5UrV8S4ceNEcHCwCA0NFQ8//LAoLS11w9VY19A1VlRUiKFDh4pWrVoJjUYj4uPjxaOPPmoxNFQIz77G+v4dr1mzRjnGlu/l2bNnxV133SUCAgJEZGSkmD17ttDr9U18NXXd6Pry8vLEHXfcISIiIoSfn5/o1KmTmDNnjiguLrY4j6denxBCPPLIIyI+Pl74+vqKVq1aidTUVCXsCNG8Pz8hGr4+b/j8rLk+8HjSZygJIYRr64yIiIiIPIvX9+EhIiIiYuAhIiIir8fAQ0RERF6PgYeIiIi8HgMPEREReT0GHiIiIvJ6DDxERETk9Rh4iIiICGfPnsXkyZORkJCAgIAAdOzYEQsXLoROp7M47ujRoxg0aBD8/f3Rtm1bLF26tM65Nm3ahK5du8Lf3x89evTAtm3bmuoy6sXAQ0QeYdKkSRg9enSTv+/atWshSRIkScLMmTNddt727dsr5y0qKnLZeYmckZKSgrVr11rd9/PPP8NkMuHdd9/Fjz/+iDfeeAPvvPMO5s+frxxTUlKCoUOHIj4+HocOHcKrr76KF154Ae+9955yzHfffYdx48Zh8uTJOHLkCEaPHo3Ro0fj+PHjjX15DXP53M1ERNdBPcskyI+FCxeKoqIice3atSYv25o1a0RoaKi4ePGiKCkpcdl5CwsLxf/+9z8BwC3XRWTN4MGDLZZfuZGlS5eKhIQE5fm//vUv0aJFC4u1J+fOnSsSExOV5/fff7+4++67Lc6TnJwsHnvssXrfZ+LEiWLUqFEW25566imLJZQ2bdokunfvLvz9/UVERIRITU0VZWVlNl8La3iIqNFdvHhReaSnpyM0NNRi29NPP42wsDCEh4e7pXySJCEmJgYhISEuO2erVq0QERHhsvMRuUNxcbHF93jfvn2444474Ovrq2wbNmwYTp48iWvXrinHpKWlWZxn2LBh2Ldvn8PluHjxIsaNG4dHHnkEP/30EzIyMnDPPfdA2LE6FgMPETW6mJgY5REWFqYEDPkRHBxcp0krJSUFM2bMwMyZM9GiRQtER0dj5cqVKC8vx8MPP4yQkBB06tQJ27dvt3iv48eP46677kJwcDCio6Px0EMP4fLly3aXuX379njllVfwyCOPICQkBO3atbOottfpdJg+fTpiY2Ph7++P+Ph4LFq0yOHfEZGnOXPmDN566y089thjyrb8/HxER0dbHCc/z8/Pb/AYeb8jLl68CIPBgHvuuQft27dHjx498MQTTyA4ONjmczDwEJHHWrduHSIjI7F//37MmDEDU6dOxZgxY3D77bfj8OHDGDp0KB566CFUVFQAAIqKinDnnXeid+/eOHjwIHbs2IGCggLcf//9Dr3/smXL0K9fPxw5cgRPPPEEpk6dipMnTwIAli9fjq1bt+K///0vTp48iQ0bNqB9+/auunQil3jllVcQHBysPL7++ms8/vjjFtvy8vLqvO78+fMYPnw4xowZg0cffdQNJbfUq1cvpKamokePHhgzZgxWrlyp1CjZioGHiDxWr1698Nxzz6Fz586YN28e/P39ERkZiUcffRSdO3fGggULcOXKFRw9ehQA8Pbbb6N379545ZVX0LVrV/Tu3RurV6/Gnj17cOrUKbvf/49//COeeOIJdOrUCXPnzkVkZCT27NkDAMjLy0Pnzp0xcOBAxMfHY+DAgRg3bpxLr5/IWY8//jiys7OVR79+/fDSSy9ZbIuLi7N4zYULFzBkyBDcfvvtFrWagLm2tqCgwGKb/DwmJqbBY+T9tjIajcrParUau3btwvbt29GtWze89dZbSExMRG5urs3nY+AhIo/Vs2dP5We1Wo2WLVuiR48eyja52rywsBAA8MMPP2DPnj0Wf7127doVAJCTk+PU+8vNcPJ7TZo0CdnZ2UhMTMSTTz6JnTt32n+BRI0sIiICnTp1Uh4BAQGIioqy2Obj46Mcf/78eaSkpKBv375Ys2YNVCrLmNC/f3/s3bsXer1e2bZr1y4kJiaiRYsWyjG7d++2eN2uXbvQv3//Bst6fUj65ZdfLJ5LkoQBAwbgxRdfxJEjR+Dr64tPPvnE5t+Fz40PISJyD41GY/FckiSLbZIkAQBMJhMAoKysDCNGjMCSJUvqnCs2NtYl7y+/V58+fZCbm4vt27fjyy+/xP3334+0tDR8/PHHdr8PkSeQw058fDxee+01XLp0Sdkn1848+OCDePHFFzF58mTMnTsXx48fx5tvvok33nhDOfapp57C4MGDsWzZMtx9993YuHEjDh48WKe26HpZWVlYuXIlUlNT8dVXX+GLL75Ax44dkZubi8LCQuzevRtDhw5FVFQUsrKycOnSJdxyyy02Xx8DDxF5jT59+uB///sf2rdvb/FXa2MJDQ3F2LFjMXbsWNx3330YPnw4rl69ytFZ1Czt2rULZ86cwZkzZ9CmTRuLffJoqLCwMOzcuRPTpk1D3759ERkZiQULFmDKlCnKsbfffjs+/PBDPPfcc5g/fz46d+6MLVu2oHv37g2+/5AhQ/C///0P06dPR1JSEtasWYNp06bh1VdfxYwZM7B3716kp6ejpKQE8fHxWLZsGe666y6br4+Bh4i8xrRp07By5UqMGzcOzzzzDCIiInDmzBls3LgRq1atglqtdtl7vf7664iNjUXv3r2hUqmwadMmxMTEuG1oPZEtMjIy6t03adIkTJo06Ybn6NmzJ77++usGjxkzZgzGjBljV9lCQkKwZcsWi20PPfSQ8vOOHTvsOt/12IeHiLxGXFwcvv32WxiNRgwdOhQ9evTAzJkzER4eXqcvgrNCQkKwdOlS9OvXD7feeivOnj2Lbdu2ufx9iMg1JGHPrD1ERF5m7dq1mDlzZqMs/5CRkYEhQ4bg2rVrrPkhasCkSZNQVFRUp4bHldikRUQ3veLiYgQHB2PatGlWOzw74ne/+12dUSZEZF1963u5Emt4iOimVlpaqgyHDQ8PR2RkpEvO++uvvypDdzt06MCmLiI3Y+AhIiIir8c/OYiIiMjrMfAQERGR12PgISIiIq/HwENERERej4GHiIiIvB4DDxEREXk9Bh4iIiLyegw8RERE5PX+PxyyfmNV1mS4AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f, ax = compiled_sched.plot_pulse_diagram(plot_kwargs={\"x_range\": (200e-6, 200.4e-6)})" ] } ], "metadata": { "file_format": "mystnb", "kernelspec": { "display_name": "python3", "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" }, "mystnb": { "execution_timeout": 120 }, "source_map": [ 8, 38, 47, 53, 60, 73, 81, 106, 129, 143, 153, 162, 167, 171, 175, 188, 197, 202, 205, 212, 215, 220, 225, 229, 235, 240, 244, 248, 251, 269, 283, 287, 295, 300, 303, 314, 345, 352, 364, 380, 385, 392 ] }, "nbformat": 4, "nbformat_minor": 5 }