some information

wangrongqian2019 · Feb 7, 2020 · 8470393 · 8470393
1 parent c9bcb74
commit 8470393
Show file tree

Hide file tree

Showing 2 changed files with 160 additions and 0 deletions.
diff --git a/.ipynb_checkpoints/srg_test-checkpoint.ipynb b/.ipynb_checkpoints/srg_test-checkpoint.ipynb
@@ -0,0 +1,6 @@
+{
+ "cells": [],
+ "metadata": {},
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
diff --git a/srg_test.ipynb b/srg_test.ipynb
@@ -0,0 +1,154 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 211,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy.matlib\n",
+    "import numpy as np"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "2 * 2  matrix"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 212,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "a = np.array([1,2,-1,2])\n",
+    "a = a.reshape(2,2)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "eigenvalue"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 213,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "g,h = numpy.linalg.eig(a)\n",
+    "ss = np.zeros(2)\n",
+    "tt = np.zeros(2)\n",
+    "ss[0] = g[0].real\n",
+    "tt[0] = g[0].imag\n",
+    "ss[1] = g[1].real\n",
+    "tt[1] = g[1].imag"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 214,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "n = 1000\n",
+    "s = np.zeros(n)\n",
+    "t = np.zeros(n)\n",
+    "p = np.zeros(n)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "SRG:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 215,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "for i in range(n):\n",
+    "    x = np.random.rand(2)\n",
+    "    x = 20*(x.reshape(2,1)-0.5)\n",
+    "    #print(x)\n",
+    "    u = np.matmul(a,x)\n",
+    "    #print(u)\n",
+    "\n",
+    "    unorm = np.linalg.norm(u)\n",
+    "    xnorm = np.linalg.norm(x)\n",
+    "    r = unorm / xnorm \n",
+    "    #print(unorm,xnorm,r,np.vdot(u,x))\n",
+    "\n",
+    "    cos_theta = np.vdot(u,x)/(unorm * xnorm)\n",
+    "    theta = np.arccos(cos_theta)\n",
+    "    angle = theta*360/2/np.pi\n",
+    "    #print(theta,angle)\n",
+    "    \n",
+    "    s[i] = r * np.cos(theta)\n",
+    "    t[i] = r * np.sin(theta)\n",
+    "    p[i] = r * np.sin(-theta)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 216,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAcF0lEQVR4nO3df5xU9X3v8ddnhl1YVgMSQCiyBbzcEJtgA9ugsbfam4dGRSPaxCaatrE+5Hpv7AND7JWoN65JrMlNQsO9pk21mh+NTdybkFXBxEvbaPq4Cs3C4lKrqYqCrqAYfkRhcZfle/+YWRxmZ3bO7Jxf35n38/Hw4e7smZkPsOc93/P9fs455pxDRET8lUm6ABERqY2CXETEcwpyERHPKchFRDynIBcR8dy4JN506tSpbs6cOUm8tYiItzZv3vy6c25a8eOJBPmcOXPo7u5O4q1FRLxlZjtKPa6pFRERzynIRUQ8pyAXEfGcglxExHMKchERzynIRUQ8l0j7oUg5XT19fOWRX9K3vz+U1zODoBf4bMrA++dO4fHn91LqKVkzhpxj1uQW5ryzZcR2LU0ZDg8eZVJLE28cHmSo6EUs//9y5UxuaaLjw7/FsvfNClawSJ4lcRnb9vZ2pz7y+jYcyK/s7+c3Jrfw5x96F8veN+u4oB4ORhnprFOn8OKv+nllfz+TJzax/9DgiA+AkyY2cdrME9m4fR9DzpE14+NLZvPFZe9NpGaJnpltds61j3hcQS7lFIbxpJYmzGDfocFRR7ktTRn+YPEp/GhzH/2DQwWPZ/mDxbNGPC7hmz+9le17DgX6kMwafO3y39ZRgCcU5HJMV08fHQ8+xf7+wVjfVyPw9DLKT/kMO2liE7derKmfJCnIG0RXTx+3PfQU+w7FG9LSWIaD/6SJTTgHB/oHj5tCk2goyOtIV08fn+ncOmIxLe00Im9sWsytnYLcQ1fe/QT/7/m9SZcxJsWH6knNkfvUtdJI5k9v5dDA0RGL4TI6BXmK1dt0SAa44ow2fvbMnpq6VszgyiVtdd+FcUvXNu7buPO4gB9r18q8aRN59rWDcZYfmneMz3Jw4Kg6cEahIE+J4tDOAEeTLSmwrMGQG32U29KU4Y7LFmp0laBburbx/U0v1c001vhxGVqaspqHJ8IgN7PZwHeBGeQy6S7n3JrRntMoQe7z1Eixs06dwn3XnJl0GRKiwqOjjMFRT3O/kbppogzymcBM59wWMzsR2Awsc879W7nn1GuQn7v6Ua8OaxtpB5Cx6erp49P3bx0xrTO8BtLanOXgQLrOC2jKwFc+Wp+98bFNrZjZA8CdzrkN5baplyBfcvsGXn1jIOkyRjV/eisbVp6TdBlSx5I6L6GS4Q+bWXU0HRNLkJvZHODnwHucc78u+tlyYDlAW1vb4h07St6xKNXSGNw6M0/SrKunj5vW9nJo8PiVoAlZ43DM/bP1sH4TeZCb2QnAY8Dtzrm1o23ry4j8lq5t3LdpZ+D2tahkgNV/qLCW+lK8hpREuPs2DRNpkJtZE7AOeMQ5t7rS9mkO8lu6tvG9jTtjfc/CnuuJTRn+wvNRg8hYldr/xmWMIzGsxI7LGF/96Omp3veiXOw04DvAXufc9UGek8YgjzvAMwZXNECPtEgY4jzXIs3TlVEG+e8C/wxs4+2W6Juccw+Xe05agjzO9sCTT2xm083nxvJeIvUujn13YlOG/sGjqepd1wlBeXGFt7pFROIX9o1JiiV9PoWCHFh460/59Vvh97zqdGKR9IlqOsaAv0yo+aBhgzyKEbiu4ibilyh63ZO4FlDDBXlXTx/X3781tNdL8lNYRMITRajHNeXSUEEe1qny9XACgYiMLsyj9qjXxhoiyLt6+rjh/zxZc8/pJ85QW6BIoym8R+24DAyO8bKkUR69132Q1/qpqvAWkUK1nlsSRabUdZDXMpWi/m4RGU0tg8SvhzwyLxfk40J7h4RcefcTYwpxjcBFJIjCRcxqB43X37+V7h17I88ar0fkYxmJh/0JKSKNZSxTLmEd+ZcbkWdqfuWELLz1p1WF+Pzprbz4paUKcRGpyReXvZcXv7SUs06dEvg5r74xwJV3PxFZTV5OrVx59xOBz9DUHLiIRGF4yiXozECUlwbxckQe9C9k/vRWhbiIRGrDynMCj84X3Fz2WoI18TLIg/jEGW26aJWIxOK+a87kE2e0VdwuqhtneBfk/+Gz6ytuc/KJzepIEZFYfXHZe6uaNw+Td0F+pMIHmqZTRCQp911zJvOnt466TRSLnt4FeSWaThGRJFXKoCgWPesqyN8xPpt0CSIisaurID950oSkSxARiV1dBXkYl64VEalFV09f7O9ZV0EOMGdV5a4WEZGofLrCDW1OPrE59Pf0LsgnZK3iNueufjT6QkREiiy5fQOVOsWj6KrzLsifuf3Cits8+9pBlty+IYZqpF6s376e8354Hgu/s5Dzfnge67fryE6C6+rpY86q9bz6xkAi7+9dkAOBmu5ffWMgstNhpb6s376ejsc72HVwFw7HroO76Hi8Q2EugdzStS3w/YEr9ZiPlZdBHvQmp4eHHHNWreeWrm0RVyQ+W7NlDYeHDh/32OGhw6zZsiahisQX565+NPAlbSdkLbLzXLwMcqjuk+17G3cyZ9V6zZ1LSbsP7q7qcRHIzYcH7ZR7x/hsoGnhsfI2yDesPKfqE4Cefe2gRugywozWGVU9Lo3r3NWPMmfV+qrmw08+sZne286PtC5vgxyg97bzA11xrNj3Nu7UYqgcs2LRCiZkjz+ZbEJ2AisWrUioIkmbrp4+5q5aX/W5KhOyFsu1n7y+1VuhBTc/PKZLRLY0ZbjjsoW6c1CDW799PWu2rGH3wd3MaJ3BikUrWDpvadJlSYLGcku3QhOyFvp0SrlbvdVNkENtd7sedtapUwIvpopI/ak1wCG6m7s3RJBD7hDoxh/18taRo2N+Dd0eTqTxhDEQjDo7GibIh4XxqQpw0sQmbr34tzT1IlJnunr6uO2hp9h3aDCU15s/vTXyy2g3XJBD7h/q0/dvrXjKbDU0WhfxVxij7lKimkop1pBBPiys0XmxiU0Z/kILpSKptuT2DaGfOm/AX/7hb8e+7zd0kA/r6unjs2t76R8c+/x5JVosFUlWVKPuYXFMoZSjIC9h7qr1oU67lKJgF4nWWFuPq5E1+Nrl8Y/Ai0Ua5GZ2L3AR8Jpz7j2Vtk9LkEM8o/RCzVnjf37k9MR/IUR8FPVou9isyS38+YfelZr9Neog/z3gTeC7vgV5oVu6tvH3m3ZyNOaDFHXGiIwUxdx2UHEtXlYr8qkVM5sDrPM5yAtFtUAalMJdGklXTx8dDz7F/v5wWgGrlTG4Ykk6w7tQ4kFuZsuB5QBtbW2Ld+zYEcr7xuWWrm3ct3Fn5HPqo9F8u/gs6bAuluSi5VglHuSFfBiRj6arp4///sMnGYh4gSWoyS1NdHxYo3dJh6SPZkfj+2BIQR6hsM8QC1trc5bbL32vgl5CcUvXNr6/6SWGEuh4q5bvwV1MQZ6Arp4+PtO5lZQM3AEwg+F/coMRU0UGXJnShR6Jnk8hXSxtHSZRiLpr5fvAOcBU4FXgVufcPeW2b5QgH5a2ucFqZM2YN20i2/ccYsg5smZ8fMlsBb1nunr6+Mojv6Rvf3/SpYQmrZ0lUdIJQSnU1dPHzT/exsGBoaRLiU29HerGyecBQa0aMbRLUZB7ohF31vnTWzk0cHTMo8XiKaIkF39HWwhvzhqX/85sfvbMHl7Z38/E5iyHBoYS7YRKIx+7SeKiIPdcmjsB0qwpA0MOjrpcr/D4cZmSZ/EOfxicNLEJ5+BA/yC/MbmF318wjfv/ZSeFT5k/vZU9bwwc92GbMThz3hQef36vgjkgjbKrpyCvY3GftiwymuEPxUZYfIxbuSAfl0QxEq5Kc86FC12VulZEqqEzkNNBQd4Alr1vVsUdraunj5vW9nKoaNqhuGtF6tvwB71G035RkAsQLOyHaSrHTy1NGe7QjVDqkoJcqjbW9sFyvcxhd60k6axTp/CLF/dF3rWiNk4ppCCX2FQz6h+r4vbNOLtWfLh6ntQnda2IiHiiXNdKJoliREQkPApyERHPKchFRDynIBcR8ZyCXETEcwpyERHPKchFRDynIBcR8ZyCXETEcwpyERHPKchFRDynIBcR8ZyCXETEcwpyERHPKchFRDynIBcR8ZyCXETEcwpyERHPKchFRDynIBcR8ZyCXETEcwpyERHPKchFRDynIBcR8ZyCXETEcwpyERHPKchFRDwXSpCb2flm9ksze87MVoXxmiIiEkzNQW5mWeAbwAXAacDHzey0Wl9XRESCCWNE/n7gOefcdufcAPAD4JIQXldERAIII8hnAS8VfP9y/rHjmNlyM+s2s+49e/aE8LYiIgLhBLmVeMyNeMC5u5xz7c659mnTpoXwtiIiAuEE+cvA7ILvTwFeCeF1RUQkgDCC/BfAfDOba2bNwMeAB0N4XRERCWBcrS/gnDtiZtcBjwBZ4F7n3FM1VyYiIoHUHOQAzrmHgYfDeC0REamOzuwUEfGcglxExHMKchERzynIRUQ8pyAXEfGcglxExHMKchERzynIRUQ8pyAXEfGcglxExHMKchERzynIRUQ8pyAXEfGcglxExHMKchERzynIRUQ8pyAXEfGcglxExHMKchERzynIRUQ8pyAXEfGcglxExHMKchERzynIRUQ8pyAXEfGcglxExHPjki5ApCrrVsLmb4MbKnjQAFd6+0wzuCPgjoJlYFwLDB4ssWH+NVqm5L7t3weTToH558GWv4OjA29vOnUBHHwN+vcWPD0Di6+CtjPggU/B0AAjZJvhkm/Awsur+iOLVGLOldkBItTe3u66u7tjf19JwHc+DC889vb32fEw9Nbx20xdkAvXAy/FW1tSWqbkPiiaJ8LAIcp+CAV5nQu+rA+GBmJmm51z7SMeV5BL1YrDWTyRP+qYNBs++Dl9AHioXJBraqXRrVsJ3feU/7ll4Z3z4VfPFk1niH/yg7YDL8Haa3L/jWZ4uuii1dGXJjVRkDeC3k74x8/ndmDL5OaLg3JD8Poz0dUm6eWO5j7kR/ugB5h7NvzJg/HUJCUpyOtBNVMd1YS4SBAvPAYdk0r/TKP6WCjIfVFpCqShedK1YlnIZEv/rF6NNqrPjodL7tRcfQi02JkmvZ3wkxuPD4h6E0bXimVh8Sf9HOUdm+Z6ufaulbqVhcu+qYAvQV0radMIoT1MbXLh6+2Eh64vc3RRrwwuu6uhf48iCXIz+yjQAbwbeL9zLlA6N1SQ+96qV9y14vNoWPz/fSyWaYZljXOSVVRB/m7gKPA3wA0KctK/oxzXtaK+Yilj3Urovhcvp33GT4LP7ky6ikhE0kfunHs6/+K1vIzf0h7cAM2tcNHXFdQS3EWrKx91lbxcQgq8deD4Lpr2q+v+CDKUOXIze5QKI3IzWw4sB2hra1u8Y8eOmt83Mb2d5TsTkqA5aEmbNA9wPA72MU+tmNk/ADNK/Ohm59wD+W0epV6nVgpPpkmaTryQepCmaRvP9qlIu1bqMsiT7NtuaoWLNRUiDaa3Ex64bmR7alw8CHVdayWIdSth87fiPftR0yIiOQsvH7kfxLlPFp6h6tn0S61dK5cC/xuYBuwHtjrnPlTpeakZkccx6s42w9Bg7ixBdYaI1C7O+feUtdvqhKBCUZ9MoVG2SHziPLku4ZG6ghxy/+Bdnzr+uhlhOmEm3KArBYokKo6zXi+7O5GBWmMH+Z1LorkUqxYlRdIvqqmYBM7PaMwgj+IfcOoCuG5TuK8pIvGJYgE1pimXxgry3k748bXhnHGm+W6R+hXmNIxl4dJor9rYOO2HYXSieNZ6JCJjVNjyWOsRvBvK3T5v58bY86N+RuRhBHgdX2xHRAIKpSkimmuq1/eIvJYQ1+hbgPXb17Nmyxp2H9zNjNYZrFi0gqXzliZdliShcJTe2wlrrwWqnaaNd3Tuf5CP9XDIg9NxJR7rt6+n4/EODg8dBmDXwV10PN4BoDBvdIWhPpbut+57crf/i3iNLRPpq0fttqnVh/gJM6HjgEJcjlmzZc2xEB92eOgwa7asSagiSaXrNuUGgNVae01uZB8hf0fkX5gBbrCKJxi0/6mmUWSE3Qd3V/W4NLDCAWA1I/S11+T+H9HI3M8R+bqVMNQfbNtMc+4srI79CnEpaUZrqas0l39cBMiN0NuvDr79cJhHwM8g3/ytYNudMBM+t0c94DKqFYtWMCE74bjHJmQnsGLRioQqEm9ctDo3VZttCbb9VxdEUoafQR7kjKy5Z+u6JxLI0nlL6fhABzNbZ2IYM1tn0vGBDi10SnD/Y3du4FjJm7sieXv/5si/EOBwVx0pUqWl85YquKU2NzwT3XWdKvBvRF5pbtyaFOIikozrNoFViNU72kJ/W/+CvJJbX0+6AhFpZJf+zeg/f+tA6G9ZX0E+NZqFBBGRwBJorqivIN+3I+kKRERiV19BHrS3XEQkKhGfxVlKfQU5RNanKSISyNrlo/+8mpOIAvIvyMdPGv3nb+5SmItIMu5cAlS4NHgEZ5j7F+RBrhf+5q7cafwiInFZtzKRHnLwMciD6r5HI3MRicedS4LdE6HSjMIY+RnkQeeY3twVSfO9iAiQG4V3TAo+Eo/oDmR+BvlFq4P3jL91IPcX3TE5kdVkEalDvZ25TKnmzmQRLHIO8zPIYQwXeXe5y0h+58ORlSQida63E74wPX9J2irudzz37Egvo+1vkEPumirV3rHjhcdyI3QthopINXo74cf/BYbequ55l90d+fWfzLkqPlVC0t7e7rq7u8N7wVquOGZNuj6LiJRWy43dIRfiIZ6yb2abnXPtxY/7PSIfVu2dOgq5wfwcukbpIpLX2wm3TakhxC30EB/13epiRF4ojOsBt1+t28KJNJreTvjHz8OBl2p7nQjvh1BuRF5/QQ65Bc0XHgvntTLNsOwbul2cSD3q7YQHrqt+3ruciAeBjRXkw766INxbK+nOQyL+W7cSuu+lqq6TSmI6ii8X5P7d6q0aNzxT+2JFoeGOl2EnzNR9QUV80NsJP7kR+veG95otU+CCL6fiaL2+R+TFwpxyKWXqgtzCq4gkK4rgHpYZB8v+OpEAb8wRebHhaZE72iK53RKvP3P8iF2LpiLRizK0CzW1wsVfT8UIvFhjjciL9XbC2muBoXjeTwunIrXr7YQfXwsu4v02hQOxSBY7zewrwMXAAPA8cJVzbn+l56UmyAuFOZceVIo/4UVSI64R97DseLjkzlTul1EF+XnAPznnjpjZlwGcczdWel4qg7zQ8Km47mgy75/CkYBI5Ho74aHrYfBgMu/vwRpX5O2HZnYp8BHn3JWVtk19kBeKezQwGrU/Sr3o7YSuT8HRgWTr8GyfiiPIHwLud859r8zPlwPLAdra2hbv2OHpHe/DOvsrbB6MJqSBpGkANCxF7YJjNeYgN7N/AGaU+NHNzrkH8tvcDLQDl7kAnwxejchHE9eiS00MLrvL619eSbE0BvawSbPhg5+rq9/9yEbkZvYnwLXAB51zh4I8p26CvNi6lbD52ykP9hLqYKQiERh1ztoI9czIsFgGFl9Vt2tMUS12ng+sBs52zu0J+ry6DfJSoj4JKS5qnawPI0bQ+UC2THKL+2PRoIOPqIL8OWA88Kv8Qxudc9dWel5DBXkpx+bZX4amicmt0o/VaDu9WirjVXgUaFlY/Mnc4z4eGZbT3AoX6XcKGvWiWT5JuvUqVBlovwqe/b+5D6tJp7w9V1m4WGzZ6sNmOKx8PXQOeoQ292zYuz3399dyEvTvY8RURnZ8eFftS4OpC3K//8W/M3KMgtxn61bC5m/5dehbPIfa1AKnXwFP/j0M9sdYRhVTBplm+M0z4YWfU3L+d/iDZ9JsmDJv5HZNrTB4KBe8hw/Uz4g4Cp61/aWFgrwehXETjTiNZQQu/kvxmZK+0UWz6tFofeNpXGRViNeR/BFXHbb4+UhBXq+CHrbG2TKpEXl6jToFZdD+p/6uSzQABXmju2h1sB30uEXKCl0rp39s5Fx4UnPk9c4y4By5FsIxdK1orrouKMglmIWXV3f43HbG2y2WhR0Ixx4fY9dKowjStTJaL7VGzw1Fi53ih7HO+aeqa8Xy5w0cUnudjIkWO8VvOvwXKSuTdAEiIlIbBbmIiOcU5CIinlOQi4h4TkEuIuI5BbmIiOcS6SM3sz1AmDftnAq8HuLrRUV1hkt1hs+XWhu1zt90zk0rfjCRIA+bmXWXapJPG9UZLtUZPl9qVZ3H09SKiIjnFOQiIp6rlyC/K+kCAlKd4VKd4fOlVtVZoC7myEVEGlm9jMhFRBqWglxExHPeBLmZnW9mvzSz58xsVYmft5nZz8ysx8x6zezChOq818xeM7N/LfNzM7P/lf9z9JrZorhrzNdRqc4r8/X1mtnjZnZ63DXm6xi1zoLtfsfMhszsI3HVVqKGirWa2TlmttXMnjKzRG6qGuDffpKZPWRmT+brvCruGvN1zM7v00/n61hRYpvE96eAdUa7PznnUv8fkAWeB+YBzcCTwGlF29wF/Nf816cBLyZU6+8Bi4B/LfPzC4GfkLt77RnAppTW+QHgpPzXF6S1zoLfj38CHgY+kkSdAf9OJwP/BrTlv5+e0jpvAr6c/3oasBdoTqDOmcCi/NcnAv9eYr9PfH8KWGek+5MvI/L3A88557Y75waAHwCXFG3jgHfkv54EvBJjfW8X4dzPyf3il3MJ8F2XsxGYbGYz46nubZXqdM497pzbl/92I3BKLIWNrKPS3yfAnwE/Al6LvqLyAtR6BbDWObczv30i9Qao0wEnmpkBJ+S3PRJHbccV4dwu59yW/NdvAE8Ds4o2S3x/ClJn1PuTL0E+C3ip4PuXGfkP2gF8wsxeJjcy+7N4SqtakD9L2lxNbtSTOmY2C7gU+GbStQTwH4GTzOxRM9tsZn+cdEFl3Am8m9xgaBuwwrmg98uLhpnNAd4HbCr6Uar2p1HqLBT6/uTLrd6sxGPFfZMfB77tnPuamZ0J/J2ZvSfpX8ASgvxZUsPMfp/cL97vJl1LGV8HbnTODeUGkKk2DlgMfBBoAZ4ws43OuX9PtqwRPgRsBf4zcCqwwcz+2Tn36ySKMbMTyB1xXV+ihtTsTxXqHN4mkv3JlyB/GZhd8P0pjJw6uRo4H8A594SZTSB3wZpED7dLCPJnSQUzWwj8LXCBc+5XSddTRjvwg3yITwUuNLMjzrmuZMsq6WXgdefcQeCgmf0cOJ3cnGqaXAV8yeUmdJ8zsxeABcC/xF2ImTWRC8f7nHNrS2ySiv0pQJ2R7k++TK38AphvZnPNrBn4GFB8N96d5EY6mNm7gQnAnlirDOZB4I/zq+1nAAecc7uSLqqYmbUBa4E/SuGI8Rjn3Fzn3Bzn3Bzgh8B/S2mIAzwA/CczG2dmE4El5OZT06ZwXzoZeBewPe4i8nP09wBPO+dWl9ks8f0pSJ1R709ejMidc0fM7DrgEXIdCvc6554ys88D3c65B4HPAHeb2afJHVp9Mj+iiJWZfR84B5ian6+/FWjK/zm+SW7+/kLgOeAQudFP7ALU+TngncBf5Ue7R1wCV5sLUGdqVKrVOfe0mf0U6AWOAn/rnBu1rTKJOoEvAN82s23kpi5udM4lccnYs4A/AraZ2db8YzcBbQW1pmF/ClJnpPuTTtEXEfGcL1MrIiJShoJcRMRzCnIREc8pyEVEPKcgFxHxnIJcRMRzCnIREc/9fybY7gpZQFwfAAAAAElFTkSuQmCC\n",
+      "text/plain": [
+       "<Figure size 432x288 with 1 Axes>"
+      ]
+     },
+     "metadata": {
+      "needs_background": "light"
+     },
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "import matplotlib.pyplot as plt\n",
+    "plt.scatter(s, t)\n",
+    "plt.scatter(s, p)\n",
+    "plt.scatter(ss, tt)\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python [conda env:.conda-devito] *",
+   "language": "python",
+   "name": "conda-env-.conda-devito-py"
+  },
+  "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.5"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 4
+}