add plot of smooth Fourier filter

meta/smooth dealiasing.ipynb

+{
+ "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
+}