{ "cells": [ { "cell_type": "code", "execution_count": 93, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 102, "metadata": {}, "outputs": [], "source": [ "def put(T,b,offset=(0,0)):\n", " ny,nx = b.shape\n", " i,j = offset\n", " T[i:i+ny,j:j+nx] = b" ] }, { "cell_type": "code", "execution_count": 103, "metadata": {}, "outputs": [], "source": [ "nx,ny = 10, 20\n", "T = np.zeros( (ny,nx),dtype=np.int)\n", "b = np.array([ [1,0],[1,0],[1,0],[1,1] ],dtype=np.int)" ] }, { "cell_type": "code", "execution_count": 104, "metadata": {}, "outputs": [], "source": [ "put(T,b,(0,2))\n", "put(T,b,(ny-1-3,1))\n", "put(T,b.T,(ny-1-3,5))\n", "put(T,b.T[::-1,:],(ny-1-2,3))" ] }, { "cell_type": "code", "execution_count": 130, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 130, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJgAAAD4CAYAAAAQNi97AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAMEUlEQVR4nO3df+hddR3H8eeruTm2FBXbdD8qiSWssBVjGhJopm1DWoXURpRZ8FVxUFDQKrD+DKKiULRVo4TSgloNWs4hQQn+2Hdjbq5t+m1M9vUrWyVtmub81rs/7vnK7e5c9905973vPff7esC455zP597zuePFOfee+/28jyICsyxvmuoB2GBzwCyVA2apHDBL5YBZqnOmegBlZuncmM3cnr/uO694edJ9n94zp+f7H1T/5l+cjFdV1taXAZvNXK7UdT1/3W3bdk+674cXLOv5/gfV4/Fw1zafIi1VrYBJWinpoKQRSRtK2iXpB0X7Hknvq7M/a57KAZM0A7gbWAUsBdZJWtrRbRWwpPg3BNxTdX/WTHWOYCuAkYg4FBEngQeANR191gD3RctjwAWSLq2xT2uYOgFbCBxpWx8ttp1pHwAkDUkaljT8Gq/WGJb1kzoBK/ta2vnL+WT6tDZGbIyI5RGxfCbn1hiW9ZM6ARsFFretLwLGKvSxAVYnYDuAJZIukzQLWAts6eizBfhM8W3yKuB4RDxfY5/WMJUvtEbEuKT1wDZgBrApIvZJuq1ovxfYCqwGRoCXgVvqD9mapNaV/IjYSitE7dvubVsO4I46+7Bm85V8S+WAWSoHzFI5YJbKAbNUDpilcsAslQNmqRwwS+WAWSoHzFI5YJbKAbNUDpilcsAslQNmqRwwS+WAWao6M7sXS/qjpP2S9kn6QkmfayQdl7S7+HdnveFa09T5m/xx4EsRsUvSecBOSdsj4i8d/f4cETfW2I81WOUjWEQ8HxG7iuUXgf10mbVt01dP6oNJejvwXuDxkub3S3qS1oTbL0fEvi6vMUSrQAqzmfrib9vGXEusF2oHTNKbgV8DX4yIEx3Nu4C3RcRLklYDv6VVaecUEbER2Ahwvi5y8f4BUbc+2Exa4fp5RPymsz0iTkTES8XyVmCmpIvr7NOapc63SAE/AfZHxHe79Lmk6IekFcX+/lF1n9Y8dU6RVwOfBvZKmvjA8jXgrfD6DO+bgNsljQOvAGvD966ZVurUpniE8vJM7X3uAu6qug9rPl/Jt1QOmKVywCyVA2apHDBL1Ze3ksnin3TOPh/BLJUDZqkcMEvlgFkqB8xSOWCWygGzVA6YpXLALJUDZqkcMEvlgFmqurOKDkvaW5QFGC5pl6QfSBqRtEfS++rsz5qnF39NcW1E/L1L2ypa8yCXAFcC9xSPNk1knyLXAPdFy2PABZIuTd6n9ZG6AQvgIUk7i6n/nRYCR9rWR+lSv0LSkKRhScOv8WrNYVm/qHuKvDoixiTNA7ZLOhARf2prL5vWVjov0qUDBlOtI1hEjBWPx4DNwIqOLqPA4rb1RbSKoNg0Uad0wNyiLhiS5gI3AE91dNsCfKb4NnkVcDwinq88WmucOqfI+cDmovTEOcAvIuJBSbfB66UDtgKrgRHgZeCWesO1pqlTOuAQ8J6S7fe2LQdwR9V9WPP5Sr6lcsAslQNmqRwwS+WAWSoHzFI5YJbKAbNUDpilcsAslQNmqRwwS+WAWSoHzFI5YJbKAbNUDpilcsAsVZ1JH5cXJQMm/p2Q9MWOPtdIOt7W5876Q7YmqfM3+QeBZQCSZgDP0Zq61unPEXFj1f1Ys/XqFHkd8NeIeLZHr2cDolcBWwvc36Xt/ZKelPQHSe/q9gIuHTCYVPcOx5Jm0Zqt/a6IONrRdj7w34h4SdJq4PsRseR0r3m+LoordV2tcdnZ83g8zIl4ofTux704gq0CdnWGCyAiTkTES8XyVmCmpIt7sE9riF4EbB1dTo+SLlEx9VvSimJ//+jBPq0halXXkTQHuB64tW1be+mAm4DbJY0DrwBro+452Rql9mewDP4M1izZn8HMunLALJUDZqkcMEvlgFkqB8xSOWCWygGzVA6YpXLALJUDZqkcMEvlgFkqB8xSOWCWygGzVA6YpXLALNVpAyZpk6Rjkp5q23aRpO2SnikeL+zy3JWSDkoakbShlwO3ZpjMEeynwMqObRuAh4s5jg8X6/+nKCdwN61pbUuBdZKW1hqtNc5pA1bcg/uFjs1rgJ8Vyz8DPlry1BXASEQcioiTwAPF82waqfoZbP7ErZGLx3klfRYCR9rWR4ttNo3Umhd5GmXTmLrOkZM0BAwBzGZO1pjsLKt6BDsq6VKA4vFYSZ9RYHHb+iJaNSxKRcTGiFgeEctncm7FYVm/qRqwLcDNxfLNwO9K+uwAlki6rCiQsrZ4nk0jk7lMcT/wKHC5pFFJnwe+BVwv6RlapQO+VfRdIGkrQESMA+uBbcB+4FcRsS/nbVi/cukAq82lA2zKOGCWygGzVA6YpXLALJUDZqkcMEvlgFkqB8xSOWCWygGzVA6YpXLALJUDZqkcMEvlgFkqB8xSOWCWygGzVFVrU3xb0gFJeyRtlnRBl+celrRX0m5Jw70cuDVD1doU24F3R8QVwNPAV9/g+ddGxLKIWF5tiNZklWpTRMRDxbQ0gMdoTao1O0UvPoN9DvhDl7YAHpK0sygN0JWkIUnDkoZf49UeDMv6Qd17dn8dGAd+3qXL1RExJmkesF3SgeKIeIqI2AhshNa8yDrjsv5R+Qgm6WbgRuBT3W70HhFjxeMxYDOtkk42jVQKmKSVwFeAj0TEy136zJV03sQycAPwVFlfG1xVa1PcBZxH67S3W9K9Rd/Xa1MA84FHJD0JPAH8PiIeTHkX1rdcm8Jqc20KmzIOmKVywCyVA2apHDBL5YBZKgfMUjlglsoBs1QOmKVywCyVA2apHDBL5YBZKgfMUjlglsoBs1QOmKWqWjrgm5KeK/4ef7ek1V2eu1LSQUkjkjb0cuDWDFVLBwB8rygJsCwitnY2SpoB3A2sApYC6yQtrTNYa55KpQMmaQUwEhGHIuIk8ACwpsLrWIPV+Qy2vqius0nShSXtC4EjbeujxbZSLh0wmKoG7B7gHcAy4HngOyV9yqYxdZ0jFxEbI2J5RCyfybkVh2X9plLAIuJoRPwnIv4L/IjykgCjwOK29UXAWJX9WXNVLR1wadvqxygvCbADWCLpMkmzgLXAlir7s+Y6bXWdonTANcDFkkaBbwDXSFpG65R3GLi16LsA+HFErI6IcUnrgW3ADGBTROxLeRfWtxpfOmDb2O5Jv+6HFyyrOqSeOJOx9oPJ/n+5dIBNGQfMUjlglsoBs1QOmKVywCyVA2apHDBL5YBZKgfMUtW604edmTP5qSrrJ7Cz/XOVj2CWygGzVA6YpXLALJUDZqkcMEvlgFmqyfxN/iZaNx49FhHvLrb9Eri86HIB8M+IOOVijKTDwIvAf4Bx3xh++pnMhdaf0ro/5H0TGyLikxPLkr4DHH+D518bEX+vOkBrttMGLCL+JOntZW2SBHwC+GBvh2WDou5PRR8AjkbEM13aA3hIUgA/LG78XkrSEDAEMJs5NYdVX9NmAE3W2Z5ZVTdg64D736D96ogYkzSP1u2XDxTFVE5RhG8jtKat1RyX9YnK3yIlnQN8HPhltz4RMVY8HgM2U15iwAZYncsUHwIORMRoWaOkuZLOm1gGbqC8xIANsMlUOLwfeBS4XNKopM8XTWvpOD1KWiBpohjdfOARSU8CTwC/j4gHezd0a4LJfItc12X7Z0u2jQGri+VDwHtqjs8azlfyLZUDZqkcMEvlgFkqB8xSTatZRU0qVjcofASzVA6YpXLALJUDZqkcMEvlgFkqB8xSOWCWygGzVA6YperLexVJ+hvwbMfmi4FBnF85CO/rbRHxlrKGvgxYGUnDgzgzfFDf1wSfIi2VA2apmhSwrrPCG25Q3xfQoM9g1kxNOoJZAzlglqrvAyZppaSDkkYkbZjq8fSSpMOS9kraLWl4qseToa8/g0maATwNXA+MAjuAdRHxlykdWI8UFSCXD3KBvn4/gq0ARiLiUEScBB4A1kzxmOwM9HvAFgJH2tZHi22DYqJA386iAN/A6fdpayrZ1r/n9DM36QJ9TdXvR7BRYHHb+iJgbIrG0nPToUBfvwdsB7BE0mWSZtGqSbZlisfUE9OlQF9fnyIjYlzSemAbMAPYFBH7pnhYvTIf2Nwq1M05wC8GsUBfX1+msObr91OkNZwDZqkcMEvlgFkqB8xSOWCWygGzVP8DV3GUyWKffGQAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.imshow(T)" ] }, { "cell_type": "code", "execution_count": 153, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([0, 0, 4, 4, 0, 0, 0, 0, 0, 0])" ] }, "execution_count": 153, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.sum(np.cumsum(T[10::-1,:],axis=0)>0,axis=0)" ] }, { "cell_type": "code", "execution_count": 152, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "6" ] }, "execution_count": 152, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.max(np.sum(np.cumsum(T[10:,:],axis=0)>0,axis=0)+\\\n", " np.sum(np.cumsum(T[10::-1,:],axis=0)>0,axis=0))" ] }, { "cell_type": "code", "execution_count": 147, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 147, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAKXklEQVR4nO3bcaidd33H8fdnubE1dcWWorRJIRVKNylsdZdZLYg0yropdn9sUKHSiZB/Nq0iSN0//dc/RPSPIVxqVbBURixYylC7aBmDEZamBdNGaanSxkbTMaZSWJLid3/c4wi3Se7NeZ57npx+3y8I95znnud5voSb932ec35JVSGprz+YegBJ0zICUnNGQGrOCEjNGQGpOSMgNbdpBJI8mORkkqNnbbs6yeNJnpt9vWp7x5S0XbZyJfAN4I4N2+4DDlbVjcDB2XNJSyhbWSyUZC/wWFXdPHv+U+D9VXUiybXAE1V103YOKml7rMy539ur6gTALARvO98Lk+wH9gPsYMef7eLKOU8paTP/y6ucrlO5mH3mjcCWVdUasAZwZa6ud2ffdp9SautQHbzofeb9dOBXs9sAZl9PznkcSRObNwKPAvfMHt8DfHeccSQt2lY+InwY+A/gpiTHk3wC+ALwwSTPAR+cPZe0hDZ9T6CqPnqeb3lzL70BuGJQas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqblBEUjymSTPJDma5OEkl481mKTFmDsCSXYDnwJWq+pmYAdw11iDSVqMobcDK8Cbk6wAu4CXh48kaZHmjkBV/QL4IvAicAL4dVX9YOPrkuxPcjjJ4TOcmn9SSdtiyO3AVcCdwA3AdcAVSe7e+LqqWquq1apa3cll808qaVsMuR34APCzqnqlqs4AjwDvHWcsSYsyJAIvArcm2ZUkwD7g2DhjSVqUIe8JHAIOAEeAH8+OtTbSXJIWZGXIzlV1P3D/SLNImoArBqXmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNTcoAgkeWuSA0l+kuRYkveMNZikxVgZuP9XgO9V1d8keROwa4SZJC3Q3BFIciXwPuDvAKrqNHB6nLEkLcqQ24F3AK8AX0/yVJIHklyx8UVJ9ic5nOTwGU4NOJ2k7TAkAivAu4CvVtUtwKvAfRtfVFVrVbVaVas7uWzA6SRthyEROA4cr6pDs+cHWI+CpCUydwSq6pfAS0lumm3aBzw7ylSSFmbopwOfBB6afTLwAvDx4SNJWqRBEaiqp4HVkWaRNAFXDErNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1NzgCSXYkeSrJY2MMJGmxxrgSuBc4NsJxJE1gUASS7AE+BDwwzjiSFm3olcCXgc8BvzvfC5LsT3I4yeEznBp4OkljmzsCST4MnKyqJy/0uqpaq6rVqlrdyWXznk7SNhlyJXAb8JEkPwe+Ddye5FujTCVpYeaOQFV9vqr2VNVe4C7gh1V192iTSVoI1wlIza2McZCqegJ4YoxjSVosrwSk5oyA1NwotwNvVN9/+empR9iyv7juT6ceQUvKKwGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAam7uCCS5PsmPkhxL8kySe8ccTNJirAzY9zXgs1V1JMkfAk8mebyqnh1pNkkLMPeVQFWdqKojs8e/BY4Bu8caTNJiDLkS+H9J9gK3AIfO8b39wH6Ay9k1xukkjWjwG4NJ3gJ8B/h0Vf1m4/eraq2qVqtqdSeXDT2dpJENikCSnawH4KGqemSckSQt0pBPBwJ8DThWVV8abyRJizTkSuA24GPA7Umenv35q5HmkrQgc78xWFX/DmTEWSRNwBWDUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0NikCSO5L8NMnzSe4bayhJizN3BJLsAP4J+EvgncBHk7xzrMEkLcaQK4E/B56vqheq6jTwbeDOccaStCgrA/bdDbx01vPjwLs3vijJfmD/7Ompf60DRwecc5Gu2XEt/zX1EFv3/DWwNPMu06ywXPPedLE7DIlAzrGtXrehag1YA0hyuKpWB5xzYZZpVliueZdpVliueZMcvth9htwOHAeuP+v5HuDlAceTNIEhEfhP4MYkNyR5E3AX8Og4Y0lalLlvB6rqtST/AHwf2AE8WFXPbLLb2rznm8AyzQrLNe8yzQrLNe9Fz5qq193GS2rEFYNSc0ZAam4hEVim5cVJrk/yoyTHkjyT5N6pZ9pMkh1Jnkry2NSzbCbJW5McSPKT2d/xe6ae6XySfGb2M3A0ycNJLp96prMleTDJySRHz9p2dZLHkzw3+3rVZsfZ9ggs4fLi14DPVtUfA7cCf3+JzwtwL3Bs6iG26CvA96rqj4A/4RKdO8lu4FPAalXdzPqb33dNO9XrfAO4Y8O2+4CDVXUjcHD2/IIWcSWwVMuLq+pEVR2ZPf4t6z+ku6ed6vyS7AE+BDww9SybSXIl8D7gawBVdbqq/mfaqS5oBXhzkhVgF5fYOpiq+jfgvzdsvhP45uzxN4G/3uw4i4jAuZYXX7L/qM6WZC9wC3Bo2kku6MvA54DfTT3IFrwDeAX4+uz25YEkV0w91LlU1S+ALwIvAieAX1fVD6adakveXlUnYP0XGvC2zXZYRAS2tLz4UpPkLcB3gE9X1W+mnudcknwYOFlVT049yxatAO8CvlpVtwCvsoXL1SnM7qXvBG4ArgOuSHL3tFNtj0VEYOmWFyfZyXoAHqqqR6ae5wJuAz6S5Oes32bdnuRb0450QceB41X1+yurA6xH4VL0AeBnVfVKVZ0BHgHeO/FMW/GrJNcCzL6e3GyHRURgqZYXJwnr96zHqupLU89zIVX1+araU1V7Wf97/WFVXbK/rarql8BLSX7/P932Ac9OONKFvAjcmmTX7GdiH5fom5gbPArcM3t8D/DdzXYY8r8It2TO5cVTug34GPDjJE/Ptv1jVf3LhDO9kXwSeGj2C+EF4OMTz3NOVXUoyQHgCOufGD3FJbZ8OMnDwPuBa5IcB+4HvgD8c5JPsB6yv930OC4blnpzxaDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNfd/co1F6rUXSR0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.imshow(np.cumsum(T[10::-1,:],axis=0)>0,extent=[0,nx,0,10])" ] }, { "cell_type": "code", "execution_count": 143, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 143, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAKeUlEQVR4nO3cb4hld33H8fenO2viRoMJoiS7gY0Q0kqgjR1qNCCSVZpWMX1QIUIkFWGe+CeKILFP8tQHIvqgCEOMCoZIWQOGUIzpaiiFsnSzWXCTVRKiJGtWN6VUJdBsQr99MNeynezuzN5z7r1z832/YJl779xzzpdh9j3n3PubSVUhqa8/WvQAkhbLCEjNGQGpOSMgNWcEpOaMgNTclhFIcl+S00mOn/XYlUkeTfL05OMVsx1T0qxs50zg28Ctmx67GzhUVdcBhyb3JS2hbGexUJL9wMNVdcPk/s+B91fVqSRXAY9V1fWzHFTSbKxMud3bq+oUwCQEbzvfE5OsAWsAu9j153u4fMpDStrKf/MSZ+rlXMw200Zg26pqHVgHuDxX1rtzYNaHlNo6XIcueptp3x34zeQygMnH01PuR9KCTRuBh4A7J7fvBH4wzjiS5m07bxE+APwbcH2Sk0k+CXwZ+GCSp4EPTu5LWkJbviZQVR87z6e8uJdeB1wxKDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNTcoAgk+XySJ5McT/JAkkvHGkzSfEwdgSR7gc8Cq1V1A7ALuH2swSTNx9DLgRXgjUlWgD3AC8NHkjRPU0egqn4FfAV4DjgF/LaqfrT5eUnWkhxJcuQVXp5+UkkzMeRy4ArgNuBa4GrgsiR3bH5eVa1X1WpVre7mkuknlTQTQy4HPgD8oqperKpXgAeB944zlqR5GRKB54CbkuxJEuAAcGKcsSTNy5DXBA4DB4GjwE8n+1ofaS5Jc7IyZOOquge4Z6RZJC2AKwal5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDU3KAIJHlLkoNJfpbkRJL3jDWYpPlYGbj914EfVtXfJnkDsGeEmSTN0dQRSHI58D7g7wCq6gxwZpyxJM3LkMuBdwAvAt9K8kSSe5NctvlJSdaSHEly5BVeHnA4SbMwJAIrwLuAb1TVjcBLwN2bn1RV61W1WlWru7lkwOEkzcKQCJwETlbV4cn9g2xEQdISmToCVfVr4Pkk108eOgA8NcpUkuZm6LsDnwHun7wz8CzwieEjSZqnQRGoqmPA6kizSFoAVwxKzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNTc4Akl2JXkiycNjDCRpvsY4E7gLODHCfiQtwKAIJNkHfAi4d5xxJM3bysDtvwZ8EXjz+Z6QZA1YA7iUPQMPd26PvHBsJvv9y6v/bCb7XSaz+tpq53x/TX0mkOTDwOmqevxCz6uq9apararV3Vwy7eEkzciQy4GbgY8k+SXwPeCWJN8dZSpJczN1BKrqS1W1r6r2A7cDP66qO0abTNJcuE5Aam7oC4MAVNVjwGNj7EvSfHkmIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOaMgNScEZCaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGrOCEjNGQGpOSMgNWcEpOamjkCSa5L8JMmJJE8muWvMwSTNx8qAbV8FvlBVR5O8GXg8yaNV9dRIs0mag6nPBKrqVFUdndz+PXAC2DvWYJLmY8iZwP9Jsh+4ETh8js+tAWsAl7JnjMMttUdeOLboEaT/Z/ALg0neBHwf+FxV/W7z56tqvapWq2p1N5cMPZykkQ2KQJLdbATg/qp6cJyRJM3TkHcHAnwTOFFVXx1vJEnzNORM4Gbg48AtSY5N/v31SHNJmpOpXxisqn8FMuIskhbAFYNSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmhvlrw2/XvmXgdWBZwJSc0ZAas4ISM0ZAak5IyA1ZwSk5oyA1JwRkJozAlJzRkBqzghIzRkBqTkjIDVnBKTmjIDUnBGQmjMCUnNGQGpuUASS3Jrk50meSXL3WENJmp+pI5BkF/APwF8B7wQ+luSdYw0maT6GnAn8BfBMVT1bVWeA7wG3jTOWpHkZ8teG9wLPn3X/JPDuzU9KsgasTe6+/M918PiAY57TrqvG3iMAbwX+YyZ7no1lmneZZoWZzfvM+LuE6y92gyERyDkeq9c8ULUOrAMkOVJVqwOOOTfLNCss17zLNCss17xJjlzsNkMuB04C15x1fx/wwoD9SVqAIRH4d+C6JNcmeQNwO/DQOGNJmpepLweq6tUknwYeAXYB91XVk1tstj7t8RZgmWaF5Zp3mWaF5Zr3omdN1Wsu4yU14opBqTkjIDU3lwgs0/LiJNck+UmSE0meTHLXomfaSpJdSZ5I8vCiZ9lKkrckOZjkZ5Ov8XsWPdP5JPn85HvgeJIHkly66JnOluS+JKeTHD/rsSuTPJrk6cnHK7baz8wjsITLi18FvlBVfwLcBHxqh88LcBdwYtFDbNPXgR9W1R8Df8oOnTvJXuCzwGpV3cDGi9+3L3aq1/g2cOumx+4GDlXVdcChyf0LmseZwFItL66qU1V1dHL792x8k+5d7FTnl2Qf8CHg3kXPspUklwPvA74JUFVnquq/FjvVBa0Ab0yyAuxhh62Dqap/Af5z08O3Ad+Z3P4O8Ddb7WceETjX8uId+5/qbEn2AzcChxc7yQV9Dfgi8D+LHmQb3gG8CHxrcvlyb5LLFj3UuVTVr4CvAM8Bp4DfVtWPFjvVtry9qk7Bxg804G1bbTCPCGxrefFOk+RNwPeBz1XV7xY9z7kk+TBwuqoeX/Qs27QCvAv4RlXdCLzENk5XF2FyLX0bcC1wNXBZkjsWO9VszCMCS7e8OMluNgJwf1U9uOh5LuBm4CNJfsnGZdYtSb672JEu6CRwsqr+cGZ1kI0o7EQfAH5RVS9W1SvAg8B7FzzTdvwmyVUAk4+nt9pgHhFYquXFScLGNeuJqvrqoue5kKr6UlXtq6r9bHxdf1xVO/anVVX9Gng+yR9+0+0A8NQCR7qQ54CbkuyZfE8cYIe+iLnJQ8Cdk9t3Aj/YaoMhv0W4LVMuL16km4GPAz9Ncmzy2N9X1T8tcKbXk88A909+IDwLfGLB85xTVR1OchA4ysY7Rk+ww5YPJ3kAeD/w1iQngXuALwP/mOSTbITso1vux2XDUm+uGJSaMwJSc0ZAas4ISM0ZAak5IyA1ZwSk5v4Xi7BJv0Qvfk8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.imshow(np.cumsum(T[10:,:],axis=0)>0,extent=[0,nx,0,10])" ] }, { "cell_type": "code", "execution_count": 124, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 124, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJoAAAD8CAYAAACPUr1QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAKtklEQVR4nO3dYajV9R3H8fdn7karNdxmmzc1F3TZsIHNyU3pidtqqQjugQx7MEMGl8TGButBbODYsz0azAxbsChh1AZbTbbbpMVGBrMyUadZ26UFXZTZbnXN6Yrrvntw/sbhdo7X+v/83vu/5/OCg/9zzs/z+wtvrud4+n9TRGB2qX1kuk/AeoNDsxQOzVI4NEvh0CyFQ7MUH63zmyV9CvgV8DngVeCbEfFmh3WvAm8D54CJiFheZ19rnro/0e4BnoqIAeCp6n43X4mIGx1Zb6ob2nrg4er4YeAbNV/PZinV+WZA0lsRMbft/psR8ckO6/4JvAkE8POIeOACrzkEDAFceYW+/IXrL/vQ5zdT/f3wFdN9CpfEf/kP78Y76vTclO/RJP0JmN/hqR9+gHO4OSKOS/oM8KSklyLi6U4LqwgfAFi+9PJ4bs+iD7BNM9x2zY3TfQqXxLPxVNfnpgwtIm7p9pykf0nqj4gTkvqBk11e43j160lJjwGDQMfQbHaq+x5tN3BHdXwH8LvJCyRdKemq88fA14EjNfe1hqkb2k+AWyX9A7i1uo+kayQNV2s+Czwj6RDwHPCHiPhjzX2tYWr9O1pEjAFf6/D4cWBtdfwKsLTOPtZ8/mbAUjg0S+HQLIVDsxQOzVI4NEvh0CyFQ7MUDs1SODRL4dAshUOzFA7NUjg0S+HQLIVDsxQOzVI4NEtRJDRJqyW9LGlE0vuuVlfL9ur5w5KWldjXmqN2aJLmAPcBa4AlwO2SlkxatgYYqG5DwM66+1qzlPiJNgiMRMQrEfEu8CitUQnt1gO7omUfMLe6DtR6RInQFgCvtd0frR77oGtsFisRWqdZC5MHelzMmtZCaUjSfkn7Xx87V/vkbGYoEdoo0D4gYyFw/EOsAVqzNyJieUQsv/rTcwqcns0EJUJ7HhiQdJ2ky4CNtEYltNsNbKo+fa4AxiPiRIG9rSFqXakOEBETku4C9gBzgAcj4qikO6vn7weGaV25PgKcATbX3deapXZoABExTCum9sfubzsOYGuJvayZ/M2ApXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqVwaJbCoVmKrNkbqySNSzpY3baV2Neao/bFKW2zN26ldf3m85J2R8SLk5bujYh1dfezZsqavWE9Lmv2BsBKSYckPSHphm4v5pEIs1PW7I0DwOKIWArcCzze7cU8EmF2Spm9ERGnIuJ0dTwM9EmaV2Bva4iU2RuS5ktSdTxY7TtWYG9riKzZGxuALZImgLPAxmpMgvWIrNkbO4AdJfayZvI3A5bCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlcGiWwqFZilIjER6UdFLSkS7PS9L2amTCYUnLSuxrzVHqJ9pDwOoLPL8GGKhuQ8DOQvtaQxQJLSKeBt64wJL1wK5o2QfMldRfYm9rhqz3aBc7NsEjEWaprNAuZmxC60GPRJiVskKbcmyCzW5Zoe0GNlWfPlcA4xFxImlvmwGKXKku6RFgFTBP0ijwI6AP3rtifRhYC4wAZ4DNJfa15ig1EuH2KZ4PYGuJvayZ/M2ApXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqVwaJYiayTCKknjkg5Wt20l9rXmKHLNAK2RCDuAXRdYszci1hXazxomaySC9bjM92grJR2S9ISkG7ot8kiE2SkrtAPA4ohYCtwLPN5toUcizE4poUXEqYg4XR0PA32S5mXsbTNDSmiS5ktSdTxY7TuWsbfNDFkjETYAWyRNAGeBjdXV69YjskYi7KD1zx/Wo/zNgKVwaJbCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlcGiWonZokhZJ+rOkY5KOSvpuhzWStF3SiKTDkpbV3deapcQ1AxPA9yPigKSrgBckPRkRL7atWQMMVLebgJ3Vr9Yjav9Ei4gTEXGgOn4bOAYsmLRsPbArWvYBcyX1193bmqPoezRJnwO+BDw76akFwGtt90d5f4znX8MjEWahYqFJ+jjwG+B7EXFq8tMdfkvH6zo9EmF2KjUfrY9WZL+MiN92WDIKLGq7vxA4XmJva4YSnzoF/AI4FhE/7bJsN7Cp+vS5AhiPiBN197bmKPGp82bgW8DfJB2sHvsBcC28NxJhGFgLjABngM0F9rUGqR1aRDxD5/dg7WsC2Fp3L2sufzNgKRyapXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqXIGomwStK4pIPVbVvdfa1ZskYiAOyNiHUF9rMGyhqJYD2uyP8Y9rwLjEQAWCnpEK0Lh++OiKNdXmMIGAK4dkHR05sx9hw/OPWiBhq87UzX57JGIhwAFkfEUuBe4PFur+ORCLNTykiEiDgVEaer42GgT9K8EntbM6SMRJA0v1qHpMFq37G6e1tzZI1E2ABskTQBnAU2VlevW4/IGomwA9hRdy9rLn8zYCkcmqVwaJbCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlKHFxyuWSnpN0qBqJ8OMOayRpu6QRSYclLau7rzVLiYtT3gG+GhGnq8vunpH0RETsa1uzBhiobjcBO6tfrUeUGIkQ56/ZBPqq2+QrnNYDu6q1+4C5kvrr7m3NUeoC4jnVpXYngScjYvJIhAXAa233R/F8jp5SJLSIOBcRNwILgUFJX5y0pNPleB2v65Q0JGm/pP2vj50rcXo2AxT91BkRbwF/AVZPemoUWNR2fyGtYS+dXsOzN2ahEp86r5Y0tzr+GHAL8NKkZbuBTdWnzxXAeEScqLu3NUeJT539wMOS5tAK99cR8XtJd8J7IxGGgbXACHAG2FxgX2uQEiMRDtOaiTb58fvbjgPYWncvay5/M2ApHJqlcGiWwqFZCodmKRyapXBolsKhWQqHZikcmqVwaJbCoVkKh2YpHJqlcGiWwqFZCodmKRyapXBoliJr9sYqSeOSDla3bXX3tWbJmr0BsDci1hXYzxqoxFVQAUw1e8N6XImfaFTXdL4AXA/c12H2BsBKSYdoXaF+d0Qc7fJaQ8BQdff0nP6Rl0uc40WYB/w7aa9MmX+uxd2eUMn/tXl1xfpjwHci4kjb458A/lf99boW+FlEDBTbuABJ+yNi+XSfR2kz5c+VMnsjIk6dH20VEcNAn6R5Jfe2mS1l9oak+ZJUHQ9W+47V3duaI2v2xgZgi6QJ4CywMUr+nV3GA9N9ApfIjPhzFX2PZtaNvxmwFA7NUvR8aJJWS3q5Gk1/z3SfTymSHpR0UtKRqVdfej0dWvUB5j5a4+mXALdLWjK9Z1XMQ7x/xOu06enQgEFgJCJeiYh3gUdpjapvvIh4Gnhjus/jvF4PzWPpk/R6aBc9lt7q6fXQLnosvdXT66E9DwxIuk7SZcBGWqPqrbCeDi0iJoC7gD3AMVpfn3X8z5eaRtIjwF+Bz0salfTtaT0ffwVlGXr6J5rlcWiWwqFZCodmKRyapXBolsKhWYr/A72Emgwj9d1zAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.imshow(np.cumsum(b[::,:],axis=0)>0)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "solution": "hidden" }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.8" }, "nbTranslate": { "displayLangs": [ "en", "pl" ], "hotkey": "alt-t", "langInMainMenu": true, "sourceLang": "pl", "targetLang": "en", "useGoogleTranslate": true }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 4 }