examples/00-mwe.ipynb

{
 "cells": [
  {
   "attachments": {
    "8a415f1d-cfcf-4c3e-9078-13167177082e.png": {
     "image/png": "iVBORw0KGgoAAAANSUhEUgAABGEAAABgCAYAAAC9g2+EAAAAAXNSR0IArs4c6QAAAIRlWElmTU0AKgAAAAgABQESAAMAAAABAAEAAAEaAAUAAAABAAAASgEbAAUAAAABAAAAUgEoAAMAAAABAAIAAIdpAAQAAAABAAAAWgAAAAAAAABIAAAAAQAAAEgAAAABAAOgAQADAAAAAQABAACgAgAEAAAAAQAABGGgAwAEAAAAAQAAAGAAAAAA1gcuCAAAAAlwSFlzAAALEwAACxMBAJqcGAAAAVlpVFh0WE1MOmNvbS5hZG9iZS54bXAAAAAAADx4OnhtcG1ldGEgeG1sbnM6eD0iYWRvYmU6bnM6bWV0YS8iIHg6eG1wdGs9IlhNUCBDb3JlIDYuMC4wIj4KICAgPHJkZjpSREYgeG1sbnM6cmRmPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5LzAyLzIyLXJkZi1zeW50YXgtbnMjIj4KICAgICAgPHJkZjpEZXNjcmlwdGlvbiByZGY6YWJvdXQ9IiIKICAgICAgICAgICAgeG1sbnM6dGlmZj0iaHR0cDovL25zLmFkb2JlLmNvbS90aWZmLzEuMC8iPgogICAgICAgICA8dGlmZjpPcmllbnRhdGlvbj4xPC90aWZmOk9yaWVudGF0aW9uPgogICAgICA8L3JkZjpEZXNjcmlwdGlvbj4KICAgPC9yZGY6UkRGPgo8L3g6eG1wbWV0YT4KGV7hBwAAJaBJREFUeAHt3QuUFNWd+PEfAgaY4SFGGEHkNeJ/iY6DuGoUoZc1ODCLiiEGXDU4HJN/HivsK/sP5n+WPCRZ12zQbMwa/8xy1BUxJEFxYMJGToFC9h8jL7M5QZlhEBREwREYHiK49/bQUNNTdauru7q7qvpb5zRdfd/3U830zK+rbomwIYAAAggggAACCCCAAAIIIIAAAggggAACCCCAAAIIIIAAAggggAACCCAQB4EupyfRTz1Xx2FCeZ7DZtV+a1ofw9Rr/WDLj4Bla5b3qQ0jx139PtbvZ/uGr10jt30n31SLidQOz64CLSpHPzLdeO9mKnW2XIva1Y9MtmpVSBvHcTP9Xy3EfIepTvSDLT8CVlqzcX4vp0214C9bVI/6Yd/wtmsEv9+imtSP9C2RnhDh104/o+M0v2IcmhbVqX6ktn5qR/9fZcu/wGbVhX5Pn9kSau9jHp4G2il9m68SsMufgd07gXVg7zXLDnt6P6GeeS8HY2CdNnV6wtjbeL4TnCEtofJw9Wcw3+CZnmWphLj66rkVc5uvOo+rbRjmlX5sLbzz9n6bn46tXlvqEYb3QVzHMF/5Om1xmq/lMME4za8Yc5mfZqoDMMUYRyn2OSxlf05qh2cEEEAAAQQQQAABBBBAAIFICLREYpQMMuwC+uwMtsIItKS6IQiTkuAZAQQQQAABBBBAAAEEEIiGAH88R+M4RWGUW6IwyIiPca19/ARh7BrsI4AAAggggAACCCCAAALhFyAIE/5jFJURtkRloBEeZwdjgjARPpIMHQEEEEAAAQQQQAABBEpOgDMXSu6Q53XCBPTyyptsvIMxQZj8g9MDAggggAACCCCAAAIIIBCUwOKgGqIdBJRAhwABInkR6GBMECYvxjSKAAIIIIAAAggggAACCORFYHleWqXRUhXoECAoVYQ8z7uDcbdMOrv/zqky4BOZlIx2mX3HRR54akWgkxj/zX8NtL2gGlv33a+5NnXT1/6vHO030DU/yIwezVtk9dOPB9JkzcxbRYb2DqStODTS+P0nA51G3c0zpapsUKBtRrmxuUt+EPjw6+puk+qqfoG3G8YG75tbX7BhzZ5WI1X9uxesvzB0lI/PM7d5hfVzLjVe0+ddqkxYn2v+z11hHVooxnXot2/K+jUd1jrMelw3jf9zmTy4Ouv6pVAx6M+9RxbWFYXN9PkThs/hzVtbpb7+FyYb/aZvMRVwy5s0YZzUVp7nll2U9DmLgvvba9Id98qxEVcUZR5h6jSLz70W0/grKy+RyunXmorkJc/0M/76iROk99UX56XfbBr1+LtLXz7Yam83oyDMZ8raZMzBvfZ6sdzf1KdCHgh4Zuuu/MuAWwyqOfcgzK8+9RmRgZcH1ZGxnfHGXJ+ZKgDTWHXYZyWKZyqgAzB3Lrkw0+KxLzc3DzPUAZjZd2/MQ8vha/K+fAC6TFMHYGYdaXLJjWdyPj7P3KTC+zmXGrH7512qRFif+UwzH5mancH9MTm43wA+48zcEvSP7WJ93pk+f8LwObzgIc/fwZd7HCrX7CH9y0P3eTjHdbT+M3QAJvyfSf7n5b9GVp97Org3wamv7dvfkO1/ooJb3U86ZectzfQzvvfA88Lzd1+bZ0hlczoSlyOli/AaAQQQQAABBBBAAAEEECiCwCu/2+XVa9ZBGK+GyS9pgU6BArvGZe/2sr9k3ybw57t72l457lrpqQRh0kV4jQACCCCAAAIIIIAAAggUWOCdvRfL6tWWqVd9WUOLqQB5CGQpYAzCXPTex1k2G/9qPXYd85pkJ1uCMF5k5COAAAIIIIAAAggggAACeRZYu+GTXj1YXgXIRyBLActYb+chY3YpZ25bt8k0/Q9UJkEYkxB5CCCAAAIIIIAAAggggEAxBDZvPeLVreVVgHwEshRoUfV0wMBxa1yy3DG91BP7f9BD9Jo5hs1yyuNMGCcV0hBAAAEEEEAAAQQQQACBAgo0rHzVqzfLqwD5COQg0OmMDXtbQ94ts79kXwlcvdNzUV7LCYogjJMKaQgggAACCCCAAAIIIIBAgQT0ejAe36jvVENpLdBw6KY0BSzTtEcdIHTQycf7Mi1HUyQ7SZKAAAIIIIBA5gL/1HWk6Ifbtq7PRfLZ3f3k9fLz3YqQjgACCCBQ4gKvbvFcD8Z4lkKJ8zH9YAQsUzMZLEBrqh7LvO2v/LdpXo7rwegKnufPmFolDwEEEEAAgVIW0MGXBU+taCe4c6r8w8mmDhzLyofL7McakmlNu0fIszcOllGH93cowwsEEEAAAQS2N5/yQrC8CpCPQI4CxkBfa8s7qvn+OXYRn+rZrgejBTgTJj7vA2aCAAIIIFBAgQ4BGNWvDsbYz4jR+7Mfbw/A6GHt2NEst//6Lc6IKeAxoisEEEAgKgKvbtrnNVTjH8helclHIAMBfbmbvg2647Z+zVqRNs7hSOGMPtA9tev2bLllEIRxkyEdAQQQQAABF4H0AEyqWCoQM/v9irNnyKQy1TOBGBsGuwgggAACZwSWLVt5Zt9lx3JJJxmBIAWMwb7LPugZZF+RbqvvWye8xm+5FSAI4yZDOgIIIIAAAg4C+hKjM5cgOeTrvGWNaxxy2pMIxLjSkIEAAgiUpMDmzZd7zVudgsCGQEEELFMvF733sSm7pPJOvn/UNF/X9WB0Jc4nMtGRhwACCERUYN3LY2TQwCNSeck2xxlsf+NSefudXjJ+3CbHfBKdBTb2qTizxotzicxSdSDmG3+skJ9flFl5SiGAAAIIxFfg1a3ZX9bgR+XgsROyaUCFnyoFKPuHAvRBFz4EjGfCiL4bUFUXH83Ft2jjkuWmyVmmTIIwJh3yEEAAgQgKNDReKbffUS/Dh4+QJxddL2OqX+swCx2AuXXGK8lLY559uk5qazZ2yOeFu8CVB/fKPLUAr+lMGPfaZ3Om10yU7wxuEzlyNo09BBBAAIHSFGhq9vwwsIKQ+dmqF+VnQTREG3EW0EEYfRZHX6dJJgMPU6c5ZZVUml6U94B5xsZgFpcjmfHIRQABBCIlkArA6EHrsy3GTXxSdFpqswdgdJoO1tjzU+V4dhfQd0DSgZhst/s+XysLKw7IoCPq2yQ2BBBAAIGSF2hY+arJQP9BbJkKkIdAwALGAIIOQJT6VvZhVy8Cy1SAIIxJhzwEEEAgQgL2AIx92DrQ8sA/Xy7pAZhUGQIxKYnMn7MNxCyYNVW+c84O6X38w8w7oyQCCCCAQGwF3tl7sWzf/oZpfsZrHkwVyUMgSwHLVC+DuwKZqsci71N7PC/JMhoShInF24BJIIBAqQssemJM8qwWN4cF33tSrrhmafLsGKcyBGKcVMxpX+m2K3nJl7nU2dy6aTXy1eNNZxPYQwABBBAoeYG1Gz7pZUAQxkuI/KAFjGfC9Nl/Kuj+ItfeOR8aFyj2XEibIEzkDjkDRgABBDoKPPxotdw39987JmbxikCMP7QHjw52DWo5tVT/y0bRt7ZmQwABBBBAICWweatxPRh9KRJBmBQWz4USMAZhTu05XKhxhLYfDwPLa+AEYbyEyEcAAQRCLjD+upO+zshwm45eyPfSkWqxWDZPgR9/YqQ8srTBs1x6Ab2gL4GYdBVeI4AAAqUp0Ha4nzz8yDOmyROAMemQly+BFtWwDgA6bh53BXKsE7fE7a/8t2lKxiCWrkgQxsRHHgIIIBABAX33oxdfuFw++9kpWY9WB2CWP/Onrre0zrrhGFZc1XuozFu8IuuZEYjJmo6KCCCAQKwE1q4f4TUfy6sA+QjkScAYSBjyblmeug1/s3phYo91nIx2eobcojr8xzn7Eb7fLJNe/DfH+qsdU9sTJzb+WLqV9elUYnv1JGmuvLFTOgmnBdrUf6eyj8wcmZQxt0Cug0DTpBPyWI+tDjntSV86ViUjV3d3zY9DxsALd8qjP+wnl466S/T6L342AjCZa71efr7M+OmqzCu4lEze4lrdYUkv8MuGAAIIIBC8wF/9nVuAw/3voxdWNsvmrZ3rTf5MP5ly08bAB/nyb4yXIun+rCw6/bMs6oS1SmtYB+ZnXDf+5qdyzlvbHaucGlwpv/70Fx3zipxoqf4nuI1h1IFzZNcFbrn+02tWdF5jxXS2ic6rkU916qhxqueCuZ3q+E24emc3aXSvpM8ganHPbs8hCOMlFOX8D9tk9dOP+57BmueWOtYZP+IKaXbMIVFUcGX4gldl9ORrpGHcMUeQ2t+Vyfsb35QN9wwQ6X7SsQyJ2Qm09Twp9c8tca185y36h3S8gzB68uXlrXL/37dKedksmffNxTopo40zYDJiSha68MQhmV4zUZY1rsm8kkvJql4qaMtdql10SEYAAQRyE6iv/4XvBlavthzrVFfVOabnmuhxa+qdqv2WLPqwsqhDlTwK6ACM299kk+64N48959S0e7RSNdvj3RMil+bUfofKfi9x0meiOJ6NMnVah3bz8mKn8Zc3o1tqPARhUhI8I5CtwImuMvH5Q7JmR7PseLRZamVmp0DMDTt6S8O/PJHsYWLvGlkzvZxATLbe1DMK6NtQP75onbFMeubS5eeq4E16Kq+dBPStpRdWHBAxBGKe+eJk2XqkmyTPdnFqRKUturdWJh/a4ZJLMgIIIOAtsPWWY7K11/uOBQec6iWTlvZ1zCMxHAIvra9Wf0QuNg0moz/mHBqwHNKimqQN5kZ18BEfd4tp/K1/fFtkXH9TkdjmeQSMMvp/SxAmtm8PJlYogYnLDsuaFWdPSmt4dEmHQMz5H/SU3f/v5TPD0WWnXPA5WXnjmSR28izwk2O/ld43B3zt6vN5HnQWzesAzK0zXvF1xx7djb58qfryOqmtCf5U6yymEfoqpkDMgllTVXClSSbrWajLjZwCMfNU+vTDXIYU+gPNABEIuYAOwMxd8gPHUdbdPFMmyVjHPBLDIfDSbzzPirayHKnrJSRZtke1PAo0/9ZS/1e/HmgPpmUnfHSkgwn6bKyhTnXWr1krMvNz3ksxOFWOcJpeC2eXefwt5uz2XIIwmShRBgEXgdqXe0jDil92yj0TiLnmhFy6dLdsUGfJ2LeV9T+T2h6dz5ixl2E/OIGf/SqEEZPgppdsKdsATGoY+vbUL6+5S/Qiv2zeAqlATOsN18mvX9qQrFA3rUa+evxscCW53ktaIEYHYFgHxtuXEggggEDcBZY++19eU7S8CpAffQHXy2rCMTVLDeMLbkO5TH3R/Psy46U5blUjm/6pPV28gjDxORNmU58KSTzmfv39BzNHR/ZAMvDoCug1Xhoefcp1AjoQM3zViE4BmFQFnX/D4LvlpeGl9cMrNX+egxPINQCTGslds9dzh6QURgbPOhDzvf/VW5p2j5CRF1XIt/vuFTnesaI9EEMApqMNrxBAAIFSFWjaPkpdivSsafp6cc+M/pgzNUIeAjkKWKq+axDm4rdPye8H5dhD1Kqb14PRs7EymdI5mRSiDAIIdBRoX+PFPQCTKr0j7QyYVHrq+aX7nxDdFhsCuQgseqqL70uQnPrT79e1v+nllEWai8Cow/vl2RsHq2BMV9FBGadNB2L0OjGcAeOkQ1qoBdSaZ8bNK99YmUwESlfg1a1qbUDzRgDG7ENuYQQsUzfHtu83Zccyz2M9mC2ZTprLkTKVimK5c8skyBW3t/cZGEWFvI1Z39bXK8iSt85pGAGbwPe+9Ud1ZyTzran1+/Xe2eONd02a9427ZPbdm2wtF213mOo54dG75ZFfsGwdiPHaJh/Sl1WzIRAdAf0FgV7P7ODXr5X9fY92HrhelF6tidZzyHmdFqPvXJgUBAovUFd3W2CdXhjwr8Dbm9WdZcybZc4mF4GCCLSoXlzXhVnz/CqRGcHcjahm5q2BTOjsKp2BNNehkQzWg8k4eBqLIEzfJX/oAJT9i6DayX4EgdY8b4Ssnv5goE3SWLuAvoTofPWL6XDFm0sg5oYHuByJ91QwAvf/vV7LxTkQowMw7beh3iyH25zL6ABMexvZjUefWn3o8Cd8VDZ+TulTX11Pfz3dSRcfnVEUAQR8COgAjD5TU2/6c04cAjFT1naRlacXpXe6K2CyMv8gUESBHz3UcT2+Ig6lU9fbXvcM3ludKpGAQHEELNWt6+9kycDEBW05j6xxavh/rQtqPRiNFYsgTM5HnQYQyEJAfzO4f95YdXvqUR3ujpRpUwRgMpWiXKYCToGYswGYbclm5nx5l2x7fYr8/OcrzzSbawBGN7Tw3z6S+vpnz7TJDgIIRFPAHoDRM9BfNKQHYmrUpRQr6588M8Ezi9GPO3YmjR0EEHAXWLbs7GewSynLJZ1kBAotYKkOXYMwycDEBYUeUnH667ovbeG/zsMwfsNoL04Qxq7BPgJ+Bco+kjXTy6Wm7FZpfGZ5xrUJwGRMFUjBz910s/T+RLC3qK5/fkkgYwu6EXsgJj0Ao/sqL2+Vf/52X9m4sf1yuiACMEHPgfYQQKA4AukBmNQo7IGYQW3dpPH77WfJpPL1M4EYuwb7CLgL6DNHPdbcXeteO7ec2eoufv/S483cGgm4dnBXNAQ8sAI1V1l5iYy4OuG7t9VPP+67TpYVjIGFDAITWXYbvmrb1nlesm9lOmqCMJlKUQ4Bg8CHbf6+/TvXef1OQw9k5SLw5R5XS9VzPXJpolPdeglnEEYPVAdiystmSe1njkrlJe1nwNgnMPDCnfLkoutl3YbxMucrxs9WezX2EUAgxgJuAZjUlFOBmMOpBIdnAjEOKCQhkCawrclzUV4rrQovYyygAzBZLR/hHYT5M8VmBUCnf1HUd+vq69RWMjBx42VOWfFKU19A6NuJG7aMF+XVbRCEMUiShYCnwOmFCdecvi7es/zpAi9+6wm5/KG75bVB3J46UzPK+RPwCq6MqX5NxlT7a5PSxRHY06tcfiX+VoXc53159iw1m4TLjBardP1gKxUB9VmmF+H12jJZA00HYs6//A7nxXy9OiAfgRIQ2N58ymuWllcB8hEosIAOxExw6lMHJvp/cJUc6OvvC2mntsKcNm5fT/H4lPT1rSZBmDAfbcYWboEsAzCpSR3+8cvJxX0d7zqRKsRzRgL6LJf35CH3ss+5Z5GDQNgF9nYrlzmPrQh6mENVg/rhtFlOiaTFWKD7yeRdkHJdbF4L1f7NndLQ1zsKGGPNUExt3R9+J1+/OeChPB9wezFv7pGfVMsbTQc7zXLduvWd0tISfP0xl1aXlyEUSJ7p4nKzlNUhHK/DkCyV5hiE0WVHH+guL8c8CNNnv2fw1Nf/W4Iw+p3DhoBfgRwDMLo7/Y3idUsvkg33DBBRvwCzIYAAAgggUCyB5BcCOd71r/YrM6XhKgIwxTqG9n71t9Mep87bi7OfBwEdgKmv/4XflvUlDa1+K1EegTwLWKr9f3Tro+9b6pbrw91y45F+ao/pYtzkHEsvCGN9aWIgR7fpVE+Z/XhDIG3RSLwFav9/d2nweQmSk8gGa53UXHuXNFZ5/sd2qk4aAggggAACgQnkEoiZOLVGGq5Rv4izIYBALgJWLpWpi0CeBIwBhtY/vi0yrn+eug5Hs41LlnsNxPIqYM+PxZkwYw7utc8p+/0+FdnXpWZJCehfNCfuqvG8NXXNDPNdk5LfGhKAKan3Tlwn+6OHmuVHD1UHNr2/+rsR2XyDGFj/NIRAqQpkE4i5LjE+eadAzuos1XcN8w5QwAqwLZpCICiBVtXQTvVwvIx5/Rp1Q68vTI/tmf1D3i2TXWZJBeBvO8dfcUojgEBSQF0+pG9Nrb/5c9v0bagbb+4iOhDjtCUDMOPivYiV07xJQwABBBAIt4C+DbWfbc/O3XL+kXP9VKEsAggggEC0BCzTcC97t5cpO9J5ow54hkw2+52gv09Zv60HVL7io8Py8Oyp7q0daXLPIweBfAmkAjHS+YwYHYB5aXj7nY8aJ3eTiW0dyxCAyddBoV0EEEAAgVwEaraWS+P3n/DVROr21aLWlGGxeV90FEYgXSChEpanJ/IagRAI6EDDF9zGcdF7H8vvB7nlRju9x7uel9pafmcYiSDMhUcOyyxhzQy/B5fyBRBwCMTYAzDJEagyW6b0l+G/H5FcjJcATAGOC10ggAACCPgTUAvOT1nbRVbWP+mv3unSBGKyYqMSAukCifQEXiMQEgHz2R471ZfPVV1CMtRgh3F0d6tXg2Ybh9qRCMI4jJskBMIjYAvEnLh2wJkzYOwD1N8M9pp7ndRuu0YauATJTsM+Agh4COh1zz6YOdqjVObZm9T6Z4nH1mRegZLxFzh9x7+VOS44nwrE7J83VqTso/i7hXiGlZWXyPjRVwU6wvrnlwTaXtwbu2RkH6mru63TNNete81056orVIV+6uH5V1+nhklAIL8Clqn55MK1U6eZikQ2b83zq0xj/0BltpgKOOURhHFSIQ0BvwI6EDOzp6rVfgmSU/VdF7TJrgucckjLVWDrLcdk4nPfdG1mzS3flarnerjmk4EAAgiUtID6DOs55LxACEZPvkZ2lLHeWSCYOTSiAzAPPq+CYQFu9UIQxg/nfV92/nL8kZ9cL9+4/w1TU9Uq0zIVIA+BIgnoBWgnuPZ9xy9ds0wZNTPVjUymFvksmizHruZlmebmlue5yoxbRdIRQAABBBBAAAEE4iGgz9LUl8uaNn0XpOHDR7gW4XJbVxoyEDgjUDnC88+vxJnC7CAQLgHnyGK4xljo0WRl4vlToNCzoD8EEEAAAQQQQKBAAvobvY8L+EgUaF5ZdWMKxOi7AW64Z4AcVIvvOgViCMBkRU6lEhS4dKTnOpeJEmRhytEQsKIxzIKOMisTLkcq6DGiMwQQKIaA6VKlYoyHPhFAAIGwCiQDMTJTGh49e+lJMsByjbo7hLpsKXn3Ix2IeVCSi83reRCACevRZFxhFBhZ+brXsNwv9/CqSX7xBW4L5tLO4k/EcQSWY2ppJ2ZlQhCmtN80zB4BBBBAAAEEEOggYA/E1P7NndJwVVuHfHsgRq8Bo8uzIRA2gUVPXBnYkMZWnZDq6tcCa2/69CmybNlKU3sJlWmZCpCHQBEEWlWfW9RDLyDNJqLXyMlqIwiTFRuVEEAAAQQQQACB+ArowMr5l98hDX07BmBSM9aBGH0XJBbhTYnwHDaB++bWBzakRxbWqSBMYM3JpaPO92osoQpYXoXIR6AIApbqkyBMO3xW68HoqqwJ0w7IvwgggAACCGQl8E9dR4p+uG2reg+Vz+7uJ6+Xe/7S7dYE6QgURSB5xoupZ25DbdIhDwFXgcoR3V3zTmcEGPLx6op8BHwJWL5Kx7vw8mynV7AzYTb1qZAtH5W5jnPWkSbXPDIQQAABBPIj8M7ei+WF1Z/01XjQp2X76jxkhXXwZcFTK5KjGlD3F3LP0eYOI9QBmBk/XZVMa9o9Qp69cbCMOry/QxleIIAAAgiUlsDYKhbnLa0jHqvZWmo233KY0SyVNtQhPZk05OE73bLkt+eeVHmhvqx1pxrgYocJWA5pGSUVLAijAzBzFrX/ouo0slkzRzslk4YAAgggkEeBPXv7it9TtoM+LTuP08tr0/YAjO5obv0LUvHFyTL5kP6sFrEHYPTrHTua5fZfC4EYjcGGAAIIlLCAXpy3svIS2b79DTeFvipjmHq0qAcbAmESaFWDme8woIRKcw3C7LrA+dJWh3bCmNSiBjU/yIEVLAgT5KBpCwEEEEAAgWIKpAdgUmP5xn9uk5HqbJemLuVnzoBJ5elnAjF2jcLsT7rj3sJ0pHppe2eXrH+xsWD90RECCERX4PO3XysPLHANwuiJJdRjsXqwIYBAzAQyCsI0neopoi4nymXbF+ngVy4zp27sBU50lSGtPUI3zV2hG1HxBjTv9v8tA7q6Xw6ZzcjmLvlBNtWoEwMBtwCMnloqyKKf3bZUGS5NchMKNn31dHUv5QJt4zf+hwhBmAJp003VkfNk4cy/dYQYcKqXYzqJ4RG44dNdvQbDujBeQiHMH//Nf3Uc1brvfs0x/XTiB+p5s6lAFnn6jBW2kApkFISZ/XhDSIfPsBAovoAOwOya81TxB8IIXAVuPH6RVD0XbKBsrmtvZMRZYKP6QmLBY+6X1uq5mwIwKRtd5ifvjZIf9mB9mJQJzwgg4E9Af65VyYX+KlE6NAI3XL/Z65IkgjChOVqZD2TdlX/pUtgYhNEBmIRLRZJjKJBRECaG82ZKCCAQI4Gyo12l7uaZrjPS+WwIBCFw5cG9smDWVJm32ByI8epres1E+XbfvSLHvUqSjwACCCCQjUDbAecYRll//feu86bXPJt990aHTKc0h2I+k2qnjJWHH3G9JGmCz+YojgACEREgCBORA8UwEUDAXWDk6u7yoIx1L0AOAgEKfPV4k+z9fK08sjS7s0R1AGZhxQHpffzDAEdFUwgggAACURMY9+leKghjHHVC5VrGEmQigEDkBM6J3IgZMAIIIIAAAkUW+HrPt0QHU/xuBGD8ilEeAQQQiK/AhOvd1w87PWvn03niS8LMECgJAYIwJXGYmSQCCCDgLFBd/ZroU7b1o67uNudCKlWfop0q53yqtmvVWGbos1i+M7hNhg8fkfH8CMBkTEVBBBBAoCQEyspbZc59M0xzTZgyyUMAgWgKEISJ5nFj1AgggAACRRbY0rV/Rovwpob56rYWeflcFtFMefCMAAIIICBSXWW8k9UtymgYTgggEC8BgjDxOp7MBgEEEECgAAKreg+VGT9d5asnfUckXeffe2Z+9oyvDiiMAAIIIBA5gQnXvec15lleBchHAIFoCWS0MK/1pYkyRt0RIpdtca+RMmdRbneTyKV/6iKQL4FdF7SJPD3Ns/mareXS+P0nTeX+WmW6L9lvquk/r9V/FWoggIAWyCYAY5ebW/+C7Ltzqnyl2y4W57XDsI8AAgiUoMDAije9blU9S7HML0EapoxAbAUyCsLEdvZMDIFwCegAjBWuITEaBBCwC2zsUyEzHvN3Boy9fmp/wVMrZO+0GvlhjzdTSTwjgAACCJSogMetqocqFr1Ab6G+qCvRo8C0ESicAJcjFc6anhBAAAEEIi5wpTordJ46iyWIbdqAU0E0QxsIIIAAAhEXGDnCuC6Mnt2tEZ8iw0cAAZtAwc6EmXWkSWbNHG3rml0EEEAAAQSiJ/APJ5tEVCBGn83itOm7IOlFePUaMG7bM1+cLOMP7nTLJh0BBBBAoIQExlad8JptwqsA+QggEB0BzoSJzrFipAgggAACIRHQgRgdbEnfUrehfmlCH8d8XV4HYCYfIgCTbsdrBBBAoFQFqqtf85r6BK8C5COAQHQECMJE51gxUgQQQACBEAksrDjQIdCSCsD0Pv5hcsFdnZ9+6RIBmBAdQIaCAAIIhEhg+vQpXqNJeBUgHwEEoiFAECYax4lRIoAAAgiETEAHW1KBGHsAJjVMna/PmFlY9xfJJAIwKRmeEUAAAQTSBcaOGZCelP5aL87LhgACMRAo2JowMbBiCggggAACCHQQSAVidKLed9ruOdosN82+RgZxCZITD2kIIIAAAkqgcoTnd+MJVWwhWAggEH0BgjDRP4bMAAEEEECgiAJuwRf7kAYdOWR/yT4CCCCAAAIdBMZe8V6H1w4vOBPGAYUkBKIo4BlyjeKkGDMCCCCAAAIIIIAAAgggEBWBgRVvSmXlJabhDlWZ/UwFyEMAgWgIEISJxnFilAgggAACCCCAAAIIIBBjgdopY71ml/AqQD4CCIRfgCBM+I8RI0QAAQQQQAABBBBAAIGYC1RX9fKaIZckeQmRj0AEBAjCROAgMUQEEEAAAQQQQAABBBCIt8CE6zzXhUnEW4DZIVAaAgRhSuM4M0sEEEAAAQQQQAABBBAIsYBeF2bSpIRphBNUJuvCmITIQyACAgRhInCQGCICCCCAAAIIIIAAAgjEX+BPrxriNcmEVwHyEUAg3AIZ3aL6P9vKZEuvkeGeSQCj29cWQCNpTYzf+B9pKby0C/Ro3mJ/mdv+zkNSI71zayOftdX4orxtbXtbnpoZ5RkEPPYlAbenmtu8tVUWPXFl8A1n2OLu3c2uJQs9NtNYXAfpkrH1wAlZ3D/+n2Eu0++UHPRnHZ9zZ4kD/UxTzdZsLT/bOHudBA6982antGwT3mrdpz7j9mRbvTTqBfy5V7zPu82ux6vQn3VOA9mzt9Up2Z52q3qx3J6Q6f6uA4dlcWXYPg//kOnwPcvpn8HjPUsVp8C64nQbml7D/nnWWCSphOr3Yx6eBtopfZuvErDLn4HdOxFzaz2/Ym4J1Tnv5WAMLMOBxDgYYztxQr3ANTiD+XZctW+pB77BGCSUpX2br15gmz8Du7Xet9QD7/wYzFe26ZulEvAO3kBHafqlY9tex8ncss0rtRun+aXm4jTP1Hyj8qznkJpP3J713ALduBwpUE4aQwABBBBAAAEEEEAAAQTyJtBXtazPhmFDAIGIChCEieiBY9gIIIAAAggggAACCCBQkgJzS3LWTBqBmAgQhInJgWQaCCCAAAIIIIAAAgggUBICV5TELJkkAjEVIAgT0wPLtBBAAAEEEEAAAQQQQCC2AtWxnRkTQyDmAgRhYn6AmR4CCCCAAAIIIIAAAgjEToAgTOwOKRMqFQGCMKVypJknAggggAACCCCAAAIIxEVgWFwmwjwQKDWBLqcn3E89E031PvqbVZHWtGLD1Gv9YMuPgGVrNu7vU6f3l236ed+Nu2/eAW0d6J8T+ng6bQmnRNJ8C1i2Grx3bRgB7LaoNvQjtVWrHW3MlruA/rnQamtmmNrXD7b8CFhpzfJeTgMJ8GWLaks/7Bvedo3g91tUk/qRviXSEyL8Wv+83Jw2/kTa6zi8dJpn1OYV5//vcTg+UXs/MV4EEEAAAQQQQAABBBBAAAEEEEAAAQQQQAABBBBAAAEEEEAAAQQQQACBSAj8D4ASOfJyk55+AAAAAElFTkSuQmCC"
    }
   },
   "cell_type": "markdown",
   "id": "3848ed8f-b092-49e1-a66e-8de9a39e6ae9",
   "metadata": {},
   "source": [
    "![superblockify_logo.png](attachment:8a415f1d-cfcf-4c3e-9078-13167177082e.png)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "abbbed56-e2a5-4930-9df7-d8db860fce20",
   "metadata": {},
   "source": [
    "More info: [superblockify.city](https://superblockify.city)"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "bf599af0-108f-4d78-bce4-0491ba7bc863",
   "metadata": {},
   "source": [
    "# Minimum working example"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6ddf99e9-c41c-4a47-878e-4f04d480fa2e",
   "metadata": {},
   "source": [
    "This notebook runs a minimum working example with `superblockify`.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5fa10ad7-f17b-4cd0-b02d-375beacf5efa",
   "metadata": {},
   "source": [
    "## Import `superblockify`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "dfa871e2-7683-428c-a4a6-df97cb62fdec",
   "metadata": {},
   "outputs": [],
   "source": [
    "import superblockify as sb"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "92a6ae87-0b1c-4960-a1b7-9c4f6b4ff673",
   "metadata": {},
   "source": [
    "## Initialize the partitioner\n",
    "\n",
    "A `.data/` folder will be created in the current working directory.\n",
    "The preprocessed graph will be saved in the `.data/graphs/{name}/` folder and the results in the `.data/results/{city_name}/` folder."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "daa39c85-bc2e-4f83-810e-d9957b4fa87c",
   "metadata": {},
   "outputs": [],
   "source": [
    "partitioner = sb.ResidentialPartitioner(\n",
    "    name=\"Scheveningen_test\",\n",
    "    city_name=\"Scheveningen\",\n",
    "    search_str=\"Scheveningen, NL\",\n",
    "    unit=\"time\",  # \"time\", \"distance\", any other edge attribute, or None to count edges\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ba2b9fa5-172e-4dab-949e-3151e6c2c360",
   "metadata": {},
   "source": [
    "If you want to select a different city, find the corresponding search string (`search_str`) at [https://nominatim.openstreetmap.org/](https://nominatim.openstreetmap.org/ui/search.html?q=Scheveningen%2C+NL).  \n",
    "The smaller the place, the quicker the partitioning. For large places sufficient memory is required."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d9580056-4c90-4b3e-b362-0cb361578eb5",
   "metadata": {},
   "source": [
    "## Run the partitioner"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "56c03f44-8603-47b4-9436-e36c7bdf649c",
   "metadata": {},
   "outputs": [],
   "source": [
    "partitioner.run(\n",
    "    calculate_metrics=True,  # set to False if you are not interested in metrics\n",
    "    make_plots=True,  # set to False if you are not interested in plots\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6aaa66e4-12aa-45cb-8307-1c756c02987c",
   "metadata": {},
   "source": [
    "## Save results"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3932594e-5aa8-4c59-9691-a8b07841cd18",
   "metadata": {},
   "source": [
    "Save the results to a GeoPackage (`.gpkg`) file:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c04e0a0c-fb4d-4068-bfec-7d881546962c",
   "metadata": {},
   "outputs": [],
   "source": [
    "sb.save_to_gpkg(partitioner, save_path=None)"
   ]
  },
  {
   "metadata": {},
   "cell_type": "markdown",
   "source": "Find the GeoPackage at `.data/results/Scheveningen/Scheveningen_test.gpkg`.",
   "id": "c602d5133ce568b2"
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "sb_env",
   "language": "python",
   "name": "sb_env"
  },
  "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.12.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}