From f5bf0fc4abb0d2ad5d6146df63fa9b4f8529f717 Mon Sep 17 00:00:00 2001 From: Chichi Lalescu <clalesc1@jhu.edu> Date: Wed, 16 Dec 2015 21:41:46 +0100 Subject: [PATCH] add plot of smooth Fourier filter --- meta/smooth dealiasing.ipynb | 804 +++++++++++++++++++++++++++++++++++ 1 file changed, 804 insertions(+) create mode 100644 meta/smooth dealiasing.ipynb diff --git a/meta/smooth dealiasing.ipynb b/meta/smooth dealiasing.ipynb new file mode 100644 index 00000000..664bc0bf --- /dev/null +++ b/meta/smooth dealiasing.ipynb @@ -0,0 +1,804 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "data": { + "application/javascript": [ + "/* Put everything inside the global mpl namespace */\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('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", + "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 = $('<div/>');\n", + " this._root_extra_style(this.root)\n", + " this.root.attr('style', 'display: inline-block');\n", + "\n", + " $(parent_element).append(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", + " 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", + " fig.waiting = false;\n", + " };\n", + "\n", + " this.imageObj.onunload = function() {\n", + " this.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 = $(\n", + " '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n", + " 'ui-helper-clearfix\"/>');\n", + " var titletext = $(\n", + " '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n", + " 'text-align: center; padding: 3px;\"/>');\n", + " titlebar.append(titletext)\n", + " this.root.append(titlebar);\n", + " this.header = titletext[0];\n", + "}\n", + "\n", + "\n", + "\n", + "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n", + "\n", + "}\n", + "\n", + "\n", + "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n", + "\n", + "}\n", + "\n", + "mpl.figure.prototype._init_canvas = function() {\n", + " var fig = this;\n", + "\n", + " var canvas_div = $('<div/>');\n", + "\n", + " canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n", + "\n", + " function canvas_keyboard_event(event) {\n", + " return fig.key_event(event, event['data']);\n", + " }\n", + "\n", + " canvas_div.keydown('key_press', canvas_keyboard_event);\n", + " canvas_div.keyup('key_release', canvas_keyboard_event);\n", + " this.canvas_div = canvas_div\n", + " this._canvas_extra_style(canvas_div)\n", + " this.root.append(canvas_div);\n", + "\n", + " var canvas = $('<canvas/>');\n", + " canvas.addClass('mpl-canvas');\n", + " canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n", + "\n", + " this.canvas = canvas[0];\n", + " this.context = canvas[0].getContext(\"2d\");\n", + "\n", + " var rubberband = $('<canvas/>');\n", + " rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n", + "\n", + " var pass_mouse_events = true;\n", + "\n", + " canvas_div.resizable({\n", + " start: function(event, ui) {\n", + " pass_mouse_events = false;\n", + " },\n", + " resize: function(event, ui) {\n", + " fig.request_resize(ui.size.width, ui.size.height);\n", + " },\n", + " stop: function(event, ui) {\n", + " pass_mouse_events = true;\n", + " fig.request_resize(ui.size.width, ui.size.height);\n", + " },\n", + " });\n", + "\n", + " function mouse_event_fn(event) {\n", + " if (pass_mouse_events)\n", + " return fig.mouse_event(event, event['data']);\n", + " }\n", + "\n", + " rubberband.mousedown('button_press', mouse_event_fn);\n", + " rubberband.mouseup('button_release', mouse_event_fn);\n", + " // Throttle sequential mouse events to 1 every 20ms.\n", + " rubberband.mousemove('motion_notify', mouse_event_fn);\n", + "\n", + " rubberband.mouseenter('figure_enter', mouse_event_fn);\n", + " rubberband.mouseleave('figure_leave', mouse_event_fn);\n", + "\n", + " canvas_div.on(\"wheel\", function (event) {\n", + " event = event.originalEvent;\n", + " event['data'] = 'scroll'\n", + " if (event.deltaY < 0) {\n", + " event.step = 1;\n", + " } else {\n", + " event.step = -1;\n", + " }\n", + " mouse_event_fn(event);\n", + " });\n", + "\n", + " canvas_div.append(canvas);\n", + " canvas_div.append(rubberband);\n", + "\n", + " this.rubberband = rubberband;\n", + " this.rubberband_canvas = rubberband[0];\n", + " this.rubberband_context = rubberband[0].getContext(\"2d\");\n", + " this.rubberband_context.strokeStyle = \"#000000\";\n", + "\n", + " this._resize_canvas = function(width, height) {\n", + " // Keep the size of the canvas, canvas container, and rubber band\n", + " // canvas in synch.\n", + " canvas_div.css('width', width)\n", + " canvas_div.css('height', height)\n", + "\n", + " canvas.attr('width', width);\n", + " canvas.attr('height', height);\n", + "\n", + " rubberband.attr('width', width);\n", + " rubberband.attr('height', height);\n", + " }\n", + "\n", + " // Set the figure to an initial 600x600px, this will subsequently be updated\n", + " // upon first draw.\n", + " this._resize_canvas(600, 600);\n", + "\n", + " // Disable right mouse context menu.\n", + " $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\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 nav_element = $('<div/>')\n", + " nav_element.attr('style', 'width: 100%');\n", + " this.root.append(nav_element);\n", + "\n", + " // Define a callback function for later on.\n", + " function toolbar_event(event) {\n", + " return fig.toolbar_button_onclick(event['data']);\n", + " }\n", + " function toolbar_mouse_event(event) {\n", + " return fig.toolbar_button_onmouseover(event['data']);\n", + " }\n", + "\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", + " // put a spacer in here.\n", + " continue;\n", + " }\n", + " var button = $('<button/>');\n", + " button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n", + " 'ui-button-icon-only');\n", + " button.attr('role', 'button');\n", + " button.attr('aria-disabled', 'false');\n", + " button.click(method_name, toolbar_event);\n", + " button.mouseover(tooltip, toolbar_mouse_event);\n", + "\n", + " var icon_img = $('<span/>');\n", + " icon_img.addClass('ui-button-icon-primary ui-icon');\n", + " icon_img.addClass(image);\n", + " icon_img.addClass('ui-corner-all');\n", + "\n", + " var tooltip_span = $('<span/>');\n", + " tooltip_span.addClass('ui-button-text');\n", + " tooltip_span.html(tooltip);\n", + "\n", + " button.append(icon_img);\n", + " button.append(tooltip_span);\n", + "\n", + " nav_element.append(button);\n", + " }\n", + "\n", + " var fmt_picker_span = $('<span/>');\n", + "\n", + " var fmt_picker = $('<select/>');\n", + " fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n", + " fmt_picker_span.append(fmt_picker);\n", + " nav_element.append(fmt_picker_span);\n", + " this.format_dropdown = fmt_picker[0];\n", + "\n", + " for (var ind in mpl.extensions) {\n", + " var fmt = mpl.extensions[ind];\n", + " var option = $(\n", + " '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n", + " fmt_picker.append(option)\n", + " }\n", + "\n", + " // Add hover states to the ui-buttons\n", + " $( \".ui-button\" ).hover(\n", + " function() { $(this).addClass(\"ui-state-hover\");},\n", + " function() { $(this).removeClass(\"ui-state-hover\");}\n", + " );\n", + "\n", + " var status_bar = $('<span class=\"mpl-message\"/>');\n", + " nav_element.append(status_bar);\n", + " this.message = status_bar[0];\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", + "\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", + "\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]);\n", + " fig.send_message(\"refresh\", {});\n", + " };\n", + "}\n", + "\n", + "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n", + " var x0 = msg['x0'];\n", + " var y0 = fig.canvas.height - msg['y0'];\n", + " var x1 = msg['x1'];\n", + " var y1 = fig.canvas.height - msg['y1'];\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, 0, fig.canvas.width, fig.canvas.height);\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", + " {\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.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", + " fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n", + " evt.data);\n", + " fig.updated_canvas_event();\n", + " return;\n", + " }\n", + " else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n", + " fig.imageObj.src = evt.data;\n", + " fig.updated_canvas_event();\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(\"No handler for the '\" + msg_type + \"' message type: \", msg);\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(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\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", + " if (e.target)\n", + " targ = e.target;\n", + " else if (e.srcElement)\n", + " targ = e.srcElement;\n", + " if (targ.nodeType == 3) // defeat Safari bug\n", + " targ = targ.parentNode;\n", + "\n", + " // jQuery normalizes the pageX and pageY\n", + " // pageX,Y are the mouse positions relative to the document\n", + " // offset() returns the position of the element relative to the document\n", + " var x = e.pageX - $(targ).offset().left;\n", + " var y = e.pageY - $(targ).offset().top;\n", + "\n", + " return {\"x\": x, \"y\": y};\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", + " {\n", + " this.canvas.focus();\n", + " this.canvas_div.focus();\n", + " }\n", + "\n", + " var x = canvas_pos.x;\n", + " var y = canvas_pos.y;\n", + "\n", + " this.send_message(name, {x: x, y: y, button: event.button,\n", + " step: event.step});\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", + "\n", + " // Prevent repeat events\n", + " if (name == 'key_press')\n", + " {\n", + " if (event.which === this._key)\n", + " return;\n", + " else\n", + " this._key = event.which;\n", + " }\n", + " if (name == 'key_release')\n", + " this._key = null;\n", + "\n", + " var value = '';\n", + " if (event.ctrlKey && event.which != 17)\n", + " value += \"ctrl+\";\n", + " if (event.altKey && event.which != 18)\n", + " value += \"alt+\";\n", + " if (event.shiftKey && event.which != 16)\n", + " value += \"shift+\";\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});\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", + "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\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"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\";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 overriden (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 = $(\"#\" + id);\n", + " var ws_proxy = comm_websocket_adapter(comm)\n", + "\n", + " function ondownload(figure, format) {\n", + " window.open(figure.imageObj.src);\n", + " }\n", + "\n", + " var fig = new mpl.figure(id, ws_proxy,\n", + " ondownload,\n", + " element.get(0));\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.get(0);\n", + " fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n", + " if (!fig.cell_info) {\n", + " console.error(\"Failed to find cell for figure\", id, fig);\n", + " return;\n", + " }\n", + "\n", + " var output_index = fig.cell_info[2]\n", + " var cell = fig.cell_info[0];\n", + "\n", + "};\n", + "\n", + "mpl.figure.prototype.handle_close = function(fig, msg) {\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).html('<img src=\"' + dataURL + '\">');\n", + " fig.send_message('closing', {});\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 dataURL = this.canvas.toDataURL();\n", + " this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\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 () { fig.push_to_output() }, 1000);\n", + "}\n", + "\n", + "mpl.figure.prototype._init_toolbar = function() {\n", + " var fig = this;\n", + "\n", + " var nav_element = $('<div/>')\n", + " nav_element.attr('style', 'width: 100%');\n", + " this.root.append(nav_element);\n", + "\n", + " // Define a callback function for later on.\n", + " function toolbar_event(event) {\n", + " return fig.toolbar_button_onclick(event['data']);\n", + " }\n", + " function toolbar_mouse_event(event) {\n", + " return fig.toolbar_button_onmouseover(event['data']);\n", + " }\n", + "\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) { continue; };\n", + "\n", + " var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n", + " button.click(method_name, toolbar_event);\n", + " button.mouseover(tooltip, toolbar_mouse_event);\n", + " nav_element.append(button);\n", + " }\n", + "\n", + " // Add the status bar.\n", + " var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n", + " nav_element.append(status_bar);\n", + " this.message = status_bar[0];\n", + "\n", + " // Add the close button to the window.\n", + " var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n", + " var button = $('<button class=\"btn btn-mini btn-danger\" href=\"#\" title=\"Close figure\"><i class=\"fa fa-times icon-remove icon-large\"></i></button>');\n", + " button.click(function (evt) { fig.handle_close(fig, {}); } );\n", + " button.mouseover('Close figure', toolbar_mouse_event);\n", + " buttongrp.append(button);\n", + " var titlebar = this.root.find($('.ui-dialog-titlebar'));\n", + " titlebar.prepend(buttongrp);\n", + "}\n", + "\n", + "\n", + "mpl.figure.prototype._canvas_extra_style = function(el){\n", + " // this is important to make the div 'focusable\n", + " el.attr('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", + " }\n", + " else {\n", + " // location in version 2\n", + " IPython.keyboard_manager.register_events(el);\n", + " }\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", + " // Check for shift+enter\n", + " if (event.shiftKey && event.which == 13) {\n", + " this.canvas_div.blur();\n", + " event.shiftKey = false;\n", + " // Send a \"J\" for go to next cell\n", + " event.which = 74;\n", + " event.keyCode = 74;\n", + " manager.command_mode();\n", + " manager.handle_keydown(event);\n", + " }\n", + "}\n", + "\n", + "mpl.figure.prototype.handle_save = function(fig, msg) {\n", + " fig.ondownload(fig, null);\n", + "}\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('matplotlib', mpl.mpl_figure_comm);\n", + "}\n" + ], + "text/plain": [ + "<IPython.core.display.Javascript object>" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/html": [ + "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAigAAAGeCAYAAAC+dvpwAAAgAElEQVR4nO3de3Bc9XXA8WNjAZIJkmlhkhgMUoCSlATjADFtCQ/bNNBCEts8ZgKB4Ac0NGni2LxKqHGUIAdCeDS23IHCDBRqE14ltFhmwkybkIDtAgFiaiRqhvBoAJmXwBbW/vrH0cZ67Gof95577u79fmZ2gJWufOZnc/X17+5eiQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABIu2ki0lHhMW0i0ikiSwaPXRL3UAAAILtmishWEVlZwTEtItItInsOea5DKo8cAACAYQ4XkS7RqOiWygJluYisHvFcs4jkZHi0AAAAVG2DVBYoPSKyuMDzORGZE8tEAAAg8yoJlBbREJlf4GO9FXwdAACAMVUSKG1SPFC2isjauIYCAADZVkmgTBMNlNkFPtYjIuvjGgoAAGQbgQIAAFKHSzwAACB1qnmRbLF38ZTzdXIiEkSkjwcPHjx48OBR8SOIfi+te5W+zbhbRK4q8HyxSz8jhVKP8eMbwoQJTTxifLCmrHk5j112qf3HuHEN7jNU+xg/fvRj3LjhD5HqH+PG6decMKEp7LprU9h996bQ1BT90dDQEMvX4TF8Tcv4fln3Kg2UDhFZM+K5/KWfcm7U1iciAclqbGz0HiFzWHMfWVn37dtDePPNEF58MYRnnw3hscdCWLcuhDvvDOH660O4/PIQFi4M4UtfCuGznw1hr71CECn8aG0N4dxzQ7j1Vv16lcrKmqdFU1NTZgKlR0YHR17+RbFLhjzXLHrPk9Yhz62SwrsqhRAoDjiBJI8198G6F/fWWyE8+WQId90Vwne/G8Ipp4Sw776jg2Xq1BCWLQth06byvi5rnqx6D5Rm0R/21yUaIAOikdI54vPaRORNGX3ppnXwcxeI7qgUek1KMQSKA04gyWPNfbDulXvppRBuvz2EBQtC2G+/4bHy+c+H8C//EsK2bcWPZ82TVe+B4olAccAJJHmsuQ/WPZpcLoT160O45JIQPv7xnaEyeXIIN94YwgcfjD6GNU8WgWKHQHHACSR5rLkP1j0+H34Ywr33hjBjxs5Q2W8/vUSUy+38PNY8WQSKHQLFASeQ5LHmPlh3G48+GsKsWTtDZdasELZs0Y+x5skiUOwQKA44gSSPNffButvJ5UK4//4Q9t9fI6W5OYQ1a1jzpBEodvoaGhq8f38zp7293XuEzGHNfbDu9vr6Qrjggp27KUcf3R527PCeKjsGA6Xf+5t5Pepramry/v0FAER0zz26iyISwqmnhvDee94TZcNgoPR5fzOvRwQKANSJZ54JYcoUjZTp00N4+23vieofgWKHQAGAOvLKKyEcdphGyp//eQjvvus9UX0jUOwQKABQZ15/PYRPf1ojZebMEPr7vSeqXwSKHQIFAOrQ//1fCAceqJFy/vnD75WC+BAodggUAKhT//M/IUyapJHyj//oPU19IlDsECgAUMd+/vMQxo8PYdddQ3jiCe9p6g+BYodAAYA6d+WVuoty8MG8aDZuBIodAgUA6tyOHSEcd5xGyje/6T1NfSFQ7BAoAJABW7aEMHFiCOPGhfDrX3tPUz8IFDsECgBkxHXX6S7KoYfy1uO4ECh2CBQAyIgdO0I46iiNlBtu8J6mPhAodggUAMiQX/9aA+WP/iiErVu9p6l9BIodAgUAMuaMMzRSLrrIe5LaR6DYIVAAIGNeeEHvi7LbbiG8/LL3NLWNQLFDoABABn3967qLsmiR9yS1jUCxQ6AAQAZt2RLChAkhNDWF8Pvfe09TuwgUOwQKAGTUeefpLsqll3pPUrsIFDsECgBk1ObN+nN6WlpCeO8972lqE4Fih0ABgAz74hd1F2XVKu9JahOBYodAAYAMe/jhnXeXzeW8p6k9BIodAgUAMiyXC+FTn9JIeeQR72lqD4Fih0ABgIxbsUID5fTTvSepPQSKHQIFADLu7bdDaGzUm7f19npPU1sIFDsECgAgfOUruouycqX3JLWFQLFDoAAAQleXBsr06d6T1BYCxQ6BAgAIO3aEsO++GinPPec9Te0gUOwQKACAEILeUVYkhMsv956kdhAodggUAEAIIYSnntJAOeQQ7olSLgLFDoECAAghaJQcdJBGyjPPeE9TGwgUOwQKAOAP8pd5li71nqQ2ECh2CBQAwB9s3Ljz1vcojUCxQ6AAAP4glwuhtZV385SLQLFDoAAAhvnOdzRQrr7ae5L0I1DsECgAgGF+/nMNlOOP954k/QgUOwQKAGCY7dtD+MhHQpgwQX9OD4ojUOwQKACAUWbP1l2Uu+/2niTdCBQ7BAoAYJSbbtJAmTfPe5J0I1DsECgAgFFeflkD5WMf466yYyFQ7BAoAICCDjtMI+Wpp7wnSS8CxQ6BAgAoKP924+uu854kvQgUOwQKAKCgBx/UQDn1VO9J0otAsUOgAAAKeuedEHbZJYTm5hB27PCeJp0IFDsECgCgqOnTdRdl/XrvSdKJQLFDoAAAirrsMg2U5cu9J0knAsUOgQIAKGrdOg2UL3zBe5J0IlDsECgAgKL6+kLYddcQJk4M4cMPvadJHwLFDoECABjT0UfrLsrGjd6TpA+BYodAAQCMafFiDZQbb/SeJH0IFDsECgBgTPfco4Fy5pnek6RPVgKlTUQ6RWSJiHQM/rOSY/PHdAw+ykGgAADG9NprGihTpnhPkj5ZCJQWEekWkT2HPFduaLQU+LwZIrKmjGMJFABASQceqJHy0kvek6RLFgJluYisHvFcs4jkZHi0FDt2aoHnN5RxLIECACjpq1/VQFm92nuSdMlCoPSIyOICz+dEZE6JY9dI4Z2WLiFQAAAxWLVKA+Wb3/SeJF3qPVBaRENkfoGP9YrIyhLHLxk8fuhrVtpEZG0ZvzaBAgAo6emnNVCOPNJ7knSp90Bpk+KBslXKC43uwa+xQUTmir7YthwECgCgpB079GZtu+4awvbt3tOkR70HyjTRuJhd4GM9IrK+jK/RLBonOSnvslAegQIAKMsxx+guyoYN3pOkB4FS2nLRHZgFsjNSytlFIVAAAGVZtEgDpbPTe5L0qPdAiXqJJx8nea2yczdlRolj+xoaGkJjY2PRR3t7u/fvPwAgBe64QwNl/nzvSZLV3t5e9HukiAQR6S/ru30Nyr9Itti7eEq9SLa3yPNdUnoXhR0UAEBZNm/WQDnsMO9J0qPed1BE9EWuVxV4vtiln7w20d2SQsq5WRuBAgAoSy4XQnNzCLvsEsL773tPkw5ZCJQOGR0T+Us/pe5lUmwHZaEUvmw0FIECACjbjBm6i/KrX3lPkg5ZCJRm0dBoHfLcKhm9q5J/Qe3Qe54skNGXctqEW90DAGJ28cUaKDfc4D1JOmQhUEQ0TjpFg6NDCr8mpU1E3pTRl33miAZJp/DDAgEARu66SwPlnHO8J0mHrASKBwIFAFC255/XQJk61XuSdCBQ7BAoAICyDQyEsMceekfZ/n7vafwRKHYIFABARY4+WndRfvMb70n8ESh2CBQAQEUuuEAD5bbbvCfxR6DYIVAAABVZuVIDZfFi70n8ESh2CBQAQEUefVQDZdYs70n8ESh2CBQAQEXefTeEceNC2HtvvbtslhEodggUAEDFDjxQd1FeecV7El8Eih0CBQBQsTlzNFD+4z+8J/FFoNghUAAAFVu2TAPlhz/0nsQXgWKHQAEAVOzuu7nlfQgEiiUCBQBQseee00A58kjvSXwRKHYIFABAxT78UG93P3Gi3v4+qwgUOwQKAKAqn/607qJs2eI9iR8CxQ6BAgCoyhlnaKA8+KD3JH4IFDsECgCgKvl38lx9tfckfggUOwQKAKAqP/2pBsq553pP4odAsUOgAACqsmmTBspRR3lP4odAsUOgAACq8uGHITQ0hLDHHtn9mTwEih0CBQBQtUMP1V2UF1/0nsQHgWKHQAEAVO3007P9M3kIFDsECgCgat/9rgbK9dd7T+KDQLFDoAAAqnbbbRooF17oPYkPAsUOgQIAqNrjj2ugzJrlPYkPAsUOgQIAqNrWrRooU6Z4T+KDQLFDoAAAItlnH42U99/3niR5BIodAgUAEMkxx2ig/OY33pMkj0CxQ6AAACKZN08D5a67vCdJHoFih0ABAESyfLkGyve/7z1J8ggUOwQKACCS++7TQDnnHO9Jkkeg2CFQAACR/Pa3GijTp3tPkjwCxQ6BAgCIZNu2EMaPD2GvvbwnSR6BYodAAQBE9olP6C7K6697T5IsAsUOgQIAiOykkzRQfvlL70mSRaDYIVAAAJH93d9poNxyi/ckySJQ7BAoAIDIVqzQQLn0Uu9JkkWg2CFQAACRrV2rgXL66d6TJItAsUOgAAAie/55DZQjjvCeJFkEih0CBQAQWX+/vtV40iTvSZJFoNghUAAAsTjgAN1F6e31niQ5BIodAgUAEIsZMzRQNmzwniQ5BIodAgUAEIuFCzVQVq/2niQ5BIodAgUAEIuODg2Uq67yniQ5BIodAgUAEIu77tJAmT/fe5LkECh2CBQAQCw2btRAOeEE70mSQ6DYIVAAALF46y0NlP33954kOQSKHQIFABCbvfbS+6Fs3+49STIIFDsECgAgNkceqbsomzd7T5IMAsUOgQIAiM0ZZ2igPPSQ9yTJIFDsECgAgNhcdpkGyk9+4j1JMggUOwQKACA2N9+sgbJokfckySBQ7BAoAIDYPPKIBsqXvuQ9STIIFDsECgAgNi++qIHymc94T5KMrARKm4h0isgSEekY/GeU4xeUcQyBAgCIzY4dIUyYEEJzs/ckychCoLSISLeI7DnkuY7BRzmmiciGIcfPFJHeEV+vEAIFABCr1lbdRdm61XsSe1kIlOUisnrEc80ikpPSkSEislVEZg/577ki8mYZxxIoAIBYHXecBsoTT3hPYi8LgdIjIosLPJ8TkTkljr1INEaqQaAAAGJ17rkaKPfd5z2JvXoPlBbREJlf4GO9IrKyxPE9Mnr3pVwECgAgVkuXaqBcd533JPbqPVDapHigbBWRtWMcm4+blaI7LfkX15b72hUCBQAQq1tu0UD59re9J7FX74EyTTQyZhf4WI+IrB/j2HzcdInI1CHPdww+VwqBAgCIVf5eKF/+svck9giU4vKBMvJzWgefn1Hi1yZQAACxeuEFDZTDD/eexF69B0qUSzwiOy/xFHr+qhLHEigAgFj194cwfnwIe+3lPYm9eg+U/OtIir2Lp9SLZMcKlFIvnu1raGgIjY2NRR/t7e3ev/8AgBqz3366i/LOO96TRNfe3l70e6SIBBHpL/G9tqZ1S+HdjmKXfobaICJrihzLDgoAIHHHHKOB8vTT3pPYqvcdFBF9UevIyMhf+il1s7UFooEzVP51LSeUOJZAAQDE7uyzNVAeeMB7EltZCJRm0XuetA55bpWM3gHJh8fIn9PTLcNfELtKyrs3CoECAIjd5ZdroNx4o/cktrIQKCIaJ52iOyIdUvg1KW2id40dedmnWXQHpnPwUejYQggUAEDsbrpJA2XxYu9JbGUlUDwQKACA2K1bp4Eyd673JLYIFDsECgAgds8/r4Fy5JHek9giUOwQKACA2G3bpoGy997ek9giUOwQKAAAEx//uEbKe+95T2KHQLFDoAAATEyfroGyaZP3JHYIFDsECgDAxGmnaaB0dXlPYodAsUOgAABMLFqkgXLzzd6T2CFQ7BAoAAAT112ngbJ0qfckdggUOwQKAMDE3XdroMyb5z2JHQLFDoECADDx+OMaKLNmeU9ih0CxQ6AAAEy8+qoGyiGHeE9ih0CxQ6AAAEwMDITQ0BDCxIkh5HLe09ggUOwQKAAAM21tuovS2+s9iQ0CxQ6BAgAwc+yxGihPPeU9iQ0CxQ6BAgAwc9ZZGigPPOA9iQ0CxQ6BAgAwc+mlGigrVnhPYoNAsUOgAADMrFihgXLppd6T2CBQ7BAoAAAzDzyggXLWWd6T2CBQ7BAoAAAzTz2lgXLssd6T2CBQ7BAoAAAzvb0aKK2t3pPYIFDsECgAADO5nN6oraFBb9xWbwgUOwQKAMDUIYfoLsqrr3pPEj8CxQ6BAgAwNWuWBsrjj3tPEj8CxQ6BAgAwNW+eBsrdd3tPEj8CxQ6BAgAwtXSpBsqPf+w9SfwIFDsECgDA1M03a6AsWuQ9SfwIFDsECgDAVFeXBsppp3lPEj8CxQ6BAgAwtWmTBsr06d6TxI9AsUOgAABMvfuuBsq++3pPEj8CxQ6BAgAw19wcwi67hLBjh/ck8SJQ7BAoAABzn/qU7qK8/LL3JPEiUOwQKAAAcyeeWJ83ayNQ7BAoAABzX/uaBsq993pPEi8CxQ6BAgAwd/nlGig33ug9SbwIFDsECgDAXGenBsoll3hPEi8CxQ6BAgAw98ADGihnn+09SbwIFDsECgDA3H//twbK8cd7TxIvAsUOgQIAMPf732ugHHyw9yTxIlDsECgAAHO5XAi77hrCxIn67/WCQLFDoAAAEnHAAbqLsnWr9yTxIVDsECgAgET8xV9ooDzzjPck8SFQ7BAoAIBEnHGGBsratd6TxIdAsUOgAAASsWiRBsrNN3tPEh8CxQ6BAgBIxLXXaqAsW+Y9SXwIFDsECgAgEatXa6Ccf773JPEhUOwQKACARPziFxoof/VX3pPEh0CxQ6AAABKxZYsGytSp3pPEh0CxQ6AAABKxfbsGyt57e08SHwLFDoECAEjMPvtopGzb5j1JPAgUOwQKACAx06ZpoLzwgvck8SBQ7BAoAIDEnHKKBsp//Zf3JPEgUOwQKACAxFxwgQbKnXd6TxIPAsUOgQIASEx7uwbKNdd4TxKPrARKm4h0isgSEekY/Ge11onInmV8HoECAEjMLbdooHz7296TxCMLgdIiIt0yPCo6Bh+VmisiOSFQAAAp09WlgXLaad6TxCMLgbJcRFaPeK5Zyg+NodaJyECZxxEoAIDE/Pa3Gih/9mfek8QjC4HSIyKLCzyfE5E5FXydJYOfzw4KACB13n5bA2X//b0niUe9B0qLaFDML/CxXhFZWebXaRORBSIyUwgUAEBK7bFHCBMmhDAw4D1JdPUeKG1SPFC2isjaMr9O/vUqBAoAILUOOUR3UV57zXuS6Oo9UKaJBsXsAh/rEZH1ZXyNmSIydci/EygAgFSaMUMDZeNG70miI1BKG/puHwIFAJBaX/2qBsr993tPEl29B0rUSzwLZHiMECgAgNS67DINlBUrvCeJrt4DJf8i2WLv4hnrRbItooEyVD5Qmsv4tQkUAECifvITDZS//3vvSaKr90AR0Zu0XVXg+WKXfvLmiMiaEY8Ng8etEb0z7Vj6GhoaQmNjY9FHe3u79+8/AKCO3HefBso553hPUp729vai3yNFJIhIf4Xf82tKh2hQDJW/9FPpjdoWVHAcOygAgERt2KCBMnOm9yTRZWEHpVn0nietQ55bJaN3VfIvqB3r5/QsFC7xAABS6tVXNVA++UnvSaLLQqCIaJx0iu6AdEjh16S0icibUviyT/74XtFb3W+Q0j/Lh0ABACRqYEBv1PaRj3hPEl1WAsUDgQIASNyUKbqL8s473pNEQ6DYIVAAAIk7+mgNlE2bvCeJhkCxQ6AAABI3d64Gyrp13pNEQ6DYIVAAAIn71rc0UG691XuSaAgUOwQKACBxV1+tgfL973tPEg2BYodAAQAk7o47NFD+5m+8J4mGQLFDoAAAEvef/6mBcuqp3pNEQ6DYIVAAAInr6dFA+exnvSeJhkCxQ6AAABL3wQcaKB/9qPck0RAodggUAICLP/7jEMaNC6G/33uS6hEodggUAICLww7TXZQXX/SepHoEih0CBQDg4uSTNVAefdR7kuoRKHYIFACAi4ULNVDWrPGepHoEih0CBQDg4sorNVB+/GPvSapHoNghUAAALm66SQNl8WLvSapHoNghUAAALh56SAPlzDO9J6kegWKHQAEAuHj6aQ2UY47xnqR6BIodAgUA4KK3VwOlrc17kuoRKHYIFACAi1wuhKamEHbbTf+9FhEodggUAICbgw7SXZTXX/eepDoEih0CBQDg5rjjNFCefNJ7kuoQKHYIFACAm7PO0kD52c+8J6kOgWKHQAEAuLn4Yg2UVau8J6kOgWKHQAEAuLnhBg2UK67wnqQ6BIodAgUA4OaeezRQ5s3znqQ6BIodAgUA4OaxxzRQ/vIvvSepDoFih0ABALj53e80UA491HuS6hAodggUAICbHTtC2GWXECZN8p6kOgSKHQIFAOBq8mTdRenr856kcgSKHQIFAODqqKM0UDZv9p6kcgSKHQIFAOBq9mwNlEce8Z6kcgSKHQIFAODqG9/QQLntNu9JKkeg2CFQAACuOjo0UDo6vCepHIFih0ABALi6/XYNlG98w3uSyhEodggUAICrRx7RQPnyl70nqRyBYodAAQC42rxZA+Woo7wnqRyBYodAAQC46uvTQJk82XuSyhEodggUAIC7SZNCGD8+hA8/9J6kMgSKHQIFAODu0EN1F+V3v/OepDIEih0CBQDg7gtf0EB57DHvSSpDoNghUAAA7ubN00C55x7vSSpDoNghUAAA7q64QgPlhhu8J6kMgWKHQAEAuPunf9JAufhi70kqQ6DYIVAAAO4efFAD5ayzvCepDIFih0ABALh78kkNlOOO856kMgSKHQIFAODujTc0UA46yHuSyhAodggUAIC7XC6E3XcPobFR/71WECh2CBQAQCocfLDuorz+uvck5SNQ7BAoAIBUmDlTA2XjRu9Jykeg2CFQAACpcN55Gij33us9SfkIFDsECgAgFZYu1UC5/nrvScpHoNghUAAAqfDP/6yB8p3veE9SPgLFDoECAEiFhx/WQDntNO9Jykeg2CFQAACpsHmzBsrnPuc9SfmyEihtItIpIktEpGPwn+WaNnjsGhHZUMGxBAoAIBXef18D5WMf856kfFkIlBYR6RaRPYc81zH4KGXaiM9rFZFe0VgphUABAKTGPvtopGzb5j1JebIQKMtFZPWI55pFJCfDo6WQQiGyZPDYw0scS6AAAFLjiCM0UHp6vCcpTxYCpUdEFhd4Picic0ocm5PRkTJt8PlCX3MoAgUAkBqzZ2ugPPKI9yTlqfdAaRGNifkFPtYrIitLHN8tIutHPEegAABqzre+pYFy663ek5Sn3gOlTYoHylYRWVvF11w4+DVPKPF5BAoAIDWuvVYDZdky70nKU++Bkt/tmF3gYz0yenekHBtF5PkyPo9AAQCkxk9/qoEyf773JOUhUCqzUETeFJEDyvhcAgUAkBqPP66BcuKJ3pOUp94DJc5LPPm3Kx9Q5ucTKACA1HjtNQ2UQw7xnqQ89R4o+RfJFnsXT6kXyQ7VJSJTK/j8voaGhtDY2Fj00d7e7v37DwDIiIGBEHbbLYSmphByOe9pVHt7e9HvkSISRKS/gu+7NadbRK4q8HyxSz+FdMroF8XOKHEMOygAgFQ58EDdRXnjDe9JSqv3HRQRvRPsyHuZ5C/9lLpRm4jIRTI6ZNpEZEGJ4wgUAECqnHCCBsrGjd6TlJaFQGkWvedJ65DnVsnoXZX8C2qH/qydmaIvpu0QvSPt8sFjNwpvMwYA1JjzztNAuftu70lKy0KgiGicdIruenRI4dektIm+Q2fobkmviAyIhsvQx4CU3n0hUAAAqfK972mgXHON9ySlZSVQPBAoAIBUuf12DZSvf917ktIIFDsECgAgVX75Sw2Uk07ynqQ0AsUOgQIASJVXXqmde6EQKHYIFABAquRyIey+uz4GBrynGRuBYodAAQCkzic/qbsoL7/sPcnYCBQ7BAoAIHVOPlkD5Re/8J5kbASKHQIFAJA6F16ogXLbbd6TjI1AsUOgAABS50c/0kBZtsx7krERKHYIFABA6txzjwbK177mPcnYCBQ7BAoAIHWeeEID5dhjvScZG4Fih0ABAKTOW29poEyZ4j3J2AgUOwQKACCV9torhPHjQ9i+3XuS4ggUOwQKACCVjjpKd1Gee857kuIIFDsECgAglb7yFQ2Uf/s370mKI1DsECgAgFS68koNlGuu8Z6kOALFDoECAEilO+/UQDn/fO9JiiNQ7BAoAIBU2rhRA+X4470nKY5AsUOgAABS6Z13NFAmT/aepDgCxQ6BAgBIrY9+VCPlvfe8JymMQLFDoAAAUuvzn9dAefJJ70kKI1DsECgAgNSaN08DZc0a70kKI1DsECgAgNRavlwDpb3de5LCCBQ7BAoAILXuvVcD5ZxzvCcpjECxQ6AAAFLrmWc0UKZP956kMALFDoECAEitbdv0Bwa2tISQy3lPMxqBYodAAQCk2p/8ie6ivPyy9ySjESh2CBQAQKrNnq2Bsnat9ySjESh2CBQAQKpdcYUGyrXXek8yGoFih0ABAKTa6tUaKOed5z3JaASKHQIFAJBq+XfyfO5z3pOMRqDYIVAAAKm2fXsIDQ0h7LFH+t7JQ6DYIVAAAKn3p3+quyhbtnhPMhyBYodAAQCk3plnaqDcf7/3JMMRKHYIFABA6v3whxoo//AP3pMMR6DYIVAAAKn38MMaKKec4j3JcASKHQIFAJB6vb0aKJMne08yHIFih0ABANSE1laNlNde855kJwLFDoECAKgJc+ZooPz7v3tPshOBYodAAQDUhB/8QAPle9/znmQnAsUOgQIAqAldXRoof/3X3pPsRKDYIVAAADXh7bdDGD8+hEmTQhgY8J5GESh2CBQAQM04/HDdRXn2We9JFIFih0ABANSMv/1bDZRVq7wnUQSKHQIFAFAz/vVfNVDOPtt7EkWg2CFQAAA146WXNFAOOMB7EkWg2CFQAAA15ROf0Ejp7vaehECxRKAAAGrKhRdqoNx4o/ckBIolAgUAUFN+9jMNlJNP9p6EQLFEoAAAakpfXwi77RZCY2MIH3zgOwuBYodAAQDUnBNPTMfP5SFQ7BAoAJqFUHIAAAbZSURBVICas3JlOt5uTKDYIVAAADXnjTdCaGgIYY89QnjvPb85CBQ7BAoAoCZ98Yu6i3LHHX4zECh2CBQAQE266y4NlBNO8JshK4HSJiKdIrJERDoG/2l9LIECAKhJ27eHMHmyRsr69T4zZCFQWkSkW0T2HPJcx+DD8ti+hoYGn9/VDGtvb/ceIXNYcx+se/KytubXXKOBMneuz68/GCj9ZX6vr0nLRWT1iOeaRSQnw8Mj7mP7RMTndzXDGhsbvUfIHNbcB+uevKyt+TvvhNDSEsK4cT67KIOBEsr7Vl+bekRkcYHncyIyx/BYAsVB1k4gacCa+2Ddk5fFNf/Rj3QX5YgjQujvT/bXrvdAaRGNifkFPtYrIiuNjhUhUFxk8QTijTX3wbonL4tr3t8fwmc+o5GyaFGyv3a9B0qbFI+MrSKy1uhYEQLFRRZPIN5Ycx+se/KyuuabNoUwcaJGytKlIQwMJPPr1nugTBONjNkFPtYjIuuNjhUhUFxk9QTiiTX3wbonL8tr3tWlP6NHJISTTgrh2Wftf00CxeZYEQLFRZZPIF5Ycx+se/KyvuaPPRbClCkaKePGhTB9egiXXBLCrbeG8NBDIfzqVyE8/bTuuDz3XAibN4fw/PMh9PSE8MILIfzv/4awZUv5j8bG+g4Uz0s8QURCU1MTjwQfrDlrnpUH686aez0aGpqCSBIPqetAyb/Qtdg7ccp5kWw1x+Y/J5R49IvutPCI78GasuZZebDurHk9PPpl7O+TOalj3SJyVYHni12+ietYAACAojpEZM2I5/KXb0rdbC3KsQAAAEU1i963pHXIc6tk9M5I/kWxS6o4FgAAoGKtoj/wb4Horkih15W0icibMvrSTTnHAgAAAAAAAAAAAEAKtIle3lwienlzydifPqZ1wgvLAQw3TfTcUok4z0t1KcoCsbjVibJu0waPXSMiGyo8NqtaRN9mPzQqOqTyk4mIyFzhnW/linp+GHn8glinq19Rz+n5Y6r9fySLZore9LTUPcWGivO8VJeiLBCLW50o6zay0FtF35k18u3jGG65iKwe8VyzVBca60RkoIrjsibq+WGaaIDnj58p+meddR9b1HP6yM+bIZxfxnK4iHSJrlu3VBYocZ6X6lKUBWJxqxNl3QqdKJYMHnt49NHqVo8Uv4vynAq+zpLBz+fPeGlRzw9bZfi7EOeKvjORdR9b1HP61ALPDw1FFLdBKguUuM5LdSvKArG41YmybjkZHSn5e97wlvHC8j/modDPoeqV8k8obaKXGGYKgVKOKH/OLxKNEVQuyrqvkcI7LV3Cn/dyVBIocZ2X6laUBWJxqxN13bpl9E+eJlDGFvUHZeblT9wESmlR/5z3yOhdAJQWdd3zu7FDX7PSJuX/P5J1lQRKXOeluuX505GzymLdFg5+zRMizFXP8gFX6OdN9cjo4Ctkpuzc+iZQSovy5zz/TXal6N/4h75gE2OL4/zSPfg1NoheVuuMbbr6V0mgxHFeqmtRFojFrY7Fum0UkeejDFXn4ljzod8cCZTSoqx5/ptslwx/PUTH4HMoLo4/682i32hzwuX6ShEoMSJQkhf3ui0UvVZ/QLSx6lrUv1UukOExQqCUFkegjPyc1sHnZ8QxYJ2K4/yyXPT/lQWyM1LYRSkPl3hixCWe5MW5bvm3Ex4Qfay6lr9kUOyFg6VeazXy3hv5QGmOZbr6FPXPebHfl5zww03HEnXd83GS1yo7d1MIw9KqeZFsNeelTIh64mZxKxfnuo3cAkdx3VL4G1uxv23mzRF9Z8PQR/6EvUb4m2UxUf+cjxUovHi2uKjr3lvk+S7hz3o5Kn2bcbXnpcyIskAsbnXiWLdOGf2iWP6GU1yHjH57dv5vm5VeqslvfXOJZ2xR/pxvkML3/GEHpbRq171NdN0L4WZt5ak0UOI8L9WlKAvE4lYn6rpdJKNPNPl7dKCwZtG/HbYOeW6VjD6R56/hj3Vr8Py7prjEM7Yof84XiH6jHSr/e8O71cYWZd2L7aAslMKXjTBcjxQPuULnlnLPS5kV5cTN4lYnyprPFP2foEP0evHywWM3CifuUlpFd54WiK5foW3wNtEXHRf6m2b++F7RW91vEN76OpaoUdgtw3cFVwmXd8oRZd0XyOhLOW3C7slYmkXXrEt0PQek8OXfYueWcs5LmRblxM3iVqfaNc9/c8yNePCzYZBGUc4tzbLzRN9Z5FgUFmXd86+76hR+thoAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALH7f9EQoe/WmJ4iAAAAAElFTkSuQmCC\">" + ], + "text/plain": [ + "<IPython.core.display.HTML object>" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "import numpy as np\n", + "%matplotlib nbagg\n", + "import matplotlib.pyplot as plt\n", + "\n", + "# the plot says a an effective kM = 0.8 is reasonable\n", + "k = np.arange(0, 2.**10)\n", + "kM = k[-1]\n", + "fig = plt.figure(figsize=(4,3))\n", + "a = fig.add_subplot(111)\n", + "a.plot(k/kM, np.exp(-36*(k / kM)**32))\n", + "fig.tight_layout()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "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.4.3" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} -- GitLab