diff --git a/benchmark/benchmark.py b/benchmark/benchmark.py index cf68c4c..53549f1 100644 --- a/benchmark/benchmark.py +++ b/benchmark/benchmark.py @@ -133,9 +133,7 @@ def source(request): @pytest.fixture(params=("overhead", "zarr", "tiff", "hdf5")) def file_type(request, source): - if request.param == "overhead": - filename, fs = source("1-byte") class Overhead(Fixture): @@ -149,7 +147,6 @@ def run(self, chunk_index): return Overhead(source) elif request.param == "zarr": - filename, fs = source(f"{BASE}.ome.zarr") class Zarr(Fixture): @@ -167,7 +164,6 @@ def run(self, chunk_index): return Zarr(source) elif request.param == "tiff": - filename, fs = source(f"{BASE}.ome.tiff") class TIFF(Fixture): @@ -189,7 +185,6 @@ def run(self, chunk_index): return TIFF(source) elif request.param == "hdf5": - filename, fs = source(f"{BASE}.ims") class HDF5(Fixture): @@ -208,12 +203,10 @@ def run(self, chunk_index): return HDF5(source) else: - raise NotImplementedError(request.param) def test_benchmark(file_type, chunk_choice, record_property): - chunk_index, chunk_meta = chunk_choice chunk_distance, chunk_record = chunk_meta diff --git a/notebooks/NumberOfChunks-100k.ipynb b/notebooks/NumberOfChunks-100k.ipynb index f942666..cf6bf4d 100644 --- a/notebooks/NumberOfChunks-100k.ipynb +++ b/notebooks/NumberOfChunks-100k.ipynb @@ -1,1059 +1,1059 @@ { - "cells": [ - { - "cell_type": "code", - "execution_count": 1, - "metadata": { - "scrolled": true - }, - "outputs": [ + "cells": [ { - "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_dpi_ratio', { dpi_ratio: fig.ratio });\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", - " // 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", - " var cursor = msg['cursor'];\n", - " switch (cursor) {\n", - " case 0:\n", - " cursor = 'pointer';\n", - " break;\n", - " case 1:\n", - " cursor = 'default';\n", - " break;\n", - " case 2:\n", - " cursor = 'crosshair';\n", - " break;\n", - " case 3:\n", - " cursor = 'move';\n", - " break;\n", - " }\n", - " fig.rubberband_canvas.style.cursor = 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", - " /* 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", - " evt.data.type = 'image/png';\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", - " evt.data\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 http://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", - " * http://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.which === this._key) {\n", - " return;\n", - " } else {\n", - " this._key = event.which;\n", - " }\n", - " }\n", - " if (name === 'key_release') {\n", - " this._key = null;\n", - " }\n", - "\n", - " var value = '';\n", - " if (event.ctrlKey && event.which !== 17) {\n", - " value += 'ctrl+';\n", - " }\n", - " if (event.altKey && event.which !== 18) {\n", - " value += 'alt+';\n", - " }\n", - " if (event.shiftKey && event.which !== 16) {\n", - " value += 'shift+';\n", - " }\n", - "\n", - " value += 'k';\n", - " value += event.which.toString();\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, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n", - "\n", - "mpl.extensions = [\"eps\", \"jpeg\", \"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.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", - " // Pass the mpl event to the overridden (by mpl) onmessage function.\n", - " ws.onmessage(msg['content']['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", - " var manager = IPython.notebook.keyboard_manager;\n", - " if (!manager) {\n", - " manager = IPython.keyboard_manager;\n", - " }\n", - "\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" + "cell_type": "code", + "execution_count": 1, + "metadata": { + "scrolled": true + }, + "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 \u2265 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_dpi_ratio', { dpi_ratio: fig.ratio });\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", + " // 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", + " var cursor = msg['cursor'];\n", + " switch (cursor) {\n", + " case 0:\n", + " cursor = 'pointer';\n", + " break;\n", + " case 1:\n", + " cursor = 'default';\n", + " break;\n", + " case 2:\n", + " cursor = 'crosshair';\n", + " break;\n", + " case 3:\n", + " cursor = 'move';\n", + " break;\n", + " }\n", + " fig.rubberband_canvas.style.cursor = 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", + " /* 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", + " evt.data.type = 'image/png';\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", + " evt.data\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 http://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", + " * http://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.which === this._key) {\n", + " return;\n", + " } else {\n", + " this._key = event.which;\n", + " }\n", + " }\n", + " if (name === 'key_release') {\n", + " this._key = null;\n", + " }\n", + "\n", + " var value = '';\n", + " if (event.ctrlKey && event.which !== 17) {\n", + " value += 'ctrl+';\n", + " }\n", + " if (event.altKey && event.which !== 18) {\n", + " value += 'alt+';\n", + " }\n", + " if (event.shiftKey && event.which !== 16) {\n", + " value += 'shift+';\n", + " }\n", + "\n", + " value += 'k';\n", + " value += event.which.toString();\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, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n", + "\n", + "mpl.extensions = [\"eps\", \"jpeg\", \"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.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", + " // Pass the mpl event to the overridden (by mpl) onmessage function.\n", + " ws.onmessage(msg['content']['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", + " var manager = IPython.notebook.keyboard_manager;\n", + " if (!manager) {\n", + " manager = IPython.keyboard_manager;\n", + " }\n", + "\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" + } ], - "text/plain": [ - "" + "source": [ + "%matplotlib notebook\n", + "from ipywidgets import *\n", + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "import math\n", + "\n", + "x = 2**15\n", + "y = 2**15\n", + "z = 1\n", + "t = 1\n", + "c = 25\n", + "chunkSizes = [32, 1024]\n", + "chunkT = 1\n", + "chunkZ = 1\n", + "\n", + "fig = plt.figure()\n", + "ax = fig.add_subplot(1, 1, 1)\n", + "plt.xlabel('Chunk size (X and Y)')\n", + "plt.ylabel('Number of chunks')\n", + "plt.title(\"XYZCT: 32k x 32k x 1 x 25 x 1\")\n", + "\n", + "def file_count(i, chunkC=1):\n", + " return math.ceil(x/i) * math.ceil(y/i) * math.ceil(z/1) * math.ceil(t/1) * math.ceil(c/chunkC)\n", + "\n", + "\n", + "for chunkC in (1, 5, 25):\n", + " numFiles = []\n", + " fileSize = []\n", + " for i in chunkSizes:\n", + " numFiles.__iadd__([file_count(i, chunkC)])\n", + " fileSize.__iadd__([i])\n", + " ax.plot(fileSize, numFiles, linewidth=4.0, label=f\"{chunkC}\")\n", + "\n", + "plt.legend(loc=\"lower left\", title=\"Chunk size (C)\", frameon=False)\n", + "\n", + "def scale(scale):\n", + " plt.yscale(scale)\n", + " plt.xscale(scale)\n", + "\n", + "scale('log')\n", + "ax.annotate('Chosen chunk size\\nC: 1, X: 256, Y: 256',\n", + " xy=((256), file_count(256)), xycoords='data',\n", + " xytext=(0, 40), textcoords='offset points',\n", + " arrowprops=dict(facecolor='black', shrink=0.05),\n", + " horizontalalignment='left', verticalalignment='bottom')\n", + "fig.savefig(\"/tmp/2d-chunks.png\", dpi=800)\n" ] - }, - "metadata": {}, - "output_type": "display_data" }, { - "data": { - "text/html": [ - "" - ], - "text/plain": [ - "" - ] - }, - "metadata": {}, - "output_type": "display_data" + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "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.7.6" } - ], - "source": [ - "%matplotlib notebook\n", - "from ipywidgets import *\n", - "import numpy as np\n", - "import matplotlib.pyplot as plt\n", - "import math\n", - "\n", - "x = 2**15\n", - "y = 2**15\n", - "z = 1\n", - "t = 1\n", - "c = 25\n", - "chunkSizes = [32, 1024]\n", - "chunkT = 1\n", - "chunkZ = 1\n", - "\n", - "fig = plt.figure()\n", - "ax = fig.add_subplot(1, 1, 1)\n", - "plt.xlabel('Chunk size (X and Y)')\n", - "plt.ylabel('Number of chunks')\n", - "plt.title(\"XYZCT: 32k x 32k x 1 x 25 x 1\")\n", - "\n", - "def file_count(i, chunkC=1):\n", - " return math.ceil(x/i) * math.ceil(y/i) * math.ceil(z/1) * math.ceil(t/1) * math.ceil(c/chunkC)\n", - "\n", - "\n", - "for chunkC in (1, 5, 25):\n", - " numFiles = []\n", - " fileSize = []\n", - " for i in chunkSizes:\n", - " numFiles.__iadd__([file_count(i, chunkC)])\n", - " fileSize.__iadd__([i])\n", - " ax.plot(fileSize, numFiles, linewidth=4.0, label=f\"{chunkC}\")\n", - "\n", - "plt.legend(loc=\"lower left\", title=\"Chunk size (C)\", frameon=False)\n", - "\n", - "def scale(scale):\n", - " plt.yscale(scale)\n", - " plt.xscale(scale)\n", - "\n", - "scale('log')\n", - "ax.annotate('Chosen chunk size\\nC: 1, X: 256, Y: 256',\n", - " xy=((256), file_count(256)), xycoords='data',\n", - " xytext=(0, 40), textcoords='offset points',\n", - " arrowprops=dict(facecolor='black', shrink=0.05),\n", - " horizontalalignment='left', verticalalignment='bottom')\n", - "fig.savefig(\"/tmp/2d-chunks.png\", dpi=800)\n" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [] - } - ], - "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.7.6" - } - }, - "nbformat": 4, - "nbformat_minor": 2 + "nbformat": 4, + "nbformat_minor": 2 } diff --git a/notebooks/NumberOfChunks-3d.ipynb b/notebooks/NumberOfChunks-3d.ipynb index 95939eb..a2d558b 100644 --- a/notebooks/NumberOfChunks-3d.ipynb +++ b/notebooks/NumberOfChunks-3d.ipynb @@ -1,96 +1,96 @@ { - "cells": [ - { - "cell_type": "code", - "execution_count": 1, - "metadata": {}, - "outputs": [ + "cells": [ { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEaCAYAAAD+E0veAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA+uklEQVR4nO3ddXxWdRvH8c+1okY+MGLECOkY3aWMbhFB9KEEQVDExKIkVJAURRAE9AFEulNakO5OGTmaAdvY9nv+uG/mGNvYWJzd2/V+vc7Lnf7ex3Ff+534HTHGoJRSSgE4WR1AKaVU0qFFQSmlVBgtCkoppcJoUVBKKRVGi4JSSqkwWhSUUkqF0aKgkg0RqSMivhbsd4OIvJnI+7Tks6rkT4tCCiQi7iJyTkReCzctvYj8IyJtRKSbiBwVkVTh5v9HRK6JSEMRWSEi/hGGABExIpLXvrybiAwUkZMict++v6ki4iUih8OtF2Jf9/H4Z8/I3ltEdolIoIhMS7CD9OQ+c4rIYhG5ZP+MXom032T1WZ/1eUTkJRE5JiIPRGS9iOQLN09E5BsRuWEfvhURia9s6l9aFFIgY4w/0B0YKyLZ7JO/BXYZY+YaYyYDvkD/cKuNAZYbY1YaYxoZY9wfD0BGYDswwxjzj335uUBz4DX7/DLAbuAlY0yJcOtuBnqH296wZ8S/BAwBpsbtKMRKKLASeDkR9wnJ77NG+XlEJCswH/gSyALsAn4Pt0h3oCW236PSQFPgrQTIqIwxOqTQAZgGzALqADeAnOHmeQG3AG+gPrZ/0Jmj2M43wAEgrX28HvAQyBODDBuAN58j+xBgWoRpdQDfcOPvAkeA3JGs/yMwN8JnWAdINPt0AQzgFdVnAHLaj8WHkaxfELgJlLOP5wKuA3Uc5bMm4OfpDvwVbjyd/XeoqH38L6B7uPldge2J/W8mJQwuT1UJlZL0xfZF4oPtS+zy4xnGmHMi0h/bX3WZgbeNMbcibkBEWmD7i62iMeaBfXI9YIcx5kJsA4lIP6CGMaZprD/Nk9v5EmgF1DbG+EWyyAfAPhHpBJzG9iXjbezfOM+5Ty9gNTDSGDMp4nxjzGkR+QT4n4iUB37B9uW44Xn3ad9von9WiPfPUwLYH27b90XktH36sYjz7T+XeN7sKmpaFFIwY8wtETkMVMPWdI/oe+B1YJ8xZmHEmSJSEFtro6sx5mS4Wf8BLkdcPoaZvn6e9Z6MJaOASkBdY8ydKPbzQERex3aq5B7wjjEmLhduiwNfAJ8aY2ZFtZAxZrKINAP+xvaXePM47NOqzxp+2/H1edyBiAXtDpA+3Pw7Eea5i4jEtbipJ+k1hRTM/kXhBazFdkrhCfZ/bEeBw5GsmxrbdYOpxph5EWbfwHYaxQqZsJ2KGB7Vl+RjxpgdwBlAgDlx3G8H4CK2Y/Isk4GSwHhjTGAc9pkJaz5rRPHxefyBDBGmZcBWxCKbnwHw14IQ/7QopFAi4gGMBrphO/3TVkRqxWITE4D7wCeRzFsLVBKR3HEOGnu3sF2E/EVEqke3oIj0AlJhu17ycRz3OxDb+fSZIuIczT7dsV20nwIMFJEscdinVZ81/Hbj6/McxnYR+fF202G7ZnE4svn2n5/6Y0XFnRaFlOt7YKExZr39WsLHwOTwt6FGRUS6YPsyamuMCY443xizFlgDLBCR8iLiYr/ltYd93edm31ZqwBlwFpHUIvLEaVD7Oe0O9v1XjmI7hbFd8HwdeAP4WES8o9lvamxfqgCp7OPhPQJewXaB9FcRierf1lhgtzHmTWAZMNEBP2tCfJ4FQEkRedm+TH/ggDHmmH3+DOB9EfEUkVzYrpNMiyaXel5WX+nWIfEHbLf2XQIyRZi+DhgaYdo0YEiEaWewfQn6RzLUtC/jBgwCTmFrUZwHfgbyRtjWBsLdfQR8BqyIJvtAbOeuww8D7fPq8OQdOU2Aq0D5CNtwAXYA/cJN6wkcBFJFsd+I+zSRfQYgNbaW0jTAKcI2WmA7xZTFPu5uPz4dHOWzJtTnsc+vh+2i8kP7MfUKN0+w3TZ90z58SzR3T+nw/IPYD7hSSimlp4+UUkr9S4uCUkqpMFoUlFJKhdGioJRSKowWBaWUUmEcupuLrFmzGi8vL6tjKKWUQ9m9e/d1Y0y2yOY5dFHw8vJi165dVsdQSimHIiLno5qnp4+UUkqF0aKglFIqjBYFpZRSYbQoKKWUCqNFQSmlVBiHLAoi0kxEJt25E+17RZRSSsWSQxYFY8wSY0z3jBkzPtf6fvcCWXf0ajynUkopx+eQRSGuJm06Tdfpu3hv9l5u3Q+yOo5SSiUZDv3w2vP6qEFR0rq5MGH9KbacusGQliVoWNKqVworpVTSkSJbCm4uTvT1Kczi3jXIniEVPX7bQ6+Ze7jhH5d3qCullONLkUXhseK5MrCwV3U+rF+Y1Yev4DN6E0v2X0LfRqeUSqlSdFEAcHV2oveLL7Ds3ZrkyZyGd2btpcdvu7l2L8DqaEoplehSfFF4rHD29MzrWY1+jYqy/rgf9UdvYsFeX201KKVSFC0K4bg4O9GjdkGWv1uTAlnT0ff3/bw5fRdX7mirQSmVMmhRiEQhD3f+6FGNL5oUY+vp6/iM3sicXRe01aCUSva0KETB2Ul4s2YBVvapRbGcGfh47gE6/rKTi7cfWh1NKaUSjBaFZ/DKmo7Z3aowqHkJdp27SYPRm/jf3+e11aCUSpa0KMSAk5PQsZoXq96rRencGfl8wSFen/I3F24+sDqaUkrFKy0KsZAnS1r+92ZlhrUqxf4Ld2gwZhMztp0jNFRbDUqp5CHJFAURcRKRoSIyXkQ6Wp0nKiLCa5XzsqpvLcrny0z/RYdpP3k7567ftzqaUkrFWYIWBRGZKiLXRORQhOkNReS4iJwSkX72yS0AT+AR4JuQueKDZ6Y0zOhSiW9fLs2Ry3dpOHYTU7acJURbDUopB5bQLYVpQMPwE0TEGZgANAKKA+1FpDhQBNhmjHkf6JnAueKFiNC2Yh7W9K1NtYJZ+WrpEdr+tI3Tfv5WR1NKqeeSoEXBGLMJuBlhciXglDHmjDEmCJiNrZXgC9yyLxMS1TZFpLuI7BKRXX5+fgkRO9ZyZEzNlI4VGP1qGU5d86fR2M1M3Hia4JBQq6MppVSsWHFNwRO4EG7c1z5tPtBARMYDm6Ja2RgzyRhTwRhTIVu2bAmbNBZEhFZlc7Pm/VrULZKNr1cc4+Uf/+LE1XtWR1NKqRizoihIJNOMMeaBMaarMeYdY8yERE8VTzzSp2bi6+UZ374sF249pOm4LXz/50keaatBKeUArCgKvkCecOO5gUux2UBSf0eziNCsTC7W9K2FT4nsjFx9gpYTtnLk0l2roymlVLSsKAo7gRdEJL+IuAHtgMWx2UBc39GcWP7jnooJr5Vj4uvluHo3gObfb2H0mhMEBWurQSmVNCX0LamzgG1AERHxFZGuxphgoDewCjgKzDHGHI7ldpN0SyGihiVzsqZvbZqWzsnYdSdp/v0WDl10jOxKqZRFHLkPnwoVKphdu3ZZHSNW1hy5yucLDnLjfhA9ahfg3ZdeIJWLs9WxlFIpiIjsNsZUiGxeknmiOaXwKZ6dNX1r06qsJxPWn6bpuC3su3Db6lhKKQVoUbBExrSujHylDL90roh/YDCtf9jK8OVHCXgU5eMZSimVKByyKDjaNYWo1C3iwaq+tXi1Yh5+2nSGxmM3s+tcxGf9lFIq8eg1hSRiy8nrfDLvAJfuPKRztfx82KAwad1crI6llEqG9JqCA6jxQlZW9a3F65XzMXXrWRqN3cz2MzesjqWUSmEcsigkl9NHEbmncuGrliWZ1a0KxkC7Sdvpv+gQ9wODrY6mlEoh9PRREvUgKJgRq44z7a9zeGZKwzcvl6Z6oaxWx1JKJQN6+sgBpXVzYUCzEvzxVlXcnJ3o8PPffDr/IPcCHj254L0rEBxoTUilVLKjRSGJq+CVheV9atK9VgF+3/kP9UdvYsPxa7aZxsCcjjCxJvgmzxaTUipxaVFwAKldnfmscTHm9ayGeyoXOv2ykw//2M+dh8FQ6yMIug9TfGD1F/DoodVxlVIOzCGLQnK90PwsZfNmZum7NehVtyAL9l7EZ/RG1j4qBW9vg3Id4a/xMLEG/LPd6qhKKQelF5od1EHfO3w0dz/HrtyjVVlPBjQrTqYrf8Hid+D2BajSE178AtzSWR1VKZXE6IXmZKhU7ows7l2DPi+9wJL9l6g3ahMrHxSFntug4puw/Qf4sTqc22J1VKWUA9Gi4MDcXJzo61OYRb2rkz1DKnr8tpte805wo/ZQ6LTMttC0JrDsAwj0tzasUsohaFFIBkrkysjCXtX5sH5hVh++gs/oTSy5UwDTcytU6QU7p8APVeH0equjKqWSOIcsCin1QnN0XJ2d6P3iCyx7tyZ5MqfhnVl76fH7Ua5V7w9dVoFLKvi1JSx+FwL0uCmlIqcXmpOh4JBQft5yllFrTpDWzZmBzUrQokRmZOPXtjuU0ueEZuPghXpWR1VKWUAvNKcwLs5O9KhdkOXv1qRA1nS89/s+us06zNXKn0HXtZAqPfzvZVj4Njy8ZXVcpVQSokUhGSvk4c4fParxRZNibDl1nXqjNjLnigem+0ao+QHsnw0TqsDxFVZHVUolEVoUkjlnJ+HNmgVY2acWxXJm4OO5B+j06wEulf8Iuv0J6bLCrHYwrxs80Bf8KJXSaVFIIbyypmN2tyoMal6CneduUn/0JmZeyILp9ifU+RQOz4cJleHIYqujKqUspEUhBXFyEjpW82LVe7UonTsjny04yOvT9nKh9LvQfSNkyAlz3oA/OsH961bHVUpZwCGLgt6SGjd5sqTlf29WZmirkuy/cIcGYzYx46w7oV3WwotfwrFlMKESHJpn64lVKZVi6C2pKdzF2w/pN+8Am09ep3L+LHzbpjT5Qv6BRb3g4m4o2hSajIL02a2OqpSKJ3pLqoqSZ6Y0zOhSiW9fLs2Ry3dpMGYTU06kJqTzKvAZDCfX2FoN+2drq0GpFECLgkJEaFsxD2v61qZawax8tfQIbSfv5HThrtBzK2QrAgvest2ldPeS1XGVUglIi4IKkyNjaqZ0rMCotmU4dc2fxmM389NhJ0I6LocGw+HMRttzDXt+1VaDUsmUFgX1BBGhdbncrOlbi9qFszF8xTFa//Q3Jwu8YWs15CgJi3vDby/b3tuglEpWtCioSHlkSM1Pb5RnXPuy/HPjPk3GbWHCAUPwG4uh8Ujb291+qAq7pmqrQalkRIuCipKI0LxMLta8Xxuf4tkZseo4LX/cxtE8r8Lbf4FnOVjaF2Y0h1vnrI6rlIoHWhTUM2V1T8WEDuX4sUM5rtwJoPn3WxizO4ig1xZA0zFwcS/8UA3+ngShoVbHVUrFgRYFFWONSuVkTd/aNC6VkzFrT9J8wlYO5WwNb2+DvFVgxUcwvSncOG11VKXUc3LIoqBPNFsnczo3xrYry+T/VuDm/SBaTNjKyL8fENhuDrSYAFcO2d4NvW0ChIZYHVcpFUv6RLN6bncePGLw0iPM2+NL4ezufNumDN4ZH9iuM5xYCbkr2QpFtsJWR1VKhROnJ5pF5FsRySAiriKyTkSui8jr8R9TOZqMaV35rm0ZfulckXsBwbT+YSvDt94hoM3/oNUkuH4CJtaALWMgJNjquEqpGIjJ6aP6xpi7QFPAFygMfJSgqZRDqVvEg1V9a/FqxTz8tPEMjcdvYXcmH+i1A17wgbUDYIoPXDtqdVSl1DPEpCi42v/bGJhljNE3sainZEjtyvDWpfm1ayUCH4XSZuI2vtp4g4etpkObX+D2efipFmwaASGPrI6rlIpCTIrCEhE5BlQA1olINiAgYWMpR1XzhWys6luLDpXzMmXLWRqN28zfaWvbWg1Fm8CfQ2Dyi3DloNVRlVKRiNGFZhHJDNw1xoSISDrA3RhzNcHTPYNeaE7a/jp9nX7zDvLPzQd0rJqPjxsWJd3p5bDsfXh4C2p+aHtXtIub1VGVSlHieqF5qjHmljHm8f2FAsyMz4AqeapWMCsr36tJ5+pezNh+ngZjNrHVrZqt1VCiNWz8GibXhUt7rY6qlLKLyemjiyLyI4S1GFYDvyVoKpVspHVzYUCzEsx5qyquzk50+PlvPlt1kXtNfoD2s+HBDZj8EqwdBMGBVsdVKsV7ZlEwxnwJ3BWRidgKwnfGmF8SPJlKVip6ZWFFn5p0r1WA2Tv+ocHoTWyUCvD2dijTHrNlFEN+qcLBw39YHVWpFC3Kawoi0jr8KPAlsANYCWCMmZ/g6Z5Bryk4pr3/3OKjuQc4dc2fthVy83mT4jw8vYjXdwzEzwk6ur/A202mkjpNZqujKpUsRXdNIbqiEF1rwBhjusRHuLjQouC4Ah6FMHbdSSZtOkNWdzeGtSpFJc8QvlvemXmBl/AKga8q9sO7VAeroyqV7DxXUXAEWhQc3wHf23z0xwGOX71Hq7KeDGhWnKNHpzJw/wQuOxk6pCvAO02mkjZtVqujKpVsxKko2J9L6AZ4AS6Pp8d3S0FE6gBfAYeB2caYDc9aR4tC8hAUHMr360/xw/pTZErrxpCWJanpBaOXdeH3gAvkCYFB5fpS0dvyxqlSyUKcbkkFFgEZgbXAsnBDTHY8VUSuicihCNMbishxETklIv3skw3gD6TG1p2GSiHcXJx436cwi3pXxyN9Knr8tptPFl+iZ5MFTC3TFwN02T+aob835oH/NavjKpWsxaSlsM8Y4/1cGxephe2LfoYxpqR9mjNwAvDB9uW/E2gPHDPGhIpIdmCUMeaZJ5O1pZD8PAoJZeKG04z78yQZUrsyuEVJ6hZyZfyyLvzv/hlyhQoDy/SmSvm3rI6qlMOKa0thqYg0fp4dG2M2ARH7SqoEnDLGnDHGBAGzgRbGmMev7LoFpHqe/SnH5+rsxDsvvcCSd2rgmTkNvWbu4YP55+nS8A+ml+uHK9Dt0PcMmt0A/3uXrY6rVLITk6LQB1theCgid0XknojcjcM+PYEL4cZ9AU8RaS0iPwG/At9HtbKIdBeRXSKyy8/PLw4xVFJWNEcG5vesxicNi7Lu2DV8Rm/kfEhd5rTbSKd0hZgfcJGWf/iwZcd4q6Mqlawk+N1HIuIFLA13+ugVoIEx5k37+BtAJWPMO7Hdtp4+ShlOXbvHR3MPsPef29Qr5sHQVqW4emER/f8eymlnQwvX7HzUZCoZM+a1OqpSDiGufR/VimyIQx5fIE+48dzApdhsQF/HmbIU8kjP3B7V+KJJMTafvI7PqI0cf1iN31/bQrf0RVkadIVW8xqzYft3VkdVyuHF5ELzknCjqbFdE9htjHkxRjt4uqXggu1C80vARWwXml8zxhyObXhtKaQ8Z6/f55O5B9hx7iZ1imRjWKtS3Lqygi+3DeKkUyhNXLLSr/FUMmXOb3VUpZKsOLUUjDHNwg0+QEkgRt1mi8gsYBtQRER8RaSrMSYY6A2sAo4Cc56nIKiUKX/WdMzuXoWBzYrz95mb1B+9iQN3KzD7ta30yFCSVY/8aLmgGeu2Drc6qlIOKdbXFEREgAPGmFIJEylGGZoBzQoVKtTt5MmTVsVQFvvnxgM+mXeAbWduUKNQVoa3LsX96+vov/VLjjqF0MA5M581nkqWLIWsjqpUkhLXJ5rHY3uwDGwtC2/gnDHm9fgM+Tz09JEKDTXM3PEPw5fb3v/cr1FR2pbLxvRVPfjx5l7SG/is4Cs0qPEl4hSTm+2USv7iWhQ6hhsNxlYQtsZjvuemRUE95nvrAZ/OP8jmk9epUiAL37xcmke3t/Llpk845BTMS04Z+KLhz2TNVszqqEpZLtl1iKenj1RkjDHM2XWBIUuPEhxq+KhBEV6vlIPfVvViwo0dpDHwiVdLmtYerK0GlaLFtaVQHRgI5MPWIZ5g6zq7QDznjDVtKajIXL7zkM/mH2T9cT8q5MvMt21Kg/8O+m/4kP1Oj6gt7vRvMAmP7JZdFlPKUnEtCseAvsBu4PF7mjHG3IjPkM9Di4KKijGG+XsuMmjJYQKDQ/mgfmE6VcnN7LXvMP7qX7gCH+VtTMu6X2urQaU4cS0KfxtjKidIsjjSoqCe5drdAD5feIg1R67inScTI9qUxi1gH/3X92WPBFKdtAysP5EcOctaHVWpRPO8b14rZ/+xLeAMzAfC3qxujNkTzzljTK8pqNgwxrB4/yUGLj7M/cAQ+tR7gW418jJ3XV/GXN6IE/CBpw9t6n2nrQaVIjxvUVgfzTZNTJ9oTkjaUlCx4XcvkAGLD7H84BVKeWZkxCulcQ8+wsB17/A3AVQmNYPqTcDTs5LVUZVKUMnu7qPHtCio57H84GW+XHiIuwGP6F33BXrWzs/C9R8y6uJaQoG+Oevwqs8YnJxdnrktpRxRXDvEGyYimcKNZxaRIfGYT6lE1bhUTta8X5tGJXMyeu0JWvywjeIlBrGgwXTKOqVl2NWNdP21Chf+2WJ1VKUSXUxOoDYyxtx+PGKMuQU810t34ov2kqriKks6N8a1L8ukN8pz3T+QFhO28r+D7oxtv5XBng05bgJova4Hvy5/i5DgIKvjKpVoYnL30QGgojEm0D6eBthljCmRCPmipaePVHy4/SCIwUuPMH/PRQpnd2dEmzLkcPuHwaveYpPxx9u4MrjOKPJ71bE6qlLxIq6v4/wNWCciXUWkC7AGmB6fAZWyUqa0boxq683UThW4+zCYVj9s5Zc9box8dRPD8jbnjAnilfW9+WVpV201qGQvRheaRaQhUA/b08yrjTGrEjpYTGhLQcW3uwGPGLr0KL/vukDBbOn4tk0Z8qa5xJBV3fkz9C6lQl34qta3FCzoY3VUpZ6b3n2kVCxtOuHHp/MPcunOQ7pWz8/79V5g499DGXZmLvcFemYpR6dGE3F1TWt1VKViLa6nj5RKcWoVzsbK92ryWqW8/LzlLI3HbyFLnndY2Gw+dV0yM+7WXjr8Vo3jJ5dZHVWpeOWQRUHvPlKJIX1qV4a2KsXMNysTYgyvTtrOuK1BDH5lPaMKtucqwbTb+gk/LnyNR4H3rY6rVLyIsiiIyDr7f79JvDgxY4xZYozpnjFjRqujqBSgWqGsrOxTi07VvJi+7TwNx24iXY7uLGy5mPquWfnhzkHazazGkWMLrY6qVJxF11LIKSK1geYiUlZEyoUfEiugUklBulQuDGxegjlvVcVZhNd+/psRG+7zxctrGFe4I7dMCK9t/4Jx89sSFHjP6rhKPbfo+j5qA3QFagARr+Zq30cqxXoYFMJ3q48zZetZcmVMw/DWpSiT/SHfLuvM4kfXKBQMg/O8Tan6Pa2OqlSk4tp19pfGmK8SJFkcaVFQVtp9/hYfz93Pab/7tK2Qm8+bFOfAwR9ZNO8nOi8O5Vz1LNQdMY+0mXNYHVWpJ8T5llQRaQ7Uso9uMMYsjcd8z02LgrJawKMQxqw9yaRNp/FIn5phrUtSPsNtNn/YgYJ77nMtC2R6vyNl2vSzOqpSYeLaId5woA9wxD70sU9TKsVL7epMv0ZFWfB2dTKkcaHLtF0M2vyQGj//hf+nrXAJBpcvprOsc1Xu37hkdVylnimmfR95G2NC7ePOwF5jTOlEyBdVJn3JjkpyAoND+P7PU/yw4TRZ0rkxtWNFvFxvsPHDdhTceY/rmSDdu+0p91p/q6OqFC4+Hl7LFO5ny+8D1VtSVVKUysWZD+oXYVGv6pTyzEi+rGlJn70ATX/dwYMv2wKQZvAslv63Mv7XzlkbVqkoxKSl0B74GliPre+jWsCnxpjZCR8venpNQTkS/xsX2PDBK+TffodbGSHV262p2HGo1bFUChSnloIxZhZQBds7mucDVZNCQVDK0bj/Jw9Np23n0aDXCXUC9+HzWdqhIncvn7I6mlJhtEM8pSzw4NYV1n/0Ml5bb3LbHVx7NKNS12+tjqVSCO0QT6kkJm3mHDT5eSuhw97kkRukH7GEpe3Kc9v3qNXRVAqnRUEpC5Vq9QGVl2/kTO2seO1/wNGWrdk2sY/VsVQKFm1REBEnETmUWGGUSonSZPSgyU+bkRE9eZhGyDRmNUtfKcut8wetjqZSoGiLgv3ZhP0ikjeR8iiVYpVo+i41VvzF2Rezk+9wAMdbt2XreO0/SSWumJw+ygkcFpF1IrL48ZDQwZRKidzcM9H4hw24fNeHB+mELBM2sKy1N9dP77Y6mkohYvKcQu3IphtjNiZIohjQJ5pVShD04A7r+rUm99pLPEgNga9Xp0afn3BydrY6mnJwcX1OYSNwDnC1/7wT2BOvCWNJn2hWKYFb2ow0GreOVOM+5l4GIdukraxsXY5rJ7ZbHU0lYzHpEK8bMBf4yT7JE1iYgJmUUuEUrteZ2it2cL5xHjxPBXG+bWc2fNOR0JAQq6OpZCgm1xR6AdWBuwDGmJOAR0KGUko9yTWNOw1HrSbdj19wO4uQ/ZcdrGxZliuHN1kdTSUzMSkKgcaYoMcjIuICOO5j0Eo5sEK1OlB3xW7+aZ4fz7OPuPjaW/w5pL22GlS8iUlR2CginwFpRMQH+ANYkrCxlFJRcUmVhgbfLifDpK+4kc2JnL/tY1Uzby7tX2t1NJUMxKQo9AP8gIPAW8By4IuEDKWUerb81drw0oq9XGj1AjkvBHPljXdYN/AVbTWoOInJ3UehwHTgK2AQMN04ci96SiUjzq5u1B++mMxTvsYvhxO5Zh9idZMyXNi9zOpoykHF5O6jJsBpYBzwPXBKRBoldDClVMzlq9SCesv2cvGVYnhcCuF6pw9Z83kLQh4FPXtlpcKJyemj74C6xpg6xpjaQF1gdMLGUkrFlrOrG/W+mo/HtNFc9XQm97wTrG1clvN/L7A6mnIgMSkK14wx4d8Ccga4lkB5lFJxlLtcQ+ov28+l10qR9Woot7p+xqpPmhIcGGB1NOUAoiwKItJaRFpj6/douYh0EpGO2O482ploCZVSsebk7MxL/eeQ87fvuZzPhbyLTvNn4/Kc3vy71dFUEhddS6GZfUgNXAVqA3Ww3YmUOcGTKaXiLFfpl2iweB9X/lueLNdD8e85kJUfNCQ48KHV0VQSlaRexyki6YBNwABjzNJnLa+v41Qq5q4e2cLufm+T/8QjLuYUvPr3o1Dd/1odS1kgTh3iiUh+ERklIvNj23W2iEwVkWsRX9QjIg1F5LiInBKRfuFmfQLMicm2lVKxk714DRou2Mu1rlXIeNvwoPdwVrxXj0cP/a2OppKQmFxoXoitl9Tx2O5EejzExDSgYfgJIuIMTAAaAcWB9iJSXETqAUewnapSSiUAJ2dnan/0C/l+/4ULL7jhtfIimxpW4tjqyVZHU0lETN6n8LcxpvJz70DEC1hqjClpH68KDDTGNLCPf2pf1B1Ih61QPARa2R+ci7i97kB3gLx585Y/f/7880ZTKsXbPKY7bjM2kzYALryYk3rfLMAtnXZJn9zF6fQRMFZEBohIVREp93iIQx5P4EK4cV/A0xjzuTHmPWAmMDmyggBgjJlkjKlgjKmQLVu2OMRQStV8bxKF5v6P88VSkX/tZbY0rMqR5ROsjqUsFJOiUAroBnzNv6eORsZhnxLJtLDmijFmWkwuMiul4sd/CpSj6bx93HznRdI+MIR+8D3LetQi8N4Nq6MpC8SkKLQCChhjahtj6tqHF+OwT18gT7jx3MCl2GxARJqJyKQ7d+7EIYZSKrzqvSZQZP4fnCuVhgIb/NjWqAaHFo2xOpZKZDEpCvuBTPG4z53AC/a7mtyAdkCM7mZ6LLav47xy5Qrt2rWjYMGCFC9enMaNG3PixAk2bNhA06ZNY/8JInHu3DlKliz5XOs2btyY27dvx0uOiBYuXMjgwYMBaNCgAd7e3mFDrly5qFz538tF27ZtI3/+/GHz3d3dKVKkCN7e3vz3v//l4MGDdOrUKUFyqqQhc76SNJ2zhzt9G+IWANLvJ5Z1q8HDu35WR1OJxRgT7QBsAG4Cq7B9eS8GFj9rPfu6s4DLwCNsLYSu9umNgRPYOtr7PCbbimwoX768eZbQ0FBTpUoV8+OPP4ZN27t3r9m0aZNZv369adKkyTO3ERNnz541JUqUiJdtxaeqVasaPz+/p6b7+/ubIkWKmNWrV4dN69+/v5k7d27YeO3atc3OnTufWO+ll14y58+fT7jAKsm4ffGYWdKunDlSpKjZUKOU2bf2d6sjqXgC7DJRfK/GpKUwANsppGHE8pZUY0x7Y0xOY4yrMSa3MWaKffpyY0xhY0xBY8zQmGwrvNicPlq/fj2urq706NEjbJq3tzc1a9YEwN/fnzZt2lC0aFE6dOjwuJjh5eXF9evXAdi1axd16tQBYODAgXTp0oU6depQoEABxo0b99Q+z5w5Q9myZdm588neQC5fvkytWrXw9vamZMmSbN68+Yl9TZw4Meyv9Pz581O3bl0AVq9eTdWqVSlXrhyvvPIK/v62+8r79etH8eLFKV26NB9++OFTOU6cOEGqVKnImjXrU/P69OlD48aN8fHxCZu2bt066tWrF+3xbNasGbNnz452GZU8ZMxVhKazdnNvWG9cgkJx6TWApe+14f69m1ZHUwkpqmrhCENMWgpjx4417733XqTz1q9fbzJkyGAuXLhgQkJCTJUqVczmzZuNMcbky5cv7C/snTt3mtq1axtjjBkwYICpWrWqCQgIMH5+fiZLliwmKCgorKVw7Ngx4+3tbfbu3fvU/kaOHGmGDBlijDEmODjY3L1796l9GWNMUFCQqVGjhlm8eLHx8/MzNWvWNP7+/sYYY77++mszaNAgc+PGDVO4cGETGhpqjDHm1q1bT+1v6tSp5v33339q+vz5802ZMmVMQEBA2DQ/Pz9Tp06dJ5aLrKWwZcsW07Rp00iPp0q+7t66apb0bGaOFClqNlYraXav/NXqSCoOiEtLQUTuichd+xAgIiEicjcR6lWiqFSpErlz58bJyQlvb2/OnTv3zHWaNGkS9he4h4cHV6/anrfz8/OjRYsW/Pbbb3h7ez+1XsWKFfnll18YOHAgBw8eJH369JFuv0+fPrz44os0a9aM7du3c+TIEapXr463tzfTp0/n/PnzZMiQgdSpU/Pmm28yf/580qZN+9R2Ll++TMTbdi9evMi7777LzJkzSZUqVdj01atXU79+/Wd+dg8PDy5ditV9ASoZSJ/Jg6Y/LObBqH44hUKaPkNZ2rsl/neuWx1NxbOYvHktvTEmg31IDbyM7WU7lonN6aMSJUqwe/fuKOeH/2J0dnYmODgYABcXF0JDbY9KBAQExGidjBkzkidPHrZu3RrpvmrVqsWmTZvw9PTkjTfeYMaMGU8tM23aNM6fP8+AAQMAW0vOx8eHffv2sW/fPo4cOcKUKVNwcXFhx44dvPzyyyxcuJCGDRs+ta00adI8kd0YQ8eOHcNOO4W3YsWKSLcRUUBAAGnSpHnmcip5Kt+4I94rN3C6fjEKrj3OvoZ12LXsF6tjqXgUk2sKTzDGLATicktqnJlY3H304osvEhgYyOTJ/z7Gv3PnTjZu3Bjtel5eXmHFZN68eTHK5ebmxsKFC5kxYwYzZ858av758+fx8PCgW7dudO3alT179jwxf/fu3YwcOZLffvsNJyfb/5oqVaqwdetWTp2yvdLiwYMHnDhxAn9/f+7cuUPjxo0ZM2YM+/bte2p/xYoVC1sPYOTIkaROnZpevXo9sZwxhgMHDkTauonoxIkTz32XlUoe3DP+h6bj5hMw7gtCnYR0H3zL0h7NuHvzitXRVDxwedYC9ncqPOYEVCDcw2ZJnYiwYMEC3nvvPb7++mtSp06Nl5cXY8aM4eLFi1GuN2DAALp27cqwYcOeuG3zWdKlS8fSpUvx8fEhXbp0tGjRImzehg0bGDFiBK6urri7uz/VUvj++++5efNm2AXmChUq8PPPPzNt2jTat29PYGAgAEOGDCF9+vS0aNGCgIAAjDGMHv30y/Bq1arFBx98gDEGEeGLL74gd+7cT3z5Z86cmREjRlC2bFlEInuu8Enr16+nSZMmMT4eKvkqW78D96s2YkP/t8i/4hAHG9XD9dN3qNTyLaujqTiISd9H4duGwdg6x5tsjLH87Wvadfaz9enTh2bNmkV7V9GQIUMoVKgQ7dq1i3ZbgYGB1K5dmy1btuDi8sy/J1QKcmDDXG59ORgPv0ecruFFja9/JlNWT6tjqShE1/dRknqfQkyJSDOgWaFChbqdPHnS6jhJ2tWrV/n7779p3rx5nLd18uRJLl68GHZ7rlLhPXxwlz8HvIXX0n3cdXeCj3tQ9ZV3rI6lIvFcRUFE+kezTWOM+So+wsWFthSUSnoOb1nEtS8GkONKIGeq5KHatz+T2SOv1bFUOM/bS+r9SAaArthehqOUUk8pUaMF1VZs4WzrCuTdcYHjjRuxddYoq2OpGIqyKBhjvns8AJOANEBnYDZQIJHyKaUcUKo07jQe9isuU0fxIIMbWQZNZtkbL3Hj8lmro6lniPaWVBHJIiJDgAPY7lQqZ4z5xOqLzNpLqlKOoViVRtRYsZWzbauQd/clTjVtwqYZ31gdS0UjyqIgIiOw9Wh6DyhljBlojLmVaMmiEZvnFJRS1nJLlZbGg3/BbfpY7mVOTbZh01jWvg7XfE9YHU1FIrqWwgdALuAL4FK4ri7uJaduLhxBly5d8PDweOqhsfHjx1OkSBFKlCjBxx9/bFE6pWKmSMX61F7+F+deq0GeA1c536wlG6d8FdZzgEoaorum4GSMSROhm4sMj8cTM2RK16lTJ1auXPnEtPXr17No0SIOHDjA4cOHI+0lVamkxtUtNY36Tybtbz9y2yMtHiNmsrJdba6eP2p1NGUX624uVOKrVasWWbJkeWLajz/+SL9+/cL6YfLw8LAimlLPpVDZOtRd+hf/dKyL55HrXGjxMut/6q+thiTAIR9LDffwWqLud9CSwxy5FL9nzornysCAZiVivd6JEyfYvHkzn3/+OalTp2bkyJFUrFgxXrMplZBcXN1o8OkPnGm6hVMf9yXP6D9YtWodZUZOJFeBUlbHS7EcsqWgF5ohODiYW7dusX37dkaMGEHbtm1xxKfTlSpQqgYvLd3GhTcbkPPkTS63asu68f201WARh2wpWOV5/qJPKLlz56Z169aICJUqVcLJyYnr168/9f4EpRyBs7ML9T8cw/kmOzj28TvknbCI1Ws3UnLEBHIXLmd1vBTFIVsKClq2bMmff/4J2E4lBQUFRfraTaUcSb5ilai3cCsXezTF4+xtrrXpwNpRHxASEmx1tBRDi4IDaN++PVWrVuX48ePkzp2bKVOm0KVLF86cOUPJkiVp164d06dPj1HX10oldc7OLtR7bwQe82ZyLX8mPCctZ22L6pw/+rfV0VIEh+wl9THtEE+p5C00NJT1P3xOpskLcQ6Fa/+tz0vvf4ezs575jovn7RBPKaUs5eTkxEu9h5NzwRwuF85CnimrWde0KmcObrE6WrLlkEVB+z5SKmXJVaAUDf7YzJX325Ll8n3ute/GqmFvE/woyOpoyY5DFgW9JVWplMfJyYm63QeRZ9E8LhbPRt4Z61nfpCqn9m6wOlqy4pBFQSmVcmXPV4yGszdw7eMOZPJ7yIPXe7JicDceBQVYHS1Z0KKglHI4Tk5O1O7yBfmWLORCmRx4zdzCpsbVOL5ztdXRHJ4WBaWUw/LIXZgmM9dz/fPOuN8KIKhjH5Z/2YmgwAdWR3NYWhQcQEBAAJUqVaJMmTKUKFGCAQMGAPDRRx9RtGhRSpcuTatWrbh9+7a1QZWySM03PqbQ0mWcr+BJ/j/+ZkvD6hzdttzqWA5Ji4IDSJUqFX/++Sf79+9n3759rFy5ku3bt+Pj48OhQ4c4cOAAhQsXZvjw4VZHVcoy/8mZn6Yz1nJzYHfS3gsiuOsHLP/0DYIe3n/2yiqMFgUHICK4u7sD8OjRIx49eoSIUL9+fVxcbA/xVKlSBV9fXytjKpUkVG/XlyIrVnK+Sl7yL9jFXw2rcXjzIqtjOQx9LDA2VvSDKwfjd5s5SkGjr5+5WEhICOXLl+fUqVP06tWLypUrPzF/6tSpvPrqq/GbTSkHlTlbHppOXcX2eRNw++YH6N6PZU1n8eKgSaRJq+8Ii462FByEs7Mz+/btw9fXlx07dnDo0KGweUOHDsXFxYUOHTpYmFCppKfKy70osWItZ2vmp8CS/extVJcHe/daHStJc8i+j8K9ZKfbyZMnrY6T6AYNGkS6dOn48MMPmT59OhMnTmTdunWkTZvW6mhKJVk7Fk/G/bsZyLUbZOnYkWx93sUpTRqrY1ki2fV9lNKeaPbz8wu7s+jhw4esXbuWokWLsnLlSr755hsWL16sBUGpZ6jUvBtFlq0kU7tXuTltGmdatuSBdqj5FIcsCinN5cuXqVu3LqVLl6ZixYr4+PjQtGlTevfuzb179/Dx8cHb25sePXpYHVWpJM3ZPR05Bwwg77RpEBLK+Tf+y5UhQwl9oM81POaQp48e066zlVLPK/TBA66NHsOt337D1dOTnEOGkK5K5WevmAwku9NHSikVV05p05Lj88/I99uviLMz/3TqxOWBAwnxT9nPNWhRUEqlaGnLlyf/wgVk6dyZ27/P4UzzZvhv2Wp1LMtoUVBKpXhOadKQ/ZOP8Zo1E6fUabjw5ptc+vxzQu7etTpaotOioJRSdmm8vcm/YD7/6daNOwsWcqZZc+5t2GB1rESlRUEppcJxSpUKjw/ex+v333HOkAHfHj259Ek/QlJIh5NaFJRSKhJpSpXEa95csr7dkzvLlnG6WTPurVtndawEp0XBQdy+fZs2bdpQtGhRihUrxrZt28LmjRw5EhHh+vXrFiZUKvlxcnMj27vvkv+PObj8Jyu+vXpz8f0PCL51y+poCUaLgoPo06cPDRs25NixY+zfv59ixYoBcOHCBdasWUPevHktTqhU8pW6WDHy/zGHrO++w901azjTpCl3V660OlaC0KLgAO7evcumTZvo2rUrAG5ubmTKlAmAvn378u233yIiFiZUKvkTV1eyvf02+efOxTVnTi6+1xffd/sQnMxa6Emm62wRKQb0AbIC64wxP1oc6Snf7PiGYzePxes2i2YpyieVPol2mTNnzpAtWzY6d+7M/v37KV++PGPHjmXdunV4enpSpkyZeM2klIpa6iKF8fp9Njem/sL18eM5s2MH2T//nAxNmySLP84StKUgIlNF5JqIHIowvaGIHBeRUyLSD8AYc9QY0wNoC0T6+HVKFRwczJ49e+jZsyd79+4lXbp0DBw4kKFDhzJ48GCr4ymV4oiLC1m7dyP/wgW45cvHpY8+wrdXbx5du2Z1tDhL0L6PRKQW4A/MMMaUtE9zBk4APoAvsBNob4w5IiLNgX7A98aYmc/afkrp++jKlStUqVKFc+fOAbB582YGDhzIwYMHw3pH9fX1JVeuXOzYsYMcOXJYmFaplMWEhHBz+gz8xo5FUqUi+6efkrFliyTdarCs7yNjzCbgZoTJlYBTxpgzxpggYDbQwr78YmNMNUDfFhNOjhw5yJMnD8ePHwdg3bp1lCtXjmvXrnHu3DnOnTtH7ty52bNnjxYEpRKZODvzny6dyb9wAaleeIHLn37KhR49eHTlitXRnosVF5o9gQvhxn0BTxGpIyLjROQnYHlUK4tIdxHZJSK7/Pz8EjprkjF+/Hg6dOhA6dKl2bdvH5999pnVkZRS4aTKn598v84g+2ef8WDHTs40bcatP/7A0XqiTvCus0XEC1ga7vTRK0ADY8yb9vE3gErGmHdiu+2UcvpIKeVYgv75h8tffMmDHTtIV60aOb8ajKunp9WxwiS1rrN9gTzhxnMDl2KzARFpJiKT7ty5E6/BlFIqPrjlzUveab+QY0B/Hu7bx5lmzbk1axYmNNTqaM9kRVHYCbwgIvlFxA1oByyOzQZS2us4lVKOR5ycyNy+PQWWLCaNdxmuDBrMP506E3ThwrNXtlBC35I6C9gGFBERXxHpaowJBnoDq4CjwBxjzOGEzKGUUlZx9fQkz5Qp5PhqMAFHjnCmeQtuzvg1ybYaHPJ1nCLSDGhWqFChbidPnrQ6jlJKxcijy5e5PGAA9zdtJk25cuQcOoRU+fMneo6kdk0hzvT0kVLKEbnmzEmen34i5/DhBJ46xdmWrbgx9RdMSIjV0cI4ZFFQSilHJSJkatWSAkuWkK56da59+y3nX+tA4OnTVkcDtCg4hC5duuDh4UHJkiXDpn300UcULVqU0qVL06pVK27bXwDy6NEjOnbsSKlSpShWrBjDhw+3KLVSKjqu2T3IPeF7co0YQdC5c5xt1ZrrkyZjgoMtzeWQRSGl3ZLaqVMnVkboptfHx4dDhw5x4MABChcuHPbl/8cffxAYGMjBgwfZvXs3P/30U1j3GEqppEVEyNisKQWWLcW9Th38Ro3i3KvtCDh+wrJMDlkUUto1hVq1apElS5YnptWvXx8XF1snt1WqVMHX1xew/ZLdv3+f4OBgHj58iJubGxkyZEj0zEqpmHPJmpXc48biOWY0jy5f5mybNvhNmIB59CjxsyT6Hh3YlWHDCDwav11npypWlBxx7LJi6tSpvPrqqwC0adOGRYsWkTNnTh48eMDo0aOfKihKqaQpQ8OGpK1UiatDhnJ9/PfcW7OWXMOGkrp48UTL4JAthZR2+ig6Q4cOxcXFhQ4dbH0I7tixA2dnZy5dusTZs2f57rvvOHPmjMUplVIx5ZIlC56jviP39+MJvn6ds21f5drYsYQGBSXO/hNlL/HMGLMEWFKhQoVuibnfuP5FH9+mT5/O0qVLWbduXVg3vTNnzqRhw4a4urri4eFB9erV2bVrFwUKFLA4rVIqNtLXq0faChW4Onw4N36ciP/ateQcNow0pUol6H4dsqWgYOXKlXzzzTcsXrw47J0KAHnz5uXPP//EGMP9+/fZvn07RYsWtTCpUup5OWfKRK5vviH3jz8Qcucu515tx7XvviM0MDDB9qlFwQG0b9+eqlWrcvz4cXLnzs2UKVPo3bs39+7dw8fHB29vb3r06AFAr1698Pf3p2TJklSsWJHOnTtTunRpiz+BUiou0tetS4GlS8jYqiU3Jv/M2VateXjgQILsyyG7uXhMu85WSqU0/pu3cHlAf3IOGox7zRrPtY3ourlwyGsK4fo+sjqKUkolKveaNSi4ciVObm4Jsn2HPH2U0p5TUEqp8BKqIICDFgWllFIJQ4uCUkqpMFoUlFJKhXHIoqBPNCulVMJwyKKgF5qVUiphOGRRUEoplTC0KCillArj0E80i4gfcD6SWRmBqC44RDUv4vTIlos4LStwPUZhn190nyW+1n3WcvFxPCObllyPZ0yWjc1xi2q6Hs+YzY/JsYvJtMQ4nlHliO918xljskU6xxiT7AZgUmznRZwe2XKRLLPLys8SX+s+a7n4OJ4xOcbJ5XjG5Zjq8Yzf4xnTYxeTaYlxPBPzmEY1JNfTR0ueY17E6ZEtF912E0pc9hnTdZ+1XHwcz8imJdfjGZNlY3PcopquxzNm82P6bzkl/ZuPkkOfPrKaiOwyUXQqpWJPj2f80uMZv1LK8UyuLYXEMsnqAMmMHs/4pcczfqWI46ktBaWUUmG0paCUUiqMFgWllFJhtCgopZQKo0UhnohISxGZLCKLRKS+1XkcnYgUE5GJIjJXRHpanSe5EJF0IrJbRJpancXRiUgdEdls/z2tY3We+KJFIRoiMlVEronIoQjTG4rIcRE5JSL9AIwxC40x3YBOwKsWxE3yYnk8jxpjegBtgWR/G+Dzis0xtfsEmJO4KR1HLI+nAfyB1IBvYmdNKFoUojcNaBh+gog4AxOARkBxoL2IFA+3yBf2+epp04jF8RSR5sAWYF3ixnQo04jhMRWResAR4Gpih3Qg04j57+hmY0wjbIV2UCLnTDBaFKJhjNkE3IwwuRJwyhhzxhgTBMwGWojNN8AKY8yexM7qCGJzPO3LLzbGVAM6JG5SxxHLY1oXqAK8BnQTEf33H0FsjqcxJtQ+/xaQKhFjJigXqwM4IE/gQrhxX6Ay8A5QD8goIoWMMROtCOeAIj2e9nO0rbH9Y1ue+LEcWqTH1BjTG0BEOgHXw32pqehF9TvaGmgAZAK+tyBXgtCiEHsSyTRjjBkHjEvsMMlAVMdzA7AhcaMkG5Ee07AfjJmWeFGShah+R+cD8xM7TELT5mPs+QJ5wo3nBi5ZlCU50OMZ//SYxq8UdTy1KMTeTuAFEckvIm5AO2CxxZkcmR7P+KfHNH6lqOOpRSEaIjIL2AYUERFfEelqjAkGegOrgKPAHGPMYStzOgo9nvFPj2n80uOpHeIppZQKR1sKSimlwmhRUEopFUaLglJKqTBaFJRSSoXRoqCUUiqMFgWllFJhtCgopZQKo0VBOSwRySEis0XktIgcEZHlIlLY/vKTpfG0D6+IfevHYt3lIpIpHjKMEZFa9p83iMg/IiLh5i8UEf9weR+KyD4R2S8if4lIEfu8UiIyLa55VPKmRUE5JPuX4gJggzGmoDGmOPAZkN3aZP8yxjQ2xtyOyzZEJAtQxd6l82O3ger2+ZmAnBFWO22M8TbGlAGmYzsuGGMOArlFJG9cMqnkTYuCclR1gUfhuyg3xuwzxmy2j7rbX+V5TET+9/gvaxE5JyJZ7T9XEJEN9p8H2t+6tUFEzojIuxF3KCIFRGSviFSMMD2niGyy/3V+SERqht+XiPSwz9snImdFZL19fn0R2SYie0TkDxFxj+RztgFWRpg2G1v/O2DrXjy6njozYOvv/7El4dZV6ilaFJSjKgnsjmZ+WeA9bG/KKoD9L+tnKIqtf/xKwAARcX08w34KZh7Q2RizM8J6rwGrjDHeQBlgX/iZxpiJ9nkVsfW4OcpemL4A6hljygG7gPcjyVQ9ks+5DqhlfyNYO+D3CPML2gvQafs2R4WbtwuoGdUBUErfp6CSqx3GGF8AEdkHeGF7tWd0lhljAoFAEbnGv6eisgGLgJej6AhtJzDVXkQWGmP2RbH9scCfxpglItIUW8Haam/EuGHriC2inIBfhGkh9s/yKpDGGHMu3CUGsJ8+AhCRV4FJ/PuKyWtArijyKaUtBeWwDgPlo5kfGO7nEP79AyiYf3/vU8dwnTvY3rwVaWvDfr6/FnAR+FVE/htxGfvbzvLx77t8BVhjP/fvbYwpbozpGsnmH0aSE2ynkMYDcyLLFM5ie7bHUtu3qVSktCgoR/UnkEpEuj2eICIVRaT2M9Y7x7/F5OUY7isIaAn8V0ReizhTRPIB14wxk4EpQLkI88sDHwKvh3sF5naguogUsi+TVkQKR7Lvo0ChSKZvBoYDs56RvQZwOtx4YeC57qZSKYMWBeWQjK3P91aAj/2W1MPAQJ79RqxBwFgR2YytNRDT/d0HmgJ9RaRFhNl1gH0ishdboRkbYX5vIAuw3n6u/2djjB/QCZglIgewFYmikex6mX37EfMYY8xIY8z1SNZ5fE1hPzAMeDPcvLr2bSoVKX2fglJJnIhsAZrGw+2tqYCNQA37i2OUeooWBaWSOBGpDDw0xhyI43ZeADyNMRviJZhKlrQoKKWUCqPXFJRSSoXRoqCUUiqMFgWllFJhtCgopZQKo0VBKaVUmP8DqsZePSh0ovoAAAAASUVORK5CYII=\n", - "text/plain": [ - "
" + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEaCAYAAAD+E0veAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA+uklEQVR4nO3ddXxWdRvH8c+1okY+MGLECOkY3aWMbhFB9KEEQVDExKIkVJAURRAE9AFEulNakO5OGTmaAdvY9nv+uG/mGNvYWJzd2/V+vc7Lnf7ex3Ff+534HTHGoJRSSgE4WR1AKaVU0qFFQSmlVBgtCkoppcJoUVBKKRVGi4JSSqkwWhSUUkqF0aKgkg0RqSMivhbsd4OIvJnI+7Tks6rkT4tCCiQi7iJyTkReCzctvYj8IyJtRKSbiBwVkVTh5v9HRK6JSEMRWSEi/hGGABExIpLXvrybiAwUkZMict++v6ki4iUih8OtF2Jf9/H4Z8/I3ltEdolIoIhMS7CD9OQ+c4rIYhG5ZP+MXom032T1WZ/1eUTkJRE5JiIPRGS9iOQLN09E5BsRuWEfvhURia9s6l9aFFIgY4w/0B0YKyLZ7JO/BXYZY+YaYyYDvkD/cKuNAZYbY1YaYxoZY9wfD0BGYDswwxjzj335uUBz4DX7/DLAbuAlY0yJcOtuBnqH296wZ8S/BAwBpsbtKMRKKLASeDkR9wnJ77NG+XlEJCswH/gSyALsAn4Pt0h3oCW236PSQFPgrQTIqIwxOqTQAZgGzALqADeAnOHmeQG3AG+gPrZ/0Jmj2M43wAEgrX28HvAQyBODDBuAN58j+xBgWoRpdQDfcOPvAkeA3JGs/yMwN8JnWAdINPt0AQzgFdVnAHLaj8WHkaxfELgJlLOP5wKuA3Uc5bMm4OfpDvwVbjyd/XeoqH38L6B7uPldge2J/W8mJQwuT1UJlZL0xfZF4oPtS+zy4xnGmHMi0h/bX3WZgbeNMbcibkBEWmD7i62iMeaBfXI9YIcx5kJsA4lIP6CGMaZprD/Nk9v5EmgF1DbG+EWyyAfAPhHpBJzG9iXjbezfOM+5Ty9gNTDSGDMp4nxjzGkR+QT4n4iUB37B9uW44Xn3ad9von9WiPfPUwLYH27b90XktH36sYjz7T+XeN7sKmpaFFIwY8wtETkMVMPWdI/oe+B1YJ8xZmHEmSJSEFtro6sx5mS4Wf8BLkdcPoaZvn6e9Z6MJaOASkBdY8ydKPbzQERex3aq5B7wjjEmLhduiwNfAJ8aY2ZFtZAxZrKINAP+xvaXePM47NOqzxp+2/H1edyBiAXtDpA+3Pw7Eea5i4jEtbipJ+k1hRTM/kXhBazFdkrhCfZ/bEeBw5GsmxrbdYOpxph5EWbfwHYaxQqZsJ2KGB7Vl+RjxpgdwBlAgDlx3G8H4CK2Y/Isk4GSwHhjTGAc9pkJaz5rRPHxefyBDBGmZcBWxCKbnwHw14IQ/7QopFAi4gGMBrphO/3TVkRqxWITE4D7wCeRzFsLVBKR3HEOGnu3sF2E/EVEqke3oIj0AlJhu17ycRz3OxDb+fSZIuIczT7dsV20nwIMFJEscdinVZ81/Hbj6/McxnYR+fF202G7ZnE4svn2n5/6Y0XFnRaFlOt7YKExZr39WsLHwOTwt6FGRUS6YPsyamuMCY443xizFlgDLBCR8iLiYr/ltYd93edm31ZqwBlwFpHUIvLEaVD7Oe0O9v1XjmI7hbFd8HwdeAP4WES8o9lvamxfqgCp7OPhPQJewXaB9FcRierf1lhgtzHmTWAZMNEBP2tCfJ4FQEkRedm+TH/ggDHmmH3+DOB9EfEUkVzYrpNMiyaXel5WX+nWIfEHbLf2XQIyRZi+DhgaYdo0YEiEaWewfQn6RzLUtC/jBgwCTmFrUZwHfgbyRtjWBsLdfQR8BqyIJvtAbOeuww8D7fPq8OQdOU2Aq0D5CNtwAXYA/cJN6wkcBFJFsd+I+zSRfQYgNbaW0jTAKcI2WmA7xZTFPu5uPz4dHOWzJtTnsc+vh+2i8kP7MfUKN0+w3TZ90z58SzR3T+nw/IPYD7hSSimlp4+UUkr9S4uCUkqpMFoUlFJKhdGioJRSKowWBaWUUmEcupuLrFmzGi8vL6tjKKWUQ9m9e/d1Y0y2yOY5dFHw8vJi165dVsdQSimHIiLno5qnp4+UUkqF0aKglFIqjBYFpZRSYbQoKKWUCqNFQSmlVBiHLAoi0kxEJt25E+17RZRSSsWSQxYFY8wSY0z3jBkzPtf6fvcCWXf0ajynUkopx+eQRSGuJm06Tdfpu3hv9l5u3Q+yOo5SSiUZDv3w2vP6qEFR0rq5MGH9KbacusGQliVoWNKqVworpVTSkSJbCm4uTvT1Kczi3jXIniEVPX7bQ6+Ze7jhH5d3qCullONLkUXhseK5MrCwV3U+rF+Y1Yev4DN6E0v2X0LfRqeUSqlSdFEAcHV2oveLL7Ds3ZrkyZyGd2btpcdvu7l2L8DqaEoplehSfFF4rHD29MzrWY1+jYqy/rgf9UdvYsFeX201KKVSFC0K4bg4O9GjdkGWv1uTAlnT0ff3/bw5fRdX7mirQSmVMmhRiEQhD3f+6FGNL5oUY+vp6/iM3sicXRe01aCUSva0KETB2Ul4s2YBVvapRbGcGfh47gE6/rKTi7cfWh1NKaUSjBaFZ/DKmo7Z3aowqHkJdp27SYPRm/jf3+e11aCUSpa0KMSAk5PQsZoXq96rRencGfl8wSFen/I3F24+sDqaUkrFKy0KsZAnS1r+92ZlhrUqxf4Ld2gwZhMztp0jNFRbDUqp5CHJFAURcRKRoSIyXkQ6Wp0nKiLCa5XzsqpvLcrny0z/RYdpP3k7567ftzqaUkrFWYIWBRGZKiLXRORQhOkNReS4iJwSkX72yS0AT+AR4JuQueKDZ6Y0zOhSiW9fLs2Ry3dpOHYTU7acJURbDUopB5bQLYVpQMPwE0TEGZgANAKKA+1FpDhQBNhmjHkf6JnAueKFiNC2Yh7W9K1NtYJZ+WrpEdr+tI3Tfv5WR1NKqeeSoEXBGLMJuBlhciXglDHmjDEmCJiNrZXgC9yyLxMS1TZFpLuI7BKRXX5+fgkRO9ZyZEzNlI4VGP1qGU5d86fR2M1M3Hia4JBQq6MppVSsWHFNwRO4EG7c1z5tPtBARMYDm6Ja2RgzyRhTwRhTIVu2bAmbNBZEhFZlc7Pm/VrULZKNr1cc4+Uf/+LE1XtWR1NKqRizoihIJNOMMeaBMaarMeYdY8yERE8VTzzSp2bi6+UZ374sF249pOm4LXz/50keaatBKeUArCgKvkCecOO5gUux2UBSf0eziNCsTC7W9K2FT4nsjFx9gpYTtnLk0l2roymlVLSsKAo7gRdEJL+IuAHtgMWx2UBc39GcWP7jnooJr5Vj4uvluHo3gObfb2H0mhMEBWurQSmVNCX0LamzgG1AERHxFZGuxphgoDewCjgKzDHGHI7ldpN0SyGihiVzsqZvbZqWzsnYdSdp/v0WDl10jOxKqZRFHLkPnwoVKphdu3ZZHSNW1hy5yucLDnLjfhA9ahfg3ZdeIJWLs9WxlFIpiIjsNsZUiGxeknmiOaXwKZ6dNX1r06qsJxPWn6bpuC3su3Db6lhKKQVoUbBExrSujHylDL90roh/YDCtf9jK8OVHCXgU5eMZSimVKByyKDjaNYWo1C3iwaq+tXi1Yh5+2nSGxmM3s+tcxGf9lFIq8eg1hSRiy8nrfDLvAJfuPKRztfx82KAwad1crI6llEqG9JqCA6jxQlZW9a3F65XzMXXrWRqN3cz2MzesjqWUSmEcsigkl9NHEbmncuGrliWZ1a0KxkC7Sdvpv+gQ9wODrY6mlEoh9PRREvUgKJgRq44z7a9zeGZKwzcvl6Z6oaxWx1JKJQN6+sgBpXVzYUCzEvzxVlXcnJ3o8PPffDr/IPcCHj254L0rEBxoTUilVLKjRSGJq+CVheV9atK9VgF+3/kP9UdvYsPxa7aZxsCcjjCxJvgmzxaTUipxaVFwAKldnfmscTHm9ayGeyoXOv2ykw//2M+dh8FQ6yMIug9TfGD1F/DoodVxlVIOzCGLQnK90PwsZfNmZum7NehVtyAL9l7EZ/RG1j4qBW9vg3Id4a/xMLEG/LPd6qhKKQelF5od1EHfO3w0dz/HrtyjVVlPBjQrTqYrf8Hid+D2BajSE178AtzSWR1VKZXE6IXmZKhU7ows7l2DPi+9wJL9l6g3ahMrHxSFntug4puw/Qf4sTqc22J1VKWUA9Gi4MDcXJzo61OYRb2rkz1DKnr8tpte805wo/ZQ6LTMttC0JrDsAwj0tzasUsohaFFIBkrkysjCXtX5sH5hVh++gs/oTSy5UwDTcytU6QU7p8APVeH0equjKqWSOIcsCin1QnN0XJ2d6P3iCyx7tyZ5MqfhnVl76fH7Ua5V7w9dVoFLKvi1JSx+FwL0uCmlIqcXmpOh4JBQft5yllFrTpDWzZmBzUrQokRmZOPXtjuU0ueEZuPghXpWR1VKWUAvNKcwLs5O9KhdkOXv1qRA1nS89/s+us06zNXKn0HXtZAqPfzvZVj4Njy8ZXVcpVQSokUhGSvk4c4fParxRZNibDl1nXqjNjLnigem+0ao+QHsnw0TqsDxFVZHVUolEVoUkjlnJ+HNmgVY2acWxXJm4OO5B+j06wEulf8Iuv0J6bLCrHYwrxs80Bf8KJXSaVFIIbyypmN2tyoMal6CneduUn/0JmZeyILp9ifU+RQOz4cJleHIYqujKqUspEUhBXFyEjpW82LVe7UonTsjny04yOvT9nKh9LvQfSNkyAlz3oA/OsH961bHVUpZwCGLgt6SGjd5sqTlf29WZmirkuy/cIcGYzYx46w7oV3WwotfwrFlMKESHJpn64lVKZVi6C2pKdzF2w/pN+8Am09ep3L+LHzbpjT5Qv6BRb3g4m4o2hSajIL02a2OqpSKJ3pLqoqSZ6Y0zOhSiW9fLs2Ry3dpMGYTU06kJqTzKvAZDCfX2FoN+2drq0GpFECLgkJEaFsxD2v61qZawax8tfQIbSfv5HThrtBzK2QrAgvest2ldPeS1XGVUglIi4IKkyNjaqZ0rMCotmU4dc2fxmM389NhJ0I6LocGw+HMRttzDXt+1VaDUsmUFgX1BBGhdbncrOlbi9qFszF8xTFa//Q3Jwu8YWs15CgJi3vDby/b3tuglEpWtCioSHlkSM1Pb5RnXPuy/HPjPk3GbWHCAUPwG4uh8Ujb291+qAq7pmqrQalkRIuCipKI0LxMLta8Xxuf4tkZseo4LX/cxtE8r8Lbf4FnOVjaF2Y0h1vnrI6rlIoHWhTUM2V1T8WEDuX4sUM5rtwJoPn3WxizO4ig1xZA0zFwcS/8UA3+ngShoVbHVUrFgRYFFWONSuVkTd/aNC6VkzFrT9J8wlYO5WwNb2+DvFVgxUcwvSncOG11VKXUc3LIoqBPNFsnczo3xrYry+T/VuDm/SBaTNjKyL8fENhuDrSYAFcO2d4NvW0ChIZYHVcpFUv6RLN6bncePGLw0iPM2+NL4ezufNumDN4ZH9iuM5xYCbkr2QpFtsJWR1VKhROnJ5pF5FsRySAiriKyTkSui8jr8R9TOZqMaV35rm0ZfulckXsBwbT+YSvDt94hoM3/oNUkuH4CJtaALWMgJNjquEqpGIjJ6aP6xpi7QFPAFygMfJSgqZRDqVvEg1V9a/FqxTz8tPEMjcdvYXcmH+i1A17wgbUDYIoPXDtqdVSl1DPEpCi42v/bGJhljNE3sainZEjtyvDWpfm1ayUCH4XSZuI2vtp4g4etpkObX+D2efipFmwaASGPrI6rlIpCTIrCEhE5BlQA1olINiAgYWMpR1XzhWys6luLDpXzMmXLWRqN28zfaWvbWg1Fm8CfQ2Dyi3DloNVRlVKRiNGFZhHJDNw1xoSISDrA3RhzNcHTPYNeaE7a/jp9nX7zDvLPzQd0rJqPjxsWJd3p5bDsfXh4C2p+aHtXtIub1VGVSlHieqF5qjHmljHm8f2FAsyMz4AqeapWMCsr36tJ5+pezNh+ngZjNrHVrZqt1VCiNWz8GibXhUt7rY6qlLKLyemjiyLyI4S1GFYDvyVoKpVspHVzYUCzEsx5qyquzk50+PlvPlt1kXtNfoD2s+HBDZj8EqwdBMGBVsdVKsV7ZlEwxnwJ3BWRidgKwnfGmF8SPJlKVip6ZWFFn5p0r1WA2Tv+ocHoTWyUCvD2dijTHrNlFEN+qcLBw39YHVWpFC3Kawoi0jr8KPAlsANYCWCMmZ/g6Z5Bryk4pr3/3OKjuQc4dc2fthVy83mT4jw8vYjXdwzEzwk6ur/A202mkjpNZqujKpUsRXdNIbqiEF1rwBhjusRHuLjQouC4Ah6FMHbdSSZtOkNWdzeGtSpFJc8QvlvemXmBl/AKga8q9sO7VAeroyqV7DxXUXAEWhQc3wHf23z0xwGOX71Hq7KeDGhWnKNHpzJw/wQuOxk6pCvAO02mkjZtVqujKpVsxKko2J9L6AZ4AS6Pp8d3S0FE6gBfAYeB2caYDc9aR4tC8hAUHMr360/xw/pTZErrxpCWJanpBaOXdeH3gAvkCYFB5fpS0dvyxqlSyUKcbkkFFgEZgbXAsnBDTHY8VUSuicihCNMbishxETklIv3skw3gD6TG1p2GSiHcXJx436cwi3pXxyN9Knr8tptPFl+iZ5MFTC3TFwN02T+aob835oH/NavjKpWsxaSlsM8Y4/1cGxephe2LfoYxpqR9mjNwAvDB9uW/E2gPHDPGhIpIdmCUMeaZJ5O1pZD8PAoJZeKG04z78yQZUrsyuEVJ6hZyZfyyLvzv/hlyhQoDy/SmSvm3rI6qlMOKa0thqYg0fp4dG2M2ARH7SqoEnDLGnDHGBAGzgRbGmMev7LoFpHqe/SnH5+rsxDsvvcCSd2rgmTkNvWbu4YP55+nS8A+ml+uHK9Dt0PcMmt0A/3uXrY6rVLITk6LQB1theCgid0XknojcjcM+PYEL4cZ9AU8RaS0iPwG/At9HtbKIdBeRXSKyy8/PLw4xVFJWNEcG5vesxicNi7Lu2DV8Rm/kfEhd5rTbSKd0hZgfcJGWf/iwZcd4q6Mqlawk+N1HIuIFLA13+ugVoIEx5k37+BtAJWPMO7Hdtp4+ShlOXbvHR3MPsPef29Qr5sHQVqW4emER/f8eymlnQwvX7HzUZCoZM+a1OqpSDiGufR/VimyIQx5fIE+48dzApdhsQF/HmbIU8kjP3B7V+KJJMTafvI7PqI0cf1iN31/bQrf0RVkadIVW8xqzYft3VkdVyuHF5ELzknCjqbFdE9htjHkxRjt4uqXggu1C80vARWwXml8zxhyObXhtKaQ8Z6/f55O5B9hx7iZ1imRjWKtS3Lqygi+3DeKkUyhNXLLSr/FUMmXOb3VUpZKsOLUUjDHNwg0+QEkgRt1mi8gsYBtQRER8RaSrMSYY6A2sAo4Cc56nIKiUKX/WdMzuXoWBzYrz95mb1B+9iQN3KzD7ta30yFCSVY/8aLmgGeu2Drc6qlIOKdbXFEREgAPGmFIJEylGGZoBzQoVKtTt5MmTVsVQFvvnxgM+mXeAbWduUKNQVoa3LsX96+vov/VLjjqF0MA5M581nkqWLIWsjqpUkhLXJ5rHY3uwDGwtC2/gnDHm9fgM+Tz09JEKDTXM3PEPw5fb3v/cr1FR2pbLxvRVPfjx5l7SG/is4Cs0qPEl4hSTm+2USv7iWhQ6hhsNxlYQtsZjvuemRUE95nvrAZ/OP8jmk9epUiAL37xcmke3t/Llpk845BTMS04Z+KLhz2TNVszqqEpZLtl1iKenj1RkjDHM2XWBIUuPEhxq+KhBEV6vlIPfVvViwo0dpDHwiVdLmtYerK0GlaLFtaVQHRgI5MPWIZ5g6zq7QDznjDVtKajIXL7zkM/mH2T9cT8q5MvMt21Kg/8O+m/4kP1Oj6gt7vRvMAmP7JZdFlPKUnEtCseAvsBu4PF7mjHG3IjPkM9Di4KKijGG+XsuMmjJYQKDQ/mgfmE6VcnN7LXvMP7qX7gCH+VtTMu6X2urQaU4cS0KfxtjKidIsjjSoqCe5drdAD5feIg1R67inScTI9qUxi1gH/3X92WPBFKdtAysP5EcOctaHVWpRPO8b14rZ/+xLeAMzAfC3qxujNkTzzljTK8pqNgwxrB4/yUGLj7M/cAQ+tR7gW418jJ3XV/GXN6IE/CBpw9t6n2nrQaVIjxvUVgfzTZNTJ9oTkjaUlCx4XcvkAGLD7H84BVKeWZkxCulcQ8+wsB17/A3AVQmNYPqTcDTs5LVUZVKUMnu7qPHtCio57H84GW+XHiIuwGP6F33BXrWzs/C9R8y6uJaQoG+Oevwqs8YnJxdnrktpRxRXDvEGyYimcKNZxaRIfGYT6lE1bhUTta8X5tGJXMyeu0JWvywjeIlBrGgwXTKOqVl2NWNdP21Chf+2WJ1VKUSXUxOoDYyxtx+PGKMuQU810t34ov2kqriKks6N8a1L8ukN8pz3T+QFhO28r+D7oxtv5XBng05bgJova4Hvy5/i5DgIKvjKpVoYnL30QGgojEm0D6eBthljCmRCPmipaePVHy4/SCIwUuPMH/PRQpnd2dEmzLkcPuHwaveYpPxx9u4MrjOKPJ71bE6qlLxIq6v4/wNWCciXUWkC7AGmB6fAZWyUqa0boxq683UThW4+zCYVj9s5Zc9box8dRPD8jbnjAnilfW9+WVpV201qGQvRheaRaQhUA/b08yrjTGrEjpYTGhLQcW3uwGPGLr0KL/vukDBbOn4tk0Z8qa5xJBV3fkz9C6lQl34qta3FCzoY3VUpZ6b3n2kVCxtOuHHp/MPcunOQ7pWz8/79V5g499DGXZmLvcFemYpR6dGE3F1TWt1VKViLa6nj5RKcWoVzsbK92ryWqW8/LzlLI3HbyFLnndY2Gw+dV0yM+7WXjr8Vo3jJ5dZHVWpeOWQRUHvPlKJIX1qV4a2KsXMNysTYgyvTtrOuK1BDH5lPaMKtucqwbTb+gk/LnyNR4H3rY6rVLyIsiiIyDr7f79JvDgxY4xZYozpnjFjRqujqBSgWqGsrOxTi07VvJi+7TwNx24iXY7uLGy5mPquWfnhzkHazazGkWMLrY6qVJxF11LIKSK1geYiUlZEyoUfEiugUklBulQuDGxegjlvVcVZhNd+/psRG+7zxctrGFe4I7dMCK9t/4Jx89sSFHjP6rhKPbfo+j5qA3QFagARr+Zq30cqxXoYFMJ3q48zZetZcmVMw/DWpSiT/SHfLuvM4kfXKBQMg/O8Tan6Pa2OqlSk4tp19pfGmK8SJFkcaVFQVtp9/hYfz93Pab/7tK2Qm8+bFOfAwR9ZNO8nOi8O5Vz1LNQdMY+0mXNYHVWpJ8T5llQRaQ7Uso9uMMYsjcd8z02LgrJawKMQxqw9yaRNp/FIn5phrUtSPsNtNn/YgYJ77nMtC2R6vyNl2vSzOqpSYeLaId5woA9wxD70sU9TKsVL7epMv0ZFWfB2dTKkcaHLtF0M2vyQGj//hf+nrXAJBpcvprOsc1Xu37hkdVylnimmfR95G2NC7ePOwF5jTOlEyBdVJn3JjkpyAoND+P7PU/yw4TRZ0rkxtWNFvFxvsPHDdhTceY/rmSDdu+0p91p/q6OqFC4+Hl7LFO5ny+8D1VtSVVKUysWZD+oXYVGv6pTyzEi+rGlJn70ATX/dwYMv2wKQZvAslv63Mv7XzlkbVqkoxKSl0B74GliPre+jWsCnxpjZCR8venpNQTkS/xsX2PDBK+TffodbGSHV262p2HGo1bFUChSnloIxZhZQBds7mucDVZNCQVDK0bj/Jw9Np23n0aDXCXUC9+HzWdqhIncvn7I6mlJhtEM8pSzw4NYV1n/0Ml5bb3LbHVx7NKNS12+tjqVSCO0QT6kkJm3mHDT5eSuhw97kkRukH7GEpe3Kc9v3qNXRVAqnRUEpC5Vq9QGVl2/kTO2seO1/wNGWrdk2sY/VsVQKFm1REBEnETmUWGGUSonSZPSgyU+bkRE9eZhGyDRmNUtfKcut8wetjqZSoGiLgv3ZhP0ikjeR8iiVYpVo+i41VvzF2Rezk+9wAMdbt2XreO0/SSWumJw+ygkcFpF1IrL48ZDQwZRKidzcM9H4hw24fNeHB+mELBM2sKy1N9dP77Y6mkohYvKcQu3IphtjNiZIohjQJ5pVShD04A7r+rUm99pLPEgNga9Xp0afn3BydrY6mnJwcX1OYSNwDnC1/7wT2BOvCWNJn2hWKYFb2ow0GreOVOM+5l4GIdukraxsXY5rJ7ZbHU0lYzHpEK8bMBf4yT7JE1iYgJmUUuEUrteZ2it2cL5xHjxPBXG+bWc2fNOR0JAQq6OpZCgm1xR6AdWBuwDGmJOAR0KGUko9yTWNOw1HrSbdj19wO4uQ/ZcdrGxZliuHN1kdTSUzMSkKgcaYoMcjIuICOO5j0Eo5sEK1OlB3xW7+aZ4fz7OPuPjaW/w5pL22GlS8iUlR2CginwFpRMQH+ANYkrCxlFJRcUmVhgbfLifDpK+4kc2JnL/tY1Uzby7tX2t1NJUMxKQo9AP8gIPAW8By4IuEDKWUerb81drw0oq9XGj1AjkvBHPljXdYN/AVbTWoOInJ3UehwHTgK2AQMN04ci96SiUjzq5u1B++mMxTvsYvhxO5Zh9idZMyXNi9zOpoykHF5O6jJsBpYBzwPXBKRBoldDClVMzlq9SCesv2cvGVYnhcCuF6pw9Z83kLQh4FPXtlpcKJyemj74C6xpg6xpjaQF1gdMLGUkrFlrOrG/W+mo/HtNFc9XQm97wTrG1clvN/L7A6mnIgMSkK14wx4d8Ccga4lkB5lFJxlLtcQ+ov28+l10qR9Woot7p+xqpPmhIcGGB1NOUAoiwKItJaRFpj6/douYh0EpGO2O482ploCZVSsebk7MxL/eeQ87fvuZzPhbyLTvNn4/Kc3vy71dFUEhddS6GZfUgNXAVqA3Ww3YmUOcGTKaXiLFfpl2iweB9X/lueLNdD8e85kJUfNCQ48KHV0VQSlaRexyki6YBNwABjzNJnLa+v41Qq5q4e2cLufm+T/8QjLuYUvPr3o1Dd/1odS1kgTh3iiUh+ERklIvNj23W2iEwVkWsRX9QjIg1F5LiInBKRfuFmfQLMicm2lVKxk714DRou2Mu1rlXIeNvwoPdwVrxXj0cP/a2OppKQmFxoXoitl9Tx2O5EejzExDSgYfgJIuIMTAAaAcWB9iJSXETqAUewnapSSiUAJ2dnan/0C/l+/4ULL7jhtfIimxpW4tjqyVZHU0lETN6n8LcxpvJz70DEC1hqjClpH68KDDTGNLCPf2pf1B1Ih61QPARa2R+ci7i97kB3gLx585Y/f/7880ZTKsXbPKY7bjM2kzYALryYk3rfLMAtnXZJn9zF6fQRMFZEBohIVREp93iIQx5P4EK4cV/A0xjzuTHmPWAmMDmyggBgjJlkjKlgjKmQLVu2OMRQStV8bxKF5v6P88VSkX/tZbY0rMqR5ROsjqUsFJOiUAroBnzNv6eORsZhnxLJtLDmijFmWkwuMiul4sd/CpSj6bx93HznRdI+MIR+8D3LetQi8N4Nq6MpC8SkKLQCChhjahtj6tqHF+OwT18gT7jx3MCl2GxARJqJyKQ7d+7EIYZSKrzqvSZQZP4fnCuVhgIb/NjWqAaHFo2xOpZKZDEpCvuBTPG4z53AC/a7mtyAdkCM7mZ6LLav47xy5Qrt2rWjYMGCFC9enMaNG3PixAk2bNhA06ZNY/8JInHu3DlKliz5XOs2btyY27dvx0uOiBYuXMjgwYMBaNCgAd7e3mFDrly5qFz538tF27ZtI3/+/GHz3d3dKVKkCN7e3vz3v//l4MGDdOrUKUFyqqQhc76SNJ2zhzt9G+IWANLvJ5Z1q8HDu35WR1OJxRgT7QBsAG4Cq7B9eS8GFj9rPfu6s4DLwCNsLYSu9umNgRPYOtr7PCbbimwoX768eZbQ0FBTpUoV8+OPP4ZN27t3r9m0aZNZv369adKkyTO3ERNnz541JUqUiJdtxaeqVasaPz+/p6b7+/ubIkWKmNWrV4dN69+/v5k7d27YeO3atc3OnTufWO+ll14y58+fT7jAKsm4ffGYWdKunDlSpKjZUKOU2bf2d6sjqXgC7DJRfK/GpKUwANsppGHE8pZUY0x7Y0xOY4yrMSa3MWaKffpyY0xhY0xBY8zQmGwrvNicPlq/fj2urq706NEjbJq3tzc1a9YEwN/fnzZt2lC0aFE6dOjwuJjh5eXF9evXAdi1axd16tQBYODAgXTp0oU6depQoEABxo0b99Q+z5w5Q9myZdm588neQC5fvkytWrXw9vamZMmSbN68+Yl9TZw4Meyv9Pz581O3bl0AVq9eTdWqVSlXrhyvvPIK/v62+8r79etH8eLFKV26NB9++OFTOU6cOEGqVKnImjXrU/P69OlD48aN8fHxCZu2bt066tWrF+3xbNasGbNnz452GZU8ZMxVhKazdnNvWG9cgkJx6TWApe+14f69m1ZHUwkpqmrhCENMWgpjx4417733XqTz1q9fbzJkyGAuXLhgQkJCTJUqVczmzZuNMcbky5cv7C/snTt3mtq1axtjjBkwYICpWrWqCQgIMH5+fiZLliwmKCgorKVw7Ngx4+3tbfbu3fvU/kaOHGmGDBlijDEmODjY3L1796l9GWNMUFCQqVGjhlm8eLHx8/MzNWvWNP7+/sYYY77++mszaNAgc+PGDVO4cGETGhpqjDHm1q1bT+1v6tSp5v33339q+vz5802ZMmVMQEBA2DQ/Pz9Tp06dJ5aLrKWwZcsW07Rp00iPp0q+7t66apb0bGaOFClqNlYraXav/NXqSCoOiEtLQUTuichd+xAgIiEicjcR6lWiqFSpErlz58bJyQlvb2/OnTv3zHWaNGkS9he4h4cHV6/anrfz8/OjRYsW/Pbbb3h7ez+1XsWKFfnll18YOHAgBw8eJH369JFuv0+fPrz44os0a9aM7du3c+TIEapXr463tzfTp0/n/PnzZMiQgdSpU/Pmm28yf/580qZN+9R2Ll++TMTbdi9evMi7777LzJkzSZUqVdj01atXU79+/Wd+dg8PDy5ditV9ASoZSJ/Jg6Y/LObBqH44hUKaPkNZ2rsl/neuWx1NxbOYvHktvTEmg31IDbyM7WU7lonN6aMSJUqwe/fuKOeH/2J0dnYmODgYABcXF0JDbY9KBAQExGidjBkzkidPHrZu3RrpvmrVqsWmTZvw9PTkjTfeYMaMGU8tM23aNM6fP8+AAQMAW0vOx8eHffv2sW/fPo4cOcKUKVNwcXFhx44dvPzyyyxcuJCGDRs+ta00adI8kd0YQ8eOHcNOO4W3YsWKSLcRUUBAAGnSpHnmcip5Kt+4I94rN3C6fjEKrj3OvoZ12LXsF6tjqXgUk2sKTzDGLATicktqnJlY3H304osvEhgYyOTJ/z7Gv3PnTjZu3Bjtel5eXmHFZN68eTHK5ebmxsKFC5kxYwYzZ858av758+fx8PCgW7dudO3alT179jwxf/fu3YwcOZLffvsNJyfb/5oqVaqwdetWTp2yvdLiwYMHnDhxAn9/f+7cuUPjxo0ZM2YM+/bte2p/xYoVC1sPYOTIkaROnZpevXo9sZwxhgMHDkTauonoxIkTz32XlUoe3DP+h6bj5hMw7gtCnYR0H3zL0h7NuHvzitXRVDxwedYC9ncqPOYEVCDcw2ZJnYiwYMEC3nvvPb7++mtSp06Nl5cXY8aM4eLFi1GuN2DAALp27cqwYcOeuG3zWdKlS8fSpUvx8fEhXbp0tGjRImzehg0bGDFiBK6urri7uz/VUvj++++5efNm2AXmChUq8PPPPzNt2jTat29PYGAgAEOGDCF9+vS0aNGCgIAAjDGMHv30y/Bq1arFBx98gDEGEeGLL74gd+7cT3z5Z86cmREjRlC2bFlEInuu8Enr16+nSZMmMT4eKvkqW78D96s2YkP/t8i/4hAHG9XD9dN3qNTyLaujqTiISd9H4duGwdg6x5tsjLH87Wvadfaz9enTh2bNmkV7V9GQIUMoVKgQ7dq1i3ZbgYGB1K5dmy1btuDi8sy/J1QKcmDDXG59ORgPv0ecruFFja9/JlNWT6tjqShE1/dRknqfQkyJSDOgWaFChbqdPHnS6jhJ2tWrV/n7779p3rx5nLd18uRJLl68GHZ7rlLhPXxwlz8HvIXX0n3cdXeCj3tQ9ZV3rI6lIvFcRUFE+kezTWOM+So+wsWFthSUSnoOb1nEtS8GkONKIGeq5KHatz+T2SOv1bFUOM/bS+r9SAaArthehqOUUk8pUaMF1VZs4WzrCuTdcYHjjRuxddYoq2OpGIqyKBhjvns8AJOANEBnYDZQIJHyKaUcUKo07jQe9isuU0fxIIMbWQZNZtkbL3Hj8lmro6lniPaWVBHJIiJDgAPY7lQqZ4z5xOqLzNpLqlKOoViVRtRYsZWzbauQd/clTjVtwqYZ31gdS0UjyqIgIiOw9Wh6DyhljBlojLmVaMmiEZvnFJRS1nJLlZbGg3/BbfpY7mVOTbZh01jWvg7XfE9YHU1FIrqWwgdALuAL4FK4ri7uJaduLhxBly5d8PDweOqhsfHjx1OkSBFKlCjBxx9/bFE6pWKmSMX61F7+F+deq0GeA1c536wlG6d8FdZzgEoaorum4GSMSROhm4sMj8cTM2RK16lTJ1auXPnEtPXr17No0SIOHDjA4cOHI+0lVamkxtUtNY36Tybtbz9y2yMtHiNmsrJdba6eP2p1NGUX624uVOKrVasWWbJkeWLajz/+SL9+/cL6YfLw8LAimlLPpVDZOtRd+hf/dKyL55HrXGjxMut/6q+thiTAIR9LDffwWqLud9CSwxy5FL9nzornysCAZiVivd6JEyfYvHkzn3/+OalTp2bkyJFUrFgxXrMplZBcXN1o8OkPnGm6hVMf9yXP6D9YtWodZUZOJFeBUlbHS7EcsqWgF5ohODiYW7dusX37dkaMGEHbtm1xxKfTlSpQqgYvLd3GhTcbkPPkTS63asu68f201WARh2wpWOV5/qJPKLlz56Z169aICJUqVcLJyYnr168/9f4EpRyBs7ML9T8cw/kmOzj28TvknbCI1Ws3UnLEBHIXLmd1vBTFIVsKClq2bMmff/4J2E4lBQUFRfraTaUcSb5ilai3cCsXezTF4+xtrrXpwNpRHxASEmx1tBRDi4IDaN++PVWrVuX48ePkzp2bKVOm0KVLF86cOUPJkiVp164d06dPj1HX10oldc7OLtR7bwQe82ZyLX8mPCctZ22L6pw/+rfV0VIEh+wl9THtEE+p5C00NJT1P3xOpskLcQ6Fa/+tz0vvf4ezs575jovn7RBPKaUs5eTkxEu9h5NzwRwuF85CnimrWde0KmcObrE6WrLlkEVB+z5SKmXJVaAUDf7YzJX325Ll8n3ute/GqmFvE/woyOpoyY5DFgW9JVWplMfJyYm63QeRZ9E8LhbPRt4Z61nfpCqn9m6wOlqy4pBFQSmVcmXPV4yGszdw7eMOZPJ7yIPXe7JicDceBQVYHS1Z0KKglHI4Tk5O1O7yBfmWLORCmRx4zdzCpsbVOL5ztdXRHJ4WBaWUw/LIXZgmM9dz/fPOuN8KIKhjH5Z/2YmgwAdWR3NYWhQcQEBAAJUqVaJMmTKUKFGCAQMGAPDRRx9RtGhRSpcuTatWrbh9+7a1QZWySM03PqbQ0mWcr+BJ/j/+ZkvD6hzdttzqWA5Ji4IDSJUqFX/++Sf79+9n3759rFy5ku3bt+Pj48OhQ4c4cOAAhQsXZvjw4VZHVcoy/8mZn6Yz1nJzYHfS3gsiuOsHLP/0DYIe3n/2yiqMFgUHICK4u7sD8OjRIx49eoSIUL9+fVxcbA/xVKlSBV9fXytjKpUkVG/XlyIrVnK+Sl7yL9jFXw2rcXjzIqtjOQx9LDA2VvSDKwfjd5s5SkGjr5+5WEhICOXLl+fUqVP06tWLypUrPzF/6tSpvPrqq/GbTSkHlTlbHppOXcX2eRNw++YH6N6PZU1n8eKgSaRJq+8Ii462FByEs7Mz+/btw9fXlx07dnDo0KGweUOHDsXFxYUOHTpYmFCppKfKy70osWItZ2vmp8CS/extVJcHe/daHStJc8i+j8K9ZKfbyZMnrY6T6AYNGkS6dOn48MMPmT59OhMnTmTdunWkTZvW6mhKJVk7Fk/G/bsZyLUbZOnYkWx93sUpTRqrY1ki2fV9lNKeaPbz8wu7s+jhw4esXbuWokWLsnLlSr755hsWL16sBUGpZ6jUvBtFlq0kU7tXuTltGmdatuSBdqj5FIcsCinN5cuXqVu3LqVLl6ZixYr4+PjQtGlTevfuzb179/Dx8cHb25sePXpYHVWpJM3ZPR05Bwwg77RpEBLK+Tf+y5UhQwl9oM81POaQp48e066zlVLPK/TBA66NHsOt337D1dOTnEOGkK5K5WevmAwku9NHSikVV05p05Lj88/I99uviLMz/3TqxOWBAwnxT9nPNWhRUEqlaGnLlyf/wgVk6dyZ27/P4UzzZvhv2Wp1LMtoUVBKpXhOadKQ/ZOP8Zo1E6fUabjw5ptc+vxzQu7etTpaotOioJRSdmm8vcm/YD7/6daNOwsWcqZZc+5t2GB1rESlRUEppcJxSpUKjw/ex+v333HOkAHfHj259Ek/QlJIh5NaFJRSKhJpSpXEa95csr7dkzvLlnG6WTPurVtndawEp0XBQdy+fZs2bdpQtGhRihUrxrZt28LmjRw5EhHh+vXrFiZUKvlxcnMj27vvkv+PObj8Jyu+vXpz8f0PCL51y+poCUaLgoPo06cPDRs25NixY+zfv59ixYoBcOHCBdasWUPevHktTqhU8pW6WDHy/zGHrO++w901azjTpCl3V660OlaC0KLgAO7evcumTZvo2rUrAG5ubmTKlAmAvn378u233yIiFiZUKvkTV1eyvf02+efOxTVnTi6+1xffd/sQnMxa6Emm62wRKQb0AbIC64wxP1oc6Snf7PiGYzePxes2i2YpyieVPol2mTNnzpAtWzY6d+7M/v37KV++PGPHjmXdunV4enpSpkyZeM2klIpa6iKF8fp9Njem/sL18eM5s2MH2T//nAxNmySLP84StKUgIlNF5JqIHIowvaGIHBeRUyLSD8AYc9QY0wNoC0T6+HVKFRwczJ49e+jZsyd79+4lXbp0DBw4kKFDhzJ48GCr4ymV4oiLC1m7dyP/wgW45cvHpY8+wrdXbx5du2Z1tDhL0L6PRKQW4A/MMMaUtE9zBk4APoAvsBNob4w5IiLNgX7A98aYmc/afkrp++jKlStUqVKFc+fOAbB582YGDhzIwYMHw3pH9fX1JVeuXOzYsYMcOXJYmFaplMWEhHBz+gz8xo5FUqUi+6efkrFliyTdarCs7yNjzCbgZoTJlYBTxpgzxpggYDbQwr78YmNMNUDfFhNOjhw5yJMnD8ePHwdg3bp1lCtXjmvXrnHu3DnOnTtH7ty52bNnjxYEpRKZODvzny6dyb9wAaleeIHLn37KhR49eHTlitXRnosVF5o9gQvhxn0BTxGpIyLjROQnYHlUK4tIdxHZJSK7/Pz8EjprkjF+/Hg6dOhA6dKl2bdvH5999pnVkZRS4aTKn598v84g+2ef8WDHTs40bcatP/7A0XqiTvCus0XEC1ga7vTRK0ADY8yb9vE3gErGmHdiu+2UcvpIKeVYgv75h8tffMmDHTtIV60aOb8ajKunp9WxwiS1rrN9gTzhxnMDl2KzARFpJiKT7ty5E6/BlFIqPrjlzUveab+QY0B/Hu7bx5lmzbk1axYmNNTqaM9kRVHYCbwgIvlFxA1oByyOzQZS2us4lVKOR5ycyNy+PQWWLCaNdxmuDBrMP506E3ThwrNXtlBC35I6C9gGFBERXxHpaowJBnoDq4CjwBxjzOGEzKGUUlZx9fQkz5Qp5PhqMAFHjnCmeQtuzvg1ybYaHPJ1nCLSDGhWqFChbidPnrQ6jlJKxcijy5e5PGAA9zdtJk25cuQcOoRU+fMneo6kdk0hzvT0kVLKEbnmzEmen34i5/DhBJ46xdmWrbgx9RdMSIjV0cI4ZFFQSilHJSJkatWSAkuWkK56da59+y3nX+tA4OnTVkcDtCg4hC5duuDh4UHJkiXDpn300UcULVqU0qVL06pVK27bXwDy6NEjOnbsSKlSpShWrBjDhw+3KLVSKjqu2T3IPeF7co0YQdC5c5xt1ZrrkyZjgoMtzeWQRSGl3ZLaqVMnVkboptfHx4dDhw5x4MABChcuHPbl/8cffxAYGMjBgwfZvXs3P/30U1j3GEqppEVEyNisKQWWLcW9Th38Ro3i3KvtCDh+wrJMDlkUUto1hVq1apElS5YnptWvXx8XF1snt1WqVMHX1xew/ZLdv3+f4OBgHj58iJubGxkyZEj0zEqpmHPJmpXc48biOWY0jy5f5mybNvhNmIB59CjxsyT6Hh3YlWHDCDwav11npypWlBxx7LJi6tSpvPrqqwC0adOGRYsWkTNnTh48eMDo0aOfKihKqaQpQ8OGpK1UiatDhnJ9/PfcW7OWXMOGkrp48UTL4JAthZR2+ig6Q4cOxcXFhQ4dbH0I7tixA2dnZy5dusTZs2f57rvvOHPmjMUplVIx5ZIlC56jviP39+MJvn6ds21f5drYsYQGBSXO/hNlL/HMGLMEWFKhQoVuibnfuP5FH9+mT5/O0qVLWbduXVg3vTNnzqRhw4a4urri4eFB9erV2bVrFwUKFLA4rVIqNtLXq0faChW4Onw4N36ciP/ateQcNow0pUol6H4dsqWgYOXKlXzzzTcsXrw47J0KAHnz5uXPP//EGMP9+/fZvn07RYsWtTCpUup5OWfKRK5vviH3jz8Qcucu515tx7XvviM0MDDB9qlFwQG0b9+eqlWrcvz4cXLnzs2UKVPo3bs39+7dw8fHB29vb3r06AFAr1698Pf3p2TJklSsWJHOnTtTunRpiz+BUiou0tetS4GlS8jYqiU3Jv/M2VateXjgQILsyyG7uXhMu85WSqU0/pu3cHlAf3IOGox7zRrPtY3ourlwyGsK4fo+sjqKUkolKveaNSi4ciVObm4Jsn2HPH2U0p5TUEqp8BKqIICDFgWllFIJQ4uCUkqpMFoUlFJKhXHIoqBPNCulVMJwyKKgF5qVUiphOGRRUEoplTC0KCillArj0E80i4gfcD6SWRmBqC44RDUv4vTIlos4LStwPUZhn190nyW+1n3WcvFxPCObllyPZ0yWjc1xi2q6Hs+YzY/JsYvJtMQ4nlHliO918xljskU6xxiT7AZgUmznRZwe2XKRLLPLys8SX+s+a7n4OJ4xOcbJ5XjG5Zjq8Yzf4xnTYxeTaYlxPBPzmEY1JNfTR0ueY17E6ZEtF912E0pc9hnTdZ+1XHwcz8imJdfjGZNlY3PcopquxzNm82P6bzkl/ZuPkkOfPrKaiOwyUXQqpWJPj2f80uMZv1LK8UyuLYXEMsnqAMmMHs/4pcczfqWI46ktBaWUUmG0paCUUiqMFgWllFJhtCgopZQKo0UhnohISxGZLCKLRKS+1XkcnYgUE5GJIjJXRHpanSe5EJF0IrJbRJpancXRiUgdEdls/z2tY3We+KJFIRoiMlVEronIoQjTG4rIcRE5JSL9AIwxC40x3YBOwKsWxE3yYnk8jxpjegBtgWR/G+Dzis0xtfsEmJO4KR1HLI+nAfyB1IBvYmdNKFoUojcNaBh+gog4AxOARkBxoL2IFA+3yBf2+epp04jF8RSR5sAWYF3ixnQo04jhMRWResAR4Gpih3Qg04j57+hmY0wjbIV2UCLnTDBaFKJhjNkE3IwwuRJwyhhzxhgTBMwGWojNN8AKY8yexM7qCGJzPO3LLzbGVAM6JG5SxxHLY1oXqAK8BnQTEf33H0FsjqcxJtQ+/xaQKhFjJigXqwM4IE/gQrhxX6Ay8A5QD8goIoWMMROtCOeAIj2e9nO0rbH9Y1ue+LEcWqTH1BjTG0BEOgHXw32pqehF9TvaGmgAZAK+tyBXgtCiEHsSyTRjjBkHjEvsMMlAVMdzA7AhcaMkG5Ee07AfjJmWeFGShah+R+cD8xM7TELT5mPs+QJ5wo3nBi5ZlCU50OMZ//SYxq8UdTy1KMTeTuAFEckvIm5AO2CxxZkcmR7P+KfHNH6lqOOpRSEaIjIL2AYUERFfEelqjAkGegOrgKPAHGPMYStzOgo9nvFPj2n80uOpHeIppZQKR1sKSimlwmhRUEopFUaLglJKqTBaFJRSSoXRoqCUUiqMFgWllFJhtCgopZQKo0VBOSwRySEis0XktIgcEZHlIlLY/vKTpfG0D6+IfevHYt3lIpIpHjKMEZFa9p83iMg/IiLh5i8UEf9weR+KyD4R2S8if4lIEfu8UiIyLa55VPKmRUE5JPuX4gJggzGmoDGmOPAZkN3aZP8yxjQ2xtyOyzZEJAtQxd6l82O3ger2+ZmAnBFWO22M8TbGlAGmYzsuGGMOArlFJG9cMqnkTYuCclR1gUfhuyg3xuwzxmy2j7rbX+V5TET+9/gvaxE5JyJZ7T9XEJEN9p8H2t+6tUFEzojIuxF3KCIFRGSviFSMMD2niGyy/3V+SERqht+XiPSwz9snImdFZL19fn0R2SYie0TkDxFxj+RztgFWRpg2G1v/O2DrXjy6njozYOvv/7El4dZV6ilaFJSjKgnsjmZ+WeA9bG/KKoD9L+tnKIqtf/xKwAARcX08w34KZh7Q2RizM8J6rwGrjDHeQBlgX/iZxpiJ9nkVsfW4OcpemL4A6hljygG7gPcjyVQ9ks+5DqhlfyNYO+D3CPML2gvQafs2R4WbtwuoGdUBUErfp6CSqx3GGF8AEdkHeGF7tWd0lhljAoFAEbnGv6eisgGLgJej6AhtJzDVXkQWGmP2RbH9scCfxpglItIUW8Haam/EuGHriC2inIBfhGkh9s/yKpDGGHMu3CUGsJ8+AhCRV4FJ/PuKyWtArijyKaUtBeWwDgPlo5kfGO7nEP79AyiYf3/vU8dwnTvY3rwVaWvDfr6/FnAR+FVE/htxGfvbzvLx77t8BVhjP/fvbYwpbozpGsnmH0aSE2ynkMYDcyLLFM5ie7bHUtu3qVSktCgoR/UnkEpEuj2eICIVRaT2M9Y7x7/F5OUY7isIaAn8V0ReizhTRPIB14wxk4EpQLkI88sDHwKvh3sF5naguogUsi+TVkQKR7Lvo0ChSKZvBoYDs56RvQZwOtx4YeC57qZSKYMWBeWQjK3P91aAj/2W1MPAQJ79RqxBwFgR2YytNRDT/d0HmgJ9RaRFhNl1gH0ishdboRkbYX5vIAuw3n6u/2djjB/QCZglIgewFYmikex6mX37EfMYY8xIY8z1SNZ5fE1hPzAMeDPcvLr2bSoVKX2fglJJnIhsAZrGw+2tqYCNQA37i2OUeooWBaWSOBGpDDw0xhyI43ZeADyNMRviJZhKlrQoKKWUCqPXFJRSSoXRoqCUUiqMFgWllFJhtCgopZQKo0VBKaVUmP8DqsZePSh0ovoAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "##### %matplotlib notebook\n", + "from ipywidgets import *\n", + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "import math\n", + "\n", + "x = 1024\n", + "y = 1024\n", + "z = 1024\n", + "t = 100\n", + "c = 1\n", + "chunkSizes = (16, 32, 64, 128)\n", + "\n", + "fig = plt.figure()\n", + "ax = fig.add_subplot(1, 1, 1)\n", + "plt.xlabel('Chunk size (XYZ)')\n", + "plt.ylabel('Number of chunks')\n", + "plt.title(\"XYZCT: 1k x 1k x 1k x 1 x 100\");\n", + "\n", + "def file_count(i, chunkT=1):\n", + " return math.ceil(x/i) * math.ceil(y/i) * math.ceil(z/i) * math.ceil(t/chunkT) * math.ceil(c/1)\n", + "\n", + "for chunkT in (1, 10, 100):\n", + " numFiles = []\n", + " fileSize = []\n", + " for i in chunkSizes:\n", + " numFiles.__iadd__([file_count(i, chunkT)])\n", + " fileSize.__iadd__([i])\n", + " line, = ax.plot(fileSize, numFiles, label=f\"{chunkT}\")\n", + "\n", + "plt.legend(loc=\"lower left\", title=\"Chunk sizes (T)\", frameon=False)\n", + "\n", + "def scale(scale):\n", + " plt.yscale(scale)\n", + " plt.xscale(scale)\n", + "\n", + "scale('log')\n", + "ax.annotate('Chosen chunk size\\nT: 1, X: 64, Y: 64, Z: 64',\n", + " xy=((64), file_count(64, 1)), xycoords='data',\n", + " xytext=(0, 40), textcoords='offset points',\n", + " arrowprops=dict(facecolor='black', shrink=0.05),\n", + " horizontalalignment='left', verticalalignment='bottom')\n", + "fig.savefig('/tmp/3d-chunks.png', dpi=800)" ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "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.7.6" } - ], - "source": [ - "##### %matplotlib notebook\n", - "from ipywidgets import *\n", - "import numpy as np\n", - "import matplotlib.pyplot as plt\n", - "import math\n", - "\n", - "x = 1024\n", - "y = 1024\n", - "z = 1024\n", - "t = 100\n", - "c = 1\n", - "chunkSizes = (16, 32, 64, 128)\n", - "\n", - "fig = plt.figure()\n", - "ax = fig.add_subplot(1, 1, 1)\n", - "plt.xlabel('Chunk size (XYZ)')\n", - "plt.ylabel('Number of chunks')\n", - "plt.title(\"XYZCT: 1k x 1k x 1k x 1 x 100\");\n", - "\n", - "def file_count(i, chunkT=1):\n", - " return math.ceil(x/i) * math.ceil(y/i) * math.ceil(z/i) * math.ceil(t/chunkT) * math.ceil(c/1)\n", - "\n", - "for chunkT in (1, 10, 100):\n", - " numFiles = []\n", - " fileSize = []\n", - " for i in chunkSizes:\n", - " numFiles.__iadd__([file_count(i, chunkT)])\n", - " fileSize.__iadd__([i])\n", - " line, = ax.plot(fileSize, numFiles, label=f\"{chunkT}\")\n", - "\n", - "plt.legend(loc=\"lower left\", title=\"Chunk sizes (T)\", frameon=False)\n", - "\n", - "def scale(scale):\n", - " plt.yscale(scale)\n", - " plt.xscale(scale)\n", - "\n", - "scale('log')\n", - "ax.annotate('Chosen chunk size\\nT: 1, X: 64, Y: 64, Z: 64',\n", - " xy=((64), file_count(64, 1)), xycoords='data',\n", - " xytext=(0, 40), textcoords='offset points',\n", - " arrowprops=dict(facecolor='black', shrink=0.05),\n", - " horizontalalignment='left', verticalalignment='bottom')\n", - "fig.savefig('/tmp/3d-chunks.png', dpi=800)" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [] - } - ], - "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.7.6" - } - }, - "nbformat": 4, - "nbformat_minor": 2 + "nbformat": 4, + "nbformat_minor": 2 } diff --git a/results/std_mean.py b/results/std_mean.py index 7c5a3b5..b58d101 100644 --- a/results/std_mean.py +++ b/results/std_mean.py @@ -5,7 +5,6 @@ # source (e.g. http, local, s3) for csv_file in ["2d_benchmark_data.csv", "3d_benchmark_data.csv"]: - print(csv_file) df = pandas.read_csv(csv_file)