diff --git a/Team 12 Notebooks/Challenge 1 Team 12.ipynb b/Team 12 Notebooks/Challenge 1 Team 12.ipynb new file mode 100644 index 0000000..27fb464 --- /dev/null +++ b/Team 12 Notebooks/Challenge 1 Team 12.ipynb @@ -0,0 +1,2635 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# The Constrained Variational Quantum Eigensolver algorithm" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit.opflow import Z, I, X, Y\n", + "from qiskit.opflow import CircuitStateFn, StateFn, CircuitSampler, PauliExpectation, ListOp\n", + "from qiskit.quantum_info.operators import Operator, Pauli\n", + "from qiskit.providers.aer import AerSimulator\n", + "from qiskit.utils import QuantumInstance\n", + "from qiskit.circuit import QuantumCircuit, QuantumRegister, Parameter, ParameterVector\n", + "from qiskit.algorithms.optimizers import ADAM, SPSA\n", + "from qiskit.circuit.library import TwoLocal\n", + "from qiskit.circuit.library.n_local import EfficientSU2, RealAmplitudes\n", + "from qiskit import Aer, execute, transpile\n", + "\n", + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "from tqdm.notebook import tqdm" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "We take the Transverse Ising Hamiltonian with periodic boundary conditions :\n", + "\\begin{equation}\n", + "H=J_z\\sum_{i=0}^{1}\\sigma_i^z\\sigma_{i+1}^z + h\\sum_{i=0}^{2}\\sigma_i^x,\n", + "\\end{equation}\n" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1.0 * ZZI\n", + "+ 1.0 * IZZ\n", + "+ 1.0 * XII\n", + "+ 1.0 * IXI\n", + "+ 1.0 * IIX\n" + ] + } + ], + "source": [ + "# implement the Transverse Ising Hamiltonian H with periodic bounday conditions and parameters N, J and h.\n", + "\n", + "N = 3\n", + "h = 1\n", + "J = 1\n", + "\n", + "\n", + "H = J* ( (Z^Z^I) + (I^Z^Z) ) + h *(( X^I^I) + ( I^X^I) + ( I^I^X))\n", + "\n", + "\n", + "\n", + "H = H.reduce() \n" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
      ┌──────────┐┌──────────┐ ░                 ░ ┌──────────┐ ┌──────────┐\n",
+       "q1_0: ┤ Ry(θ[0]) ├┤ Rz(θ[3]) ├─░───■────■────────░─┤ Ry(θ[6]) ├─┤ Rz(θ[9]) ├\n",
+       "      ├──────────┤├──────────┤ ░ ┌─┴─┐  │        ░ ├──────────┤┌┴──────────┤\n",
+       "q1_1: ┤ Ry(θ[1]) ├┤ Rz(θ[4]) ├─░─┤ X ├──┼────■───░─┤ Ry(θ[7]) ├┤ Rz(θ[10]) ├\n",
+       "      ├──────────┤├──────────┤ ░ └───┘┌─┴─┐┌─┴─┐ ░ ├──────────┤├───────────┤\n",
+       "q1_2: ┤ Ry(θ[2]) ├┤ Rz(θ[5]) ├─░──────┤ X ├┤ X ├─░─┤ Ry(θ[8]) ├┤ Rz(θ[11]) ├\n",
+       "      └──────────┘└──────────┘ ░      └───┘└───┘ ░ └──────────┘└───────────┘
" + ], + "text/plain": [ + " ┌──────────┐┌──────────┐ ░ ░ ┌──────────┐ ┌──────────┐\n", + "q1_0: ┤ Ry(θ[0]) ├┤ Rz(θ[3]) ├─░───■────■────────░─┤ Ry(θ[6]) ├─┤ Rz(θ[9]) ├\n", + " ├──────────┤├──────────┤ ░ ┌─┴─┐ │ ░ ├──────────┤┌┴──────────┤\n", + "q1_1: ┤ Ry(θ[1]) ├┤ Rz(θ[4]) ├─░─┤ X ├──┼────■───░─┤ Ry(θ[7]) ├┤ Rz(θ[10]) ├\n", + " ├──────────┤├──────────┤ ░ └───┘┌─┴─┐┌─┴─┐ ░ ├──────────┤├───────────┤\n", + "q1_2: ┤ Ry(θ[2]) ├┤ Rz(θ[5]) ├─░──────┤ X ├┤ X ├─░─┤ Ry(θ[8]) ├┤ Rz(θ[11]) ├\n", + " └──────────┘└──────────┘ ░ └───┘└───┘ ░ └──────────┘└───────────┘" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "\n", + "\n", + "n_qubits = 3\n", + "params = ParameterVector('θ',(n_qubits*6))\n", + "\n", + "qr = QuantumRegister(n_qubits)\n", + "ansatz = QuantumCircuit(qr)\n", + "\n", + "\n", + " \n", + "ansatz.compose(EfficientSU2(3, reps=1, insert_barriers=True).decompose(), inplace=True)\n", + "ansatz.draw()\n", + "\n", + "\n", + "# use SU2 Ansatz\n", + " " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### The VQE algorithm" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "C:\\Users\\vince\\.conda\\envs\\q_env\\lib\\site-packages\\sympy\\core\\expr.py:3949: SymPyDeprecationWarning: \n", + "\n", + "expr_free_symbols method has been deprecated since SymPy 1.9. See\n", + "https://github.com/sympy/sympy/issues/21494 for more info.\n", + "\n", + " SymPyDeprecationWarning(feature=\"expr_free_symbols method\",\n" + ] + } + ], + "source": [ + "### Implementation with Qiskit modules (written by Julien to demonstrate the advantage of Qiskit builtin functions)\n", + "\n", + "\n", + "from qiskit.algorithms import VQE as QiskitVQE\n", + "from qiskit.algorithms.optimizers import SPSA\n", + "\n", + "vqe = QiskitVQE(ansatz, SPSA(), quantum_instance= Aer.get_backend('qasm_simulator'))\n", + "result = vqe.compute_minimum_eigenvalue(H)" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [], + "source": [ + "opt_circuit = ansatz.bind_parameters(result.optimal_parameters)\n", + "sampler = CircuitSampler(Aer.get_backend('qasm_simulator'))\n", + "\n", + "def overlap(theta):\n", + " expr = StateFn(opt_circuit).adjoint() @ StateFn(ansatz.bind_parameters(theta))\n", + " exp = PauliExpectation().convert(expr)\n", + " return np.abs(sampler.convert(exp).eval()) ** 2" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": {}, + "outputs": [], + "source": [ + "energy = vqe.get_energy_evaluation(H)\n", + "\n", + "def objective(theta):\n", + " return energy(theta) + 2 * overlap(theta)\n", + "\n", + "res = SPSA().optimize(1, objective, initial_point = np.pi * np.random.normal(0,1,size = len(ansatz.parameters)))\n" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(array([ 0.11694058, -3.0815028 , 4.16561148, -2.3187025 , 1.19466155,\n", + " 0.95748962, -2.54543286, 2.66601361, 3.45751739, -5.95875447,\n", + " 0.11767189, 2.03868092]), -2.888671875, 200)\n" + ] + } + ], + "source": [ + "print(res)\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [], + "source": [ + "## code we used during the Hackathon \n", + "\n", + "\n", + "noise_model = None # the noise model used to simulate the hardware, for now, leave it to None (no noise)\n", + "\n", + "\n", + "class VQE:\n", + " def __init__(self,ansatz, qr, H):\n", + " self.H = H\n", + " self.qr = qr\n", + " self.backend = Aer.get_backend('qasm_simulator',noise = noise_model)\n", + " self.ansatz = transpile(ansatz, self.backend)\n", + " self.shots = 2**12 # number of sample used to approximate the expectation value\n", + " #choose some classical optimizer. You may take a gradient based or gradient free optimizer \n", + " # and play with the hyperparameters as well.\n", + " self.optimizer = ADAM(maxiter=1, tol=1e-06, lr=0.02, beta_1=0.9, beta_2=0.99, amsgrad=False)\n", + " \n", + " self.lesser_eigenstate_parameters = [] # list to contain previously computed parameters\n", + " self.constraint_multiplier = 10 \n", + " \n", + "\n", + " def overlap(self, parameters, other_state_parameters):\n", + " \"\"\" \n", + " return the overlap between a state given by specific parameters\n", + " \"\"\"\n", + " \n", + " # to find the overlap < phi_0 | phi_1 > we need to construct a circuit that measures\n", + " # < \"000\" | U_0 * U_1 | \"000\" > where U are the unitary evolution operators\n", + " \n", + " total_overlap = 0\n", + " \n", + " \n", + " for i in range(len(self.H)): #loop over the pauli terms\n", + " \n", + " \n", + " # construct U_0 * \n", + " \n", + " second_circuit = self.ansatz.copy()\n", + " \n", + " second_circuit = second_circuit.bind_parameters(other_state_parameters) \n", + "\n", + " second_circuit = second_circuit.inverse()\n", + "\n", + " # compose U_0 * with U_1\n", + " \n", + " \n", + " additional_circuit = self.ansatz.copy()\n", + " additional_circuit = additional_circuit.bind_parameters(parameters)\n", + "\n", + " second_circuit.compose( self.ansatz.copy() ) \n", + "\n", + " \n", + " coeff = self.H[i].primitive.coeffs # get the coefficients\n", + " pauli = self.H[i].primitive.table.to_labels()[0] #get the pauli string\n", + " \n", + " observable, measure_which = self.basis_change(pauli) #get which qubits to measure and the basis change\n", + " \n", + " second_circuit.append(observable,self.qr) # append the basis change\n", + " second_circuit.measure_all() # measure all the qubits\n", + " \n", + " job = execute(second_circuit, self.backend, shots=self.shots)\n", + " \n", + " counts = job.result().get_counts()\n", + " \n", + " # count \"000\" ( aka project onto the |\"000\"> state)\n", + " # that probability gives the total overlap\n", + " \n", + " if \"000\" in counts: \n", + " total_overlap += coeff * counts[\"000\"] / np.sum([counts[key] for key in counts])\n", + " \n", + " # return the total overlap \n", + " return total_overlap.real\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " def expectation_value(self,parameters):\n", + " \n", + " '''\n", + " return the expectation value of the quantum circuit wrt to the observable H\n", + " \n", + " '''\n", + " \n", + " result = 0\n", + " \n", + " for i in range(len(self.H)): #loop over the pauli terms\n", + " ansatz = self.ansatz.bind_parameters(parameters) #insert the parameters\n", + " coeff = self.H[i].primitive.coeffs # get the coefficient\n", + " pauli = self.H[i].primitive.table.to_labels()[0] #get the pauli string\n", + " \n", + " observable, measure_which = self.basis_change(pauli) #get which qubits to measure and the basis change\n", + " \n", + " ansatz.append(observable,self.qr) # append the basis change\n", + " ansatz.measure_all() # measure all the qubits\n", + " \n", + " job = execute(ansatz, self.backend, shots=self.shots)\n", + " \n", + " counts = job.result().get_counts()\n", + " \n", + " result = result + coeff * self.expectation_from_counts(counts,measure_which) \n", + " \n", + " for i, other_state_parameters in enumerate(self.lesser_eigenstate_parameters) :\n", + " result += self.constraint_multiplier * self.overlap(parameters, other_state_parameters)\n", + " \n", + " \n", + " \n", + " \n", + " return result.real\n", + " \n", + " def basis_change(self,pauli):\n", + " '''transform the measurement basis as function of the observable\n", + " return: -quantum circuit to append to the ansatz\n", + " -qubit to measure\n", + " '''\n", + " observable = QuantumCircuit(len(pauli)) #the basis transformation circuit U to append to the ansatz\n", + " \n", + " measure_which = [] #which qubit we will have to measure\n", + " for i in range(len(pauli)):\n", + " if pauli[i] == \"I\":\n", + " continue \n", + " else:\n", + " measure_which.append(i)\n", + " if pauli[i] == \"X\":\n", + " observable.h(i)\n", + " \n", + " \n", + " return observable, measure_which\n", + " \n", + " def expectation_from_counts(self,counts,measure_which):\n", + " '''compute the expectation value from counts\n", + " Tipp: look at the parity of the state\n", + " remember: in qiskit we read the qubits from right to left\n", + "\n", + " '''\n", + " \n", + " total_counts = np.sum([counts[key] for key in counts])\n", + " expectation = 0\n", + " \n", + " if len(measure_which)==0:\n", + " return 1\n", + " # return the expectation value when the pauli is the identity\n", + " \n", + " for i,state in enumerate(counts):\n", + " \n", + " parity = (-1)**np.sum([int(state[k]) for k in range(len(state)) if (len(state)- 1 - k) in measure_which ])\n", + " # your code goes here, compute the parity of the state\n", + " expectation = expectation + parity * counts[state]\n", + " \n", + " return expectation / total_counts\n", + " \n", + " def gradient(self,parameters):\n", + " '''\n", + " return the gradient of the quantum circuit \n", + " \n", + " '''\n", + " \n", + " gradients = np.zeros_like(parameters)\n", + " for i,p in enumerate(parameters):\n", + " \n", + " shift = np.zeros_like(parameters)\n", + " shift[i] = np.pi/2\n", + " gradients[i] = (1/2) * (self.expectation_value( parameters + shift ) - (self.expectation_value(parameters-shift)))\n", + " \n", + " return gradients\n", + " \n", + " def update(self,parameters):\n", + " '''\n", + " update the parameters with the classical optimizer\n", + " \n", + " '''\n", + " parameters, loss, it = self.optimizer.optimize(parameters.size,\n", + " lambda param: self.expectation_value(param),\n", + " gradient_function= lambda param: self.gradient(param),initial_point=parameters)\n", + " \n", + " return loss, parameters\n", + " \n", + "\n", + " " + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [], + "source": [ + "%matplotlib notebook\n", + "\n", + "\n", + "# The following functions compute the excited states. \n", + "# Sadly they are very long computations. \n", + "# To converge succesfully, about 130 iterations are needed which takes over an hour\n", + "\n", + "# Below you can see plots which show the results after 50 iterations\n", + "\n", + "\n", + "def constrained_Hamiltonian(H, order, a):\n", + " \"\"\"\n", + " Function that uses the constrained VQE to compute different orders of excited states\n", + " \n", + " Inputs :\n", + " H : base Hamiltonian\n", + " order : which eigenvalue we're looking for\n", + " a : multiplier coefficient\n", + " \"\"\"\n", + " plt.ion()\n", + " fig = plt.figure()\n", + " fig.show()\n", + " fig.canvas.draw()\n", + " \n", + " vqe = VQE(ansatz, qr, H)\n", + " \n", + " vqe.constraint_multiplier = a\n", + " \n", + " losses = []\n", + " \n", + " eigen = np.linalg.eig(H.to_matrix())\n", + " exact_eigenvalues = np.sort(np.real(eigen[0]))\n", + " \n", + " \n", + " \n", + " for i in range(order):\n", + " \n", + " \n", + " print(f\"Constructing Hamiltonian and computing ground state of order {i}.\")\n", + " print(f\"Exact eigenvalue computed : {exact_eigenvalues[i]} eV\")\n", + " \n", + " \n", + " ax = fig.add_subplot(111)\n", + " \n", + " vec , l = iterate_to_eigenvector(H, vqe, ax, fig, exact_eigenvalues[i])\n", + " losses.append(l)\n", + " vqe.lesser_eigenstate_parameters.append(vec)\n", + " \n", + " print(f\"The quantum algorithm converged to {l[-1]}\")\n", + " \n", + " return vqe, losses, exact_eigenvalues\n", + " \n", + "\n", + "def iterate_to_eigenvector(Hm, vqe, ax, fig , exact_energy, maxiter = 100):\n", + " \"\"\"\n", + " Iterates to ground state of a Hamiltonian, returns the eigenvector\n", + " \n", + " Inputs:\n", + " H : Hamiltonian\n", + " ax : axis on which to draw\n", + " fig : figure on which to draw\n", + " maxiter : maximum number of iterations\n", + " \n", + " \"\"\"\n", + " \n", + " # select a starting vector according to a random initial vector \n", + " \n", + " v0 = np.pi * np.random.normal(0,1,size = len(ansatz.parameters))\n", + " \n", + " vqe.optimizer.set_options(kwargs = { \"lr\" , 0.08}) # set learning rate to 0.08\n", + " \n", + " \n", + " \n", + " print(\"Computing Eigenvector for new Hamiltonian\")\n", + " \n", + " loss = []\n", + " for i in tqdm(range(maxiter)):\n", + " l , v0 = vqe.update(v0)\n", + " loss.append(l)\n", + " if i == 40 :\n", + " vqe.optimizer.set_options(kwargs = { \"lr\" , 0.02})\n", + " if i == 60:\n", + " vqe.optimizer.set_options(kwargs = { \"lr\" , 0.001})\n", + " \n", + " ax.clear()\n", + " ax.plot(loss,'b.',label='VQE')\n", + " \n", + " ax.plot(exact_energy*np.ones_like(range(maxiter+1)),'k--',label='exact')\n", + " ax.legend()\n", + " ax.set_xlim([-0.2, maxiter])\n", + " ax.set_ylim([-4, max(loss) + 1])\n", + " fig.canvas.draw()\n", + " \n", + " return v0 , loss" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": { + "scrolled": false + }, + "outputs": [ + { + "data": { + "application/javascript": [ + "/* Put everything inside the global mpl namespace */\n", + "/* global mpl */\n", + "window.mpl = {};\n", + "\n", + "mpl.get_websocket_type = function () {\n", + " if (typeof WebSocket !== 'undefined') {\n", + " return WebSocket;\n", + " } else if (typeof MozWebSocket !== 'undefined') {\n", + " return MozWebSocket;\n", + " } else {\n", + " alert(\n", + " 'Your browser does not have WebSocket support. ' +\n", + " 'Please try Chrome, Safari or Firefox ≥ 6. ' +\n", + " 'Firefox 4 and 5 are also supported but you ' +\n", + " 'have to enable WebSockets in about:config.'\n", + " );\n", + " }\n", + "};\n", + "\n", + "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n", + " this.id = figure_id;\n", + "\n", + " this.ws = websocket;\n", + "\n", + " this.supports_binary = this.ws.binaryType !== undefined;\n", + "\n", + " if (!this.supports_binary) {\n", + " var warnings = document.getElementById('mpl-warnings');\n", + " if (warnings) {\n", + " warnings.style.display = 'block';\n", + " warnings.textContent =\n", + " 'This browser does not support binary websocket messages. ' +\n", + " 'Performance may be slow.';\n", + " }\n", + " }\n", + "\n", + " this.imageObj = new Image();\n", + "\n", + " this.context = undefined;\n", + " this.message = undefined;\n", + " this.canvas = undefined;\n", + " this.rubberband_canvas = undefined;\n", + " this.rubberband_context = undefined;\n", + " this.format_dropdown = undefined;\n", + "\n", + " this.image_mode = 'full';\n", + "\n", + " this.root = document.createElement('div');\n", + " this.root.setAttribute('style', 'display: inline-block');\n", + " this._root_extra_style(this.root);\n", + "\n", + " parent_element.appendChild(this.root);\n", + "\n", + " this._init_header(this);\n", + " this._init_canvas(this);\n", + " this._init_toolbar(this);\n", + "\n", + " var fig = this;\n", + "\n", + " this.waiting = false;\n", + "\n", + " this.ws.onopen = function () {\n", + " fig.send_message('supports_binary', { value: fig.supports_binary });\n", + " fig.send_message('send_image_mode', {});\n", + " if (fig.ratio !== 1) {\n", + " fig.send_message('set_device_pixel_ratio', {\n", + " device_pixel_ratio: fig.ratio,\n", + " });\n", + " }\n", + " fig.send_message('refresh', {});\n", + " };\n", + "\n", + " this.imageObj.onload = function () {\n", + " if (fig.image_mode === 'full') {\n", + " // Full images could contain transparency (where diff images\n", + " // almost always do), so we need to clear the canvas so that\n", + " // there is no ghosting.\n", + " fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n", + " }\n", + " fig.context.drawImage(fig.imageObj, 0, 0);\n", + " };\n", + "\n", + " this.imageObj.onunload = function () {\n", + " fig.ws.close();\n", + " };\n", + "\n", + " this.ws.onmessage = this._make_on_message_function(this);\n", + "\n", + " this.ondownload = ondownload;\n", + "};\n", + "\n", + "mpl.figure.prototype._init_header = function () {\n", + " var titlebar = document.createElement('div');\n", + " titlebar.classList =\n", + " 'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n", + " var titletext = document.createElement('div');\n", + " titletext.classList = 'ui-dialog-title';\n", + " titletext.setAttribute(\n", + " 'style',\n", + " 'width: 100%; text-align: center; padding: 3px;'\n", + " );\n", + " titlebar.appendChild(titletext);\n", + " this.root.appendChild(titlebar);\n", + " this.header = titletext;\n", + "};\n", + "\n", + "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n", + "\n", + "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n", + "\n", + "mpl.figure.prototype._init_canvas = function () {\n", + " var fig = this;\n", + "\n", + " var canvas_div = (this.canvas_div = document.createElement('div'));\n", + " canvas_div.setAttribute(\n", + " 'style',\n", + " 'border: 1px solid #ddd;' +\n", + " 'box-sizing: content-box;' +\n", + " 'clear: both;' +\n", + " 'min-height: 1px;' +\n", + " 'min-width: 1px;' +\n", + " 'outline: 0;' +\n", + " 'overflow: hidden;' +\n", + " 'position: relative;' +\n", + " 'resize: both;'\n", + " );\n", + "\n", + " function on_keyboard_event_closure(name) {\n", + " return function (event) {\n", + " return fig.key_event(event, name);\n", + " };\n", + " }\n", + "\n", + " canvas_div.addEventListener(\n", + " 'keydown',\n", + " on_keyboard_event_closure('key_press')\n", + " );\n", + " canvas_div.addEventListener(\n", + " 'keyup',\n", + " on_keyboard_event_closure('key_release')\n", + " );\n", + "\n", + " this._canvas_extra_style(canvas_div);\n", + " this.root.appendChild(canvas_div);\n", + "\n", + " var canvas = (this.canvas = document.createElement('canvas'));\n", + " canvas.classList.add('mpl-canvas');\n", + " canvas.setAttribute('style', 'box-sizing: content-box;');\n", + "\n", + " this.context = canvas.getContext('2d');\n", + "\n", + " var backingStore =\n", + " this.context.backingStorePixelRatio ||\n", + " this.context.webkitBackingStorePixelRatio ||\n", + " this.context.mozBackingStorePixelRatio ||\n", + " this.context.msBackingStorePixelRatio ||\n", + " this.context.oBackingStorePixelRatio ||\n", + " this.context.backingStorePixelRatio ||\n", + " 1;\n", + "\n", + " this.ratio = (window.devicePixelRatio || 1) / backingStore;\n", + "\n", + " var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n", + " 'canvas'\n", + " ));\n", + " rubberband_canvas.setAttribute(\n", + " 'style',\n", + " 'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n", + " );\n", + "\n", + " // Apply a ponyfill if ResizeObserver is not implemented by browser.\n", + " if (this.ResizeObserver === undefined) {\n", + " if (window.ResizeObserver !== undefined) {\n", + " this.ResizeObserver = window.ResizeObserver;\n", + " } else {\n", + " var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n", + " this.ResizeObserver = obs.ResizeObserver;\n", + " }\n", + " }\n", + "\n", + " this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n", + " var nentries = entries.length;\n", + " for (var i = 0; i < nentries; i++) {\n", + " var entry = entries[i];\n", + " var width, height;\n", + " if (entry.contentBoxSize) {\n", + " if (entry.contentBoxSize instanceof Array) {\n", + " // Chrome 84 implements new version of spec.\n", + " width = entry.contentBoxSize[0].inlineSize;\n", + " height = entry.contentBoxSize[0].blockSize;\n", + " } else {\n", + " // Firefox implements old version of spec.\n", + " width = entry.contentBoxSize.inlineSize;\n", + " height = entry.contentBoxSize.blockSize;\n", + " }\n", + " } else {\n", + " // Chrome <84 implements even older version of spec.\n", + " width = entry.contentRect.width;\n", + " height = entry.contentRect.height;\n", + " }\n", + "\n", + " // Keep the size of the canvas and rubber band canvas in sync with\n", + " // the canvas container.\n", + " if (entry.devicePixelContentBoxSize) {\n", + " // Chrome 84 implements new version of spec.\n", + " canvas.setAttribute(\n", + " 'width',\n", + " entry.devicePixelContentBoxSize[0].inlineSize\n", + " );\n", + " canvas.setAttribute(\n", + " 'height',\n", + " entry.devicePixelContentBoxSize[0].blockSize\n", + " );\n", + " } else {\n", + " canvas.setAttribute('width', width * fig.ratio);\n", + " canvas.setAttribute('height', height * fig.ratio);\n", + " }\n", + " canvas.setAttribute(\n", + " 'style',\n", + " 'width: ' + width + 'px; height: ' + height + 'px;'\n", + " );\n", + "\n", + " rubberband_canvas.setAttribute('width', width);\n", + " rubberband_canvas.setAttribute('height', height);\n", + "\n", + " // And update the size in Python. We ignore the initial 0/0 size\n", + " // that occurs as the element is placed into the DOM, which should\n", + " // otherwise not happen due to the minimum size styling.\n", + " if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n", + " fig.request_resize(width, height);\n", + " }\n", + " }\n", + " });\n", + " this.resizeObserverInstance.observe(canvas_div);\n", + "\n", + " function on_mouse_event_closure(name) {\n", + " return function (event) {\n", + " return fig.mouse_event(event, name);\n", + " };\n", + " }\n", + "\n", + " rubberband_canvas.addEventListener(\n", + " 'mousedown',\n", + " on_mouse_event_closure('button_press')\n", + " );\n", + " rubberband_canvas.addEventListener(\n", + " 'mouseup',\n", + " on_mouse_event_closure('button_release')\n", + " );\n", + " rubberband_canvas.addEventListener(\n", + " 'dblclick',\n", + " on_mouse_event_closure('dblclick')\n", + " );\n", + " // Throttle sequential mouse events to 1 every 20ms.\n", + " rubberband_canvas.addEventListener(\n", + " 'mousemove',\n", + " on_mouse_event_closure('motion_notify')\n", + " );\n", + "\n", + " rubberband_canvas.addEventListener(\n", + " 'mouseenter',\n", + " on_mouse_event_closure('figure_enter')\n", + " );\n", + " rubberband_canvas.addEventListener(\n", + " 'mouseleave',\n", + " on_mouse_event_closure('figure_leave')\n", + " );\n", + "\n", + " canvas_div.addEventListener('wheel', function (event) {\n", + " if (event.deltaY < 0) {\n", + " event.step = 1;\n", + " } else {\n", + " event.step = -1;\n", + " }\n", + " on_mouse_event_closure('scroll')(event);\n", + " });\n", + "\n", + " canvas_div.appendChild(canvas);\n", + " canvas_div.appendChild(rubberband_canvas);\n", + "\n", + " this.rubberband_context = rubberband_canvas.getContext('2d');\n", + " this.rubberband_context.strokeStyle = '#000000';\n", + "\n", + " this._resize_canvas = function (width, height, forward) {\n", + " if (forward) {\n", + " canvas_div.style.width = width + 'px';\n", + " canvas_div.style.height = height + 'px';\n", + " }\n", + " };\n", + "\n", + " // Disable right mouse context menu.\n", + " this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n", + " event.preventDefault();\n", + " return false;\n", + " });\n", + "\n", + " function set_focus() {\n", + " canvas.focus();\n", + " canvas_div.focus();\n", + " }\n", + "\n", + " window.setTimeout(set_focus, 100);\n", + "};\n", + "\n", + "mpl.figure.prototype._init_toolbar = function () {\n", + " var fig = this;\n", + "\n", + " var toolbar = document.createElement('div');\n", + " toolbar.classList = 'mpl-toolbar';\n", + " this.root.appendChild(toolbar);\n", + "\n", + " function on_click_closure(name) {\n", + " return function (_event) {\n", + " return fig.toolbar_button_onclick(name);\n", + " };\n", + " }\n", + "\n", + " function on_mouseover_closure(tooltip) {\n", + " return function (event) {\n", + " if (!event.currentTarget.disabled) {\n", + " return fig.toolbar_button_onmouseover(tooltip);\n", + " }\n", + " };\n", + " }\n", + "\n", + " fig.buttons = {};\n", + " var buttonGroup = document.createElement('div');\n", + " buttonGroup.classList = 'mpl-button-group';\n", + " for (var toolbar_ind in mpl.toolbar_items) {\n", + " var name = mpl.toolbar_items[toolbar_ind][0];\n", + " var tooltip = mpl.toolbar_items[toolbar_ind][1];\n", + " var image = mpl.toolbar_items[toolbar_ind][2];\n", + " var method_name = mpl.toolbar_items[toolbar_ind][3];\n", + "\n", + " if (!name) {\n", + " /* Instead of a spacer, we start a new button group. */\n", + " if (buttonGroup.hasChildNodes()) {\n", + " toolbar.appendChild(buttonGroup);\n", + " }\n", + " buttonGroup = document.createElement('div');\n", + " buttonGroup.classList = 'mpl-button-group';\n", + " continue;\n", + " }\n", + "\n", + " var button = (fig.buttons[name] = document.createElement('button'));\n", + " button.classList = 'mpl-widget';\n", + " button.setAttribute('role', 'button');\n", + " button.setAttribute('aria-disabled', 'false');\n", + " button.addEventListener('click', on_click_closure(method_name));\n", + " button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n", + "\n", + " var icon_img = document.createElement('img');\n", + " icon_img.src = '_images/' + image + '.png';\n", + " icon_img.srcset = '_images/' + image + '_large.png 2x';\n", + " icon_img.alt = tooltip;\n", + " button.appendChild(icon_img);\n", + "\n", + " buttonGroup.appendChild(button);\n", + " }\n", + "\n", + " if (buttonGroup.hasChildNodes()) {\n", + " toolbar.appendChild(buttonGroup);\n", + " }\n", + "\n", + " var fmt_picker = document.createElement('select');\n", + " fmt_picker.classList = 'mpl-widget';\n", + " toolbar.appendChild(fmt_picker);\n", + " this.format_dropdown = fmt_picker;\n", + "\n", + " for (var ind in mpl.extensions) {\n", + " var fmt = mpl.extensions[ind];\n", + " var option = document.createElement('option');\n", + " option.selected = fmt === mpl.default_extension;\n", + " option.innerHTML = fmt;\n", + " fmt_picker.appendChild(option);\n", + " }\n", + "\n", + " var status_bar = document.createElement('span');\n", + " status_bar.classList = 'mpl-message';\n", + " toolbar.appendChild(status_bar);\n", + " this.message = status_bar;\n", + "};\n", + "\n", + "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n", + " // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n", + " // which will in turn request a refresh of the image.\n", + " this.send_message('resize', { width: x_pixels, height: y_pixels });\n", + "};\n", + "\n", + "mpl.figure.prototype.send_message = function (type, properties) {\n", + " properties['type'] = type;\n", + " properties['figure_id'] = this.id;\n", + " this.ws.send(JSON.stringify(properties));\n", + "};\n", + "\n", + "mpl.figure.prototype.send_draw_message = function () {\n", + " if (!this.waiting) {\n", + " this.waiting = true;\n", + " this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_save = function (fig, _msg) {\n", + " var format_dropdown = fig.format_dropdown;\n", + " var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n", + " fig.ondownload(fig, format);\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_resize = function (fig, msg) {\n", + " var size = msg['size'];\n", + " if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n", + " fig._resize_canvas(size[0], size[1], msg['forward']);\n", + " fig.send_message('refresh', {});\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n", + " var x0 = msg['x0'] / fig.ratio;\n", + " var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n", + " var x1 = msg['x1'] / fig.ratio;\n", + " var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n", + " x0 = Math.floor(x0) + 0.5;\n", + " y0 = Math.floor(y0) + 0.5;\n", + " x1 = Math.floor(x1) + 0.5;\n", + " y1 = Math.floor(y1) + 0.5;\n", + " var min_x = Math.min(x0, x1);\n", + " var min_y = Math.min(y0, y1);\n", + " var width = Math.abs(x1 - x0);\n", + " var height = Math.abs(y1 - y0);\n", + "\n", + " fig.rubberband_context.clearRect(\n", + " 0,\n", + " 0,\n", + " fig.canvas.width / fig.ratio,\n", + " fig.canvas.height / fig.ratio\n", + " );\n", + "\n", + " fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n", + " // Updates the figure title.\n", + " fig.header.textContent = msg['label'];\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n", + " fig.rubberband_canvas.style.cursor = msg['cursor'];\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_message = function (fig, msg) {\n", + " fig.message.textContent = msg['message'];\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n", + " // Request the server to send over a new figure.\n", + " fig.send_draw_message();\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n", + " fig.image_mode = msg['mode'];\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n", + " for (var key in msg) {\n", + " if (!(key in fig.buttons)) {\n", + " continue;\n", + " }\n", + " fig.buttons[key].disabled = !msg[key];\n", + " fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n", + " if (msg['mode'] === 'PAN') {\n", + " fig.buttons['Pan'].classList.add('active');\n", + " fig.buttons['Zoom'].classList.remove('active');\n", + " } else if (msg['mode'] === 'ZOOM') {\n", + " fig.buttons['Pan'].classList.remove('active');\n", + " fig.buttons['Zoom'].classList.add('active');\n", + " } else {\n", + " fig.buttons['Pan'].classList.remove('active');\n", + " fig.buttons['Zoom'].classList.remove('active');\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.updated_canvas_event = function () {\n", + " // Called whenever the canvas gets updated.\n", + " this.send_message('ack', {});\n", + "};\n", + "\n", + "// A function to construct a web socket function for onmessage handling.\n", + "// Called in the figure constructor.\n", + "mpl.figure.prototype._make_on_message_function = function (fig) {\n", + " return function socket_on_message(evt) {\n", + " if (evt.data instanceof Blob) {\n", + " var img = evt.data;\n", + " if (img.type !== 'image/png') {\n", + " /* FIXME: We get \"Resource interpreted as Image but\n", + " * transferred with MIME type text/plain:\" errors on\n", + " * Chrome. But how to set the MIME type? It doesn't seem\n", + " * to be part of the websocket stream */\n", + " img.type = 'image/png';\n", + " }\n", + "\n", + " /* Free the memory for the previous frames */\n", + " if (fig.imageObj.src) {\n", + " (window.URL || window.webkitURL).revokeObjectURL(\n", + " fig.imageObj.src\n", + " );\n", + " }\n", + "\n", + " fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n", + " img\n", + " );\n", + " fig.updated_canvas_event();\n", + " fig.waiting = false;\n", + " return;\n", + " } else if (\n", + " typeof evt.data === 'string' &&\n", + " evt.data.slice(0, 21) === 'data:image/png;base64'\n", + " ) {\n", + " fig.imageObj.src = evt.data;\n", + " fig.updated_canvas_event();\n", + " fig.waiting = false;\n", + " return;\n", + " }\n", + "\n", + " var msg = JSON.parse(evt.data);\n", + " var msg_type = msg['type'];\n", + "\n", + " // Call the \"handle_{type}\" callback, which takes\n", + " // the figure and JSON message as its only arguments.\n", + " try {\n", + " var callback = fig['handle_' + msg_type];\n", + " } catch (e) {\n", + " console.log(\n", + " \"No handler for the '\" + msg_type + \"' message type: \",\n", + " msg\n", + " );\n", + " return;\n", + " }\n", + "\n", + " if (callback) {\n", + " try {\n", + " // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n", + " callback(fig, msg);\n", + " } catch (e) {\n", + " console.log(\n", + " \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n", + " e,\n", + " e.stack,\n", + " msg\n", + " );\n", + " }\n", + " }\n", + " };\n", + "};\n", + "\n", + "// from https://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n", + "mpl.findpos = function (e) {\n", + " //this section is from http://www.quirksmode.org/js/events_properties.html\n", + " var targ;\n", + " if (!e) {\n", + " e = window.event;\n", + " }\n", + " if (e.target) {\n", + " targ = e.target;\n", + " } else if (e.srcElement) {\n", + " targ = e.srcElement;\n", + " }\n", + " if (targ.nodeType === 3) {\n", + " // defeat Safari bug\n", + " targ = targ.parentNode;\n", + " }\n", + "\n", + " // pageX,Y are the mouse positions relative to the document\n", + " var boundingRect = targ.getBoundingClientRect();\n", + " var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n", + " var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n", + "\n", + " return { x: x, y: y };\n", + "};\n", + "\n", + "/*\n", + " * return a copy of an object with only non-object keys\n", + " * we need this to avoid circular references\n", + " * https://stackoverflow.com/a/24161582/3208463\n", + " */\n", + "function simpleKeys(original) {\n", + " return Object.keys(original).reduce(function (obj, key) {\n", + " if (typeof original[key] !== 'object') {\n", + " obj[key] = original[key];\n", + " }\n", + " return obj;\n", + " }, {});\n", + "}\n", + "\n", + "mpl.figure.prototype.mouse_event = function (event, name) {\n", + " var canvas_pos = mpl.findpos(event);\n", + "\n", + " if (name === 'button_press') {\n", + " this.canvas.focus();\n", + " this.canvas_div.focus();\n", + " }\n", + "\n", + " var x = canvas_pos.x * this.ratio;\n", + " var y = canvas_pos.y * this.ratio;\n", + "\n", + " this.send_message(name, {\n", + " x: x,\n", + " y: y,\n", + " button: event.button,\n", + " step: event.step,\n", + " guiEvent: simpleKeys(event),\n", + " });\n", + "\n", + " /* This prevents the web browser from automatically changing to\n", + " * the text insertion cursor when the button is pressed. We want\n", + " * to control all of the cursor setting manually through the\n", + " * 'cursor' event from matplotlib */\n", + " event.preventDefault();\n", + " return false;\n", + "};\n", + "\n", + "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n", + " // Handle any extra behaviour associated with a key event\n", + "};\n", + "\n", + "mpl.figure.prototype.key_event = function (event, name) {\n", + " // Prevent repeat events\n", + " if (name === 'key_press') {\n", + " if (event.key === this._key) {\n", + " return;\n", + " } else {\n", + " this._key = event.key;\n", + " }\n", + " }\n", + " if (name === 'key_release') {\n", + " this._key = null;\n", + " }\n", + "\n", + " var value = '';\n", + " if (event.ctrlKey && event.key !== 'Control') {\n", + " value += 'ctrl+';\n", + " }\n", + " else if (event.altKey && event.key !== 'Alt') {\n", + " value += 'alt+';\n", + " }\n", + " else if (event.shiftKey && event.key !== 'Shift') {\n", + " value += 'shift+';\n", + " }\n", + "\n", + " value += 'k' + event.key;\n", + "\n", + " this._key_event_extra(event, name);\n", + "\n", + " this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n", + " return false;\n", + "};\n", + "\n", + "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n", + " if (name === 'download') {\n", + " this.handle_save(this, null);\n", + " } else {\n", + " this.send_message('toolbar_button', { name: name });\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n", + " this.message.textContent = tooltip;\n", + "};\n", + "\n", + "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n", + "// prettier-ignore\n", + "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n", + "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n", + "\n", + "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n", + "\n", + "mpl.default_extension = \"png\";/* global mpl */\n", + "\n", + "var comm_websocket_adapter = function (comm) {\n", + " // Create a \"websocket\"-like object which calls the given IPython comm\n", + " // object with the appropriate methods. Currently this is a non binary\n", + " // socket, so there is still some room for performance tuning.\n", + " var ws = {};\n", + "\n", + " ws.binaryType = comm.kernel.ws.binaryType;\n", + " ws.readyState = comm.kernel.ws.readyState;\n", + " function updateReadyState(_event) {\n", + " if (comm.kernel.ws) {\n", + " ws.readyState = comm.kernel.ws.readyState;\n", + " } else {\n", + " ws.readyState = 3; // Closed state.\n", + " }\n", + " }\n", + " comm.kernel.ws.addEventListener('open', updateReadyState);\n", + " comm.kernel.ws.addEventListener('close', updateReadyState);\n", + " comm.kernel.ws.addEventListener('error', updateReadyState);\n", + "\n", + " ws.close = function () {\n", + " comm.close();\n", + " };\n", + " ws.send = function (m) {\n", + " //console.log('sending', m);\n", + " comm.send(m);\n", + " };\n", + " // Register the callback with on_msg.\n", + " comm.on_msg(function (msg) {\n", + " //console.log('receiving', msg['content']['data'], msg);\n", + " var data = msg['content']['data'];\n", + " if (data['blob'] !== undefined) {\n", + " data = {\n", + " data: new Blob(msg['buffers'], { type: data['blob'] }),\n", + " };\n", + " }\n", + " // Pass the mpl event to the overridden (by mpl) onmessage function.\n", + " ws.onmessage(data);\n", + " });\n", + " return ws;\n", + "};\n", + "\n", + "mpl.mpl_figure_comm = function (comm, msg) {\n", + " // This is the function which gets called when the mpl process\n", + " // starts-up an IPython Comm through the \"matplotlib\" channel.\n", + "\n", + " var id = msg.content.data.id;\n", + " // Get hold of the div created by the display call when the Comm\n", + " // socket was opened in Python.\n", + " var element = document.getElementById(id);\n", + " var ws_proxy = comm_websocket_adapter(comm);\n", + "\n", + " function ondownload(figure, _format) {\n", + " window.open(figure.canvas.toDataURL());\n", + " }\n", + "\n", + " var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n", + "\n", + " // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n", + " // web socket which is closed, not our websocket->open comm proxy.\n", + " ws_proxy.onopen();\n", + "\n", + " fig.parent_element = element;\n", + " fig.cell_info = mpl.find_output_cell(\"
\");\n", + " if (!fig.cell_info) {\n", + " console.error('Failed to find cell for figure', id, fig);\n", + " return;\n", + " }\n", + " fig.cell_info[0].output_area.element.on(\n", + " 'cleared',\n", + " { fig: fig },\n", + " fig._remove_fig_handler\n", + " );\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_close = function (fig, msg) {\n", + " var width = fig.canvas.width / fig.ratio;\n", + " fig.cell_info[0].output_area.element.off(\n", + " 'cleared',\n", + " fig._remove_fig_handler\n", + " );\n", + " fig.resizeObserverInstance.unobserve(fig.canvas_div);\n", + "\n", + " // Update the output cell to use the data from the current canvas.\n", + " fig.push_to_output();\n", + " var dataURL = fig.canvas.toDataURL();\n", + " // Re-enable the keyboard manager in IPython - without this line, in FF,\n", + " // the notebook keyboard shortcuts fail.\n", + " IPython.keyboard_manager.enable();\n", + " fig.parent_element.innerHTML =\n", + " '';\n", + " fig.close_ws(fig, msg);\n", + "};\n", + "\n", + "mpl.figure.prototype.close_ws = function (fig, msg) {\n", + " fig.send_message('closing', msg);\n", + " // fig.ws.close()\n", + "};\n", + "\n", + "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n", + " // Turn the data on the canvas into data in the output cell.\n", + " var width = this.canvas.width / this.ratio;\n", + " var dataURL = this.canvas.toDataURL();\n", + " this.cell_info[1]['text/html'] =\n", + " '';\n", + "};\n", + "\n", + "mpl.figure.prototype.updated_canvas_event = function () {\n", + " // Tell IPython that the notebook contents must change.\n", + " IPython.notebook.set_dirty(true);\n", + " this.send_message('ack', {});\n", + " var fig = this;\n", + " // Wait a second, then push the new image to the DOM so\n", + " // that it is saved nicely (might be nice to debounce this).\n", + " setTimeout(function () {\n", + " fig.push_to_output();\n", + " }, 1000);\n", + "};\n", + "\n", + "mpl.figure.prototype._init_toolbar = function () {\n", + " var fig = this;\n", + "\n", + " var toolbar = document.createElement('div');\n", + " toolbar.classList = 'btn-toolbar';\n", + " this.root.appendChild(toolbar);\n", + "\n", + " function on_click_closure(name) {\n", + " return function (_event) {\n", + " return fig.toolbar_button_onclick(name);\n", + " };\n", + " }\n", + "\n", + " function on_mouseover_closure(tooltip) {\n", + " return function (event) {\n", + " if (!event.currentTarget.disabled) {\n", + " return fig.toolbar_button_onmouseover(tooltip);\n", + " }\n", + " };\n", + " }\n", + "\n", + " fig.buttons = {};\n", + " var buttonGroup = document.createElement('div');\n", + " buttonGroup.classList = 'btn-group';\n", + " var button;\n", + " for (var toolbar_ind in mpl.toolbar_items) {\n", + " var name = mpl.toolbar_items[toolbar_ind][0];\n", + " var tooltip = mpl.toolbar_items[toolbar_ind][1];\n", + " var image = mpl.toolbar_items[toolbar_ind][2];\n", + " var method_name = mpl.toolbar_items[toolbar_ind][3];\n", + "\n", + " if (!name) {\n", + " /* Instead of a spacer, we start a new button group. */\n", + " if (buttonGroup.hasChildNodes()) {\n", + " toolbar.appendChild(buttonGroup);\n", + " }\n", + " buttonGroup = document.createElement('div');\n", + " buttonGroup.classList = 'btn-group';\n", + " continue;\n", + " }\n", + "\n", + " button = fig.buttons[name] = document.createElement('button');\n", + " button.classList = 'btn btn-default';\n", + " button.href = '#';\n", + " button.title = name;\n", + " button.innerHTML = '';\n", + " button.addEventListener('click', on_click_closure(method_name));\n", + " button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n", + " buttonGroup.appendChild(button);\n", + " }\n", + "\n", + " if (buttonGroup.hasChildNodes()) {\n", + " toolbar.appendChild(buttonGroup);\n", + " }\n", + "\n", + " // Add the status bar.\n", + " var status_bar = document.createElement('span');\n", + " status_bar.classList = 'mpl-message pull-right';\n", + " toolbar.appendChild(status_bar);\n", + " this.message = status_bar;\n", + "\n", + " // Add the close button to the window.\n", + " var buttongrp = document.createElement('div');\n", + " buttongrp.classList = 'btn-group inline pull-right';\n", + " button = document.createElement('button');\n", + " button.classList = 'btn btn-mini btn-primary';\n", + " button.href = '#';\n", + " button.title = 'Stop Interaction';\n", + " button.innerHTML = '';\n", + " button.addEventListener('click', function (_evt) {\n", + " fig.handle_close(fig, {});\n", + " });\n", + " button.addEventListener(\n", + " 'mouseover',\n", + " on_mouseover_closure('Stop Interaction')\n", + " );\n", + " buttongrp.appendChild(button);\n", + " var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n", + " titlebar.insertBefore(buttongrp, titlebar.firstChild);\n", + "};\n", + "\n", + "mpl.figure.prototype._remove_fig_handler = function (event) {\n", + " var fig = event.data.fig;\n", + " if (event.target !== this) {\n", + " // Ignore bubbled events from children.\n", + " return;\n", + " }\n", + " fig.close_ws(fig, {});\n", + "};\n", + "\n", + "mpl.figure.prototype._root_extra_style = function (el) {\n", + " el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n", + "};\n", + "\n", + "mpl.figure.prototype._canvas_extra_style = function (el) {\n", + " // this is important to make the div 'focusable\n", + " el.setAttribute('tabindex', 0);\n", + " // reach out to IPython and tell the keyboard manager to turn it's self\n", + " // off when our div gets focus\n", + "\n", + " // location in version 3\n", + " if (IPython.notebook.keyboard_manager) {\n", + " IPython.notebook.keyboard_manager.register_events(el);\n", + " } else {\n", + " // location in version 2\n", + " IPython.keyboard_manager.register_events(el);\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype._key_event_extra = function (event, _name) {\n", + " // Check for shift+enter\n", + " if (event.shiftKey && event.which === 13) {\n", + " this.canvas_div.blur();\n", + " // select the cell after this one\n", + " var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n", + " IPython.notebook.select(index + 1);\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_save = function (fig, _msg) {\n", + " fig.ondownload(fig, null);\n", + "};\n", + "\n", + "mpl.find_output_cell = function (html_output) {\n", + " // Return the cell and output element which can be found *uniquely* in the notebook.\n", + " // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n", + " // IPython event is triggered only after the cells have been serialised, which for\n", + " // our purposes (turning an active figure into a static one), is too late.\n", + " var cells = IPython.notebook.get_cells();\n", + " var ncells = cells.length;\n", + " for (var i = 0; i < ncells; i++) {\n", + " var cell = cells[i];\n", + " if (cell.cell_type === 'code') {\n", + " for (var j = 0; j < cell.output_area.outputs.length; j++) {\n", + " var data = cell.output_area.outputs[j];\n", + " if (data.data) {\n", + " // IPython >= 3 moved mimebundle to data attribute of output\n", + " data = data.data;\n", + " }\n", + " if (data['text/html'] === html_output) {\n", + " return [cell, data, j];\n", + " }\n", + " }\n", + " }\n", + " }\n", + "};\n", + "\n", + "// Register the function which deals with the matplotlib target/channel.\n", + "// The kernel may be null if the page has been refreshed.\n", + "if (IPython.notebook.kernel !== null) {\n", + " IPython.notebook.kernel.comm_manager.register_target(\n", + " 'matplotlib',\n", + " mpl.mpl_figure_comm\n", + " );\n", + "}\n" + ], + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/html": [ + "" + ], + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Constructing Hamiltonian and computing ground state of order 0.\n", + "Exact eigenvalue computed : -3.493959207434937 eV\n", + "Computing Eigenvector for new Hamiltonian\n" + ] + }, + { + "data": { + "application/vnd.jupyter.widget-view+json": { + "model_id": "40dad3d8a98c4e1fb4a0a5375eb9d2bf", + "version_major": 2, + "version_minor": 0 + }, + "text/plain": [ + " 0%| | 0/50 [00:00open comm proxy.\n", + " ws_proxy.onopen();\n", + "\n", + " fig.parent_element = element;\n", + " fig.cell_info = mpl.find_output_cell(\"
\");\n", + " if (!fig.cell_info) {\n", + " console.error('Failed to find cell for figure', id, fig);\n", + " return;\n", + " }\n", + " fig.cell_info[0].output_area.element.on(\n", + " 'cleared',\n", + " { fig: fig },\n", + " fig._remove_fig_handler\n", + " );\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_close = function (fig, msg) {\n", + " var width = fig.canvas.width / fig.ratio;\n", + " fig.cell_info[0].output_area.element.off(\n", + " 'cleared',\n", + " fig._remove_fig_handler\n", + " );\n", + " fig.resizeObserverInstance.unobserve(fig.canvas_div);\n", + "\n", + " // Update the output cell to use the data from the current canvas.\n", + " fig.push_to_output();\n", + " var dataURL = fig.canvas.toDataURL();\n", + " // Re-enable the keyboard manager in IPython - without this line, in FF,\n", + " // the notebook keyboard shortcuts fail.\n", + " IPython.keyboard_manager.enable();\n", + " fig.parent_element.innerHTML =\n", + " '';\n", + " fig.close_ws(fig, msg);\n", + "};\n", + "\n", + "mpl.figure.prototype.close_ws = function (fig, msg) {\n", + " fig.send_message('closing', msg);\n", + " // fig.ws.close()\n", + "};\n", + "\n", + "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n", + " // Turn the data on the canvas into data in the output cell.\n", + " var width = this.canvas.width / this.ratio;\n", + " var dataURL = this.canvas.toDataURL();\n", + " this.cell_info[1]['text/html'] =\n", + " '';\n", + "};\n", + "\n", + "mpl.figure.prototype.updated_canvas_event = function () {\n", + " // Tell IPython that the notebook contents must change.\n", + " IPython.notebook.set_dirty(true);\n", + " this.send_message('ack', {});\n", + " var fig = this;\n", + " // Wait a second, then push the new image to the DOM so\n", + " // that it is saved nicely (might be nice to debounce this).\n", + " setTimeout(function () {\n", + " fig.push_to_output();\n", + " }, 1000);\n", + "};\n", + "\n", + "mpl.figure.prototype._init_toolbar = function () {\n", + " var fig = this;\n", + "\n", + " var toolbar = document.createElement('div');\n", + " toolbar.classList = 'btn-toolbar';\n", + " this.root.appendChild(toolbar);\n", + "\n", + " function on_click_closure(name) {\n", + " return function (_event) {\n", + " return fig.toolbar_button_onclick(name);\n", + " };\n", + " }\n", + "\n", + " function on_mouseover_closure(tooltip) {\n", + " return function (event) {\n", + " if (!event.currentTarget.disabled) {\n", + " return fig.toolbar_button_onmouseover(tooltip);\n", + " }\n", + " };\n", + " }\n", + "\n", + " fig.buttons = {};\n", + " var buttonGroup = document.createElement('div');\n", + " buttonGroup.classList = 'btn-group';\n", + " var button;\n", + " for (var toolbar_ind in mpl.toolbar_items) {\n", + " var name = mpl.toolbar_items[toolbar_ind][0];\n", + " var tooltip = mpl.toolbar_items[toolbar_ind][1];\n", + " var image = mpl.toolbar_items[toolbar_ind][2];\n", + " var method_name = mpl.toolbar_items[toolbar_ind][3];\n", + "\n", + " if (!name) {\n", + " /* Instead of a spacer, we start a new button group. */\n", + " if (buttonGroup.hasChildNodes()) {\n", + " toolbar.appendChild(buttonGroup);\n", + " }\n", + " buttonGroup = document.createElement('div');\n", + " buttonGroup.classList = 'btn-group';\n", + " continue;\n", + " }\n", + "\n", + " button = fig.buttons[name] = document.createElement('button');\n", + " button.classList = 'btn btn-default';\n", + " button.href = '#';\n", + " button.title = name;\n", + " button.innerHTML = '';\n", + " button.addEventListener('click', on_click_closure(method_name));\n", + " button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n", + " buttonGroup.appendChild(button);\n", + " }\n", + "\n", + " if (buttonGroup.hasChildNodes()) {\n", + " toolbar.appendChild(buttonGroup);\n", + " }\n", + "\n", + " // Add the status bar.\n", + " var status_bar = document.createElement('span');\n", + " status_bar.classList = 'mpl-message pull-right';\n", + " toolbar.appendChild(status_bar);\n", + " this.message = status_bar;\n", + "\n", + " // Add the close button to the window.\n", + " var buttongrp = document.createElement('div');\n", + " buttongrp.classList = 'btn-group inline pull-right';\n", + " button = document.createElement('button');\n", + " button.classList = 'btn btn-mini btn-primary';\n", + " button.href = '#';\n", + " button.title = 'Stop Interaction';\n", + " button.innerHTML = '';\n", + " button.addEventListener('click', function (_evt) {\n", + " fig.handle_close(fig, {});\n", + " });\n", + " button.addEventListener(\n", + " 'mouseover',\n", + " on_mouseover_closure('Stop Interaction')\n", + " );\n", + " buttongrp.appendChild(button);\n", + " var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n", + " titlebar.insertBefore(buttongrp, titlebar.firstChild);\n", + "};\n", + "\n", + "mpl.figure.prototype._remove_fig_handler = function (event) {\n", + " var fig = event.data.fig;\n", + " if (event.target !== this) {\n", + " // Ignore bubbled events from children.\n", + " return;\n", + " }\n", + " fig.close_ws(fig, {});\n", + "};\n", + "\n", + "mpl.figure.prototype._root_extra_style = function (el) {\n", + " el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n", + "};\n", + "\n", + "mpl.figure.prototype._canvas_extra_style = function (el) {\n", + " // this is important to make the div 'focusable\n", + " el.setAttribute('tabindex', 0);\n", + " // reach out to IPython and tell the keyboard manager to turn it's self\n", + " // off when our div gets focus\n", + "\n", + " // location in version 3\n", + " if (IPython.notebook.keyboard_manager) {\n", + " IPython.notebook.keyboard_manager.register_events(el);\n", + " } else {\n", + " // location in version 2\n", + " IPython.keyboard_manager.register_events(el);\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype._key_event_extra = function (event, _name) {\n", + " // Check for shift+enter\n", + " if (event.shiftKey && event.which === 13) {\n", + " this.canvas_div.blur();\n", + " // select the cell after this one\n", + " var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n", + " IPython.notebook.select(index + 1);\n", + " }\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_save = function (fig, _msg) {\n", + " fig.ondownload(fig, null);\n", + "};\n", + "\n", + "mpl.find_output_cell = function (html_output) {\n", + " // Return the cell and output element which can be found *uniquely* in the notebook.\n", + " // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n", + " // IPython event is triggered only after the cells have been serialised, which for\n", + " // our purposes (turning an active figure into a static one), is too late.\n", + " var cells = IPython.notebook.get_cells();\n", + " var ncells = cells.length;\n", + " for (var i = 0; i < ncells; i++) {\n", + " var cell = cells[i];\n", + " if (cell.cell_type === 'code') {\n", + " for (var j = 0; j < cell.output_area.outputs.length; j++) {\n", + " var data = cell.output_area.outputs[j];\n", + " if (data.data) {\n", + " // IPython >= 3 moved mimebundle to data attribute of output\n", + " data = data.data;\n", + " }\n", + " if (data['text/html'] === html_output) {\n", + " return [cell, data, j];\n", + " }\n", + " }\n", + " }\n", + " }\n", + "};\n", + "\n", + "// Register the function which deals with the matplotlib target/channel.\n", + "// The kernel may be null if the page has been refreshed.\n", + "if (IPython.notebook.kernel !== null) {\n", + " IPython.notebook.kernel.comm_manager.register_target(\n", + " 'matplotlib',\n", + " mpl.mpl_figure_comm\n", + " );\n", + "}\n" + ], + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/html": [ + "" + ], + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 23, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "\n", + "\n", + "## The figure below shows the convergence\n", + "# Although converging to the ground state was relatively accurate, the |1> and |2> states are very inaccurate\n", + "# this is likely because the error from the ground state carries over, \n", + "# and because the lagrange multiplier needs more iterations to converge\n", + "\n", + "fig = plt.figure()\n", + "for i, loss in enumerate(losses):\n", + " plt.plot(range(50), loss,'b.',label=f'VQE {i}')\n", + " plt.plot(exact_eigenvalues[i]*np.ones_like(range(50)),'k--',label=f'exact eigenvalue {i}')\n", + " \n", + "plt.legend(loc= \"best\")" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.9.7" + } + }, + "nbformat": 4, + "nbformat_minor": 4 +} diff --git a/Team 12 Notebooks/Challenge 2 Team 12.ipynb b/Team 12 Notebooks/Challenge 2 Team 12.ipynb new file mode 100644 index 0000000..bcf4ef7 --- /dev/null +++ b/Team 12 Notebooks/Challenge 2 Team 12.ipynb @@ -0,0 +1,282 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "id": "5b143860", + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit.opflow import Z, I, X, Y\n", + "from qiskit.opflow import CircuitStateFn, StateFn, CircuitSampler, PauliExpectation, ListOp, TensoredOp\n", + "\n", + "from qiskit.providers.aer import AerSimulator\n", + "from qiskit.utils import QuantumInstance\n", + "from qiskit.circuit import QuantumCircuit, QuantumRegister, Parameter, ParameterVector\n", + "from qiskit.algorithms.optimizers import ADAM, SPSA\n", + "from qiskit.circuit.library import TwoLocal\n", + "from qiskit.circuit.library.n_local import EfficientSU2, RealAmplitudes\n", + "from qiskit import Aer, execute\n", + "\n", + "import numpy as np\n", + "import matplotlib.pyplot as plt" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "id": "e92f08f2", + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit_nature.operators.second_quantization import FermionicOp\n", + "from qiskit_nature.mappers.second_quantization import JordanWignerMapper" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "id": "b2512b5c", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Fermionic Operator\n", + "register length=4, number terms=20\n", + " (1+0j) * ( --II )\n", + "+ (1+0j) * ( +-II )\n", + "+ (1+0j) * ( -+II )\n", + "+ (1+0j) * ( ++II )\n", + "+ (-2+0j) * ( --II )\n", + "+ (-2+0j) * ( -+II )\n", + "+ (1+0j) * ( I-+I )\n", + "+ (-1+0j) * ( I++I )\n", + "+ (-1+0j) * ( I--I )\n", + "+ (1+0j) * ( ...\n", + "1.0 * ZZII\n", + "+ 1.0 * IYYI\n", + "+ 1.0 * IIXX\n" + ] + }, + { + "name": "stderr", + "output_type": "stream", + "text": [ + ":1: UserWarning: The default value for `display_format` will be changed from 'dense' to 'sparse' in version 0.3.0. Once that happens, you must specify display_format='dense' directly.\n", + " a0 = FermionicOp('-III')\n" + ] + } + ], + "source": [ + "a0 = FermionicOp('-III')\n", + "a0_dagger = FermionicOp('+III')\n", + "a1 = FermionicOp('I-II')\n", + "a1_dagger = FermionicOp('I+II')\n", + "a2 = FermionicOp('II-I')\n", + "a2_dagger = FermionicOp('II+I')\n", + "a3 = FermionicOp('III-')\n", + "a3_dagger = FermionicOp('III+')\n", + "\n", + "H = (FermionicOp('-III') + FermionicOp('+III')) @ (FermionicOp('I-II') + FermionicOp('I+II')) @ (FermionicOp('IIII') - 2*(FermionicOp('+III') @ FermionicOp('-III'))) \n", + "H += (FermionicOp('I-II') - FermionicOp('I+II')) @ (FermionicOp('II+I') - FermionicOp('II-I')) @ (FermionicOp('IIII') - 2*(FermionicOp('I+II') @ FermionicOp('I-II'))) \n", + "H += (FermionicOp('IIII') - 2*(FermionicOp('II+I') @ FermionicOp('II-I'))) @ (FermionicOp('IIII') - 2*(FermionicOp('III+') @ FermionicOp('III-'))) \n", + "\n", + "print(H)\n", + "\n", + "mapper = JordanWignerMapper()\n", + "\n", + "print(mapper.map(H))\n", + "H = mapper.map(H)" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "id": "73b1f9d9", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
      ┌──────────┐┌──────────┐ ░                                ░  ┌──────────┐»\n",
+       "q0_0: ┤ Ry(θ[0]) ├┤ Rz(θ[4]) ├─░───■────■────■──────────────────░──┤ Ry(θ[8]) ├»\n",
+       "      ├──────────┤├──────────┤ ░ ┌─┴─┐  │    │                  ░  ├──────────┤»\n",
+       "q0_1: ┤ Ry(θ[1]) ├┤ Rz(θ[5]) ├─░─┤ X ├──┼────┼────■────■────────░──┤ Ry(θ[9]) ├»\n",
+       "      ├──────────┤├──────────┤ ░ └───┘┌─┴─┐  │  ┌─┴─┐  │        ░ ┌┴──────────┤»\n",
+       "q0_2: ┤ Ry(θ[2]) ├┤ Rz(θ[6]) ├─░──────┤ X ├──┼──┤ X ├──┼────■───░─┤ Ry(θ[10]) ├»\n",
+       "      ├──────────┤├──────────┤ ░      └───┘┌─┴─┐└───┘┌─┴─┐┌─┴─┐ ░ ├───────────┤»\n",
+       "q0_3: ┤ Ry(θ[3]) ├┤ Rz(θ[7]) ├─░───────────┤ X ├─────┤ X ├┤ X ├─░─┤ Ry(θ[11]) ├»\n",
+       "      └──────────┘└──────────┘ ░           └───┘     └───┘└───┘ ░ └───────────┘»\n",
+       "«      ┌───────────┐\n",
+       "«q0_0: ┤ Rz(θ[12]) ├\n",
+       "«      ├───────────┤\n",
+       "«q0_1: ┤ Rz(θ[13]) ├\n",
+       "«      ├───────────┤\n",
+       "«q0_2: ┤ Rz(θ[14]) ├\n",
+       "«      ├───────────┤\n",
+       "«q0_3: ┤ Rz(θ[15]) ├\n",
+       "«      └───────────┘
" + ], + "text/plain": [ + " ┌──────────┐┌──────────┐ ░ ░ ┌──────────┐»\n", + "q0_0: ┤ Ry(θ[0]) ├┤ Rz(θ[4]) ├─░───■────■────■──────────────────░──┤ Ry(θ[8]) ├»\n", + " ├──────────┤├──────────┤ ░ ┌─┴─┐ │ │ ░ ├──────────┤»\n", + "q0_1: ┤ Ry(θ[1]) ├┤ Rz(θ[5]) ├─░─┤ X ├──┼────┼────■────■────────░──┤ Ry(θ[9]) ├»\n", + " ├──────────┤├──────────┤ ░ └───┘┌─┴─┐ │ ┌─┴─┐ │ ░ ┌┴──────────┤»\n", + "q0_2: ┤ Ry(θ[2]) ├┤ Rz(θ[6]) ├─░──────┤ X ├──┼──┤ X ├──┼────■───░─┤ Ry(θ[10]) ├»\n", + " ├──────────┤├──────────┤ ░ └───┘┌─┴─┐└───┘┌─┴─┐┌─┴─┐ ░ ├───────────┤»\n", + "q0_3: ┤ Ry(θ[3]) ├┤ Rz(θ[7]) ├─░───────────┤ X ├─────┤ X ├┤ X ├─░─┤ Ry(θ[11]) ├»\n", + " └──────────┘└──────────┘ ░ └───┘ └───┘└───┘ ░ └───────────┘»\n", + "« ┌───────────┐\n", + "«q0_0: ┤ Rz(θ[12]) ├\n", + "« ├───────────┤\n", + "«q0_1: ┤ Rz(θ[13]) ├\n", + "« ├───────────┤\n", + "«q0_2: ┤ Rz(θ[14]) ├\n", + "« ├───────────┤\n", + "«q0_3: ┤ Rz(θ[15]) ├\n", + "« └───────────┘" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "n_qubits = 4\n", + "params = ParameterVector('θ',(n_qubits*4))\n", + "\n", + "qr = QuantumRegister(n_qubits)\n", + "ansatz = QuantumCircuit(qr)\n", + "\n", + " \n", + "ansatz.compose(EfficientSU2(4, reps=1, insert_barriers=True).decompose(), inplace=True)\n", + "ansatz.draw()" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "id": "1cf897d0", + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit.providers.aer import AerSimulator\n", + "\n", + "sim = AerSimulator()\n", + "\n", + "from qiskit.algorithms.optimizers import SPSA\n", + "\n", + "optimizer = SPSA(maxiter=300)\n", + "\n", + "from qiskit.algorithms import VQE as VQE_q\n", + "\n", + "vqe = VQE_q(\n", + " ansatz=ansatz,\n", + " optimizer=optimizer,\n", + " quantum_instance=sim # use qiskit.utils.QuantumInstance for error mitigation and other features\n", + ")\n", + "\n", + "vqe_result = vqe.compute_minimum_eigenvalue(H)" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "id": "05af4242", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "{ 'aux_operator_eigenvalues': None,\n", + " 'cost_function_evals': 600,\n", + " 'eigenstate': { '0001': 0.15625,\n", + " '0010': 0.16237976320958225,\n", + " '0100': 0.6952686081652184,\n", + " '0111': 0.6817945071647321,\n", + " '1000': 0.03125},\n", + " 'eigenvalue': (-2.23828125+0j),\n", + " 'optimal_parameters': { ParameterVectorElement(θ[7]): -6.395952848173277,\n", + " ParameterVectorElement(θ[4]): 2.250509332217714,\n", + " ParameterVectorElement(θ[1]): 0.4324056782064047,\n", + " ParameterVectorElement(θ[6]): -3.4016515185029284,\n", + " ParameterVectorElement(θ[11]): 6.226143981483694,\n", + " ParameterVectorElement(θ[3]): -6.218365322324229,\n", + " ParameterVectorElement(θ[13]): -4.321134059757356,\n", + " ParameterVectorElement(θ[15]): 1.7460062016947542,\n", + " ParameterVectorElement(θ[10]): -3.144056546098471,\n", + " ParameterVectorElement(θ[0]): 7.846190848261606,\n", + " ParameterVectorElement(θ[14]): -3.4620044347295287,\n", + " ParameterVectorElement(θ[2]): 0.07498512390436408,\n", + " ParameterVectorElement(θ[12]): 5.114535339252314,\n", + " ParameterVectorElement(θ[5]): 7.495171057359882,\n", + " ParameterVectorElement(θ[9]): -6.273228533679424,\n", + " ParameterVectorElement(θ[8]): -0.015414042392154728},\n", + " 'optimal_point': array([ 7.84619085, 0.43240568, 0.07498512, -6.21836532, 2.25050933,\n", + " 7.49517106, -3.40165152, -6.39595285, -0.01541404, -6.27322853,\n", + " -3.14405655, 6.22614398, 5.11453534, -4.32113406, -3.46200443,\n", + " 1.7460062 ]),\n", + " 'optimal_value': -2.23828125,\n", + " 'optimizer_evals': None,\n", + " 'optimizer_time': 8.168923377990723}\n" + ] + } + ], + "source": [ + "print(vqe_result)" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "id": "895a21f4", + "metadata": {}, + "outputs": [ + { + "ename": "TypeError", + "evalue": "Can't instantiate abstract class Eigensolver with abstract methods compute_eigenvalues", + "output_type": "error", + "traceback": [ + "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", + "\u001b[1;31mTypeError\u001b[0m Traceback (most recent call last)", + "\u001b[1;32m\u001b[0m in \u001b[0;36m\u001b[1;34m\u001b[0m\n\u001b[0;32m 1\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0mqiskit\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0malgorithms\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mEigensolver\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[0mes\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mEigensolver\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 4\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[0mes\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcompute_eigenvalues\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mH\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;31mTypeError\u001b[0m: Can't instantiate abstract class Eigensolver with abstract methods compute_eigenvalues" + ] + } + ], + "source": [ + "from qiskit.algorithms import Eigensolver\n", + "\n", + "es = Eigensolver()\n", + "\n", + "es.compute_eigenvalues(H)" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.8.11" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}