Streams-TimeseriesForecast.ipynb

{
 "cells": [
  {
   "attachments": {
    "how-it-works.jpg": {
     "image/jpeg": "/9j/4AAQSkZJRgABAQAAkACQAAD/4QCARXhpZgAATU0AKgAAAAgABQESAAMAAAABAAEAAAEaAAUAAAABAAAASgEbAAUAAAABAAAAUgEoAAMAAAABAAIAAIdpAAQAAAABAAAAWgAAAAAAAACQAAAAAQAAAJAAAAABAAKgAgAEAAAAAQAAAlKgAwAEAAAAAQAAAXIAAAAA/+EJIWh0dHA6Ly9ucy5hZG9iZS5jb20veGFwLzEuMC8APD94cGFja2V0IGJlZ2luPSLvu78iIGlkPSJXNU0wTXBDZWhpSHpyZVN6TlRjemtjOWQiPz4gPHg6eG1wbWV0YSB4bWxuczp4PSJhZG9iZTpuczptZXRhLyIgeDp4bXB0az0iWE1QIENvcmUgNS40LjAiPiA8cmRmOlJERiB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiPiA8cmRmOkRlc2NyaXB0aW9uIHJkZjphYm91dD0iIi8+IDwvcmRmOlJERj4gPC94OnhtcG1ldGE+ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgPD94cGFja2V0IGVuZD0idyI/PgD/7QA4UGhvdG9zaG9wIDMuMAA4QklNBAQAAAAAAAA4QklNBCUAAAAAABDUHYzZjwCyBOmACZjs+EJ+/+IPYElDQ19QUk9GSUxFAAEBAAAPUGFwcGwCEAAAbW50clJHQiBYWVogB+MAAQAEAA0AOAAnYWNzcEFQUEwAAAAAQVBQTAAAAAAAAAAAAAAAAAAAAAAAAPbWAAEAAAAA0y1hcHBsAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAARZGVzYwAAAVAAAABiZHNjbQAAAbQAAAQ2Y3BydAAABewAAAAjd3RwdAAABhAAAAAUclhZWgAABiQAAAAUZ1hZWgAABjgAAAAUYlhZWgAABkwAAAAUclRSQwAABmAAAAgMYWFyZwAADmwAAAAgdmNndAAADowAAAAwbmRpbgAADrwAAAA+Y2hhZAAADvwAAAAsbW1vZAAADygAAAAoYlRSQwAABmAAAAgMZ1RSQwAABmAAAAgMYWFiZwAADmwAAAAgYWFnZwAADmwAAAAgZGVzYwAAAAAAAAAIRGlzcGxheQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAG1sdWMAAAAAAAAAIwAAAAxockhSAAAAFAAAAbRrb0tSAAAADAAAAchuYk5PAAAAEgAAAdRpZAAAAAAAEgAAAeZodUhVAAAAFAAAAfhjc0NaAAAAFgAAAgxkYURLAAAAHAAAAiJubE5MAAAAFgAAAj5maUZJAAAAEAAAAlRpdElUAAAAFAAAAmRyb1JPAAAAEgAAAnhlc0VTAAAAEgAAAnhhcgAAAAAAFAAAAop1a1VBAAAAHAAAAp5oZUlMAAAAFgAAArp6aFRXAAAADAAAAtB2aVZOAAAADgAAAtxza1NLAAAAFgAAAup6aENOAAAADAAAAtBydVJVAAAAJAAAAwBmckZSAAAAFgAAAyRtcwAAAAAAEgAAAzpoaUlOAAAAEgAAA0x0aFRIAAAADAAAA15jYUVTAAAAGAAAA2plc1hMAAAAEgAAAnhkZURFAAAAEAAAA4JlblVTAAAAEgAAA5JwdEJSAAAAGAAAA6RwbFBMAAAAEgAAA7xlbEdSAAAAIgAAA85zdlNFAAAAEAAAA/B0clRSAAAAFAAABABwdFBUAAAAFgAABBRqYUpQAAAADAAABCoATABDAEQAIAB1ACAAYgBvAGoAac7st+wAIABMAEMARABGAGEAcgBnAGUALQBMAEMARABMAEMARAAgAFcAYQByAG4AYQBTAHoA7QBuAGUAcwAgAEwAQwBEAEIAYQByAGUAdgBuAP0AIABMAEMARABMAEMARAAtAGYAYQByAHYAZQBzAGsA5gByAG0ASwBsAGUAdQByAGUAbgAtAEwAQwBEAFYA5AByAGkALQBMAEMARABMAEMARAAgAGMAbwBsAG8AcgBpAEwAQwBEACAAYwBvAGwAbwByIA8ATABDAEQAIAZFBkQGSAZGBikEGgQ+BDsETAQ+BEAEPgQyBDgEOQAgAEwAQwBEIA8ATABDAEQAIAXmBdEF4gXVBeAF2V9pgnIAIABMAEMARABMAEMARAAgAE0A4AB1AEYAYQByAGUAYgBuAP0AIABMAEMARAQmBDIENQRCBD0EPgQ5ACAEFgQaAC0ENAQ4BEEEPwQ7BDUEOQBMAEMARAAgAGMAbwB1AGwAZQB1AHIAVwBhAHIAbgBhACAATABDAEQJMAkCCRcJQAkoACAATABDAEQATABDAEQAIA4qDjUATABDAEQAIABlAG4AIABjAG8AbABvAHIARgBhAHIAYgAtAEwAQwBEAEMAbwBsAG8AcgAgAEwAQwBEAEwAQwBEACAAQwBvAGwAbwByAGkAZABvAEsAbwBsAG8AcgAgAEwAQwBEA4gDswPHA8EDyQO8A7cAIAO/A7gDzAO9A7cAIABMAEMARABGAOQAcgBnAC0ATABDAEQAUgBlAG4AawBsAGkAIABMAEMARABMAEMARAAgAGEAIABDAG8AcgBlAHMwqzDpMPwATABDAEQAAHRleHQAAAAAQ29weXJpZ2h0IEFwcGxlIEluYy4sIDIwMTkAAFhZWiAAAAAAAADzFgABAAAAARbKWFlaIAAAAAAAAHHAAAA5igAAAWdYWVogAAAAAAAAYSMAALnmAAAT9lhZWiAAAAAAAAAj8gAADJAAAL3QY3VydgAAAAAAAAQAAAAABQAKAA8AFAAZAB4AIwAoAC0AMgA2ADsAQABFAEoATwBUAFkAXgBjAGgAbQByAHcAfACBAIYAiwCQAJUAmgCfAKMAqACtALIAtwC8AMEAxgDLANAA1QDbAOAA5QDrAPAA9gD7AQEBBwENARMBGQEfASUBKwEyATgBPgFFAUwBUgFZAWABZwFuAXUBfAGDAYsBkgGaAaEBqQGxAbkBwQHJAdEB2QHhAekB8gH6AgMCDAIUAh0CJgIvAjgCQQJLAlQCXQJnAnECegKEAo4CmAKiAqwCtgLBAssC1QLgAusC9QMAAwsDFgMhAy0DOANDA08DWgNmA3IDfgOKA5YDogOuA7oDxwPTA+AD7AP5BAYEEwQgBC0EOwRIBFUEYwRxBH4EjASaBKgEtgTEBNME4QTwBP4FDQUcBSsFOgVJBVgFZwV3BYYFlgWmBbUFxQXVBeUF9gYGBhYGJwY3BkgGWQZqBnsGjAadBq8GwAbRBuMG9QcHBxkHKwc9B08HYQd0B4YHmQesB78H0gflB/gICwgfCDIIRghaCG4IggiWCKoIvgjSCOcI+wkQCSUJOglPCWQJeQmPCaQJugnPCeUJ+woRCicKPQpUCmoKgQqYCq4KxQrcCvMLCwsiCzkLUQtpC4ALmAuwC8gL4Qv5DBIMKgxDDFwMdQyODKcMwAzZDPMNDQ0mDUANWg10DY4NqQ3DDd4N+A4TDi4OSQ5kDn8Omw62DtIO7g8JDyUPQQ9eD3oPlg+zD88P7BAJECYQQxBhEH4QmxC5ENcQ9RETETERTxFtEYwRqhHJEegSBxImEkUSZBKEEqMSwxLjEwMTIxNDE2MTgxOkE8UT5RQGFCcUSRRqFIsUrRTOFPAVEhU0FVYVeBWbFb0V4BYDFiYWSRZsFo8WshbWFvoXHRdBF2UXiReuF9IX9xgbGEAYZRiKGK8Y1Rj6GSAZRRlrGZEZtxndGgQaKhpRGncanhrFGuwbFBs7G2MbihuyG9ocAhwqHFIcexyjHMwc9R0eHUcdcB2ZHcMd7B4WHkAeah6UHr4e6R8THz4faR+UH78f6iAVIEEgbCCYIMQg8CEcIUghdSGhIc4h+yInIlUigiKvIt0jCiM4I2YjlCPCI/AkHyRNJHwkqyTaJQklOCVoJZclxyX3JicmVyaHJrcm6CcYJ0kneierJ9woDSg/KHEooijUKQYpOClrKZ0p0CoCKjUqaCqbKs8rAis2K2krnSvRLAUsOSxuLKIs1y0MLUEtdi2rLeEuFi5MLoIuty7uLyQvWi+RL8cv/jA1MGwwpDDbMRIxSjGCMbox8jIqMmMymzLUMw0zRjN/M7gz8TQrNGU0njTYNRM1TTWHNcI1/TY3NnI2rjbpNyQ3YDecN9c4FDhQOIw4yDkFOUI5fzm8Ofk6Njp0OrI67zstO2s7qjvoPCc8ZTykPOM9Ij1hPaE94D4gPmA+oD7gPyE/YT+iP+JAI0BkQKZA50EpQWpBrEHuQjBCckK1QvdDOkN9Q8BEA0RHRIpEzkUSRVVFmkXeRiJGZ0arRvBHNUd7R8BIBUhLSJFI10kdSWNJqUnwSjdKfUrESwxLU0uaS+JMKkxyTLpNAk1KTZNN3E4lTm5Ot08AT0lPk0/dUCdQcVC7UQZRUFGbUeZSMVJ8UsdTE1NfU6pT9lRCVI9U21UoVXVVwlYPVlxWqVb3V0RXklfgWC9YfVjLWRpZaVm4WgdaVlqmWvVbRVuVW+VcNVyGXNZdJ114XcleGl5sXr1fD19hX7NgBWBXYKpg/GFPYaJh9WJJYpxi8GNDY5dj62RAZJRk6WU9ZZJl52Y9ZpJm6Gc9Z5Nn6Wg/aJZo7GlDaZpp8WpIap9q92tPa6dr/2xXbK9tCG1gbbluEm5rbsRvHm94b9FwK3CGcOBxOnGVcfByS3KmcwFzXXO4dBR0cHTMdSh1hXXhdj52m3b4d1Z3s3gReG54zHkqeYl553pGeqV7BHtje8J8IXyBfOF9QX2hfgF+Yn7CfyN/hH/lgEeAqIEKgWuBzYIwgpKC9INXg7qEHYSAhOOFR4Wrhg6GcobXhzuHn4gEiGmIzokziZmJ/opkisqLMIuWi/yMY4zKjTGNmI3/jmaOzo82j56QBpBukNaRP5GokhGSepLjk02TtpQglIqU9JVflcmWNJaflwqXdZfgmEyYuJkkmZCZ/JpomtWbQpuvnByciZz3nWSd0p5Anq6fHZ+Ln/qgaaDYoUehtqImopajBqN2o+akVqTHpTilqaYapoum/adup+CoUqjEqTepqaocqo+rAqt1q+msXKzQrUStuK4trqGvFq+LsACwdbDqsWCx1rJLssKzOLOutCW0nLUTtYq2AbZ5tvC3aLfguFm40blKucK6O7q1uy67p7whvJu9Fb2Pvgq+hL7/v3q/9cBwwOzBZ8Hjwl/C28NYw9TEUcTOxUvFyMZGxsPHQce/yD3IvMk6ybnKOMq3yzbLtsw1zLXNNc21zjbOts83z7jQOdC60TzRvtI/0sHTRNPG1EnUy9VO1dHWVdbY11zX4Nhk2OjZbNnx2nba+9uA3AXcit0Q3ZbeHN6i3ynfr+A24L3hROHM4lPi2+Nj4+vkc+T85YTmDeaW5x/nqegy6LzpRunQ6lvq5etw6/vshu0R7ZzuKO6070DvzPBY8OXxcvH/8ozzGfOn9DT0wvVQ9d72bfb794r4Gfio+Tj5x/pX+uf7d/wH/Jj9Kf26/kv+3P9t//9wYXJhAAAAAAADAAAAAmZmAADypwAADVkAABPQAAAKW3ZjZ3QAAAAAAAAAAQABAAAAAAAAAAEAAAABAAAAAAAAAAEAAAABAAAAAAAAAAEAAG5kaW4AAAAAAAAANgAAp0AAAFWAAABMwAAAnsAAACWAAAAMwAAAUAAAAFRAAAIzMwACMzMAAjMzAAAAAAAAAABzZjMyAAAAAAABDHIAAAX4///zHQAAB7oAAP1y///7nf///aQAAAPZAADAcW1tb2QAAAAAAAAGEAAAoC8AAAAA0OXuAAAAAAAAAAAAAAAAAAAAAAD/wAARCAFyAlIDASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9sAQwACAgICAgIDAgIDBAMDAwQFBAQEBAUHBQUFBQUHCAcHBwcHBwgICAgICAgICgoKCgoKCwsLCwsNDQ0NDQ0NDQ0N/9sAQwECAgIDAwMGAwMGDQkHCQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0NDQ0N/90ABAAm/9oADAMBAAIRAxEAPwD9qKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigD/9D9qKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigD/9H9qKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigD/9L9qKKKKACiiigAooooAKK8N+KX7Qnw7+FLPp+rXL6jrQXcuk6eBLcjOdplJIjgVscGRhnsK+MNU/bA+M3jS/fTvhzoVtZYcp5dpZz63drk8eYyBI4W55BUj3wK7aGX1qq5oqy7s+2yLw9zvNaX1ijT5KW/PN8qt3XVrzSa8z9QMj1oyPWvy2Gq/t06mxuYY/FESv8AMEW30mBRn0WQbh9CTT/M/bv9PFX/AHzov+Fb/wBl/wDT2H/gX/APcfhfJaSzPCp/9ff/ALU/UbI9aMj1r8ufM/bv9PFX/fOi/wCFHmft3+nir/vnRf8ACj+y/wDp7D/wL/gB/wAQwf8A0NML/wCDf/tT9Rsj1oyPWvy58z9u/wBPFX/fOi/4UeZ+3f6eKv8AvnRf8KP7L/6ew/8AAv8AgB/xDB/9DTC/+Df/ALU/UbI9aMj1r8ufM/bv9PFX/fOi/wCFHmft3+nir/vnRf8ACj+y/wDp7D/wL/gB/wAQwf8A0NML/wCDf/tT9Rsj1oyPWvy58z9u/wD6mr/vnRf8KPM/bv8ATxV/3zov+FH9l/8AT2H/AIF/wA/4hg/+hphf/Bv/ANqfqNketGR61+XPmft3+nir/vnRf8KPM/bv9PFX/fOi/wCFH9l/9PYf+Bf8AP8AiGD/AOhphf8Awb/9qfqNketGR61+XPmft3+nir/vnRf8KPM/bv8ATxV/3zov+FH9l/8AT2H/AIF/wA/4hg/+hphf/Bv/ANqfqNketGR61+XPmft3+nir/vnRf8KPM/bv9PFX/fOi/wCFH9l/9PYf+Bf8AP8AiGD/AOhphf8Awb/9qfqNketGR61+XPmft3+nir/vnRf8Khbx7+2v4QR73VrXXp4U5Y32k2l/EAO22xAk/Gj+y30qwfz/AOAC8LastKOY4aUuyq3b/wDJT9TaK/OnwR+3Ffw3X9m/Erw+jCIqs91pG+O4hJ6tLYz/ADqP9lXZz2U190+DPHXhP4g6Mmv+D9Th1Kzc7WaM4eJ+6Sxth43HdWAIrlxGDrUf4i07nyvEHBmb5Lrj6LUXtJaxfbVbN9nZ+R1tFFFcp8uFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRQSAMntzQAUV82/EL9p3wJ4Nlm0zRA3iPU4iyPHaOFtYnXqJLg5XIPUIGI9K+dZP2ivjv45naPwPpiwxM21RpmnyXpH1uJP3fHuBX0WD4Wx+Ih7VxUI95O3/B/A+HzTxCybBVfq8ZurU/lprmfpfSN/K9/I/RzIoyK/OhLr9si4HmKPEAH/XLT4/0IFP3/ALZX/Uwf986dXV/qo/8AoKpf+B/8A4l4iJq6y7E/+Cv/ALY/RTIoyK/Ovf8Atlf9TB/3zp1G/wDbK/6mD/vnTqP9VP8AqLpf+B/8AP8AiIf/AFLsT/4K/wDtj9FMijIr869/7ZX/AFMH/fOnUb/2yv8AqYP++dOo/wBVP+oul/4H/wAAf/EQ1/0LsT/4K/8Atj9FMijIr869/wC2V/1MH/fOnUb/ANsr/qYP++dOo/1U/wCoul/4H/wBf8RD/wCpdif/AAV/9sfopkUZFfnXv/bK/wCpg/7506jf+2V/1MH/AHzp1H+qn/UXS/8AA/8AgD/4iGv+hdif/BX/ANsfopkUZFfnXv8A2yv+pg/7506jf+2V/wBTB/3zp1H+qn/UXS/8D/4Af8RDX/QuxP8A4K/+2P0UyKMivzr3/tlf9TB/3zp1G/8AbK/6mD/vnTqP9VP+oul/4H/wBf8AEQ/+pdif/BX/ANsfopkUZFfnXv8A2yv+pg/7506jf+2V/wBTB/3zp1H+qn/UXS/8D/4Af8RD/wCpdif/AAV/9sfopkUZFfnXv/bK/wCpg/7506gy/tlKNx/4SDj/AGdPP9KP9VP+oul/4H/wA/4iGv8AoXYn/wAFf/bH6KUV+cD/ABR/an8F/v8AXrS/mgHJOoaSs8f4vahcD6kV6X4I/bD0m+kS18d6UdPBO032nM1zAp7+ZEQJU99u/FZ1+EcfGHtKPLUX9x3/AMr/ACuaYTxKyipVVDFqdCb6VI8v4ptL1bSPtSisvRtb0jxFpsOsaFeQ39lcDdFPbuJI2/EdCO4PIrUr5mUXFuMlZo+/hOM4qcHdPZoKKKKkoKKKKACiiigD/9P9qKKKKACiiigAr4h/ab/aTuPB0s3w7+H1wF10hV1LUIx5h08SgbYYVw267kDDAwfLBBwWZVPv/wAdPiV/wqr4b6l4ntgr6k+yz0yN+Va9uDtjLDqUj5kfH8Kmvjr9kP4Qr4n1K4+MXi4NeR215Ouli4AY3Oobybm+k9SkjMqcAeZvbBGwj0sFRhGDxNZXS2Xd/wBf1ofpnBWTYHD4OrxLnEeajSfLCH/Pypul6Lta27d1Fp3/AIL/ALIT6tEni74xicC8IuU0TzW8+VnwTJqE4YuztxmMNk8+YzZwPv7QvD+h+GNOi0jw5p9tpdlAoSO3tIlhjVR0ACgfrWvRXNicXUru838uh8xxHxbmWd1nUxtT3ekVpGPov1d2+rEwKULnoM/hWL4i16x8M6PcazqAkeOAALFCpkmmkY4SKNBks7twoHevmrV/DXx9+JkjXF9dReE9Kk5hsBcMsqp283yPnZsdcyKO22vnM1zh4RqlRpSq1HqoxX4yk9IrtffWydnaMl4f+vJ1sRXjRpJ25pvd9oxWsmuttFpdq6v9YbPQZ+nNNwPSvjaH9nz4n6RIL3RfGKC5TlcS3cOT9TJIp/FSPavWvA3izx/o+oweEvirY7Zrj93YazAQ9tcyD/ljKyAKkrD7pKoHPGM4z5+B4ixEqipZjhJUbuyd1KPknKPwt7K+nnsenmPCmGhRlWyvGwr8qu42cJ2W7UZfElu7O9tbWue34HpRgelLS4YHBBB9CK+qPihuB6UYHpTiCDgjB96T2FACYHpRgelLUc00NvGZp5EijXq8jBFH1LEAUAPwPSjA9KFZWUMpDAjIIOQQfQiloATA9KMD0pahe4t4pEhkljSSXIRGdVZ8f3QTlvwoAlwPSjA9KWigBMD0pRkdCR9OKKKAPMPiH8HPh18UbTyPF+jwz3KKRDfwfuL6Dd1Mc8eHGcDIJIPQgivzo8bfDr4ofsreLIPGfhLUnudInkWCPURH+6lBPy2mpwLhTuz8jrgEkhDGxCt+s1ZOu6HpHibR7zQNetY73T9QhaC5t5RlJI3GCD/Q9jzXdhcdOl7ktY9j7nhXjnGZU/quI/e4WWkqctVZ78t9n5bPr3Xn3wd+LWhfF/wmniDS1FtewMINSsGfc9pc4BK543RuDujfA3KexyK9Xr8lbB9Y/ZP+Pgtbu4km0GbZHNI5IF5oly/ySuTwZrN8lmycbXOB5uB+s8bpKiyRsHRwGVhyGUjII9iOaWNw0aclKn8MtULjnhuhlmJp4nAS5sNXXPTfl1j3urrfWzV9bj6KKK4j4cKKKKACiiigAooooAKKKKACiiigAooooAZJJHDG00rKiIpZmY4VVUZJJPAAHJNfnh8WvjT4m+LGuj4c/C9LiTTLiQwA2pKT6mw+8S+R5dqOvUbh8zHaQG9E/ay+Js2mafB8NdHkZZ9UiE+ptHneLQnbHAMc5nYHI6lFI716T+z98H4Phz4cXWNWhB8R6vEr3bNgm1ibDLbIe23rIR95/YAD7HLKFDLcIs0xUeacv4cX/wClP+u3dNfmGf4zF57mcuHsvm4UoJe2mt9fsL1699b6Ralyfww/ZX8L+HYINT8erHrmqYVvsn/MPtiP4QmB5xH95/l/uqK+q7a3t7OBba0iSCFAAscShEAHYKoAFTUV87j8zxONqe0xM23+C9Fsj7jJshwGVUVQwNNRXV9X6vd/P5aBgUYFFFcB64YFGBRRQAYFGBRRQAYFGBRRQAYFGBRRQAYFGBRRQAYFGBRRQAYFGBRRQAYFGBRRQAc4wDwa8e+IPwM+H3xEjkm1GwWw1MjCalYqsNwD23gDbKv+y4New0V0YbFVsPNVKEnGXdHHj8vw2Nouhi6anB9Gr/8ADPz3PzEntfij+y74sjmhkW70m+kJyAwsNSVeWV1yTDcBe/LDqC6ggfoR4B8daD8RfDVt4n8PyEwzZSWGTAlt50+/FIB0ZT+DDBGQav8Ai3wnonjbw/d+GvEEAns7xMHs8bjlJIz/AAuh5UjvX56/D/Wtb/Z4+MU/hHxFMTpV5NHa3jAYikgmOLS9UdivRsdAWHRRX2EnDPcNKdksTBX0+2v8/wDgLZ6fmEIVeD8fTpKTlgKrsru/spPzf2Xq35Xb1V5fpdRR+v0or4c/WwooooAKKKKAP//U/aiiiigAooooA/Nb9uXxHeX/AIp8LeBtPYM1rayX4jPQ3t9ILW0yB2x5g/Gvv3wH4VsfA/gzRPCGnJst9IsYLVR3JjQbmPqWbJJ7mvze/aPUXH7V+iwS8o0vhCPH+y1+xIr9T2++3+8f516uN93DUoLbf8n+p+qcaN4fhzKMDT+FwlN+blZ6+nNJLyY2iiivKPys858aeNYfC+pW1tPYpc7bT7aszyBPKzeW1mTypxhbktu7Bcd81zup/Fe5hvrqy0qwsp/sa6/K8txqCwRtHoEsMcgUhGy8nmn/AGYyvzHmur8VweH9HmuPGPiSOS6s/sX9m3QkUTW9rZSuHldodp3I7hPOY7sKoOAoNLa+FfAPiCTT9bsLTT762s4JorVII4pLQC4eORnCBSu/dEuG69fWsVWp+0dJNc29utu/p5nRLDVvYqu4vkva/S/a/fW9vO/U4LX/AIu6kmkeI4/D+nwx6xomj3erst7IRFDb+THLYyyKq7sXO9ht4CtDIMnbV+/+Jd6l3q2nT6Aby20WN11CaN2MKTxWIv8AeGdQpgJKxq3+s3ndt2gmvTbXw3p1te6rqMiNc3GsbFunuAJMwRKVjgAwAIUDNhTnl2JJJNctqF74B03W7HRILLT7nXLiH7FbWlvDG86Wirhg5UHyrdEJ3bjtwdoBLAErVqdON6rST016t7Lzb7dQw2GrVp8tCLbWunRLVt9kt29kdXp15f3/AIdg1DUoI7W6urPz3hhlMyRmRCwUSFULEAjJwOc44r8Kf2TvjV8RvAGhaz4Y+JmrXl34V+KumeKI/C2o3l9cXUllrmhedDJA08xZ4vtEQXAMuAyKyjc7V+9kkSm3a3jwqmMxrgYABXaMD0HpX5+v+wbpOq/srn9nXxPr9tf6laa1d6/pmvpp5RLS9ublpx/ozyuzRlHaGVRIC6McFTgjY5zhfh1+0Z458GfA39nD4a+BtKtvFXxC+J+jyG2l16+lgs4YbCIzXFxdTqk0zsc4UY5wTngAyfEf9pT4u+Lv2afjZZRaPYeFPiB8N55tH8QNZalM0cFlLDv/ALQ0ydFjl88KQYlYp8wyWHQel3f7GutWXgf4NxeDPGcOk+Pfgtby22k67LpxubG7huk8u4jnsjMrbZFxgiQleeoNdB4V/ZEa2+FPxT8G+PfFsuveJ/i/LdXHiDX7ezS0SJ5o/LgS2ttz7YoFA4d2LHPIGAADyrUf2qPiF8Ff2e/hXd+ObPw7N4w8cyWun6Zc6lrMkGlR6eLZJBqOqXkse9JNpHmqARvbgnv438V/2o5Pj3+yv+0H4J1y30Jtd8C2umrNfeF9ROq6HqNrfTRtFPaXDIjEo6Mrgjggc5yB9Q+JP2SfHHir4T+AfDGpeP7OPxt8MtTivfD3iCHRF+wi3hjEC2lzYSTOZUMIwzebuLYbpkHT1/8AZm+KHjz4LfEL4afEP4gaVfap43+zR2t3p3h2LTrDSLeAxs0UcEcvnTCRkZiXlwGbIAwcgHGfDv8AaP8Aix4T+J3wo+D3xY8F6Ro+i/EfRHPhy90zU3vru3Om2qSFb9WjjjVmTG4RFlUsArMATXHa3+3N8UdP8I638e7LwFo83wd0DxN/wj890+qSrr88KXS2b3kNuITBt8xvlQuCehIHzD1Dwr+yX8QR8Z/AHxO+KPxKXxXpvwu06ey8NabFpMdhMZLmAQPJeTI7LL8oydiIWKp90Bg3Dav+wd4r1LR9a+EUfxKEfwb1/wATHxJd+HW0kHVU3zrdPZxagJgqwmddwbytwGOp5IB0/jb9qn4u6l8RviH4O+A/grRtf0/4V6XZalr17rmpS2Ul697bm6WCwSKKRSwiGMyMvz5zgAFvnn45ftCeAdc+KP7IX7RGuTDQvDt3a+JtYuTckM1qBawq8JYcM4lxGuPvMRX0v8Qf2SfHk3j/AMc+Nfgh8QrfwRF8TNJtdK8S2d7pP9qD/Q4TbxT2bieHynEJK4IIySe/Emt/sOeDNUvPgXp/2u2vPC/wZgvbeXStWsRdtrQu0iG+Rg6RxsJY/Nb92wLYxtxyAeyfs2fFLxx8a/h8nxQ8UaJYeHtI1+V7jw3ZW9wbq8bStxEU964JiSaUDPlR52D7xydq/QlfMv7NvwD1r9niz8TeELXxLHrHgu+1e41Tw7pbWjQ3GipdyNJNbed5rJNESQVwiFW3Elt3y/TVABRRRQAUUUUAfD/7cXhC21DwRo3jTyw0uj6gtjcHAJez1P8AdlD7CYRt+B9a9q/Zp8U3Hi74KeGdQvZDNd2cEmmXMjdXl0+RoCx/3ggNZ37VkMcvwD8VNIobyIradM9niuI2U/ga4f8AYldm+Dt0hJITxBqYX2BcH+Zr1G+bAJvpK34XP1OpL61wBGVTejX5Y/4ZRcmvvd/uPr+iiivLPywKKKKACiiigAooooAKKKKACiiigAo47nA7n2orJ16Z7bQtSuI+Gis7l1+qxMRTirtImUuVOXY/OzwLB/wuH9pW51q/AmsrW9udSZTll+zacwhtU9AC4V8eoP1r9KeTyetfn1+xVbJJrPiS9YAyRaZp0YPfEskrN+e0V+gtfWcZT5cesLH4acYxX3X/AF/A/OfC2lz5O8wn8dec5yfnzOP5pv5sKKKK+SP0g/PX4uftEfHO+/aC1b4B/A2TwPoVz4b0e31S8vvG08o/tGW7UtHBaRxMmFXgOxJODnAwN3RfGX9qPx/+z3+z1ofjv4oeHdKX4havcxWiaPYT3FxpLP5qiRmvEjbyT9nPmKp3Av8Au0Mhxnz/APan8O6X4r+IuoWXxE/ZkvviTpX9lQ2mi+KtBuIf7RkkYOXhnYSQy28cTv8AuySxUlmHXjzS+/Z4+O8P/BO62+Fd/pF3qXiuz1yz1a10AXKXF5Z6VDqKTx2aySSBHkggBIQPgD5F6AUAfbGvftgfAHwlo/hTXPF/iKXRLPxlY3moaVLf6fd26tDp5CzmYPEDAQ5CIsoVpGZQgJYA1PFv7Z/7PHgqHSJdb8QXbNrOhReJoYLLSb69uINHmVXS7u4oIXe1jKsCfOCkdxxXi3xD8IeN/it+0f8Asy/E8+B9VstD0O28QSa9BqkELPpL3FuI7dbtVeSMO0iArtLY4NeUfFj4X/EvwN+098Q/H8Wg/EPXfCfxF0fS7e2l8AT2AeOS1hNvLYXkd5ykBA3AoFA3kgk52gH3f43/AGlvgx4A8LeGvGOt681zpnjEj+wW0mzudUm1EeUZi0MNpHJIVWMZZiAF6Hnisq0/aw+A9/8AB3U/jxY+I2uPB2i3QstSuY7K5a6s7nzUh8qez8v7TG4aRCQU+4wf7hzXx54p8P8Ax1+D/gD4G+DPhP4c8c6N4HtU1CfxVbaB/ZeteLtJnuAXtrSGe9VrYRiSVhvC8RqVLDAVvFdN+Bvx00z9lz9pH4Wah4I8Q3Gv+JvFdtr2jvObWabVYL6aAsd1u6xNcxLGZLry1WMFsRg/dAB9/wAP7dX7L9xf32lw+L3a7tLOO+ihGnXhk1CCT7raevk7r3IycQB/lDN91SR1L/tb/AJPhRp3xq/4Sbd4X1a//smxdLO5e9uNR3mP7JHZLGblpwVJKBM7Ru+7zXlEPw28Ur+1p8LPGw8PTroujfDS60q81HyFENrfN5YjgduqyABwBjjJA6mvlux+Bvii3/ZnTQfGPw78aXOpab8VNV8R2Z8K3FrZa5pMJkkeDUbaO5by5wykRiPB+9u+UqDQB+nXwp+Mvw++NWjX2t/D2+nu4dMvW06+iu7K4sLm1u0VXMUsFzHHIrbWVunQivljx/8AGL9p3Rf2s9B+B/hlfBB8NeJLd9Zs5r22vWv49KsjEt2ksiTCP7Sx83ycRlMbN3eur/Y2Hx3/ALD8XT/F86+dHl1dH8Kf8JbHZxeIWsjCvmteCyzHzJwu75+D/DivMPiTc/E3/huLwh4+0v4YeKNT8M+FtJuvDlxq1qlsbad9TaNhdRFpgxt4dx8wlQ3ynaDxkA9y8b/trfs5fD3W/EXhzxN4juodS8KXa2mrWsGl3l1Lal0EgkYQxP8AuApGZvuBiFzuIFdX47/ag+Cnw8svDt7reuy3o8WWX9p6NFo1jc6vNd2G1X+1LFZxyuIMMP3hAXnFfIGr/CD4iXMH7ZYh8MXrzeO4PK8Nv5K51bFiYgtuc5YCX1wATn1rgPFfww8ZeHvh78C9QPgn4m2Pivwn4Bt9Im17wFcWH22yuiiK2nXdldtho9wLvIAMcAFucAH0p8bf25fAXw58G/Dnxt4MH/CRad8QNStvIme1vI1XRzII7q4RRFuNzEzKEtmAkk52g4xXpHjf9sn9n/4ejTB4l1rUI5NW0aLxDFBBouoXNxBpU/8Aq7i7iigZ7VW6YmCsCCCAQa+PPiL8N/2o/GH7MPwT1bx/oNz4j8d+B/H2neJtb0uy+zDUW0q0ln8pQiGKB7hYmQOFKjuec10vxp1D9pfx98TNT0q00Px/o/hHxR4PtItBsfDltoihNRvEdrq28QXl6sxhSJn2yrDIVAU+VliCwB9l+L/2lfgn4F8K+FfHPiXxNDb+H/GkqQ6LqUcUk1tcM8RmBZ41PlrsBOXwARjrxVX4f/tQfBD4maH4p8Q+GvEYhtPBI3+IBqltPpc2nRbGkE08N2kcixMqttcrglSOoIr4C0D4R/Fa7+CH7LnhbVfAurW+o+APHUb69aXMMUjWtramTF42HKeRudWRsk8ZAzit3xt+zx8V/G3xB/bCt9L0e5soPiFpOhR+Gb2XbFb6nNYwAywpITgBmXYS2Bk9e9AH2d8Lv2tfgN8YtYl0DwT4hla/TTzq0cWpafd6X9o09c7rm3N3FEJolwSzJnA5NU/An7Yv7PXxI8YWHgfwl4lkn1HWHu49JkuNPu7Sy1N7Ftk62V3PEkFyUPaNzkcivE/Cus/FD4x+Ab74MTfCDWPh8B8Pr3w9c+J9Ye1tjaakbdYIrewERlklgkYs+/KBduSOa+dvgP8ACHxnqDfCX4d/E/wH8Vra/wDh5f2t093dajpkPhnSJ7EMyz2jRZmnjlZVjeLc2UdsE9wD9kqKCckn1OaKACvi79sXwfDd6Ho/jWFP3tnOdMu2HU290CY8/wC7KMD03GvtGvD/ANo+zjvPgr4mEn/LvBFcofR4ZkYH869rh3EyoZlRnH+ZL5PR/mfLcbYCGMyLFUpq9oOS9Y+8vxRs/BDxVL4x+Fnh/Wbpy90tt9kumPUz2jGFyfrtz+Ner18q/shXMkvw11C2c5W21u5VPYSJHIf1Y19VVjneHjQzCtShspO3p0OnhPGTxeTYXEVHeThG/m0rN/egoooryz6AKKKKAP/V/aiiiigBrkqjMqliASFHUkDOB7npX4Z+Iv8Agpj+09puvajp0XwpsLJLa6lhW2u7HVZLiEIxASV408tpFxhinyk9Mjmv3OpCATk8/UZoA/l4+If7Z/xp8XfFix8e6v4PsLTVbaTSHjtEtL9UY6dP5sI2yASfvGODgcj7vNfUR/4Kf/tREkn4aaOMknJ07WMD/wAc7V9M/tEAf8Na6HwP+Pnwf2/6fjX6msq724H3j2HrXfi6co06blK9192x99xhgcTQy/Laleu6kZ07xTSXIrR0XfdfceXfBPxx4g+JXwo8L+PPFOhyeG9W1vT47u60uUMGtpHzxhwHAYDcoYBgDzXqNFFcB8CNdElRopVDo4KsrDIIPBBB4INfL/if9nN4r6XVPhprMmgtKxdrPzJY4Ax5/dvEyso/2WDgdBgcV9RUV5ObZHgsygoYuF7bPZr0a1/R9T3Mk4izDKakp4Gpbm3TScX6p6P13XRnxrH8D/jPqTi21vxaYbbPLC+upyR7IDHn8Wr6A+Hfwt8N/Di2lOmB7rULoD7Vfz4M0uOdoxwqZ/hHU8kk16VRXFlvCmXYKt9ZhFyqLaUm5Nel9F2ulex6Wccb5rmNB4WrJRpveMIqKfrbV+jdutgooor6M+RCiiigAooooAKKKKACiiigAooooAKKKKACiiigD59/ao/5ID4w/wCvWL/0clcB+xJ/yR+8/wCxh1L/ANCWu/8A2qP+SA+MP+vWL/0clcB+xJ/yR+8/7GHUv/Qlr04/7g/8X6H6jhf+Tf1v+wlf+kI+waKKK8w/LgooooAKKKKACiiigAooooAKKKKACsTxN/yLWr/9eF1/6KetusTxN/yLWr/9eF1/6KetKXxr1M638OXoz4e/Yn/4/fFX/XhpP/oU9fflfAf7E/8Ax++Kv+vDSf8A0KevvyvpeM/+RxV/7d/9JR8D4V/8kxhv+3//AE5MKKKK+WP0IPpRiiigAooooAKMUUUAFFFFABRRRQAUe9FFABgUe1FFABRRRQAUf0oooAKKKKACvHP2g/8Aki3i7/sHn/0NK9jrxz9oP/ki3i7/ALB5/wDQ0r0Mq/32j/ij+aPH4h/5FWK/69z/APSWeYfse/8AJPtZ/wCw5N/6Iir6yr5N/Y9/5J9rP/Ycm/8AREVfWVdvEv8AyNK3qeVwJ/yT+E/w/qwooorwj60KKKKAP//W/aiiiigAp21z0Vj+B/wqNwxRgh2sQQpIzg44OO+DzX4P+JP2OP8AgpBf+INSvY/iVqF2txdzyi4t/Gl7YxTB5GYOlqqbbdSCMRL8sYwo4GSAfTH7RKt/w1toeVP/AB8+D+x/5/jX6nMj72+RvvHsfX6V/LL4/wD2d/2vdA+Kdl4W8VeK7668SzSaSsN0/ie5uXVrubZaYuWUOuyTkHH7vqK+nB+xd/wUhSQF/iDqrhW+ZT4+vvmweRkJnmu/F1XKnTTjay+/Y+94vzGticvy6lVoSgoU7JvaatHVeWn4n79UV5Z8EfD/AI/8KfCXwr4c+KerJrnivT9Nig1W/RzKJrhc5/eMFaQqMKXIBcjJHNep1wHwQUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAfPv7VH/JAfGH/XrF/6OSuA/Yk/5I/ef9jDqX/oS13/AO1R/wAkB8Yf9esX/o5K4D9iT/kj95/2MOpf+hLXpx/3B/4v0P1HC/8AJv63/YSv/SEfYNFFFeYflwUUUUAFFFFABRRRQAUUUUAFFFFABWJ4m/5FrV/+vC6/9FPW3WJ4m/5FrV/+vC6/9FPWlL416mdb+HL0Z8PfsT/8fvir/rw0n/0KevvyvgP9if8A4/fFX/XhpP8A6FPX35X0vGf/ACOKv/bv/pKPgfCv/kmMN/2//wCnJhRRRXyx+hBRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAV45+0H/yRbxd/wBg8/8AoaV7HXjn7Qf/ACRbxd/2Dz/6GlehlX++0f8AFH80ePxD/wAirFf9e5/+ks8w/Y9/5J9rP/Ycm/8AREVfWVfJv7Hv/JPtZ/7Dk3/oiKvrKu3iX/kaVvU8rgT/AJJ/Cf4f1YUUUV4R9aFFFFAH/9f9qKKKKACiivnzUv2sv2ZNI1C50nU/il4Utryyme3uIJNUhDxSxMVdGGeGVgQR2IoA+Nv2iP8Ak7XQ/wDr58H/APpca/U5vvt/vH+dfip8c/j78Etc/aa0fxNo/jnQr3SYbjws0l7BexvAotbwvMS4OB5a8t6Cv0ab9sP9lXcx/wCFs+EepP8AyFIfX616eOknRpJPp+iP0zjvEUamVZTGnNNxpWdmnZ2ho+x9H0Vl6JrejeJdHsvEHh6+t9T0zUYEubO8tJVmt7iGQZSSORSVZWHIINaleYfmYUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAfPv7VH/JAfGH/XrF/6OSuA/Yk/5I/ef9jDqX/oS13/AO1R/wAkB8Yf9esX/o5K4D9iT/kj95/2MOpf+hLXpx/3B/4v0P1HDf8AJv63/YSv/SEfYNFFFeYflwUUUUAFFFFABRRRQAUUUUAFFFFABWJ4m/5FrV/+vC6/9FPW3WJ4m/5FrV/+vC6/9FPWlL416mdb+HL0Z8PfsT/8fvir/rw0n/0KevvyvgP9if8A4/fFX/XhpP8A6FPX35X0vGf/ACOKv/bv/pKPgfCv/kmMN/2//wCnJBRRRXyx+hBRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAV45+0H/yRfxb/wBg8/8AoaV7HXjn7Qf/ACRbxd/2Dz/6GlehlX++0f8AFH80ePxD/wAirFf9e5/+ks8w/Y9/5J9rP/Ycm/8AREVfWVfJv7Hv/JPtZ/7Dk3/oiKvrKu3iX/kaVvU8rgT/AJJ/Cf4f1YUUUV4R9aFFFFAH/9D9qKKKKAGuiyIyOMqwKkeoIwf0r8q7r/gkn8D2up30vxj4t06zeWR4LONrGVLeORy4jV5bZpGVN2AXYsR1JPNfqtRQB/Oj8Vv2Efh/4L+Oem/Dmx8T6/c2d3NoKNczi089RqdyYZMbIQnyAZXK9eua+1P+HSfwfjk3J4+8YAq2Qdum9Qf+vWr/AO0R/wAna6H/ANfPg/8A9LjX6mt99v8AeP8AOvRxlKEKVOUVq1r9yP0TjTK8Lhctyyrh4csqlK8nrq7R1/FnnXwn+Gnh/wCDvw58P/DLws9xLpfh2zSzt5btw9xKFyWkkKqq73YljtVVGcAAcV6HRRXnH52FFFFABRRRQAUUUUAFFFFABRVe5vLSyQSXk8cCnoZXCA/TJp/nQ7EfzF2yY2HcAGzyMHvkVPPG/LfUV1sS0UUVQwooooAKKKKACiiigAooooA+ff2qP+SA+MP+vWL/ANHJXAfsSf8AJILz/sYdS/8AQlrv/wBqj/kgPjD/AK9Yv/RyVwH7En/JH7z/ALGHUv8A0Ja9OP8AuD/xfofqOF/5N/W/7CV/6Qj7BooorzD8uCiiigAooooAKKKKACiiigAooooAKxPE3/Itav8A9eF1/wCinrbrE8Tf8i1q/wD14XX/AKKetKXxr1M638OXoz4e/Yn/AOP3xV/14aT/AOhT19+V8B/sT/8AH74q/wCvDSf/AEKevvyvpeM/+RxV/wC3f/SUfA+Ff/JMYb1n/wCnJBRRRXyx+hBRRRQAUUUUAFFFFABRVa8vLTT7Wa+v5o7e2t0Mks0rBI40XkszHgAepr4l+I/7W/l3D6P8MLRJzv8AKGqXiMySMen2a3BDSZ/hLkZ6qrCvTyzJ8Vj58mGje276L1f6bng59xLl+T0lVx1S19ktZP0X6uyXVo+5ACegJ+gpBycAZNfmnB4M/ag+KAXUNSudUgtZBuRtQvf7LhIb+7BCFfH+8masr+zD8cIgbi21axd/+mWs3YYn/e8vGa9x8NYSm+StjIKXZa/jdfkfJLjzMa37zC5XVlB7N3i7d7cr/Nn6RnI4PH1or8z5ta/ab+CxFxqz6jJp0ROWvMarp5Uf3pATLGPdmTFfSXwo/aa8N+PLmHQvEkKaFrE+FhJk32V0x6LFKcFXPZHAz/CWrlx3C+Jo0vrFCSqU+8Xe3qv8r26noZV4gYDE4hYLFwlQrP7NRWv6P8r2v0TPp6iiivmj7sKKKKACiiigAooooAK8c/aD/wCSL+Lv+wef/Q0r2OvHP2g/+SLeLv8AsHn/ANDSvQyr/faP+KP5o8fiH/kVYr/r3P8A9JZ5h+x7/wAk+1n/ALDk3/oiKvrKvk39j3/kn2s/9hyb/wBERV9ZV28S/wDI0rep5XAn/JP4T/D+rCiiivCPrQooooA//9H9qKKKKAGu2xGfBbaCcDknAzge57V+IGrf8FgNas9VvbO3+FtlbR291PAkOp65Ja30axSMgW4hFqwim+X50DHaeMnrX7g1jS+HPD08rTz6VYSyOSzu9rCzMT1JJQkk+poA/mf+Iv7f+q+OPi5YfEyTwdp9pJZy6O4tI9WeWJv7Ln84Aym3UjzM4Py/L15r6rH/AAWI195Mf8K10NdzdT4klwMnv/ofSvcP2hNF0aL9rDRIItPtEiNz4QzGtvGqHdfEH5QuOR19a/UNvC/hje3/ABJtO+8f+XOD1/3K78XCap03KV01p5bH3vF+ExtHL8tniq/PGVO8VypcitHS636avscN8EPiYfjJ8JPCnxRbS5dFPiXTYb82Er+YYDJnIWQqhkQkZR9q7lIOBmvVKaiJGixxqFVQAqqMAAdAAOAB6U6uA+CCiiigAooooAKKKKAK13d29jbvdXT7I06nqSTwAAOSxPAA5JqvdWp1CKLfLcW8ZG6SKNvLdsj7rMPmXHcKR9aW1uY795z5QMdtOY43bB3PGPmYccbWJUH1Brx7xf4un1KeTTdOkMdlGSrMpwZmHXn+56Dv1NdmV5VUzOo6cH7nV/P8b9F19DxM6zqhl9D2tXW+iXVv/Lu/zukekah4q8MWLeVd3UUjxfwopmZcfQHB/Wue8U/Ev4a+Htci8IeMdTtbC9mtoLxYb5CsXlz7hHl2Xy1Y7TwSCK4vwBoH/CQ+KbOydc28LfaZx28uIggf8CbAr4k/aY1x9d+OXiuffujs7iHTowOgW0hQEf8AfbvX2dDhLBVsV9Tcn8Lk3p3SStb1PzzNuP8AHYTLf7Q9nG0pqCTT25W3fVdkr/gfqGCW02GXw89vJEFVovm3wyx/3VkUnGezfNj0rTjZ3jV5EMbEAspIO0nqMjg4r8k/hJ8ZfEXwq1OMQPJeaDI/+maXnKFSfmktweI5h1AGFfowydw/VCxurLxDa6Z4l0S7EltcRLNFKvMc9tMM4Iz16EHqpGOma+az/h2vlVRNPmg9F0X3dGl0vZ9D6rg/jLC59RlKlHlqRtzRve1+qfVfc11Wze/RR9KK8U+0CiiigAooooAKKKKAPn39qj/kgPjD/r1i/wDRyVwH7En/ACSC8/7GHUv/AEJa7/8Aao/5ID4w/wCvWL/0clcB+xJ/ySC8/wCxh1L/ANCWvTj/ALg/8X6H6jhf+Tf1v+wlf+kI+waKKK8w/LgooooAKKKKACiiigAooooAKKKKACsTxN/yLWr/APXhdf8Aop626xPE3/Itav8A9eF1/wCinrSl8a9TOt/Dl6M+Hv2J/wDj98Vf9eGk/wDoU9fflfAf7E//AB++Kv8Arw0n/wBCnr78r6XjP/kcVf8At3/0lHwPhX/yTGG/7f8A/Tkwooor5Y/QgooooAKKKKACgkAZPH14orwX9pDxpN4N+F98LCUxX+tSLpds6nDL54Jldcc5SIMR711YLCTxVeGHp7yaX9ehwZrmNPAYOrja3wwTb87dPV7LzPlj4zfEjXvjX40g+GngINc6QlyYIooztXULmI/PPK3T7PDglcjbwXIJ2V9QfCT4IeF/AOli/tXi1HxBJ8s2rSw+YsZB+dbVH4RAeA3JY8kmvLv2SvAtlY6Nc+Mr+3D3esRSQWmRxBp0T+Xxnn/SJFOMdUQdq+jvGHiv+w0XSdK2rclBlgBtgjxgADpuI6DoBX1uaV6lat/YWV6UoXUujcurbXRde79Ej8x4bwVJUZcV5771WprHryx+yoLo2tV2T6Pmb6u8n0qwtBb6zdRujD5jdspaTnPIwM/gKxdU8V+BdB0yxbWdRsNNsNWklhszcOIIrh48b1TOBwTz0zXh9vBeaxqMNsjNNd3kqxKzksxZzgEk+nX6V5N+2rc2tj4k8I+CLNsw6No0szJ23XMqopI9SsLfmaVLg/DfWKWFnNuUk7uy0il53erstX3OrH8fV6WAr5hTox5YcqSbd25SStfTZXbVux90W0dsdNP9kvHcRTKWiLyGaFw3bdlsqfQZAr4++NX7N1nrenXPjDwHp4sNTUPLe6OmDDdbeWeAD5Ul43ADCyeitg18ofDX4s+LvhdqKXGgzmbT2cG50uZj9lnTPzbV6RSYzh1xz94MK/Tu0+I/hDXPBum+NoNZttKsrxlMEt9MsCrOuQ9vLkgbwcqy9QRkdqxxWV5jw7iI18HJyg9EktH1s137a+mq0jA55kvGOAqYXHU+ScFdpte705oSsttnot7STTV/Av2ZPjVd+JI/+FeeLbkz6naxM+m3UpPm3NvFgPDITy00XYn5mTr8ysa+xq/LT403vhrwv8V7fxv8LdXsbzzpF1Ro7CZZFtr+F/3qNs4VbgHOM939a+1bP9pH4SXiWccWrvNe3ixAWdraT3EwmkA/dBI4ySwbIwM9PStOIclnUlDHYKk+WoruKT919dOi/wCDbSxfBXFVKhCrlOa4iPPQfLGbkrTi9rNvVr77NJ63PdqKarblDYIyAcMMEZ7EdjTq+LP1MKKKKACiiigArxz9oP8A5It4u/7B5/8AQ0r2OvHP2g/+SLeLv+wef/Q0r0Mq/wB9o/4o/mjx+If+RViv+vc//SWeYfse/wDJPtZ/7Dk3/oiKvrKvk39j3/kn2s/9hyb/ANERV9ZV28S/8jSt6nlcCf8AJP4T/D+rCiiivCPrQooooA//0v2oooooAKTIpH37G8vAbB27umccZ9s9a/CTxFoX/BXj+3tR+zya7NH9rn2SaZqOiRWLr5jbWtklPmrCVwUEmHC8NyDQB9I/tEH/AIy10P8A6+fB/wD6XGv1NYje3+8f51/Lb8QNM/b4j+KtlF43/tv/AITEyaT9n8690x597TYscPEfJH73O3PQ/e4r6eGg/wDBXzzMSjxYF3fNt1XQM4zzj5utd+Lq89OnGzVl9+2x97xfmrxeX5dRdGcPZ07XkrKWkdYvqtN/NH750V5Z8Ebf4m2vwl8K2/xlmjuPGsemxDWpYihDXXOdxiAjL4xvKfKWzt4r1OuA+CCiiigAooooAKhuZ1traW5flYY3kPbhAT/SpqZIAUYMu8YOVxnIx0x70pXs7AzmvEmovaeFrm9iHlSSQqFA6q02B+Y3V87gYGB2r37xRt1bwdPcWgO14Y7hARztQhiMeoGa8APIx6+lfovBfI8HKUXdt/hZW/U/I+P5TeNpp/Dy6fe7/ofTnwY0H7Hok+uzLiXUXxGSORBESB/302T9MV8T+I/2M/i3r3iLVtdfVtBB1LULu8O+W4JAuJnkGcRDorAEdsV9J6d8X/EOl2EGnW1jYmO2iWKMbXXhRgdG6muv/aI+I83w7+EOoapFIsGsanEum2IU/durpSGdc84iTfJ9FrB1s1wmYOVKylWdl10W33J6mlfB8P5llCjiuZww0XJ2vHVq79W2nZH4zajZtpuo3enNLHO1nczWxmgbdFIYHaMvGSASjFcqSBkV+in7IfiCfU/h7f6DdHeui6i0cOe0F0gmCfRWLADsK/N5VVFCJ91QFGfQDA/Sv0S/Y50i4tfBeuazMpWLUdUCQnH3ltYljcj/AIGSPwr6rjWMf7Klz7pq3rf/ACufmfhTOa4hj7LZxlf0t/8AJcv4H1lp1o9hYw2Tv5nkrsVsY+UE7evouB+FXao6fbz20Ui3Db3knml4JICu5Kjn0XHFXq/FKKtTikrabH9RR2QUUUVoUFFFFABRRRQB8+/tUf8AJAfGH/XrF/6OSuA/Yk/5JBef9jDqX/oS13/7VH/JAfGH/XrF/wCjkrgP2JP+SQXn/Yw6l/6EtenH/cH/AIv0P1HC/wDJv63/AGEr/wBIR9g0UUV5h+XBRRRQAUUUUAFFFFABRRRQAUUUUAFYnib/AJFrV/8Arwuv/RT1t1ieJv8AkWtX/wCvC6/9FPWlL416mdb+HL0Z8PfsT/8AH74q/wCvDSf/AEKevvyvgP8AYn/4/fFX/XhpP/oU9fflfS8Z/wDI4q/9u/8ApKPgfCv/AJJjDf8Ab/8A6cmFFFFfLH6EFFFFABRRRQAV8I/to3WoSSeF9NtY5ZI0h1C7Plxs/wC9HlxpnaDjhm6193V8afGH9qAeG9aTw98PEtNRmspv+JheXCmS3JQ4a3i2kbmzw8gOEPChmzj6PhWGI/tCFXD0+dxv1stU1q7O258N4iVcF/YtShja/slNxSaXM3aSdlG6vtrroj6S+HVlZaP4F0XTbWPy49PsILYkjBbyIxuPr1J689a8SvryTUL2e+mOXnkZz9CeB+AwK6L4YfH/AME/EaJdGmK6FrlwGX7DcMBHPIwO428uAsmTztOH9RXLSRPbyPbygq8TFGB6gqcGvouGsHXoV66xcbVNN/Nu7XdN9UfM8U5jhcVgcK8vqKVJXWnRpRsmt00ujs0evfBrQft+vTa3MuYtNTbHnoZ5QQP++UyfxFeK/H39nH4s/E34n6l4r0VdK/s2SC0trT7RevFL5cEZzuQQOFzI7/xHjBr1Dwp8SdQ8JaZ/ZdnYW0yGVpWkdnV2Z/XGRwAAPavdv+E/g074b3PxE8RRpaW9pZT30qIxYeXFu2gZGSXwAB6mssxxGY4PHvG04qz9yN9fPZO+rR0ZdgMmzTKlluIm1y3qStptpq2mrJP8Ln4n+OPBes/D7xReeEfEDWrahYiIzfZJvPiUyoHC79qHcFILAqMZHrX1H+yRqFrrKeKPh3rtvDf6ZcQQ6ittcxiWPczGGYbWyMMAh+uTXyP4g13UvE+u6h4k1ht19ql1LeXHOcSTNu2g+iDCD2UV9W/sbaTNN4u8Ra4FPk2umwWu7HHmTyl9uemQqZ+hFfX8RqSyeo6z95KLuv5rrb5n5JwNKD4ooxwifs5SmrP+Rxlo++m/mjrfiH+yLDdNc6p8N7pbaQuXj0m6YmDYQPlimO5oznOA25MYHFesfBP4DaZ8MrVNZ1rZfeJpo8POOYrNW6xW+Rn2eTAZ+2FwKxf2otb1Lw/8PZGsr6ezuNS1exhtWtp3hlCRKzygMhVgCFOQDyOtdX+zdLq1z8INGvtZubi7uLuW8nWW5leaQxPO/l5dyWI2gYyeK/Pq2PzKrkiq1qvuObjZ/E1ru+q0/wA2fsmX5RkWH4plh8NhrVY0+e9/dT5raR6S1TVtEtktz3WiiivkT9OCiiigAooooAK8c/aD/wCSLeLv+wef/Q0r2OvHP2g/+SLeLv8AsHn/ANDSvQyr/faP+KP5o8fiH/kVYr/r3P8A9JZ5h+x7/wAk+1n/ALDk3/oiKvrKvk39j3/kn2s/9hyb/wBERV9ZV28S/wDI0rep5XAn/JP4T/D+rCiiivCPrQooooA//9P9qKKKKACkwPSlr5z1H9rz9l7SNQudK1L4peF7e7spnt7iFtQTdHLGSrocZGVYEH3oA+P/ANogD/hrXQ/+vnwf/wClxr9TWA3tx/Ef51+KHxy/aE+B2u/tMaP4o0bxxot7pENx4XaS9hug0Ci0vC85LYx+7XlvQV+i7ftm/sobmP8Awtjwt1J/4/09fpXp46SdKkk+n6I/TOO8TRq5XlMac03GlZ2adnaGj7H0vRWToWu6L4n0ay8ReHL631PS9SgS5s7y0kWaCeGQZV43UkMrDoRWtXmH5mFFFFABRRRQAUUUUAZqSWMM/wDYyxiPdE0iIRhHQsd4X1IJ+Yds5rwjxR4bn8P3p2qWspWJhk7DP8DejD9RX0K0UTujyIrNGSUJHKkjBIPbI4rPQyXvn2GqWY2c88SQSoTxgnkMO6kcdiRXpZLnFTLq2qvGXTyW3kmtbdGtD57P8hp5lRVOTtJfC+3k/J/15/P/AIYsf7T8SaXYYyJruIN/uq25v0BrwX9sz4hf8JP8R4vB1lKWsfC0RjkAJ2tf3KhpT6Zji2J9WYV90aR4T0jRdah13TldXhD7Iy26IM6lQwzz8ueADXzun7JfhK81q51zxPruratNeXEt1cLmK2Ekszl3JeNd/JPQMMDjpivtaHEWXPGxxdVu0I6K2t3v5befU/Oc04OzmWVyy3CqN6k05PmsuWKuvPWWuib021PhLwF4B8Q/EbxDF4d8OxZlOGuLhhmKziPWWU9gP4V6ueBxkj9a9A8J23g/wbaeEfDGIo7G2FvDI/BLH78rYHLkln92pfCmg+DvB9sPDHhGzt7GKLLyQ2y5O4Y+aV+SXP8AttuNdBaWb28s1xPM0007ck/KqoudqKuSAADyerHk9gPmOJeJZZtJUqStSV+ut7bvz7Lprr3+n4I4IpZFSlOcuetPSUlsl/LHrvu9L2WiskXI0EcaxglgihQWOScDGSe5Pen0UV8+kfoAUUUUAFFFFABRRRQB8+/tUf8AJAfGH/XrF/6OSuA/Yk/5I/ef9jDqX/oS13/7VH/JAfGH/XrF/wCjkrgP2JP+SP3n/Yw6l/6EtenH/cH/AIv0P1HC/wDJv63/AGEr/wBIR9g0UV8+XmgftUvdzvYeOfAMNq0rmCObwjqUsiRFjsV3XW0V2C4BYKoJ5AHSvMPy4+g6K+c/+Ee/a0/6H34ef+Edqf8A8vqP+Ee/a0/6H34ef+Edqf8A8vqAPoyivnP/AIR79rT/AKH34ef+Edqf/wAvqP8AhHv2tP8Aoffh5/4R2p//AC+oA+jKK+c/+Ee/a0/6H34ef+Edqf8A8vqP+Ed/a0/6H34ef+Edqf8A8vqAPoyivnP/AIR79rT/AKH34ef+Edqf/wAvqP8AhHv2tP8Aoffh5/4R2p//AC+oA+jKPft0r5z/AOEe/a0/6H34ef8AhHan/wDL6vmH4reAf+CgGofGTwbqfgDxl4fW2tdPnTU76y06fStDSFrmJjDe2V1eag17M6qxiMXkui7h5i7s0AfpVWJ4m/5FrV/+vC6/9FPUugxa5BotjD4mubW81ZII1vbixge1tZbgAb3ihkkmeNC3Kq0jkDjcai8Tf8i1q/8A14XX/op60pfGvUzrfw5ejPh79if/AI/fFX/XhpP/AKFPX35XwH+xP/x++Kv+vDSf/Qp6+/K+l4z/AORxV/7d/wDSUfA+Ff8AyTGG/wC3/wD05MKKKK+WP0IKKKKACiiigDE8R6FB4l0a50O5ubq0gu18uWSylMExQ/eUSAEqGHBIwcdCK+Pvi5+zR4M8PfDzUdb8C2l4NS0tUuSsl1LceZaxH98ixsSu4R5K4GcjFfblNZVdSjgMrAgqeQQeCCO4Ir1MtzfE4KcZUZtRum1fR+p4Ge8NYDNaMoYqmnJxcVJpNxv1i3s09dD88P2cPDPwe8baVdaX4t0O2vNds7kMJpXkKz2txkwvgMFGxgY8gddpP3hX1F4w8HLpAS+0eJhZKipJHlnaIrxuJYlip7kk4PWvj34oeA/EXwC8fQeM/CC/8SSe4Z7NmXfBH5hzJY3CjrG38GeoC4+ZBn7N+F/xm8PfFK2UWERt71Yx9qsncNNA+Pmypxvhb+CRcg9G2txX12bYvFUK8c4w03Uw8+mr5W9097K+z2T07X/MuGMNgqlCfDeY040sXB6NRSc0r2kn9rTe7u17297eZHO07eTjj3PaoP2wvGieHfAvhr4Sae+Jr+GK81BQelraYEaHn/lrPg9OQhr6GuvA2gXV2l7FG1uySJIywnEbFWDYKnIGcc4xXknxC/Z00v4m+O77xt4m8QaiPtflIlnbRwokMEK7UiWRlZscsSeCSxPFXHiHL62Jo1qzajC7tb7Wltu2r+R14zhHN6GX4nC4RKU6vLG90koauW9t9I28/I/M7R9H1bxDqdto2h2kt9f3jhILeEZdye/oqjqzn5VHJNfrP8HvhtF8LfBMOhqUuNSnY3eoTLwst04A2qevloAEX2GaveCvh78PvhmBp/hexjt7y7AWSZy097OB/fkbL7B+CCvEfjl8f7Dwa93oPhC9W+8QNG0JljINrpasPndznbLcH+FeiDBbHfjzXO8Rn9WOAy6D5b6/5y7JdFrd2+U8PcN4Lg/DTzbNqqdW1lbp5RvZyk+9lZdldvw/9o/xRqHjfxzpXw40WU3s+kMLRmQ5WXV70gMBgfdgUhe+AW/umv0I8KeH7bwp4Z0rwzZ/6nS7OG1U+vlKAT+Jya+Qf2Yvg3eW06fFLxdHILiVWbSYbjJl/e533km75t7gkJnnBZjy2B9u15/EWIoUqdLKsK7wpKzfeXX+u7aPe4Ey/FVJ18+x8eWpXa5U/swXwr5q3qknuwooor5Y/RQooooAKKKKACvHP2g/+SLeLv8AsHn/ANDSvY68c/aD/wCSLeLv+wef/Q0r0Mq/32j/AIo/mjx+If8AkVYr/r3P/wBJZ5h+x7/yT7Wf+w5N/wCiIq+sq+Tf2Pf+Sfaz/wBhyb/0RFX1lXbxL/yNK3qeVwJ/yT+E/wAP6sKKKK8I+tCiiigD/9T9qKKKKAGuiyI0bjKsCpHqCMEflX5Z33/BJb4ET3k82neLvGOm2sksjw2cE9i8VtG7lhEjSWrSFEztUuzNgDJJ5r9TqKAP50/iv+wh8OfBXx0034c2HibxHc2V3NoKNc3D2huVGp3JhkwVt1T5FGU+Xr1yK+1x/wAEk/gtHJuXx542yjZB8zTjyDx/y51a/aI/5O10P/r58H/+lxr9TW++3+8f516OMpQhSpyit1+iP0TjTK8Lhctyyrh4KMqlK8n3do6/izz34V/Dbw78H/h3oHwz8KGdtK8PWaWdu904knkC5LPIyhVLuxLHaqrk8ADivQKKK84/OwooooAKKKKACiiigD5o+Ifi/wAX+HPHN3Jret6n4S8L28GntpepQaNFqei3EkrMt2NUmAe4t/LYKoJ+zxKjB/Nc7gmnbfFvxnZr4ifVvCf2uW18ZReGNItbC/t/MnS4WIo8ry+WiCMOZHYkkodqqzKA3a+KvhRo/i+7v5NU1TXIrDWI4YdU0q31B00++ih4CPEwYxq6/LKIWjEq8Sbqvx/DfQU1K91Lz75vtuu2/iIwNcZgiv7aJYg0S7cpG6qu9M7SRkYJOQDyvSfjp4YXx8PCVppUUJl1r/hHp3g1CGS5jv1i8x5DYj5/say/uDOMMX+Yx+V+8qrbfHfSrQ+K7x9JkvL7w/o0utizsNcttReWJJzbm2k2uILW68wDCszR4PEnytj1a1+GGi2PiK616z1DVoIL2+l1O40mO9ZdNkvZofJklMONwDr8zRBxC0hMpTzCWrjdF/Z88B6Jo8Xhe0v9Yl06x0H/AIRq1s5b/clro5aNordFVAf3ZhHlyvulPO52rNUaataK020J5V2Kr/GjxBpmq6hpuseBrizi0fXNJ0nVLiLUrSVFGuvClrdQKNjTohnX7SGEbptbyxLhd3T/AA9+K0Hj/W9Y0iCygtDpMk8UkR1CGW/haCd4Nt5YgLNbGXYZIid6tGQSwJxXR6p8PPDmsTatPei4La1e6Vf3WyYqPO0d43ttmPuruiXevRxkHrVTRfhpoukeKx40lvdS1TU4LS70+yl1K4FwbOzvpo7i4hifYsjK8kMf+teQqEVVIArQo9FooooAKKKKACiiigAooooA+ff2qP8AkgPjD/r1i/8ARyVwH7En/JH7z/sYdS/9CWu//ao/5ID4w/69Yv8A0clcB+xJ/wAkfvP+xh1L/wBCWvTj/uD/AMX6H6jhf+Tf1v8AsJX/AKQj7BooorzD8uCiiigAooooAKKKKACiiigAooooAKxPE3/Itav/ANeF1/6KetusTxN/yLWr/wDXhdf+inrSl8a9TOt/Dl6M+Hv2J/8Aj98Vf9eGk/8AoU9fflfAf7E//H74q/68NJ/9Cnr78r6XjP8A5HFX/t3/ANJR8D4V/wDJMYb/ALf/APTkwooor5Y/QgooooAKKKKACiiigDO1fSNL1/TbjR9atYr2yu0Mc0Ey7kdT2IP6HqD0r4R8e/sp+JdCvzr/AMKr6S6SFmlhtJJzb6hbk54guQVD8cDcVYj7zGvv+ivWyvOsVl8m6EtHunqn/XlqfOcQcK5dnMFHGQ96O0lpJej/AEaa8j8y7f44/G74dXQj8X6abiWBtjSazZSwykAYANzCVjcY6Ehvrmte4/bE8Rzyi4i0PTFnGNo/tC5aMEcD90oAP071+jTokq7JVDr/AHWAYfkeKpppelo/mJZ2yv8A3hCgP57a9V5zlNT3q+AjzXvo+VX72S/zPllwbn1FezwubSUP70FN/wDgTlf8j82rrxj+0f8AGYnT9Ltb6GynAWRdNtm0y1K54EtzIQ7KO4D9M/Kc4r3T4T/sq6V4ang13x/JBqt9CRJDp0A/0GBwcguSAZ2U8jIVM87Sea+wOwXsOg7D8KKjF8U1pUnh8HBUoPfl3fq/+AmdmXeH2Fp4hY3M6ssRVWzm/dXpHX7m2r6pIMAcCiiivlz9ACiiigAooooAKKKKACvHP2g/+SLeLv8AsHn/ANDSvY68c/aD/wCSLeLv+wef/Q0r0Mq/32j/AIo/mjx+If8AkVYr/r3P/wBJZ5h+x7/yT7Wf+w5N/wCiIq+sq+Tf2Pf+Sfaz/wBhyb/0RFX1lXbxL/yNK3qeVwJ/yT+E/wAP6sKKKK8I+tCiiigD/9X9qKKKKAGu4jRpGzhQWOBk4AyeK/IbVP8Agr38M7PUrq0s/AesXEEMzxxTSajZW7yIrEBmidt8ZOMlGwy9CAc1+vdeS6h8A/gZq19canqnw88L3d5dyNNPPNpFq8ssjnLO7GPLMx5JPWgD8Dvih+354V8cfGrTviXa+FL22trSXQ3a1e+tndv7LuDM3zqSg3g4GTx3r7DP/BYT4esxx8PdUyScZ1ewHX/gVQ/Hr4R/CrS/2otG0XTfB2hWmnyXHhQPaQ6dBHAwnvSsoKKm0+YvDcfMODX6Uv8As5fs+lnB+GfhLBJGP7FtOn/fuu/FxqKnTc3dW0/A+94voY+nl+XSxlVSg6d4JK3LG0dH36a+Ru/CH4maP8ZPhp4d+J+g21zZ2HiOyS8ht7sBZ4gxIKPtJUkMCMg4I5FekVVsrKy02zg0/TreK1tbaNYoYIEEcUUaDCqiKAFUDoAMVargPggooooAKKKKACoLm3gu7aW1uUEkM0bRyI3RkcEMD7EHFT0UAfF50/w38LrP4zeOfCmnadZal4X1WKLT7q7Mk1ppUM2k6d5jvGZlCWsPmtcTRq0YYbmypYtW7qfxE8caTruqeCYfFenXUdprPhS3j8TS2VuUji11p1ubOSOOVbc3KCBXhfK7VuYt0b4/efVQsrMfaMW8Q+1ndcfIv747QmZOPnOwBfmz8ox0rjdQ+G/g6+0qx0KHTLSx0yy1KDVPsVpbQx2000G4hZIgmxlbd83GTgc0AeDaT8RfGuq69p/gCXxZp1vH/bHiSzfxZFZ2+y9TRkt3htY45JDardD7S63LLuDfZJtkcZY+VwWifGDXJNN1jx7b3WlNruq+HvAEU19EgNjFBqeqX9i2rLB5rt9j2M1zEGkKbeDIVBavs6Xwd4RuNEh8NXGh6bLpFtt8mwezha0j29NkJTYuMnGB3q/LomizNI82n2jtLbCykLQRkvajkQtleYhn7h+X2oA+Tdb+K3jXStbufBieIFv7ey8Ry6Q+t6bY2k+oOh0UajHD9nMgtfNjmYCaRV4jKDYpYyLb0f4j/Ezw/D4X+Hfiqa1HjLxZpWgXOmmaWO7YT79viEF0CrcCxjAlVwAv71E6AZ9i8T/B3wvr40VLBY9Fh0I3RtrS1sbOawP2wKJC9ncwyQeYCoKyhRIuWAba7hjwh8HvDngy60Gewury5j8N2F9ZadBceV5UT6nN593OojjTa8hARY12wwxjZEirxQB6ycZOOmePpRRRQAUUUUAFFFFABRRRQB8+/tUf8kB8Yf8AXrF/6OSuA/Yk/wCSP3n/AGMOpf8AoS13/wC1R/yQHxh/16xf+jkrgP2JP+SP3n/Yw6l/6EtenH/cH/i/Q/UcL/yb+t/2Er/0hH2DRRRXmH5cFFFFABRRRQAUUUUAFFFFABRRRQAVieJv+Ra1f/rwuv8A0U9bdYnib/kWtX/68Lr/ANFPWlL416mdb+HL0Z8PfsT/APH74q/68NJ/9Cnr78r4D/Yn/wCP3xV/14aT/wChT19+V9Lxn/yOKv8A27/6Sj4Hwr/5JjDf9v8A/pyYUUUV8sfoQUUUUAFFFFABRRRQAVS1G1uL2wubO0u5rGaaJ0jurdUaWBiOJEEqvGWXqAysvqCKu1h+I/D9h4p0a40LUnuY4LnYTJaXElrcRvE4kR45oiroyuoIIPbByMggHgHhn4ieNNO1jVPAEZbxjr//AAkuqWOnS6jNBpiJpunWlncSy3MtvAV+WS6WNBHCWYuucKGYUl+MXjLVNRl8UeD/AA9Lq9jF4SOp3uiXGoQ2c1tdWV/cwXcULFHimuiYXiTMiQsUGZFVtw9Mi+C3hG3tbRbW71q31G01K91ddYj1OYanJeaihjuWlnJIljkTC+S6mJAkexF8qMrRvPgJ4DuNMt9HtJtb021h0U+HpFsNXuoHudOeRpZI7iUOZZZZXdme4LeeSzESAsxIB63pWpWus6XZ6xZEm2v7eG6hLDBMc6CRcjscMMir9QW1tb2dtDZ2kawwW8aRRRoMKkcYCqoHYBQAKnoAKKKKACiiigAooooAKKKKACvHP2g/+SLeLv8AsHn/ANDSvY68c/aD/wCSLeLv+wef/Q0r0Mq/32j/AIo/mjx+If8AkVYr/r3P/wBJZ5h+x7/yT7Wf+w5N/wCiIq+sq+Tf2Pf+Sfaz/wBhyb/0RFX1lXbxL/yNK3qeVwJ/yT+E/wAP6sKKKK8I+tCiiigD/9b9qKKKKACjIprlwjGMAtg7QeATjjPtmvww8R/FT/grBFr+pRW/h/U7eJbqYJDYaFYXFpGoY4WCV5N8kQH3XblhycHigD6F/aI/5O10P/r58H/+lxr9TWI3t/vH+dfy5/ELxz+3Te/Fix1Hxpp+qR+Llk0g28cmkWkUheKbNliJSUOZc7efmPB4r6hb4sf8FZ8ndo2vDJOSPDenZHrgb+td+LrKdOnFJqy/yPveL82ji8vy2jGnKPs6fLeSspaR1jrqtN/Q/eqivLvgpqHxJ1X4UeF9R+MFlFp3jK406J9YtYQqrHcnOQVUsqsVwWVSQGyBXqNcB8EFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAHz7+1R/yQHxh/wBesX/o5K4D9iT/AJI/ef8AYw6l/wChLXf/ALVH/JAfGH/XrF/6OSuA/Yk/5I/ef9jDqX/oS16cf9wf+L9D9Rwv/Jv63/YSv/SEfYNFFFeYflwUUUUAFFFFABRRRQAUUUUAFFFFABWJ4m/5FrV/+vC6/wDRT1t1ieJv+Ra1f/rwuv8A0U9aUvjXqZ1v4cvRnw9+xP8A8fvir/rw0n/0KevvyvgP9if/AI/fFX/XhpP/AKFPX35X0vGf/I4q/wDbv/pKPgfCv/kmMN/2/wD+nJhRRRXyx+hBRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAV45+0H/AMkW8Xf9g8/+hpXsdeOftB/8kW8Xf9g8/wDoaV6GVf77R/xR/NHj8Q/8irFf9e5/+ks8w/Y9/wCSfaz/ANhyb/0RFX1lXyb+x7/yT7Wf+w5N/wCiIq+sq7eJf+RpW9TyuBP+Sfwn+H9WFFFFeEfWhRRRQB//1/2oooooAKKK82vfjL8IdNu5tP1Dxv4dtrq2kaKaGXVLZJI5FOGVlMmQwPBB5FAH5+/tEf8AJ2uh/wDXz4P/APS41+prffbn+I/zr8e/j78S/h1qH7U2i6rYeKdGubJLjwmWuYr+B4VEN6Wky4faNi8tzwOTX6Yv8cvgruY/8J/4Zxkn/kL2vr/10r1Me06NK3b9Efp/HtWEsqyhRabVLXy0geo0VXtLu1v7WK9sZo7i3nRZIpoXDxyIwyGVlJDAjoQcVYryz8wCiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigD59/ao/5ID4w/wCvWL/0clcB+xJ/ySC8/wCxh1L/ANCWu/8A2qP+SA+MP+vWL/0clcB+xJ/yR+8/7GHUv/Qlr04/7g/8X6H6jhf+Tf1v+wlf+kI+waKKK8w/LgooooAKKKKACiiigAooooAKKKKACsTxN/yLWr/9eF1/6KetusTxN/yLWr/9eF1/6KetKXxr1M638OXoz4e/Yn/4/fFX/XhpP/oU9fflfAf7E/8Ax++Kv+vDSf8A0KevvyvpeM/+RxV/7d/9JR8D4V/8kxhv+3//AE5MKKKK+WP0IKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACvHP2g/wDki3i7/sHn/wBDSvY68c/aD/5It4u/7B5/9DSvQyr/AH2j/ij+aPH4h/5FWK/69z/9JZ5h+x7/AMk+1n/sOTf+iIq+sq+Tf2Pf+Sfaz/2HJv8A0RFX1lXbxL/yNK3qeVwJ/wAk/hP8P6sKKKK8I+tCiiigD//Q/aiiiigBroJEaNs4YFTg4OCMHmvx71b/AIJBeBL7VLu8s/iFqcEE80kkUc+l2lzKiOxIV5m+aVhnl2+Zup5r9hqKAP5rPid+wF4f8D/GjT/hpb+L7m6gu5dEjN02nW8bL/alwYW/dqdp2AZH97oa+vz/AMEePBqscfEe9yCcZ0Syxx7V6R+1qj+Gv2hND8WXHy2z2uh34Y8YXSr4tOeeMKpU/jX6kpKk6CeM7kkAdSO6uMj9DXoYujCNGnOC3X+R+g8X5RhMNlWWYvCwt7WD5ndu8ly33fRt7HnPwf8AhnpXwb+GXhz4YaJd3N9ZeHLJLKK5uyDNKFJJZtvAyScKOFHAr0miivPPz4KKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAPn39qj/kgPjD/AK9Iv/RyVwH7En/JH7z/ALGHUv8A0Ja6H9sHWE0z4HanZ5Bk1a90+wRe7B51eT8o1Y/hVP8AY00yWw+CVvdSgr/aWrapdpkdYzcMisPZguRXprTL3f8Am/Q/UqK5PD6o5faxKt52h/wD6sooorzD8tCiiigAooooAKKKKACiiigAooooAKxPE3/It6v/ANeF1/6Katuq15bLe2k9m/3biKSI/SRSv9aqDtJNkVIuUHFHwZ+xOR9t8Uj/AKcNJ/8AQp6+/K/Oj9ky9bw38Tta8H32FluLGa1GeMzaXOwIx6lWJH0r9F6+p40i/wC1Zz6SUWvSyX6H574VTT4co0vtQc012fM5fk0FFFFfKH6KFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABXjf7QhA+C3i3P/AD4Ef+PpXslfN/7VWtR6Z8I7vTi22XWby1sox3ZQ/myY+iIa9PJabqZhQiv5o/meBxVXjRybF1JdKc/xi0vvehzv7Hv/ACT3Wf8AsOTf+iIq+sq+b/2VdJfTvhHa3silG1W+vL0Z7xl/LQ/iqZFfSFb8RTUszrtfzNfdocvBFKVPIMJGW/In9+v6hRRRXin1IUUUUAf/0f2oooooAKKKKAPiX9tzwL/bXgbTfHEEXmN4duGgvQFyTp1/hHJ9EjlEbseygnpXp37L3xHT4gfCywtryfzNZ8Oquk6kG++zQqPJmI9JotrccA5XqDXvmraVp2uaXd6Lq9ul1Y38EltcwSDKSxSqVdSPQg1+UV9a+NP2RPi+l3ZrJf6Leh1gLsVj1XTAxYwOxO1bu2zkFj1+bIWR9vrYa2Iw/wBX+0tV/l/X5I/XOG40+I+H5cOuSWJotzo3+0ndyh+LfzT2iz9baK43wJ4+8K/Ejw9B4m8I3q3lnL8rr92aCUfeimjPzRyKeCpA9q7KvLlFxfLJWZ+U4jDVcPVlQrxcZRdmmrNPs0FFFFSYhRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFfMH7Qn7Q2mfCzTZfDvhyWK88X3cf7mHh49PR+PtFzjoRz5cWQzsOyhmGtGjOrNQgtT1MmybF5pi44LBQ5py+5Lq2+iXVnzP+2L42ufGnj3Rvhb4WxdTaNIA6L8wbWb8CKGI7ef3UT/AD8cCQHsa/QrwF4TtPAvgvRPB1jzDo1jBZg/3mjUBmPuzZNfDv7JHwZvtQ1FfjZ4yEkvmGWXRVuPmluZrjPm6jJkDltzCEgAMGZxwy4/Q+u7HzjFRw0HpHf1/r87H3XiBjsNhaOH4awMuaGHT53/ADVH8X3O/o249AooorzD8xCiiigAooooAKKKKACiiigAooooAKKKKAPze+Nmnah8IvjpZ/EHSYiba/uE1aBQMB5VHlXkA93Q5H+/noK/QrQdc0zxNotl4g0aYT2OoQpcQSL3RxnB9COhHYiuL+K/w30/4oeEZ/D904t7uM/aLC6Iybe5UEKSO6OMq47qTXxP8IPiprHwO8SXvw++IMM0Gk/aD56EF20+dus0Y/jt5R8zbev3xzuA+2lS/tnLoOlrXoqzXWUelvT8790flEK/+q2d1ViNMHiZcyl0hUe6fZP8rdIya/SOiqen6hYatZQalpdxFd2lygkhngcSRyIwyGVhwQauV8U007M/VoyUkpRd0wooopDCiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKPpQAV+cP7RPim6+KPxN034eeE2+0R6bN9giZDuSTUbg7Zn46i3QYJ9n9K9o+P37QVn4WtLrwb4Kukm1yVTHd3kZDR6chHzANyGuCOg6R/ebnCml+zH8Gbjw7bj4i+KrdotUvI2XTraUHzLa2l5aaQNyJpvQ/MqdfmZq+1ybDf2Zh3m2KVpWtTi9231t2/TXtf8AKuKsc8/xkeGsud4Jp1pLaMU/hT/mv+KS6St9UeG9Cs/C/h/TfDmnDFtplrFaxe6xKFz+J5/GtuiivjJzlOTlJ3bP1KnTjTgqcFZJWS7JBRRRUlhRRRQB/9L9qKKKKACiiigArjfHngHwt8SfDs/hfxbZi7s5iHQg7JYJV+7LDIPmjkXsw+h4rsqKqMnF80XZm2HxNXD1Y16EnGcXdNOzTXVM/KnxN8GPjf8As769N4u+HN5d6jpq9b/TofNkMCklY9QsRnzABxvjUjklfKya9K8EftzxvCtp4+8OtLNEfLmvdDlR49y9TJbTMrxt/sK7ketfoZ0ORwa8z8W/Br4WeOpftPirwvpl/cjpctAsdwPpLHtfP416X1+nVVsTC77rR/1+HkfqC49yzNaap8T4NVJrT2kPdn80rJ/eorpE8jg/bK+BkqBpr3Vrcn+GTR7wkf8AfEbCp/8AhsX4Df8AQW1L/wAE19/8Zpsv7HPwJlkZ103U4gxJ2x6veKoz2A8w4HtUX/DGvwL/AOfLV/8Awc3f/wAXT/4Tv7/4Ef8AGu3r/tS/8Fk//DYvwG/6C2pf+Ca+/wDjNH/DYvwG/wCgtqX/AIJr7/4zUH/DGvwL/wCfLV//AAc3f/xdH/DGvwL/AOfLV/8Awc3f/wAXR/wnf3/wC3h33xX/AJS/yJ/+GxfgN/0FtS/8E19/8Zo/4bF+A3/QW1L/AME19/8AGah/4Y1+Bf8Az5av/wCDm7/+LpP+GNfgX/z5av8A+Dm7/wDi6P8AhO/v/gFvDvviv/KX+RP/AMNi/Ab/AKC2pf8Agmvv/jNH/DYvwG/6C2pf+Ca+/wDjNQ/8Ma/Av/ny1f8A8HN3/wDF0n/DGvwL/wCfLV//AAc3f/xdH/Cd/f8AwC3h33xX/lL/ACJ/+GxfgN/0FtS/8E19/wDGaP8AhsX4Df8AQW1L/wAE19/8ZqH/AIY1+Bf/AD5av/4Obv8A+LpP+GNfgX/z5av/AODm7/8Ai6P+E7+/+AW8O++K/wDKX+RP/wANi/Ab/oLal/4Jr7/4zR/w2L8Bv+gtqX/gmvv/AIzUP/DGvwL/AOfLV/8Awc3f/wAXR/wxr8C/+fLV/wDwc3f/AMXR/wAJ39/8At4d98V/5S/yJv8AhsX4Df8AQW1L/wAE19/8Zo/4bF+A3/QW1L/wTX3/AMZqH/hjX4F/8+Wr/wDg5u//AIuj/hjX4F/8+Wr/APg5u/8A4uj/AITv7/4Bbw774r/yl/kTf8Ni/Ab/AKC2pf8Agmvv/jNH/DYvwG/6C2pf+Ca+/wDjNQf8Ma/Av/ny1f8A8HN3/wDF0f8ADGvwL/58tX/8HN3/APF0f8J39/8AALeHffFf+Uv8if8A4bF+A3/QW1L/AME19/8AGayNT/bT+DtrCzaXFreqSgHbFFp0lvuI/wBu58tB+JrQ/wCGNfgX/wA+Wr/+Dm7/APi609N/ZG+A2ny+bJoE1/0+TUNQubqPj1R5Np/KlfLltz/gHN4eQ95LEy8n7Oz+6z/E+UfF37XXxQ+Id43hX4W6PJo8tx+7VbMf2jrLhxjoqmG25/jIdccll613Pwa/ZCvLi8Txb8ah5ryP9oGiNMbh5pTg79Rny3mk4G6IMwPR3cYA+6PDnhHwt4Psxp/hXSLLSLYDHlWUCQrj32gE/ia6KlPMFGPJho8q79f6/q5lj/EOGHw0sBw3hlhqct5XvUl/29uvvbXSSGoiRIscahEUBVVRgADgAAcAAdBTqKK8w/MWwooooAKKKKACiiigAooooAKKKKACiiigAooooAK8i+K3wb8M/FXT1F/mx1a3QraalCoMkY67JF48yInqpOR1Ug8167RXRhcVVw9VVqErSXU48fl+HxtCWFxUFKEt0/60fZrVdD8yTb/HX9m28ke3DS6IzlmZVa70mYf3mUYe2c9z8vuWr2nwz+2P4cu4UHirQryzk6NPp7peQE+oB2SD6YOPWvsxlV1KOAysMMCMgj0IPBrynxD8DvhN4ole41XwzZefIcvPbKbWUn1LQlTmvqHn2X43XM8P7/8ANDRv1Wn4t+h+ex4PzrKnbIMZ+76U6qvFekkm0vJJebbOLj/ap+Dbrl7+/jPo2m3Gf0UipP8Ahqb4Mf8AQTvv/Bbdf/EVE37KvwcJytjqC+w1KfH6tTf+GVPg7/z56l/4Mp/8aztw3/09/wDJToT46/6hv/KhP/w1N8GP+gnff+C26/8AiKP+Gpvgx/0E77/wW3X/AMRUH/DKnwd/589S/wDBlP8A40f8MqfB3/nz1L/wZT/40W4b71f/ACUd+Ou2G/8AKhP/AMNTfBj/AKCd9/4Lbr/4ij/hqb4Mf9BO+/8ABbdf/EVB/wAMqfB3/nz1L/wZT/40f8MqfB3/AJ89S/8ABlP/AI0W4b71f/JQvx12w3/lQn/4am+DH/QTvv8AwW3X/wARR/w1N8GP+gnff+C26/8AiKg/4ZU+Dv8Az56l/wCDKf8Axo/4ZU+Dv/PnqX/gyn/xotw33q/+Shfjrthv/KhP/wANTfBj/oJ33/gtuv8A4ij/AIam+DH/AEE77/wW3X/xFQf8MqfB3/nz1L/wZT/40f8ADKnwd/589S/8GU/+NFuG+9X/AMlC/HXbDf8AlQn/AOGpvgx/0E77/wAFt1/8RR/w1N8GP+gnff8Agtuv/iKg/wCGVPg7/wA+epf+DKf/ABo/4ZU+Dv8Az56l/wCDKf8Axotw3/09/wDJQvx12w3/AJUJ/wDhqb4Mf9BO+/8ABbdf/EUf8NTfBj/oJ33/AILbr/4ioP8AhlT4O/8APnqX/gyn/wAaP+GVPg7/AM+epf8Agyn/AMaLcN96v/koX467Yb/yoT/8NTfBj/oJ33/gtuv/AIij/hqb4Mf9BO+/8Ft1/wDEVB/wyp8Hf+fPUv8AwZT/AONH/DKnwd/589S/8GU/+NFuG+9X/wAlC/HXbDf+VCf/AIam+DH/AEE77/wW3X/xFIf2pvgzjjUr4/8AcNuv/iKh/wCGVPg7/wA+epf+DKf/ABo/4ZU+Dv8Az56j/wCDKf8Axotw33q/+Shfjrthv/Khhax+138ObOM/2NYatqsuDhRALVPxaYr+gJrwPXPjp8YvjJcv4b8B2E1jbTfI9vpGZbhlOcie8O1YlI64Kf7x6V9Z6V+zb8GNKlEy+HkvXBz/AKfPLdD/AL5kYr+ley6bpemaParY6TaQWVugwsVvGsSDH+yoArSOa5NhPeweHc5dHN7fJXT/AAfmc8+HeKMyXs80xsadN7xop3a/xNJr/wAmXdHyj8Hf2YLDwxPb+JfiAYdQ1WFhLb2EZ32lq4OQ7kj9/KDyCRsU8gE/NX17RRXz2Y5niMdV9tiZXf4L0R9rkeQ4HKcMsNgYcsevdvu31f5bKyCiiiuA9gKKKKACiiigD//T/aiiiigAooooAKKKKACkJC8kgfU4/nS14H8dtH0zxFN8O9B1u0jv9Ov/ABrZxXVpMu+GeMWd44WROjLuUNg8ZAPagD3zsD2PIPrUc00NvE89w6xRRqWd3YKqqOpLEgADuTXxT4l1yD4AeLPHth8OYQtifDOhapHoRMk9jpusalqUmmQTxWwcSRx3S4Z7eJkWRrclAskjs3N/FjVvjB4h+Hfjrwb4gJfR1j8OT2ut6l4cOnrdrd6klvqGmTWMlyRIFCpKswIV4pvJdGKM7AH3rb3tlec2lxDPhUf91IsnySDKN8pPDDlT0I6VZr42+IfxA8W/DeXxrD4Vt9Ai1fTbz4daNDfNpXlxzf29qC2EzTxQzIzxwxysbeIOBFyMsDUvjL4vfE74c2vjTw3ez2PiDW9O1jwpp+jajDpcqLHF4suRaK1xYwzM1y9lIkkgSOSLzl2Rkod0lAH2JRXyRc+PPj3baJ4jtLHSb3Um0/UtENlrT+HjbXb6TeyFdS26U12pubqxEbOGQojpKmyKV42R/Y/hV4ovvGngqW+Gtx6lqMN7f2DXb6LcaPLbywSsqR3Wm3b+ak8S7fMGVSbh0Co6gAHqnI4PFFfGHw+8R/FuH4efDXwtH4hsLnXPGN7qX2jWrrTmP2Ox09Z55xDbG4cT3EzqFiMkoSJGJIfywr6MPxR+K+p+ItI+G9ne6Va6xb+NNX8K61rDaZI8NxbWWjDVoLmytjNhJiJYY5Fd2jV1k64AoA+vqieeGOWOF5EWSXd5aFgGfaMttBOTgcnHTvXyrpPxG+KOtarZfC3+0dLtvEkWva5puoa/Hp7PBJZaNBb3Cvb2LzkJczJdxI4eRo0KSMAQVWs3RfFviLW/iJ4LHjSO0m1rwrqXjfR72506N4ra8+w2cDpcQxyFmiMsLJ5kW9xHLuQO4UMQD7Bor478D/Fv4q3enfCTx34kuNIv9I+LzGNNGtLU20ujveaZc6tYiG8M8n2wRx2xguC0a+Yz+enlIhjPF+E/id8ffHEena94fulu9Rm8DazrV1o/9mvZ6TpHiOJ4BZ6XdzsrPeCTM8SvFICDE02GSRFUA++Kiingn3iGRJPLcxvsYNtcdVbBOGHcHkV438PfifN8TfEtxdeHFjPhS28P6RefaGXM0mp6urXPkZJBT7Na+WZFK5LTLyNpB8l1L4i+M7DxBYeGvB6aLo83ib4s6n4XurttNMv+iQaLcagbho1miEt4z2yqZXO0rgFTigD7DpAytwpBx6EH+VeU/CLxJ4o13TNd0zxlcWt9qvhrxBf6LLfWdubSK7jgEcsMvkF5PLcxTKsgDFS6llCqQo+P/hX8PNc8T+Efh/q3g/wOPCusQa3DqN/43kvbVZbmwtbyZrnMMMks9yNQjHleTOqhUfexV41FAH6L1E80MUkcMkipJKSI0ZgGcqMkKDycDk46DmvmDwn8VvHnjN/DHhG2Nrp/iSO38Rf8Jdctp7tb2s2gynTh5EDzhoxeXrpPb73cPbI+Cchx5N4Mg+I+t+D/ANmy5ufEVpLrmoRX1y+qXOnNM8MM2hyltsbTnzZ9x/1kj7SCcqeKAPv2ivkPwl8VPif4ytPC/gkXumad4lvtW8Xadqmtxae0tu8HhG8Nm8trZSTfu5rxnjcK8kiRDeMv8pM2k/EP4veIdY8N+AP7T0XTNZOseJ9G1zVYdNkuYLkaHHE0FxZW8lwogeYTo0kcjzJGweMF8B6APraivIPgX44174g/DuHW/FC239r2Wq63od7LZI0VtcT6HqVzpzXEcTs7RLP9n80Rl38vdt3NjJ9foAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigArzDVPippeieNtH8Eaxo+s2cviDUZNK0vUJLZP7Pu7qO0mvWVJVlZwPIglOWRQShHpn0/v6V8jeKPFVz8Sfix8MtG0bw9r2n6j4M8aavf6yNR02aC2tbK30bUrGK7F5tNpNFdS3lu0AileRldtyK0UqoAfWK3dq8stvHNG80IDSRK6mRARkblBLLntkCuQ+HfjzSfiR4O0TxlpUclpHrmnwalHZXLxG7giuBlRMsTuAe2QSpPQmviv4W+GLS1ufhNYweF9QtPipo15cP481z+y57Uyxi0nj1R73UXjiivob66aA28avJ5jCOVI1WAmPH+EHhjRbez+B9p4e8IarpfjrRhcTeLtQXR7jT510xtPu47pLu9ljjS4S5u3tzBBvcvKqSKgETMgB+jkN1a3DyR280UrQtslWORXMbf3WCklT7HmnmWJTtZ0BzjBYA5xnpnPTn6V8Kfsw6LbeEPGMHhfw9oMV7pkPhiRbrxNc+F7vw3rsEsN2vkWWsTT/utUvLoyTzyzRBWSRGkcYuEJ9Wh+G+la38WfiX4t1HQUutXis9Nt9B1C6hY+S76VPBIbN3+RHJmaOSSPDENsY4wAAfS3mxAAl1wQWB3DkDqR6j3qOK6tp3kjgmjlaIhZFjdXKEjIDBSSpI5Ga/PzT9cvfFegeCdK0rw/4it7zQvhV4r0vUkvNHu7VbPU47PT7cWTPLGqSXDyRv5QiLrKqlo2ZSCfRvAHw5svA3xh+G134S8O/2JYX3wz1S21+WztWghnvbWbRjYrfOAA9zGHuvJ80mQAy7eN9AH17LdWsMsUE00cck5IiR3VXkI6hFJBYj2BoN1bLcLZmaMXDKXWEuokKDqwTO4geuMV8d+NNI8BR/FDxjdfGHwdf+KJ9QbRT4Smt9HudTkW0jSFTb2VxBGVsriPUVaaRjJCVBSZpNiZTh/E3hnTpNT8Q2N/4U1G7+Lt545ivNB16HTJ3lisRdJJZXEerKqwwWFpp4kiuITMin97CUkecCUA++/tdoJ1tfPi85wxWLzF8xgv3iFzuOO+Bx3rkZPHOm3EWtr4dhl1298P6hDpt9ZWjxRypPKsEjYad442WOG4WRiG6AquXG2vjzxL8N1j0TUPHOk+GHHjGD402N3a6pFZOdTj02XXrS3uZIpQplW0fT2lWQr+6MBYtxk0us+CNH0i3+N+keHfCb2Oq3/jTRdShez0iSIXmnzy6VI0kM8cQSdPtK3EkqoxKPvdwNxYgH3i91aJOts88QlkLKkZkUO5Xk7VJ3HA5IA4FSebHgkuuAdpORwfTr19utfAfin4XW2p3Xxv8AGVz4VebxFbfEzRdT8OaotnJ9vEdtZaEvn2UyjzPKDCeOQxHY6iRHyu8VpeO9UvIbXxj8Pf7C16TV7j4raBq1uLfS7qe2uNLm1bS7g3iXUaNB5ESJIJSXDRlG3LjBIB9zG6tluFs2mjFwyl1hLqJCg6sEzuI98YqTzY8Z3rjO3O4Y3en19utfAHiTwzYSap4hsdR8LahefF278fC90LXoNMnaSLT/ALcstjcR6sqLDDYWum+ZFcQmZAR5kBSR5gJdbXb26Gn6l8Nxo2ttrUfxg03WNsemXbWraZNrtrqC3i3Yj+ztbiA/ORIdrBkYAggAH3E17arJLAJo2ngjMrwh1MqoATkpncAccEjFcn8N/Hel/E7wD4c+IeiwzW1j4m0q01e2t7nZ9oigvIxIiyhGZQwB5wSM9Ca+PPCHh+28N/HfHh7QRr8+q+JfEs2pahrHhi6sPEGiR30E0z3P9uk/Zr6wZ1gtLWHbuMDRBWP2dgOJ8Nf8Jr8NfgP8JvFPhfT7+08TXPg23+FeoWVxYSRvba3OUh025uo5RFK0Gn3yzjeh2NFO0iFl2tQB+lnmx7Q+9dp6NuGDnjrnHWgyRgkM6jGM5IGM9M+me3rXzX8X/B3/AAh/wB0zwx4R0281Kz8H3vhOYWtkj3N22n6HqdlNcyJGu6WeRIIZJSihpJSCqhmIB8i+KAm+I1r8Rb7QtH1e80vVNb+G1vH5um3dub21g1WFrx445I0keCOF2859oCqHLYUE0AfaWt+LPD3h7RpPEGq30MenxSxQNOjCRfNmlWFFGwnLeY4Ugcjv0Ncz4x+I0fgm31TUNU0a+l07Tm0iOO9gktmiu5dWuxabIl83zAbZmR5S6KrI48ouwYL8k/E/4WeHdLk+LOj6V4JiPhhH8D+I4tPsdHM9rLfQ3sg1K5tbaKJhLci0gjEywq0jKFG0lgDb+JuiSjWvifrHhvQb/wDsf/hHfhpPZLZaZcYe00vV725mit4Ej3u9tB8zW8aGVAyjZllBAPvRhtYqexI/KkpqTR3Ki5iyY5h5iblKna/zDIIBBwehGR3p1ABRRRQAUUUUAFFFFAH/1P2oooooAKKKKACiiigArjPG/wAP/DHxF0210nxRDcyxWV7DqFs1neXNhcQ3UG4JIk9rJFKpAZgQGwykgggkV2deB/tBaH471vw9oQ8H2+q6nptnrUdz4l0jQdRGk6vqekrbzqILS8aW3EbLdtBK6+fAZYo2TfhirAHVaX8FPhho/hzX/CcWgx3On+KQy66NRlmv7jUw0Yi/0u4unknm2xgIm5zsUALgAVBp3wR+HWnaXqelGyvL+PWBYi+n1PU73UbuWPTXElpGbm5mkmEUMg3KgYLuLMQS7E+UfDz4geHvDkGmeGPB8viG9W+8cHw5qemeL7m5k1Pw7NLpst+YUe4Es00bCNJYt88qFJiY5PLVEHOfEj4h+Mtc+I/gvR/Dlva2t1onxYuvD0Qlu50tr63HhK7vjJdrGm4hHmJEK7gzRIdyk/KAfS+ufDPwR4kudSu9a01bmXV7nSLu8YyyL5k+gzi5sGwrAL5EyhwBgMR82RxSeIfhh4F8V/29/wAJDpUd6fEsVjDqReSRWkGmM0lo6FWUwy28jF4pIiro+GDBlBHg3/C/PHi3+n+BToOkHxc3j6TwJfSm7uF0pSNEk12K/gIjM7K9qIw0DAFZSyCRlUSNJ/wuz4paNa6frni3w/oMWlw+NbbwPqq2F7cS3Tz316lhb39p5kSxiLzZoTJbyHeEZyJGKKsgB62vwU+H6aLdaJJbahIL29g1Ga9l1a+fU2u7YBYJRfNP9pUxKNqBXChSVxgkV2XhPwjoXgnSv7H8PRSxwtPLdSyXNxLd3E9xMd0ks087vLLIxxlnYnAAHAAr5jufG/jjz/D8XgG10zS4r34qa5oWpx3l3dzi6W1W8bfu2sVWdoQ7xrhYyFVPlBz6t8a9c8d6FpPhmXwPNYwT3vifR7G7a8aVQ0NzMFMa+UrfJI3yyZGdhO35qAL7/A34bNoEPhtbC7js7PVZNbsGj1K8W506+lL75LG5E3nWikSSJ5cLLH5bsm3YxU8prX7PPha4fwJpWgQf2donhPVdV1a5CXd2mozXOpWc8DXK3qSfaGu3mnMks8khd8tkkmsTxh8bPHHh+28eeL7DQ9Kl8KfDGUQ64bi6mTUL5ba1gvNRmsVSN40S1gm/crN81xIjIfKXbI1r49ar4lt/FPwVPg6SL7TqPjiWIR3U8sFrNC3h/VpcXHlBmdE2iUR4+Z0UZU/MADuF+A/wvi8LaV4RttLmtrXRL+XVbC6gv7qLU4dRuN/n3f8AaCyi7a4uPNkE8jSFpg7CTcGIOnpnwd+HmjW/hq30nS5LUeEbu6v9Mkiu7kSi5vg4u5LmTzN921yZGef7QZPNkPmPlwCPMPD/AMdPEUviGw8DeI9HsY9ch8dy+CdXmsriVrPH9hya5b3lpvTzD5sHlI8MuPLdnAdwqs9zVfibHfeM/Cmlaho8UrD4nXHha0mW7miMRt9Hur0XRRAFlchWjMMmYwG353KooA7HT/gj4H8O3M+qeFbWWw1GO2vYtLNxdXN5Y6VJeoyu9nYzStb24Yt8yxKmUygwpIrxD4W/s+eNvhxrvhe80WaPQrHw3bSQXdnD4o1vV7HWYfskkMVsNPv/APRdNhW4Mc4eBWeJY/JjHlk1L8NfiF4ttfDXgfwB4QgtL3XvEl54v1CS91ua4a1tNN0fVpEmY+XvlllaS5giii3oqoWbOIwjeifs+6rrUnwy1vU/EKJFqEHivxkZo5b1p4EMGr3YVVuZQCsACjZlQI48DaAuKAOn+CPwwj+FPgh9BkW1F/qeq6prupmyDfZhe6vdSXUkcJcK5ih3iKMsAxRASBnFdGfhp4JOpWurnTQbuy16fxNBKZZMpq1zbPZy3AG7GWt5HTYRsAOQM8180p+1BrNh4W8Y+IL2Hwx4hbw/4Du/HFpJ4a1C5ns5DZhjJZSTTQp5iMdpiu4gVlUsTGmFD9nqfxv8YeBF8QS/EXQ9OmWx8GS+MtPj0G5lkklS2dYp7GU3KRr5gkli8uYFVcO2UTZlgD6D0bw5o/h+XU59It/IfV7+TUr1tzN5t3KiI7/MTtysajAwOOnWjw54c0bwnoln4b0C3+y6dYJ5VvDuZ9iFixG5yWPLE8mvn7wze/Eyb46aBD8RYdKtZ5PA+qXCxaHeXU1p5r39h5kckdwke97fKrHPtzIGchIuVPVTfEbxvL458Z6Jp+jaUfD/AIMtYJri+mvZhe3Es9g92Io7ZYDGArBAXMvKscLleQC38OPhg/hbxr8QfiJrFvpkWs+OtTtpJP7OWQgWGnWyW1ss0koBkmfDyyFVVVLhBu2b33/DHwk8BeDrXQ7LQLCWG38M3F3caPHLeXFwti17G0UqQiaR9sIjYpHCP3cS4WNVAAHAa38bNT0jTfDuoDSLaQa14E1vxfMDcOPJn0m2s50tlGz5o5DdMGckMoQYU7jtd4G+Kfjq+8e6H4J8faRpNm3ivwnL4q0yXR7qe5Fv9iktIru0uTPHFvKm9haGZAocbwyJtBcA7C/+Cvw31DS4dKbS5LZLXWr3xFa3Fld3FneWuqalJLJd3EFzDIk0TTtNIJFVgrI7IQVJFa3h34X+BfCaaEnh/S1tD4cF9/Z7+bLJIr6mQ15LK7szTzXLjfLLKXkdyWLFiSe/ooA57wx4V0Hwbpsmj+G7UWdnLe32ovEGZ83WpXEl3cyZck5knldyOgJwMDAroaKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKd8+3+IoPrtH9KbxkZ6ZGfpXxP4E1PwPN431C/8AiR4j1G0+JsfjnVrez0mLVbtLh7CKaePTLePTEleB7CbTVSZiIvKZ90zlZQxUA+2jv2jdu29s5x+HahvMx8+7H+1nH61+avwo+I+jap4+8I6zoV6NMj8Y6L4rN9YTeJtQ1nUftyotzbW2qQXEYtrTUoSs2IEdZIRHJDGrxRMydB8IbHXNE8K/Bjxf4f1HWbzxF4w+Gl9Lqi3+pXd/HqOoRaXa3dtNLFcyMgnS4JCOAh2uUORjAB+hTb8Dfux23Zx+Ga5vxb4q0jwRocniDxE80NlFNa27NHE8r+ZeTx20QCKCcGWVQT0UZY8AmvgvTde8GW3wX1fVfg/4y1zV/iVL4EubvVUTUL7Ub5LsNGL25vbNknWz1GCfzkhiEEcilXiigZYvLXE1W+8OW3gPxdqPhLxXpV/4RsrHwpe3sGm6zqOv2ts8Wr20o1CW/vtyxO9orNcrvLFVE8oGSzAH6H6N4u0fxBrXiHw/ps00l74WvodO1JHjdFjuJ7aK7RY2YYkHkzIdyZAJK/eBA6MhwdpDZ645yfwr87fiN4nl0+6+OsGlat/Z8c/jrwde6vPFdT2UkXha907TILm8NxbI9xBZSLHLE95CB5arKRJGY2dKuoahaW1jbQS+LbYfCy58dafbeIW0TVtQnsNLtRptzI0B1eZo5EsrjUVsxMsMhjjZmichJHAAPvnRPFeka/qev6NpcspufDeoJpmpI8TxKlzJbQ3aqpYASL5NxGdy5UElc7lIHR5ONuTj0zx+VfL/AOzePD66v8WB4TvZtQ0X/hNIfsM8k8t2jRf2JpmRDcSlnlhV9wRtzAY2qcKAPp+gBcn1Pp17Ubm6AnjpzSUUALk+p46UbmA2gkD0ycUlFAC5ONuTj0zx+VG5gMAnHpnikooAXJxtycDtnj8q5rWPCWg69rejeINXgkubzw/LJPp26eUQwzyoYzKYQwieRUZlR3VmQM20jcc9JRQAdOnGKUsxOSST9TSUUAKCRjBIx09qNzZzk569e9JRQAUUUUAFFFFABRRRQAUUUUAf/9X9qKKKKACiiigAooooAK4vxr4Lh8Z21ig1fVtDvNMuvtlnfaPdfZ5o5djRkOjLJBcRsjsDFPHImcMAGVWHaUUAeESfs/8AhuXRrq2fW9dOvXeuweJm8TG5iOrJq9tAtpFPHmL7KiJbL5HkCDyTGWDISzEsf9nvwu4t7r+2tfGq2viWbxdFqv2xGvBq0+ntprSZMRiMf2diBD5fl8427eK96ooA8W0T4FeEtIl0XUbi71LVNX0jxHN4sm1W9mjN3qOrz2EumNNd+XHHGVFpKY0jjSNEVECgBQK3dU+FHhjV9Kk0e8e7NvL4otfFp2ShW/tCzu4r2IAgf6nzYVyndcjPNemUUAeSXXwc8PTWCWdlqGqadNB4oufF0F5azR/aIdRvHkaYL5kciGF1mkj2MpwrcEMAw6nxx4KsfHeiRaNeXt7prW17aajbXmnyJHcwXVlIJYnQyJIhww5V0ZWGQRXZUUAeHeIvgH4V8SahrU1zqWr22leKp4LnxNoltPEmna5NbxRQbrtTE0wEsUMUc6xSRpOiBJAylgev8efDjSPH7+H7m+vL/Tbzwvqbavpd1psqwywXZtZ7Pd86OjqIbiQbGUqc8gjivQqKAPBm/Z98MHQ4rRdY1pfEEWvjxT/wlImgOsNrX2c2ZumJgNqc2h+zeV5HlCH5AgAGNfSPgl4W0iXSro32rX11pPiq58Yx3N7dmaWXVLu2mtZDKdoDQ+XO+2IAKnyhAFVVHsVFAHgafs9eGrK10hdC13X9I1Hw/qeq6jpGq21zC95Zxa3KZr6xUzQyRy2Mznd5MySbWWNlIaOMr2GhfCfwpofgHVfhv/peoaNrcmryX4vpzNPOdcllmvN0uAfneZyP7oOBwBXplFAHz1d/s5+HtY06/wBN8UeJvEuupqXhi/8AB85vbuBR/Y9+qq0Sxw28UayR7crMF8wn77MAAO58TfCTwZ4ymnbxLBNfW934Zu/CVzavKVgm0y+aN5gwXDeYTEuHBBXnFemUUAeTeEfhLa+GPEdt4tv/ABHr/iTVbPSptFgn1i5ikEdjLNHNsEcEMKFw0SgykGRwPnZjgjrLPwVotnq3iXVwJZn8V/Z/7QilfMRW3tvsqqgABUNF97k5PPFdbRQB89aT+zj4Y060t7G+8Q+I9Yt7Dw9qnhXTk1C9jf7Fo+qpFHJBGY4YzI8aQRrHPMZJsL87tzn0qD4d+H7fxR4f8Xobg6h4a0K78O2RMn7s2V69o8vmLj5pC1nFhuMfNxzXd0UAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABUflReb5+xfNC7RJtG/b6bsZx7ZxUlFAESwQIWKRopZt7EIoLN03Hjk479aeEQbcKBt4XgcfT0/CnUUARpFFGzvGio0h3OyqFLH1YgAk+5pEghjRo440RGJLKqKFJPUkAYJPfNS0UAM8tMk7V5G08DlfTp09ulNEEAh+zCJBDgr5YRQmD1G3G3B+lS0UAMREjG1FVQMDCgAccDoB0FPoooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAP/W/aiiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooA//Z"
    }
   },
   "cell_type": "markdown",
   "metadata": {
    "id": "8ab7a830-78fb-45c8-aa35-d76919f41a9b"
   },
   "source": [
    "# Timeseries Forecasting Sample\n",
    "\n",
    "In this notebook, you'll see examples of how to:\n",
    "\n",
    "- define structured schemas by using `NamedTuple`, also with containers and nested named tuples as members\n",
    "- create a data source as a composite\n",
    "- insert window punctuations into a stream by using a callable\n",
    "- aggregate data within a tumbling punctuation based window\n",
    "- use the [`Stream.flat_map`](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.topology.html#streamsx.topology.topology.Stream.flat_map) function to create additional tuples\n",
    "- invoke an SPL operator\n",
    "\n",
    "Familiarity with [Python](https://www.python.org/about/gettingstarted/) is highly recommended.\n",
    "\n",
    "If you are new to the  Streams Python API, [start with this tutorial](https://github.com/IBMStreams/sample.starter_notebooks/blob/latest/Streams-RollingAverageSample.ipynb).\n",
    "\n",
    "The notebook is divided into following sections, from which _Create the application_ is the most extensive one.\n",
    "\n",
    " 1. [Setup your Streams instance](#setup)\n",
    " 1. [Create the application](#create)\n",
    " 1. [Submit the application](#launch)\n",
    " 1. [Connect to the running application to view data](#view)\n",
    " 1. [Stop the application](#cancel)\n",
    "\n",
    "# Overview\n",
    "\n",
    "**About the sample**\n",
    "\n",
    "This sample simulates a call center, which creates an event when a call is terminated. For simplicity, these records contain only following data: \n",
    "\n",
    "- the call center agent\n",
    "- the timestamp when the call has been terminated\n",
    "- the duration of the call\n",
    "\n",
    "The duration of the calls as well as the time between two calls is randomly distributed and follows a log-normal distribution. During the night hours and at weekends, the calls shall happen less frequently, so that the call termination events form an _irregular complex seasonal timeseries_. The timeseries is irregular as the events are not equally spaced in time. The data has two seasons: The inner season is a day, while the outer season is a week.\n",
    "\n",
    "These call events are aggregated for every hour, in such a ways that we get one sampling point for every hour. We are only interested in the hourly call volume, i.e. the sum of the durations of all calls. For hours, in which we see no activity in the call center, the hourly aggregates are missing. To create a _regular timeseries_, i.e. at equally spaced points in time, we create artificial aggregates. Here we use the [`Stream.flat_map`](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.topology.html#streamsx.topology.topology.Stream.flat_map) function with a suitable callable class.\n",
    "\n",
    "Finally, we use the [`BATS`](https://www.ibm.com/support/knowledgecenter/SSCRJU_5.5/com.ibm.streams.toolkits.doc/toolkits/dita/tk$com.ibm.streams.timeseries/op$com.ibm.streams.timeseries.modeling$BATS.html) SPL operator from the timeseries toolkit to forecast the hourly call volumes for the next hours. The outcome of the forecast operator is visualized in a chart with [matplotlib](https://matplotlib.org/).\n",
    "\n",
    "The data generation can be speeded up to make the timestamps increase significantly faster than realtime. The default value in this notebook is 7200, defined in variable `time_compression`. \n",
    "\n",
    "**How it works**\n",
    "   \n",
    "The Python application created in this notebook is submitted to the IBM Streams service for execution.\n",
    "\n",
    "![how-it-works.jpg](attachment:how-it-works.jpg)\n",
    "\n",
    "\n",
    "## Documentation\n",
    "    \n",
    "- [Streams Python development guide](https://ibmstreams.github.io/streamsx.documentation/docs/latest/python/)\n",
    "- [Streams Python API](https://streamsxtopology.readthedocs.io/)\n",
    "\n",
    "\n",
    "\n",
    "# Prerequisites\n",
    "\n",
    "This notebook can be used as-is from within an IBM Cloud Pak for Data project and with a local Streams installation.\n",
    "To use this notebook, you must install `matplotlib` within your Python environment, for example:\n",
    "\n",
    "    pip install matplotlib\n",
    "\n",
    "If you are not running this notebook from within IBM Cloud Pak for Data, [follow these steps to make sure you have installed all the prerequisites](https://ibmstreams.github.io/streamsx.documentation/docs/python/1.6/python-appapi-devguide-2/).\n",
    "\n",
    "<a name=\"setup\"></a>\n",
    "\n",
    "\n",
    "# 1. Set up a connection to the Streams instance\n",
    "\n",
    "\n",
    "To submit the application for execution, you have to connect to the Streams instance. The information required to connect to the instance depends on the target installation of Streams. \n",
    "\n",
    "Choose the option that matches your development environment.\n",
    "\n",
    "- **Option 1**: [I'm running the notebook from an IBM Cloud for Data project](#cpd)\n",
    "- **Option 2**: [I'm using IBM Watson Studio, Jupyter Notebooks, or any other development environment](#notcpd)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "16da67f8-5688-439a-a4a4-d1dbf267e697"
   },
   "outputs": [],
   "source": [
    "try:\n",
    "    import matplotlib.pyplot as plt\n",
    "except ImportError:\n",
    "    print(\"Please install matplotlib into your Python environment\")\n",
    "    raise"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "bdb29bf6-cb2e-4b7a-b98a-72d1ad7d4eb8"
   },
   "source": [
    "<a name=\"cpd\"></a>\n",
    "### Option 1: Connect to a Streams instance from an IBM Cloud Pak for Data  project\n",
    "\n",
    "In order to submit a Streams application you need to provide the name of the Streams instance.\n",
    "\n",
    "1. From the navigation menu, click **Services > Instances**.\n",
    "2. Update the value of `streams_instance_name` in the cell below according to your Streams instance name\n",
    "3. Run the cell and skip to section 1.2\n",
    "\n",
    "The cell below defines a function called `submit_topology` that will be used later on to submit the `Topology` once it is defined.\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "8c58cbde-707d-43df-88cb-7f8d4dab7692"
   },
   "outputs": [],
   "source": [
    "from streamsx.topology import context\n",
    "\n",
    "streams_instance_name = ## Change this to Streams instance\n",
    "try:\n",
    "    from icpd_core import icpd_util\n",
    "    cfg = icpd_util.get_service_instance_details(name=streams_instance_name, instance_type=\"streams\")\n",
    "except TypeError:\n",
    "    cfg = icpd_util.get_service_instance_details(name=streams_instance_name)\n",
    "except ImportError:\n",
    "    # not within a Cloud Pak for Data project\n",
    "    cfg = dict()\n",
    "\n",
    "def submit_topology(topo, context_type=context.ContextTypes.DISTRIBUTED):\n",
    "    global cfg\n",
    "    # Disable SSL certificate verification if necessary\n",
    "    cfg[context.ConfigParams.SSL_VERIFY] = False\n",
    "    return context.submit (context_type, topo, config=cfg)\n",
    "\n",
    "print(\"Setup complete, continue to section 1.2\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "451bffa7-1393-465e-8817-d57994f28ef8"
   },
   "source": [
    "<a name=\"notcpd\"></a>\n",
    "### Option 2: Connect to a Streams instance from IBM Watson Studio and other environments\n",
    "\n",
    "*Skip this section if you are running the notebook from a Cloud Pak for Data project.*\n",
    "\n",
    "The code for each scenario is available in the development guide.  \n",
    "Each snippet will define a function called `submit_topology` that will be used later on to submit the `Topology` once it is defined.\n",
    "\n",
    "- Choose the tab that best matches your environment. \n",
    "- Copy the code under the heading **Copy this code snippet**.\n",
    "- Paste it in the cell below.\n",
    "\n",
    "    [Connection instructions from the development guide](https://ibmstreams.github.io/streamsx.documentation/docs/python/1.6/python-appapi-devguide-2/#connect)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "eb546fcd-8059-4869-9b3c-7450947f5e6a"
   },
   "outputs": [],
   "source": [
    "## Paste snippet from development guide\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "ddb8082b-9129-4464-814b-c0232fc48ae8"
   },
   "source": [
    "### 1.2 Import  the `streamsx` package and verify the package version\n",
    "We need at least version 1.16 of the streamsx package."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "a5b564db-6ee3-4852-9db4-ae176c08bbb2"
   },
   "outputs": [],
   "source": [
    "import streamsx.topology.context\n",
    "version = streamsx.topology.context.__version__\n",
    "v = version.split('.')\n",
    "vv = int(v[0] + v[1])\n",
    "\n",
    "print(\"INFO: streamsx package version: \" + version)\n",
    "\n",
    "if vv < 116:\n",
    "    raise ValueError (\"streamsx version too low. Please upgrade streamsx to at least version 1.16\")\n",
    "\n",
    "#For more details uncomment line below.\n",
    "#!pip show streamsx"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "1ee8cf71-26f8-40c2-bfb4-230f43ca5a40"
   },
   "source": [
    "<a name=\"create\"></a>\n",
    "# 2. Create the application\n",
    "Creating an application means to define schemas, callables, and finally a `Topology`, which is basically a data structure that defines the streaming application. We are going to define the schemas and callables when they are needed to define a stream or transform.\n",
    "\n",
    "## 2.1 Define the schema for the data source\n",
    "The data source creates a stream of tuples with structured data with three attributes: \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "857446b9-5451-42c1-8609-f8af8be860d1"
   },
   "outputs": [],
   "source": [
    "import typing\n",
    "\n",
    "# schema of the source stream\n",
    "class CallTerminationEventSchema(typing.NamedTuple):\n",
    "    callcenter_agent: str\n",
    "    call_termination_time: float  # seconds since Epoch\n",
    "    call_duration_secs: float\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "f09ca8a2-9e7a-4118-ac95-b545b7cec7d3"
   },
   "source": [
    "## 2.2 Define a callable to punctuate the source stream\n",
    "Punctuations, or window markers, are signals within the stream. They flow together with the data tuples through the application. The meaning of these signals fully depends on the application design. In this application, a window marker separates the data according the full hour of their timestamps. Whenever the hour of the timestamp of two sequential tuples changes, a window marker is to be inserted.\n",
    "\n",
    "The callable below returns `True`, when the hour of the call termination timestamp is higher than the previous one, `False` otherwise. A callable for the [`Stream.punctor`](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.topology.html#streamsx.topology.topology.Stream.punctor) function must return a boolean value.\n",
    "\n",
    "We will use this callable class when we define the source of our data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "95227709-143f-4104-9f86-4e3321c65cf7"
   },
   "outputs": [],
   "source": [
    "# Callable to punctuate the source stream\n",
    "class TimeStampHourChanged:\n",
    "    def __init__(self):\n",
    "        self.previous = None\n",
    "\n",
    "    def __call__(self, call_termination: CallTerminationEventSchema):\n",
    "        if self.previous is None:\n",
    "            self.previous = call_termination\n",
    "            return False\n",
    "        prev_hour = self.previous.call_termination_time//3600\n",
    "        this_hour = call_termination.call_termination_time//3600\n",
    "        self.previous = call_termination\n",
    "        if this_hour > prev_hour:\n",
    "            return True\n",
    "        return False\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "f628d34f-0999-413c-bd70-e8d9d8b635ad"
   },
   "source": [
    "## 2.3 Define a data source as a composite\n",
    "[Composite transforms](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.composite.html#module-streamsx.topology.composite) encapsulate multiple stream transforms in a single class for hiding compexity and re-use. In this notebook we will define a [Source](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.composite.html#streamsx.topology.composite.Source) that emits a punctuated stream. The subgraph in our composite consists of the actual source followed by a punctuation inserter pseudo transform.\n",
    "\n",
    "As we create a _source composite_, we derive our class from `streamsx.topology.composite.Source`. The transforms of the source must be implemented in the overwritten `populate` method. This method must return a `Stream` object:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "a1cbfbeb-69ac-47f3-9d7c-78ca6eb1869b"
   },
   "outputs": [],
   "source": [
    "from streamsx.topology.composite import Source\n",
    "from datetime import datetime\n",
    "import time\n",
    "import random\n",
    "\n",
    "class CallsSeasonal(Source):\n",
    "    def __init__(self, speedup = 3600.0, n_tuples_no_delay = 0):\n",
    "        self.starttime = datetime.now().timestamp()  # seconds since Epoch as float\n",
    "        self.speedup = speedup                       # time scaling compared with a wall clock\n",
    "        self.n_no_delay = n_tuples_no_delay\n",
    "        \n",
    "    # generator method\n",
    "    def create_call_terminations(self) -> typing.Iterable[CallTerminationEventSchema]:\n",
    "        next_timestamp = self.starttime;\n",
    "        delta_t = 0.0\n",
    "        i = 0\n",
    "        while True:\n",
    "            if i > self.n_no_delay:\n",
    "                time.sleep(delta_t/self.speedup)\n",
    "            i += 1\n",
    "            \n",
    "            # we assume that call durations are distributed according a log-normal distribution\n",
    "            # and do not depend on week day and hour of day\n",
    "            duration = random.lognormvariate(0, 0.7) * 120\n",
    "            ct = CallTerminationEventSchema(\"Superagent\", next_timestamp, duration)\n",
    "\n",
    "            # distance to next call termination\n",
    "            # decompose the timestamp into hour and weekday, so that we can model seasonal data\n",
    "            dt = datetime.fromtimestamp(next_timestamp)\n",
    "            weekday = dt.isoweekday()  # Monday = 1, Sunday = 7\n",
    "            hour = dt.hour             # 0..23\n",
    "\n",
    "            # the time difference between the call terminations shall depend on hour and weekday\n",
    "            # and follows also a lognormal distribution, which we scale dependent on hour and day.\n",
    "            if weekday in [6, 7]:  # Sat, Sun\n",
    "                w = 3\n",
    "            else:\n",
    "                w = 1\n",
    "            if hour in [0, 1, 2, 3, 4]:\n",
    "                h = 600\n",
    "            elif hour in [5, 23]:\n",
    "                h = 300\n",
    "            elif hour in [6, 7, 22]:\n",
    "                h = 90\n",
    "            else:\n",
    "                h = 30\n",
    "            delta_t = random.lognormvariate(0, 1) * h * w\n",
    "            next_timestamp += delta_t\n",
    "            yield (ct)\n",
    "\n",
    "    # define the subgraph; here we have a source followed by a punctuation inserter\n",
    "    def populate(self, topology, name, **options) -> streamsx.topology.topology.Stream:\n",
    "        if name is None:\n",
    "            name = 'Calls'\n",
    "        \n",
    "        call_terminations = topology.source(self.create_call_terminations, name=name + '_gen')\n",
    "        # insert a window marker when the hour of the timestamp has changed, before the new hour\n",
    "        # The callable class will keep the previous tuple for detecting the change\n",
    "        return call_terminations.punctor(TimeStampHourChanged(), before=True, replace=False, name=name + '_punct_insert')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "2f10e32a-d009-4ef5-b0ee-6b22af3f0fc2"
   },
   "source": [
    "## 2.4 Define the topology and add the source"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "382367ec-0fd4-4019-be46-ab9e96ea31ef"
   },
   "outputs": [],
   "source": [
    "from streamsx.topology.topology import Topology\n",
    "\n",
    "# time stamps grow 7200 times faster than in realtime. With 7200, one hour elapses in 0.5 seconds\n",
    "time_compression = 7200\n",
    "\n",
    "topo = Topology(name=\"CallcenterWorkloadForecast\", namespace=\"sample\")\n",
    "source_composite = CallsSeasonal(speedup=time_compression, n_tuples_no_delay=14*24*38)\n",
    "terminated_calls = topo.source(source_composite, name=\"CallGenerator\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "7db310a3-a79f-4f1e-bb96-f6338142682b"
   },
   "source": [
    "## 2.5 Aggregate the call termination events for every hour\n",
    "The `terminated_calls` stream contains window markers such that onyl data from one hour are delimited by window markers. For the outcome of the aggregation, we define the new schema `CallAggregateSchema`. The actual aggregation is performed by the callable class `CallDurationSum`, to which an `Iterable` with the tuples within the window is passed. The schema of the transformed stream is derived from the type hint of the `__call__` function. The punctuation based tumbling window is defined by [`Stream.batch('punct')`](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.topology.html#streamsx.topology.topology.Stream.batch)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "88445319-fa96-4626-9b91-d632059723f6"
   },
   "outputs": [],
   "source": [
    "# Schema definition\n",
    "class CallAggregateSchema(typing.NamedTuple):\n",
    "    num_calls: int = 0\n",
    "    n_distinct_callcenter_agents: int = 0\n",
    "    sample_timestamp: int = 0\n",
    "    aggregated_call_duration: float = 0.0\n",
    "    min_call_duration: float = 0\n",
    "    max_call_duration: float = 0\n",
    "\n",
    "# Aggregation Callable\n",
    "class CallDurationSum:\n",
    "    def __call__(self, window_content) -> CallAggregateSchema:\n",
    "        distinct_agents = set()\n",
    "        min = 0.0\n",
    "        max = 0.0\n",
    "        sum = 0.0\n",
    "        for tuple in window_content:\n",
    "            if tuple.call_duration_secs > max:\n",
    "                max = tuple.call_duration_secs\n",
    "            if min == 0.0 or tuple.call_duration_secs < min:\n",
    "                min = tuple.call_duration_secs\n",
    "            sum += tuple.call_duration_secs\n",
    "            distinct_agents.add(tuple.callcenter_agent)\n",
    "        # calculate the sample timestamp. The sample timestamp is the timestamp of \n",
    "        # the last item in the window aligned to the next full hour, for example, if the\n",
    "        # items in the window have timestamps 09:03:45, 09:45:34, the sample timestamp \n",
    "        # of the aggregate will be 10:00:00.\n",
    "        agg_sample_time = 0\n",
    "        if len(window_content) > 0:\n",
    "            div = int(window_content[-1].call_termination_time) // 3600\n",
    "            mod = int(window_content[-1].call_termination_time) % 3600\n",
    "            if (mod == 0):\n",
    "                agg_sample_time = div * 3600\n",
    "            else:\n",
    "                agg_sample_time = (div + 1) * 3600\n",
    "        \n",
    "        return CallAggregateSchema(len(window_content), len(distinct_agents), agg_sample_time, sum, min, max)\n",
    "\n",
    "# define a window that contains the data between two window markers, and aggregate these data\n",
    "hourly_call_volumes = terminated_calls.batch('punct').aggregate(CallDurationSum())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "c4480268-c105-40cc-986a-514a960b3a80"
   },
   "source": [
    "\n",
    "## 2.6 Regularize the time series\n",
    "During night hours, the simulated call center can be idle. For these hours, the aggration does not create tuples as the window marker separates time periods, in which the hour differs by more than one. To have a regular time series with equally spaced timestamps, we create artificial tuples, in which all the data is zero. We can use the [`Stream.flat_map`](https://streamsxtopology.readthedocs.io/en/stable/streamsx.topology.topology.html#streamsx.topology.topology.Stream.flat_map) transform to create addional tuples.\n",
    "\n",
    "The callable returns an `Iterable`, that contains the potentially missing tuples and the current tuple. The `flat_map` transfrom emits a tuple for each item in the Iterable returned by `__call__`. The schema of the transformed stream is derived from the type hint."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "46b547b0-ef58-40f0-a4a7-c8caa635508f"
   },
   "outputs": [],
   "source": [
    "# create missing sample points for hours, in which no call happened. We create a flat_map callable that\n",
    "# creates an Iterable, which contains missing tuples + current tuple\n",
    "class Regularize:\n",
    "    def __init__(self):\n",
    "        self.previous = None\n",
    "\n",
    "    def __call__(self, tupl: CallAggregateSchema) -> typing.Iterable[CallAggregateSchema]:\n",
    "        if self.previous is None:\n",
    "            self.previous = tupl\n",
    "            return [tupl]\n",
    "        \n",
    "        next_sample_time = self.previous.sample_timestamp + 3600\n",
    "        timestamp_diff = tupl.sample_timestamp - next_sample_time\n",
    "        \n",
    "        result = []\n",
    "        while timestamp_diff > 0:\n",
    "            # add artificial equally spaced sample points\n",
    "            result.append(CallAggregateSchema(0, 0, next_sample_time, 0.0, 0.0, 0.0))\n",
    "            timestamp_diff -= 3600\n",
    "            next_sample_time += 3600\n",
    "        # finally add the received tuple\n",
    "        result.append(tupl)\n",
    "        self.previous = tupl\n",
    "        return result\n",
    "\n",
    "hourly_call_volumes_regularized = hourly_call_volumes.flat_map(Regularize())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "1dafb190-41f6-42c6-803f-73d2afca2c86"
   },
   "source": [
    "## 2.7 Invoke an SPL forecasting operator\n",
    "To forecast the hourly call volume for the next hours, we use the [BATS operator](https://www.ibm.com/support/knowledgecenter/SSCRJU_5.5/com.ibm.streams.toolkits.doc/toolkits/dita/tk$com.ibm.streams.timeseries/op$com.ibm.streams.timeseries.modeling$BATS.html) from the timeseries toolkit. The schema of the stream is defined by the NamedTuple `Prediction`. As a speciality we have here\n",
    "- a list of primitive types (`predicted_values` attribute)\n",
    "- a list of another structured type (`prediction_intervals` attribute)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "0e0a4f7a-8ffb-4c0d-b113-02e73880ae63"
   },
   "outputs": [],
   "source": [
    "import streamsx.spl.op as op\n",
    "\n",
    "# The BATS operator specifies the structure of a prediction interval.\n",
    "# The prediction interval MUST have the attributes 'lower' and 'upper' as float in this sequence. \n",
    "class PredictionInterval(typing.NamedTuple):\n",
    "    lower: float\n",
    "    upper: float\n",
    "\n",
    "# outgoing schema of the BATS operator\n",
    "class Prediction(typing.NamedTuple):\n",
    "    sample_timestamp: int\n",
    "    aggregated_call_duration: float\n",
    "    predicted_values: typing.List[float]\n",
    "    prediction_intervals: typing.List[PredictionInterval]\n",
    "\n",
    "# invoke the operator com.ibm.streams.timeseries.modeling::BATS\n",
    "from streamsx.spl.types import uint64\n",
    "\n",
    "bats_invoke = op.Invoke(topo,\n",
    "                        inputs = hourly_call_volumes_regularized,\n",
    "                        kind = \"com.ibm.streams.timeseries.modeling::BATS\",\n",
    "                        schemas = Prediction,\n",
    "                        params = {'initSamples': 14 * 24,  # 14 days * 24 hours\n",
    "                                  'stepAhead': 8,          # predict up to next 8 hours\n",
    "                                  'confidenceLevel': 0.9,\n",
    "                                  'useBoxCox': False,      # cannot be used with zero values\n",
    "                                  'forecastedAllTimeseriesSteps': 'predicted_values',\n",
    "                                  'allConfidenceIntervals': 'prediction_intervals'}\n",
    "                       )\n",
    "\n",
    "# parameter, that references an attribute of the input stream:\n",
    "bats_invoke.params['inputTimeSeries'] = bats_invoke.attribute(hourly_call_volumes_regularized,\n",
    "                                                              'aggregated_call_duration')\n",
    "\n",
    "# output stream:\n",
    "hourly_forecasts = bats_invoke.outputs[0]\n",
    "\n",
    "# Create a view on this stream for data visualization:\n",
    "forecast_view = hourly_forecasts.view(name=\"ForecastedData\", description=\"hourly forecasts\")\n",
    "\n",
    "# Terminate stream by publishing the data to be picked up by another streaming application in this service instance.\n",
    "# Other Stream terminations like publishing to Kafka, writing into data bases, or whatever are possible, of course.\n",
    "hourly_forecasts.publish(topic='callcenter/hourlyforcasts')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "67039add-63a7-4490-bcc5-1b4e04396198"
   },
   "source": [
    "<a name=\"launch\"></a>\n",
    "# 3. Submit the application\n",
    "A running Streams application is called a *job*. By submitting the topology we create a job."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "0215aeac-830a-47dd-b44e-ba66c5b4de0a"
   },
   "outputs": [],
   "source": [
    "submission_result = submit_topology(topo)\n",
    "if submission_result.job:\n",
    "    streams_job = submission_result.job\n",
    "    print(\"Job-ID: \", streams_job.id , \"\\nJob name: \", streams_job.name)\n",
    "else:\n",
    "    print(\"Submission failed: \"   + str(submission_result))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "aa80f9cb-63cb-44ce-855c-883cb652822a"
   },
   "source": [
    "<a name=\"view\"></a>\n",
    "# 4. Visualize the data from the job with `matplotlib`\n",
    "Now that the job is started, use the `View` object to start retrieving data from a `Stream`. We will be using the view object `forecast_view`, which we have created from the `hourly_forecasts` stream.\n",
    "The forecasting operator will suppress initial input tuples to calculate the BATS algorithm parameters from. Therefore it takes some time until the first data is visualized."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "88be23b7-b2bd-436f-813e-641468103bf0"
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "\n",
    "plot_window = 120\n",
    "def create_plot(timeticks, current_vals, next_forecasts, forecasts_lower, forecasts_upper):\n",
    "    figure, axes = plt.subplots(figsize=(10, 5), dpi=96)\n",
    "    plt.title('Actual and forecasted hourly call volumes')\n",
    "    axes.set_xlabel('timeticks [hour]')\n",
    "    axes.set_ylabel('hourly call volume [seconds]')\n",
    "    axes.grid(True, axis='both') # axis = {x,y,both}\n",
    "    axes.plot(timeticks, current_vals, linewidth=3, label='actual data')\n",
    "    axes.plot(timeticks, next_forecasts, linewidth=1, label='next hour forecast')\n",
    "    axes.plot(timeticks, forecasts_lower, linewidth=1, linestyle='dashed', label='forecast interval lower')\n",
    "    axes.plot(timeticks, forecasts_upper, linewidth=1, linestyle='dashed', label='forecast interval upper')\n",
    "    axes.set_xlim(0, plot_window)\n",
    "    axes.set_ylim(-2000, 20000)\n",
    "    plt.legend(loc='upper center', ncol=4, borderaxespad=0.0, fontsize='small')\n",
    "    return figure, axes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "ac54d868c7444de2b440abd02ef439c1"
   },
   "source": [
    "**Be patient, the BATS forecasting algorithm is being trained. It can take some minutes until data is visualized.**"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "d7fedbb9b0ee476d8b9d89640d7896c4"
   },
   "source": [
    "%matplotlib notebook\n",
    "\n",
    "from collections import deque\n",
    "from queue import Empty\n",
    "n_tuples_to_visualize = 360\n",
    "next_forecast = 1     # 1..8 forecast horizon = next 1 hour\n",
    "\n",
    "timeticks = deque(maxlen=plot_window)\n",
    "aggregated_call_durations = deque(maxlen=plot_window)\n",
    "next_forecasts = deque(maxlen=plot_window)\n",
    "forecasts_lower = deque(maxlen=plot_window)\n",
    "forecasts_upper = deque(maxlen=plot_window)\n",
    "figure, axes = create_plot([], [], [], [], [])\n",
    "\n",
    "try:\n",
    "    tuple_queue = forecast_view.start_data_fetch()\n",
    "    timetick_offset = -1\n",
    "    i = 0\n",
    "    timeouts = 0\n",
    "    batch = []\n",
    "    while (i < n_tuples_to_visualize) and (timeouts < 360):\n",
    "        # visualize the data from the view in batches as we want the figure not be drawn too frequently\n",
    "        batch.clear()\n",
    "        while len(batch) < max(1, time_compression/3600):\n",
    "            try:\n",
    "                tuple = tuple_queue.get(timeout=1)\n",
    "                batch.append(tuple)                \n",
    "                timeouts = 0\n",
    "                i += 1\n",
    "            except Empty:\n",
    "                timeouts += 1\n",
    "                break\n",
    "        \n",
    "        if len(batch) > 0:\n",
    "            for tuple in batch:\n",
    "                if timetick_offset == -1:\n",
    "                    timetick_offset = tuple['sample_timestamp']\n",
    "                timeticks.append((tuple['sample_timestamp'] - timetick_offset)/3600)\n",
    "                aggregated_call_durations.append(tuple['aggregated_call_duration'])\n",
    "                next_forecasts.append(tuple['predicted_values'][next_forecast-1])\n",
    "                forecasts_lower.append(tuple['prediction_intervals'][next_forecast-1]['lower'])\n",
    "                forecasts_upper.append(tuple['prediction_intervals'][next_forecast-1]['upper'])\n",
    "            axes.lines[0].set_xdata(timeticks)\n",
    "            axes.lines[0].set_ydata(aggregated_call_durations)\n",
    "            axes.lines[1].set_xdata(timeticks)\n",
    "            axes.lines[1].set_ydata(next_forecasts)\n",
    "            axes.lines[2].set_xdata(timeticks)\n",
    "            axes.lines[2].set_ydata(forecasts_lower)\n",
    "            axes.lines[3].set_xdata(timeticks)\n",
    "            axes.lines[3].set_ydata(forecasts_upper)\n",
    "            axes.set_xlim(timeticks[0], timeticks[0] + plot_window)\n",
    "            figure.canvas.draw()\n",
    "finally:\n",
    "    forecast_view.stop_data_fetch()\n",
    "    print('data fetch from view ended')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "a8bfe565-b9bb-49a0-9f5c-24012c802375"
   },
   "source": [
    "\n",
    "<a name=\"cancel\"></a>\n",
    "\n",
    "# 5. Cancel the job\n",
    "\n",
    "The Streams job is running in the Streams service. You can cancel it within the notebook or delete it from **Projects** > **Jobs**."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "23a56879-63e5-4882-a625-c819abc27913"
   },
   "outputs": [],
   "source": [
    "# cancel the job directly using the Job object\n",
    "submission_result.cancel_job_button()\n",
    "# or, use \n",
    "# streams_job.cancel()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "id": "80b861b4-bc21-4333-9f0a-dd249ba95796"
   },
   "source": [
    "## Learn more\n",
    "\n",
    "Learn more about the [Python API from the documentation](https://streamsxtopology.readthedocs.io/en/stable/index.html/). You can also visit the [Streams community for more resources](https://ibm.biz/streams-articles).\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "23620cf7-de97-4a31-8f1d-839dd41072a8"
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.6",
   "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.6.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}