2019-11-11 13:36:11 +00:00
|
|
|
|
{
|
|
|
|
|
"cells": [
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 1,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAADFcAAAA/CAYAAACfdcacAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAgAElEQVR4Ae2djdXURrKGgeMAvL4Z4AzAG8FCBrAbgXEG5jgCDs4AHMFeyMBsBDZkYDJYXzLwrZpPGvVILan/u7rq1Tn6PqnVP1VPvd0tzUij+/fu3XtC66+0rpfPf/3117frxPv370flX5fHfhkCFIePVNMjT23PKW7vPelIAgGTBDBmjRd2ieMbdCRbRxI1I5uYHevQd/vEWmKfhBb6aMFtVaIuXPuw3Z4A+mV75jEtSuuz0EtM9NrnlaYXJgDNtNdBbosSdZTrE8qDQA0CGN9qUI2vE2NWPLNRSqCPjRKpWzvRJ295YA8EQAAELBHA3K072pjjdccX3oEACJQngHmxPNMRapQ2X0KHI6hGr42x/eEo/1cOpqfONm/+udpf78bmX5fHfh6B76n4N6sqfA/JrLJgFwTMEsCYNU7oJY9v0JFMHUnWjExi9qxC320bc8l9ElpoqwW3Ncm6cO3EdnsC6JftmYe0KLXPQi8h0WufR6pemAQ0014PqS1K1lGqTygHAjUJYHyrSfe8boxZ54xGz4E+NlYE0SfHihesBQEQAIEaBDB316Dav07M8f1jAAtAAATGJIB5ccy4pVotdb6EDlMjinI5BGL7w27+68MV9LaDDzEWxeYPqZueAvnrLB+1e/8sj8Tj5Ns7suvZmW2h/lG+T+u6qI0v6zTsgwAI3BGoMWbNbEv377le6f9rjdmSx7fSOtKunVb+SdaM9H5sxb7SfXfNrZXW1+3W2C8xtkvuk6W1oCn2KXqK8V+yLlJ8R5lyBEr3S9eyGI265TRux7KQ2mdr6CWWjTZ9kP/Zn8VJ1QvHqoZmZg1Y187Mgf+XYCFZR66v2AYBKQRqjm97PpaYM/bq7p0eO45hzOodsfrt9+hjrlexmnTLjrydOs6gT44cddgOAiAAAmUI9Jy7rc7bc+Rq+o85fqaM/yAAAiAQRwDzYhwvCblz5lOp82VPHUqIKWzoQyC2Pxzlf9DHhW2rNEC8oNQfyNj7R+u2pPyUybdfj/yaj8n3BhaCAAi4BKz2b81jthvfmtvataPdv5raQN1jEdCkdYztcdrTFPs4z+9yW/c/hRnKtCUAjS68wWJhsd6yzmbyX+VncetYl963rh2XJ1i4NLANAnoJaJ4zMI7p1e2onlnVpOZxZlQtwm4QAAEQAIFzAlbn7ZmMdf9nDvgPAiAAAiBwRwDzQpoSwC2NG0rFESCdfU3rQ18pSn+2d8yXv2Xa2m7a/7FV+9c3V7Rq8KCd5/SAwfpVMAfZhzqk2behAgFjQaACAav926rfJSWknWG0f3xCRIC/ofOBzyVBoy4QqEwgWuuV7cmpPsoX9Nl7UbxyAiO0rHX/hYYFZjkEoNEFBlgsLNZb1tlE+4/5/yqhaHbXkvo2ollAR+1FAObnzMHolFF0Xz+tUU4Gzb7JoQxLYghY1aRVv2O0gbzKCeB8JD/AYJjPEDVEE7A+f0X7j34arbFuBRCrbuivDSMGVxTYGIdA9LwwjmtVLQW3qnjrVz7IeP0TkfgvrT+7RMh2finC59x79SYG/6S6+EfdHrttZG7f2E11/0xtvaH/P2TWe1pcxMMV5Cw/TfKGrZ0gMxBe5idlvicYX+6Sxvrr+saWa/PPjcbkm5t02R4hdtJsl2bPJqgHCSPbfuCW9xD5eh27OMPku5rxy+v0nZ9Xv634vMciNV27djL8uzkhSuXbu9zULzZmjDAfboxOSLDkf4bWE8jWLeL6MsUwZD5T0WdTyLq8uHwEs5TmxJXR4v8Utw1freO1JX+1aHQjzoQEsNiHZp2N6z9t80POIXM/AzU7/89qctlxWiS/uRoV/zNYDKGjKbabWK3PFaZ8Nb402LSdkdCV+ciMQnWQERvxRYmB2s8CXd84EFO8Q+dEMbEbXaej219SCFo0GcvE9XvSw3D9cO3zyLqWZnspe6Z6cM62FmvE/sgMS+koAtdhVmn2HBq7Ojiy7StXsneJxfU8mSub2Aw/h4WCyfC/6/VpqH++fCX1P9WFeckHulHayDEoqcUSuKXZU8Kn0Dos+75mRCxMz4trHqH74HZOaupnmDPPUR3moO81Xq4zTGyf0rHn62Mx+1TPI8r/ZCrD3zmeLlQm6AEJn91U+Tsq/3rn2GnboRlEPFxBxnKA5rdWsNPXp0oYIh3/SOu3oU4Jy+f6xqZp88/FzXGaH4iZ09/SxjWec6LA/+/IprmDz+Z9oI1Zl3Naq/9g2Yp0XjuW+rdLyvVb85jm+lx622XIdWvjmORf7ZOe0kE8qG/kMfzAreBD0ubUYMMTMiZpPaGdFkVcX4LGJEV9NoWvy4vLBzFLaUhoGS3+WxuvMT5PHUrB5wwxQ4OW/hrjc2he62xc/4PnMePz/6wtlx2nBfObK1D0P4nFQDo6PVegOSX6S4Me8e/JXAEjS+dQe/J0+7q2Mc/1jf0f1b/T8WovuELS0c+WQGjR5OJR2Jbr96j9cO3pyP1SWp/MZqngfGStr+b7Chiq03VzESwNSmO5WNZ+y52/uHUtc1goyST/e16fhjp2kC97TuK6FYypB4jGOKQgBtLG4iJ9Ywz1bKyUFouNgQ0TkuaFhvZJbQrcDiKjYLw+8E7EoddkBd+fn7XQ+d0nquATxetZREXr+8yDi1J7H6it17R+TdvVXtrwINiiShnJwZun1qiZF5Tm3uTOAXxIafxF1VCLxze2X41/O8HgAf++s47wYMU9svfGbvJNgt03NpGNEmzaCfuSLJTlYmChLaP9my/01Y7ZhaRxWo2HIZdRMzdo9+80wEuGIcfwxfz0LcwD453HevqtmjEpXcn7JT28OLMZZgr9NzNeY3web3zeH4nCjijsr2GOB+Syzsbjv5l5LEAeh1k87Di/SX6GWByeK9D8+olWfo3250PxGD44OiMr51B7EvX0dTVjnsc3xjCyf4fj1V6MJaRb72dzDBRqcnbt8L/H75H74drXIful0D6ZxXL085G1sHrsj85Qo6576IDbFMqyOQ7P/MU2aJrDDpka9z9rTmKwo4+ph+IY5ODoMRA6Fmf3jUHkc2Om0Fjc2Nhix/i8kIwY3M7RjT5en3tYPwfp7AlrjVZ+GGy9PCHG/AP04pYTu9lettt9zqC4D13fXEEA+BUgPLm6bwfgG8h/L+5p4wp3fGMrVPjXGCeaAwFRBKz27x2/MaZFqHOHIdeggmOqf1SOT3b4Icq/0zlB1qvGGCYWEKhNIFXrte1KqX/Hl8MxyXKf3eHF6A+ZpcRGYhnr/kuMCWy6JQCNLjzAYmGx3rLOZsf/03mMypk/Z99hxxI75bfW4ej7qSygo/aRB/Nz5mC0z2inr6sY83Z8Yxgq/NuPKo5IJWBVkzt+ox9KFSrsqkaA+oL5661cuGCYSxDlYwjszF9chYk5LNV/9NMYlfXNi1j15c+tIwb9YwALwgmkzgvhLejMCW464jrIeM335r+i9eZND2Q7vzmi2lsfCkTYa7dT72+0zXneO2lFN6u/uYKDQOuLHat/onR+M8V1oZsq39LqBo0vQD5TGr86RNQS6xsbP5J/omDDGBBoTMBq/471G2PaVpixDLmGkThW8o9Pdt7S+ozq5wcvsYBAdwKVtN7Fr1hfAsYk1X02lhcHNYBZl9inNGrd/xRmKNOWADS68AaLhcV6yzqbWP8D5zHV8/+soVh2XC6Q39zEMP8rsTChI2FBBvPzgJhmFNvXRxrzYn1jqYzk37m0kUMaAauajPUb/VCacmFPIwKmz0cKMQbDQiBRzR2B2PmLS2mawyr5j346TgdDrPrHCjHoHwNY4BCoNC84LejcBDedcV15JXq8Zg2SvfxgBT/Qz/fgu/fl87E/aRW3nNg928u+sA/VlqoPV5CTHBR+OOINbd/cLDntP6KA7b5WhPLwr1g/o/UxraKWXN/YGcn+1YTNftP6B603mqjZ5sSbn77ilV9zw5qs2rlq+zP5BJYVQJM2ssauOTb0X+T4tYcs128qL3bM3vO5dHouQ7ZHMsca/lGdPBbvnci
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2}}{2}\\\\\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2}}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}1 & 0\\\\0 & i\\end{matrix}\\right], \\ \\left[\\begin{matrix}1 & 0\\\\0 & 1\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2}}{2}\\\\\\frac{\\sqrt{2} i}{2} & - \\frac{\\sqrt{2} i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2} i}{2}\\\\\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2} i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}1 & 0\\\\0 & -1\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2} i}{2}\\\\\\frac{\\sqrt{2} i}{2} & \\frac{\\sqrt{2}}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2}}{2}\\\\\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2}}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}1 & 0\\\\0 & - i\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2}}{2}\\\\\\frac{\\sqrt{2} i}{2} & \\frac{\\sqrt{2} i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2} i}{2}\\\\\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2} i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2} i}{2}\\\\\\frac{\\sqrt{2} i}{2} & - \\frac{\\sqrt{2}}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{1}{2} + \\frac{i}{2} & \\frac{1}{2} - \\frac{i}{2}\\\\\\frac{1}{2} - \\frac{i}{2} & \\frac{1}{2} + \\frac{i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2}}{2}\\\\- \\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2}}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}0 & 1\\\\1 & 0\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2}}{2}\\\\- \\frac{\\sqrt{2} i}{2} & \\frac{\\sqrt{2} i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}0 & 1\\\\i & 0\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{1}{2} - \\frac{i}{2} & \\frac{1}{2} + \\frac{i}{2}\\\\- \\frac{1}{2} + \\frac{i}{2} & \\frac{1}{2} + \\frac{i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}0 & i\\\\1 & 0\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & \\frac{\\sqrt{2} i}{2}\\\\- \\frac{\\sqrt{2} i}{2} & - \\frac{\\sqrt{2}}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{1}{2} - \\frac{i}{2} & - \\frac{1}{2} + \\frac{i}{2}\\\\- \\frac{1}{2} + \\frac{i}{2} & - \\frac{1}{2} + \\frac{i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}0 & -1\\\\1 & 0\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2} & - \\frac{\\sqrt{2}}{2}\\\\- \\frac{\\sqrt{2} i}{2} & - \\frac{\\sqrt{2} i}{2}\\end{matrix}\\right], \\ \\left[\\begin{matrix}\\frac{1}{2} - \\frac{i}{2} & \\frac{i \\left(-1 + i\\right)}{2}\\\\- \\frac{1}{2} + \\frac{i}{2} & \\frac{i \\left(-1 + i\\right)}{2}\\end{matrix}\\right]\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡⎡√2 √2 ⎤ ⎡ √2 √2 ⎤ ⎡√2 √2⋅ⅈ ⎤ ⎡ √2 √2⋅\n",
|
|
|
|
|
"⎢⎢── ── ⎥ ⎢ ── ── ⎥ ⎢── ──── ⎥ ⎢ ── ───\n",
|
|
|
|
|
"⎢⎢2 2 ⎥ ⎡1 0⎤ ⎡1 0⎤ ⎢ 2 2 ⎥ ⎢2 2 ⎥ ⎡1 0 ⎤ ⎢ 2 2 \n",
|
|
|
|
|
"⎢⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ \n",
|
|
|
|
|
"⎢⎢√2 -√2 ⎥ ⎣0 ⅈ⎦ ⎣0 1⎦ ⎢√2⋅ⅈ -√2⋅ⅈ ⎥ ⎢√2 -√2⋅ⅈ ⎥ ⎣0 -1⎦ ⎢√2⋅ⅈ √2\n",
|
|
|
|
|
"⎢⎢── ────⎥ ⎢──── ──────⎥ ⎢── ──────⎥ ⎢──── ──\n",
|
|
|
|
|
"⎣⎣2 2 ⎦ ⎣ 2 2 ⎦ ⎣2 2 ⎦ ⎣ 2 2 \n",
|
|
|
|
|
"\n",
|
|
|
|
|
"ⅈ⎤ ⎡√2 -√2 ⎤ ⎡ √2 -√2 ⎤ ⎡√2 -√2⋅ⅈ ⎤ ⎡ √2 -√2⋅ⅈ ⎤ ⎡1 ⅈ 1\n",
|
|
|
|
|
"─⎥ ⎢── ────⎥ ⎢ ── ────⎥ ⎢── ──────⎥ ⎢ ── ──────⎥ ⎢─ + ─ ─\n",
|
|
|
|
|
" ⎥ ⎢2 2 ⎥ ⎡1 0 ⎤ ⎢ 2 2 ⎥ ⎢2 2 ⎥ ⎢ 2 2 ⎥ ⎢2 2 2\n",
|
|
|
|
|
" ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ \n",
|
|
|
|
|
" ⎥ ⎢√2 √2 ⎥ ⎣0 -ⅈ⎦ ⎢√2⋅ⅈ √2⋅ⅈ⎥ ⎢√2 √2⋅ⅈ ⎥ ⎢√2⋅ⅈ -√2 ⎥ ⎢1 ⅈ 1\n",
|
|
|
|
|
" ⎥ ⎢── ── ⎥ ⎢──── ────⎥ ⎢── ──── ⎥ ⎢──── ──── ⎥ ⎢─ - ─ ─\n",
|
|
|
|
|
" ⎦ ⎣2 2 ⎦ ⎣ 2 2 ⎦ ⎣2 2 ⎦ ⎣ 2 2 ⎦ ⎣2 2 2\n",
|
|
|
|
|
"\n",
|
|
|
|
|
" ⅈ⎤ ⎡ √2 √2⎤ ⎡ √2 √2 ⎤ ⎡ 1 ⅈ 1 ⅈ⎤ ⎡\n",
|
|
|
|
|
" - ─⎥ ⎢ ── ──⎥ ⎢ ── ── ⎥ ⎢ ─ - ─ ─ + ─⎥ ⎢\n",
|
|
|
|
|
" 2⎥ ⎢ 2 2 ⎥ ⎡0 1⎤ ⎢ 2 2 ⎥ ⎡0 1⎤ ⎢ 2 2 2 2⎥ ⎡0 ⅈ⎤ ⎢\n",
|
|
|
|
|
" ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢\n",
|
|
|
|
|
" ⅈ⎥ ⎢-√2 √2⎥ ⎣1 0⎦ ⎢-√2⋅ⅈ √2⋅ⅈ⎥ ⎣ⅈ 0⎦ ⎢ 1 ⅈ 1 ⅈ⎥ ⎣1 0⎦ ⎢\n",
|
|
|
|
|
" + ─⎥ ⎢──── ──⎥ ⎢────── ────⎥ ⎢- ─ + ─ ─ + ─⎥ ⎢\n",
|
|
|
|
|
" 2⎦ ⎣ 2 2 ⎦ ⎣ 2 2 ⎦ ⎣ 2 2 2 2⎦ ⎣\n",
|
|
|
|
|
"\n",
|
|
|
|
|
" √2 √2⋅ⅈ⎤ ⎡ 1 ⅈ 1 ⅈ⎤ ⎡ √2 -√2 ⎤ ⎡ 1 ⅈ ⅈ⋅(-1 \n",
|
|
|
|
|
" ── ────⎥ ⎢ ─ - ─ - ─ + ─⎥ ⎢ ── ──── ⎥ ⎢ ─ - ─ ──────\n",
|
|
|
|
|
" 2 2 ⎥ ⎢ 2 2 2 2⎥ ⎡0 -1⎤ ⎢ 2 2 ⎥ ⎢ 2 2 2 \n",
|
|
|
|
|
" ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ ⎥, ⎢ \n",
|
|
|
|
|
"-√2⋅ⅈ -√2 ⎥ ⎢ 1 ⅈ 1 ⅈ⎥ ⎣1 0 ⎦ ⎢-√2⋅ⅈ -√2⋅ⅈ ⎥ ⎢ 1 ⅈ ⅈ⋅(-1 \n",
|
|
|
|
|
"────── ────⎥ ⎢- ─ + ─ - ─ + ─⎥ ⎢────── ──────⎥ ⎢- ─ + ─ ──────\n",
|
|
|
|
|
" 2 2 ⎦ ⎣ 2 2 2 2⎦ ⎣ 2 2 ⎦ ⎣ 2 2 2 \n",
|
|
|
|
|
"\n",
|
|
|
|
|
"+ ⅈ)⎤⎤\n",
|
|
|
|
|
"────⎥⎥\n",
|
|
|
|
|
" ⎥⎥\n",
|
|
|
|
|
" ⎥⎥\n",
|
|
|
|
|
"+ ⅈ)⎥⎥\n",
|
|
|
|
|
"────⎥⎥\n",
|
|
|
|
|
" ⎦⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 1,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"import sympy\n",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"import numpy as np\n",
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"\n",
|
|
|
|
|
"M = sympy.Matrix\n",
|
|
|
|
|
"simplify = sympy.simplify\n",
|
|
|
|
|
"sqrt = sympy.sqrt\n",
|
|
|
|
|
"i = sympy.I\n",
|
|
|
|
|
"_i = i\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"from sympy.physics.quantum import Dagger\n",
|
|
|
|
|
"from sympy.physics.quantum import tensorproduct\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"sympy.init_printing()\n",
|
|
|
|
|
"H = simplify(M([[1, 1], [1, -1]]) / sqrt(2))\n",
|
|
|
|
|
"S = M([[1, 0], [0, i]])\n",
|
|
|
|
|
"C_L = [H, S]\n",
|
|
|
|
|
"basis_vectors = [M([1, 0]), M([0, 1])]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"def getitem(m, i, j):\n",
|
|
|
|
|
" # Note that sympy uses an extremely weird way to store the matrices,\n",
|
|
|
|
|
" # in particular the internal representation is vastly different from what\n",
|
|
|
|
|
" # is printed. For instance the H matrix is stored as a Mul object\n",
|
|
|
|
|
" # (even after using simplify) instead of a matrix. \n",
|
|
|
|
|
" # Therefore m[i][j] will NOT work. This is a workaround.\n",
|
|
|
|
|
" return (Dagger(basis_vectors[i]) * m * basis_vectors[j])[0]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"def is_known(C_L, a):\n",
|
|
|
|
|
" for c in C_L:\n",
|
|
|
|
|
" test = c * Dagger(a)\n",
|
|
|
|
|
" if(simplify(getitem(test, 0, 1)) != 0):\n",
|
|
|
|
|
" continue\n",
|
|
|
|
|
" if(simplify(getitem(test, 0, 0) / getitem(test, 1, 1)) == 1):\n",
|
|
|
|
|
" return True\n",
|
|
|
|
|
" return False\n",
|
|
|
|
|
"for i in range(5):\n",
|
|
|
|
|
" for m in (H, S):\n",
|
|
|
|
|
" for c in C_L:\n",
|
|
|
|
|
" c = simplify(c*m)\n",
|
|
|
|
|
" if(is_known(C_L, c)):\n",
|
|
|
|
|
" continue\n",
|
|
|
|
|
" C_L.append(c)\n",
|
|
|
|
|
"C_L"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 2,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAALQAAABlCAYAAADkr8m4AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAIeklEQVR4Ae1d7Y3cNhTcDVyA08K6geDcQIBsOojdQeIaAuR/kHRgF5A/SQdOOohdgIG7AvIjOaSAXGZuJULeD5JvSZHv0U+AfFpRFGeGI+otpX3ePjw8bCTLdrt9iuNfYn2Fus8ldVseawVnjiajcGnB4wkauYGo7y4I+xtM+2Ium47dT59pbJWLFZw54o3CpRYPnOcWuu3OaQevbp8sCn7GNg9eLnfLD6jwHp/f46TfLPdr27aCM0e3UbhU5PHTGd2+xr5HTy4N/RqNfmTgMxV9lyvQVQF49M0xAAyw3PVo6M+OC/2zK2BZATe05d5z7CcKuKFPJPEdlhVwQ1vuPcd+ooAb+kQS32FZgeUsR5IHvk3ucNArrJyL3uHza/y9xTdPTvmpWazgzBFsFC6teHC+Y36w8syn7XIs5sdoUwAXy3fAxGnnrYcc2nrH8RQp4IYuks8ra1MgO4bGsJ58i4lDfm+CVnDm6DQKl5Y8sg2twaw5JrCC81Pi0rJPPOTIcZYfY0aB7BGajHDr4Cuj30/sOIXH5VtcgfeHTR3/WsGZo9YoXFry4LQd4+MdjLmJrTiGUyPhGH7GynnosE/DthWcOVqNwmVNHjg3p+0g5+FrnsTQrLGfOwLbHKW572bep+HvhEk9zhytRuGyJg+cOxhaGkPzKeGfWLUvVnDm6DgKlyY8RDE0RpTjl6sJ8g77+UsWNYsVnDmCjcKlFQ/pCB36AEE+QxX+SkDtD2UJ1gpOYk0to3BZk8dVhgYgxs78bddzXHmqZjiWprCCc4n50vYoXNbmIQo5KPYMCEbmDxPnEfAen1X9HtEKTmqYWkbh0oKHaISeAHGq7kds33DFNuPov1Od0rLcCs4cTUbh0pKHZNruH3QCp+k+WjE6a5uHTuIEBz4k4nTPO234l3iAL8llebzW7RweOOaqPkG9MG0nDTm+QuU9Vi7MnPTssKnrX+D6PIZourPMPCii5sWE5ikBW/WJyNAAZSLRTIa4ZniMonmrPhHF0ClQXu4K9FbADd27B7z9qgq4oavK6SfrrYAbuncPePtVFXBDV5XTT9ZbATd07x7w9qsqIJq2m5728Mkg53DVJppJKWSJhyWsKd1j5bV4eqKZmMpeZkIBXAx8UuiJZkz0loMUKeAxtEguP1i7Aoyh+abcL1j/i4HFsM4XkqILHtOummgGGPjexR9YJe9fvJgeHz9i18AjKuKiUDvWGv1BuhV4fsBp6OGNx9BUwRfTCuCC8BjadA86+IsKSKfteKtXn2jmItupYLpVmuBhCWtK91h5TZ6SF/xNJJpBzBz9wQGENcPDEtaU7rHyEp6oy5ADpz98zZMYmjXMJ3AheSs8LGGNGTZVVsITdYOhpdN2TZKFAODaiyUelrCW9FsVnqIYGleZiUQzKVUt8bCENaV7rLwWT+kIHTAhiGeooj7RTAB8YcMSD0tYL8idtbuE51WGRoM7IFOfaCalniUelrCmdI+Vl/IUhRwEMjeIW4TqRDMx0azxGEXzFn0iGqEnYdUnmskUzgSPUTRv2SdfojFOY32BUTc1f2sm6Qn48CEQp3NOEslgnyUeSawxrqk+1VJe0ieo+wNWUHnYMOT4FyuX+e/h0/l/TSQ9wajGL6z7icLJi0wgHk1Ec556n70prCmufVBf1WqJt/6aWxTF0BDXRIIWKzjnTij5OwrXWjxEMXSJ8F7XFWihgBu6hcreRjMF3NDNpPaGWijghm6hsrfRTAE3dDOpvaEWCrihW6jsbTRTQDRtNz214mt+nONVm2jGCs4avTwK11o8/EeyNVzl5+iqAC4G/5Fs1x7wxldTwGPo1aT1E/dQIDuGxrB++AViBCUeX66aaCbSdCiygBMY+X5JUcIcEu7NVSOPbENrMGtwbWTDAk5gvAeF4v9SujdXjTw85IhcHF5kT4HsEZrUpluM+gQtVnDWsMsoXGvy4LvDjI93uIWkXvA3kaAFXEzgTOmdUz4K1xIeqBvycmB7IzG0iQQt4GQCZ45hU8eMwrWEB+oGQ0tj6CrJQABg7cUKzho6jMK1Cg9RDI3RwkSiGSs4a7h5FK61eEhH6NAHCOIZqqhPNGMFZxC2YGMUriU8rjI0GtxBd/WJZqzgLPBwqDoK11IeopCD6s0N4hahOtGMFZzBkQUbo3CtwUM0Qk8Nqk/QYgVngYdD1VG41uQhmbZLJj3ByB2dy25Rjt5O4sQxFxPRtMBYq40crrXaKj1PTPMSHqgbpu2kIUdJMhC022aB8NFEMhgNeBHvJzQniWjaoKzWiok+SWme6rNctUSGRqMmEs2kyI/CgzytcGmFUxRDp4zi5a5AbwXc0L17wNuvqoAbuqqcfrLeCrihe/eAt19VATd0VTn9ZL0VcEP37gFvv6oComm76WkOX/PjHK7aRDMphUbhQZ5WuLTC6YlmUu73cvUK4GLxRDPqe8kBXqWAx9BXyeaVtCqQHUNjWO+eaAYYihO0jMKDhtLAJcfYLXFmGxrP4rtnRQKG4gQto/CgkTRwyTF0S5wecuT0iB9jRoHsEZqMplu++kQzKfVH4WGpT1pqLnnBf4gELjDCEDxwK6enTXBZEyfOHV7wpyASQw+RwAWch+AxGdoElzU1x7mDoaUxdJVkIADQexmFB3W0wqUJTlEMjRHBRKKZ1NUyCg/ytMKlFU7pCB28giCfoYr6RDMB8IWNUXiQnhUua+K8ytAAtIN+6hPNXPBw2D0KDxKywmVtnKKQYykcbiGqE80E117YmIW1zsNSn7TSfJ7l4Lfl4/VXdHjIs4FyjsxvsbLOvDLxzNPlcdq3gXcIHtTZCpdaOHGeW6zHPn38TD34OJvvR7zEem65w0G/zwW4wpjA5SSPBY7p/lh8xpjzdxQe5GqFSy2cOA+n6M4u8OGb/wFmtDR28OD1wgAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}\\frac{1}{2} & \\frac{1}{2} & \\frac{1}{2} & \\frac{1}{2}\\\\\\frac{1}{2} & - \\frac{1}{2} & \\frac{1}{2} & - \\frac{1}{2}\\\\\\frac{1}{2} & \\frac{1}{2} & - \\frac{1}{2} & - \\frac{1}{2}\\\\\\frac{1}{2} & - \\frac{1}{2} & - \\frac{1}{2} & \\frac{1}{2}\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡1/2 1/2 1/2 1/2 ⎤\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2 -1/2 1/2 -1/2⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2 1/2 -1/2 -1/2⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣1/2 -1/2 -1/2 1/2 ⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 2,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"H2 = tensorproduct.matrix_tensor_product(H, H)\n",
|
|
|
|
|
"H2"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 3,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"zero_state = M([1, 0, 0, 0])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 4,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAB0AAABlCAYAAACx4uijAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAC/klEQVRoBe1ay03DQBB1EGeUEwWYDkgJhA4QHRBKQFSAQgdJAVxSApRAKAApKYADIArAvOfsrhxn7czYgIOYkSb78ey89dvfZJMky7JEo0mS9KEj6FzTrmi71+v1jqFZhc7gPAhtUSAggalRgd2iwl/GBvuFVrfILwplZpfFMnr7hPITHJ4V6yP5caTuFHV5uyLoBE7XQCINRVXwMy0boqOsykH3yg9/o2ygP8pyJ/QWZ+/Wt8MMTGF0CR1CU5QnSBeYrVxuYlGBuiV1JfZeYdgJvZ2AiunF+OX7ZgVjeTXoz7edOhs+E4NKHW4D5PPdppc9BMU8zq6Zh3D5UC7AwvsqK/sU0+vcjQHAdZqLW6dzFI5clSjR0jsC0LDgmecmNwke7mLRgvItH8XeKwxV9ILa8uHMTixRz4hCLNo3DY4dpYwEBqFSmGkECsAU/jmeA+3MZb9U9LKBBwQYAy2WOYneURbHV6o3dYA8zm4I5gA5rq9QsWjflGuSGwTTIHjLsHZDZU1GC3oCX0Pn7xJgqk3B90MFChBpsO39R1PVmEY9NKg00AakyZv8H3pVS8btSNwIuFYt2N46ozqZSBxTnhB30M+6LmI82wbbzw4nh+F5SIcp9lbV9Y7GHv55K4MmmQXbICIuFmzHeamoVW2DmHkWbFcQGa9W0UsX7qThLP57wfaBI8GncU5WtW2C7UPvmPR+uIJP/bNY2ibYfvEOVWOKcbRg2zMnSjuJHAxUNDZNjVRLxoJtLc2dzF7xmH5DsB0IEYNi3xVdlQfPNZndppcdB8V2s80vYWLRjim/+tvN9o/Rmzv2sS8KdrNdS7V4G3Re2gTboSNa0DbBdjNQC7YDb8KMdhsUuq03M9B6flo+Va1TC7a1bHcye8VjasG2djxpL6aXxhZsgwELtjkVNqTx5uAotb+RbFBaruDM429tMZ0hAgy/v8Emhd5D2cYr/1bSL9nxz8kxfzDLEt4j8Nv1OTQm/K/Rg3+AcXxDnvZrApu1+wjY8ce8qMB2+gXddflBWpRqZgAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}\\frac{1}{2}\\\\\\frac{1}{2}\\\\\\frac{1}{2}\\\\\\frac{1}{2}\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡1/2⎤\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣1/2⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 4,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"plus_state = H2 * zero_state\n",
|
|
|
|
|
"plus_state"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 5,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"I = M([[1, 0], [0, 1]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 6,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"C_L_0 = [tensorproduct.matrix_tensor_product(I, c) for c in C_L]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 7,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"C_L_1 = [tensorproduct.matrix_tensor_product(c, I) for c in C_L]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 8,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAIgAAABkCAYAAABdPHirAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAHJElEQVR4Ae1d7Y0cNwy9Dfw7MGIgBaw7OCcVxOnAlw6S1JASnA6cDoJcCUkFhq+DSwEBYizSwOXxMDQGa2oo3S5JrUQBws5QGpF8fKP50PBud3V1dY36AVUqtw8PDzdSQ8rGQGC3293Dk73kDWK/e7Zq+BXb1Hld/l7v5PaQCLwVvPoesjckXxPkHRiThBDQGlmEmP927B9mFRJ9RpDjfuI+DqZL0h+orzD4QexkIIReZvq/GP4l6lsvQkfq1qD0sI0C/oC6B+BXUkXbc1QixTtUul+h/s+lvhayRecbHpt0oz5eO1lm9RupW/PJyjaM+xMq1D9cfYENtaDjAfUG9Wd0/l094IwdcIaQsUTGWx6W7ME27RNhzUqkbs0pL9uqCKIZa9xOT1F3go73kL0GUDSbWJVI3ZpPLrZdAkFeA6mPAlo0i1ChdqsSqVvzycW2rglSOTt8pSH5lPZI3Zq9nrZ1TRAAxcHn2WKNHc8qVpeYSN1rP6VtN9t6J4gEzrHsxbHAcT9St+bmWWzrnSA8S0hg8FlE70UsSqRuzR8327omyPI4S2BJlxGWmbz9jdStscPTtq4JsgD1J373Amg8g1C7VYnUrfnkYtslEITe4H4joPUKsrvV2SR0OVkUqVsz3sW27gkCAtBi0kc82j0uHhFqy2PeD9j8kfatSqRuzScv255phnA7gkKMpes+n80fILvD/nsYS58KWBaaLX6BPr7UfIv976CX9FuXSN2ab+a20bruNSotwL0E4CY3fJqX2d4XAjgRaf2LPv/YdX+J6Qu6+axJgswX8yaPkyBNcM3XOQkyX8ybPE6CNME1X+ckyHwxb/I4CdIE13ydkyDzxbzJ4yRIE1zzdU6CzBfzJo+r12JoVLyCDUteWvTTsgCtCU2TtEV+bxWPmBDom4lTeCdPNtJ6jXvyEnTSAiGRYrqkLcJ9qwITk5hg3LbEqWXxJiR5CQBNmbSFIG0Wr5jU3oO4JOlsIhLT2LPfLrbVEsQlSSeGA5tae/bbxTaVIJjK+OPgLST5+9CtPhfV1rPfnrapBEFUOfgHIcL8+X0NiYTDuxb17LebbTUEqYnii5pOA/bp2e+z2EYE+XIJHP8ex5FniWM57TOTrZKXJJ1esp79trbtawZZnUHoMXPpLF1GWDbct6w9++1pGxHkv4UA/MvkWf+6JOmsFXay3bPflrb9w/irM8jS0SVJh43q6Ldnv11sqyIIprSw5KVIsvTst5dtLYt15kk6JTLguZ/OFrrfmS1pqwQJy81jkolTDHX+fkIAJ2QmTn1CIzc2Eai6B9kcIRuHRiAJMnR4T3cuCXI6hkOPkAQZOrynO5cEOR3DoUcggtA6Cv0N9q0FoKFBSOc+Q4Be4xMnHv+Y/x6/9EEwr8ySPMvcCNDXao//KCEvMXMTQfU+CaJCNHeHlrWY0MQpvP6dMmlLo6cHLtcwotvEKaxaEkYmCUI09laFXlogpIXCkKStLdsscYG/l5M4tSwcTZe0hSBtFi9cau9BXJJ0CohE6i6Y1IXYBZdagrgk6RRgj9RdMKkLsQsuKkEwlfGHyVuomLxDidS95Wx0mycuKkEABgf/IADDb19rSCQcrooidavGBXZww6WGIDU4nCVJp0aR0CdSt2BON6Kz4FJDEJ4lJM+ZyVaJU5G6JX97kbnhohIEz9t8aZEuIywzSZyK1N0LEyQ7PHFRCbIYaJmkI2GwlkXqXtvR27YLLrUEcUnSKUQgUnfBpC7ELrhUEQRTWljiVKTuLmhQMMILl5bFOvMknQIWJA7TjXcOdKbSvVZE0hb5vlXMccnEqS34J23DSZGJU5PGvtntqnuQ5lHzgGEQSIIME0obR5IgNrgOM2oSZJhQ2jiSBLHBdZhRkyDDhNLGkSSIDa7DjJoEGSaUNo4kQWxwHWbUlrWYaROn8Oo5NGnrVLbB/muMQWtKT/pPXXRwJk4VkqiATUjSFlZrN5O6tHbYTQuMRIrmpC8ck4lTAEEty6JVSNKWapzSAQQ6y3/qqr0HcUnSKfg8q+4CHL7iWoK4JOkUXJ9VdwEOX7FKEEyz/GHylmX8dftWn+a2WXU3A2V4gEoQ6ObgHwQ7+PP7GhIJh6uiWXWrwHh1qCFIjS1nSdKpUST0mVW3AMX5RTXvQXiWkLTzGT5i4lSk3/TOiWblv1BbZucbPL3cSYF6qkwlCD0uwVgaXzKUZWaJUzPqJrAJd/zQR8mhpfYS45KkU0BiVt0FOHzFtQRxSdIpuD6r7gIcvuIqgmC6mzJxKtJvXxqUtan3IKtDzZN0VrqON2fVfYxD0z7u32j2pfvEJyd9ZeJUE+RzdAaxMnFqjlCf7mXVPcjpanKES0UgCXKpkXOyOwniBPSlqkmCXGrknOxeP+beL6+116pv8S6APtjJMigCiPk9XNuX3COC8H+ckvqYrLFIilIWhgB/kC0a8D9jejx3xeetigAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}1 & 0 & 0 & 0\\\\0 & 1 & 0 & 0\\\\0 & 0 & 1 & 0\\\\0 & 0 & 0 & -1\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡1 0 0 0 ⎤\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢0 1 0 0 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢0 0 1 0 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣0 0 0 -1⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 8,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"CZ = M([[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, -1]])\n",
|
|
|
|
|
"CZ"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 9,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAC0AAABlCAYAAADZEWqbAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAD30lEQVR4Ae1cXW7UMBjcVDyj0gcOsEgcgB6B7Q2AB95pj4A4ASo3aA/Qlx4BuAHtAZC2B+ABVhyAZSZ8jtxgb7Lx2CjIllz/xPt5Mpn97HxNttlutwtlaprmEPZeIZ/B9rHStrN1gEmeIW8j+doNHFPSFsadIhM486QEO+sInpbhB57VD6ivvTard732ziaYvcWAW0z4YufA4YPngSEn6Gvt+qAvMOleIAOGJV3Acdk3BCLY1YI+6B+cQ7uCLnWVKtOlmPa9R/Kc+IYvYeQMeYW8RPsC5RregO5UlqSgzWW+laGLGKqajhAj75bJA/od3HlBPu2ylnoWMtAqQGNOaJaaljFNhiARbkffGVt0f0xvcBU2f6qav1LQgHQOgPTTbTI/fYPGE+uSFGp5nALoykPGfTEXGd4cyJIaNFn+IkMXMSSVB6TR37zzJO7QzzsaWVIz3QEzSfBOQ35zmwU0AC8Blno+VnsOsiKVBw06wADLG1G2+SXcoC27/5QybYC5HX1PsAaYuv6OLEtqpumTucCw7BJY7nx315lQUYN+Diwrw8MIk3RRcecpBQ2QqmCNwxcspZoOzpChs4LOQGrQZGU6SEuGzlkyLXV5tiJyIaGvrsEaX2WzlMcsQVPT3IFdIf/yL8G+deg5d7Dmq+FckOkj5NfISaxj39EMZcyRkp7iw8SZBjQFQcpn1S6vBmtiVyNJxwGjNVgTIKXtkmoa3qMGa2JMqzXdzmMbpxqs8VmXMm0MFwnWPLSzcKV/UvvWcwZrHjsw9B4/reFKd2xKmTNY880BUru8GqxxzPZL6RexbzxXu4LOxWzfbmW6z0iuttTl1WDNjstUNb2DHOkhmaYLBGu6E5eBZqCms5q5MktNy5gmuZBIDdbEVKaWRw3WxJiWaroGa2I0o1+t6XYq2zjVYI1PvJRpYzh7sIYnwGeM+E+eJb5Ii5QMGz/MFu112beJfi5AfOvoxu8fqttnMGyrffAKBh/BeDThSpCglQ0g+ElJ6vKGEOCkJMEcqaaHQKuOV9AqJofsVKaHGFIdnyXTRV2erZjJT94UBQ0/zSd7k1+TmqU8Zgl6lDygRe4TPiPvs194act26zRgQ/bkzSjQmHyDmZOe7YcNWTDn/5VHe30Ff0xmya9JjZKHAK8zIXlNqrQ8JMGc0qC5Gia/JlVUHvAgkidvSjPttO1e2pn0mtQ/AW0bp8nBnKLyIM0OMKSS9JqUi3t0cQrYdvVrGE+Khfifh90l8kdkzukygzuHvXH8cRSH4V7JcVxauZ/gb8yEEt8l/BQ6MKUPLDOY89f+BXPcW+IxjsGcYMLYy9/xqBy2/jxphQAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}\\frac{1}{2}\\\\\\frac{1}{2}\\\\\\frac{1}{2}\\\\- \\frac{1}{2}\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡1/2 ⎤\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣-1/2⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 9,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"CZ * plus_state"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 10,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAHgAAABkCAYAAABNcPQyAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGmklEQVR4Ae1c/83bNhC1i/5dBCnQAZwNnGaCphskK7QbFOgI6QbpBkW/EbJBkG+DdIAAaT90Afc9R2foU2nxKJI+UjoCsqTTUffjiZJ1fNB+t9sdsXzAEmp3p9PpdeiAy9rIwH6//whPDiFvgN3+69GB37BN5XH7a7zj201m4E3Aqx8he0X5GOC3QNwBDWSrZREw+33qH0Y1Rf8DeKoX3Edn3tL/xPIcJ38IKlUQwq5cqZ9x+mdY3tzqguzdNgE7YTkgYbvQgmNPsBDUt1j4vKb+k5BuDdlg85Wcm7axnJ89Iqu17tE2fP4JC1Jy2n2FjWiD4gOW11h+hvIf0Q4FFTB66Cwvpjs5Lf3BNvd5wVVra7CtArhaBnUn5r/4+4Dqe8heAgSO5lqte9s9APwS6P0dQJCjmI3Ha7XubTcNsHJ0Pq2B7lpsNw0wgBPwZLSOsZRRXesWvQrbrQM8BvTa9rfXDtxA3rzt1gGWURrCSkYY34trtFXYbhrg4XWI4IVuwyKrUn1bi+2mAR6G5TusD8P2eCUjmMdrte5t9wAwK2jfBxB8Dtn9aKQFVLJFq7B9RBpmS5VI4qWECd1fBv1blipZlpyWKv+B7Dj2rcY2bHRnGz5fSpXj2STIrze8F/Jq5nNPRtMHyO6x/x6J5VRjzcbR+ivsya36BfZ/gF3ar926ts15JY5gTiA8Q8Kq/GGpjYCf/3EGMBA4gjn9u+/hGfzYe99LyoADnJSu/pQd4P4wS/LYAU5KV3/KDnB/mCV57AAnpas/ZQe4P8ySPHaAk9LVn7ID3B9mSR47wEnp6k9ZXYtmaCiBmZHPB/ssq7Im7qR7JkTZVLNJOBfr1dMZnerkc9jkBAdBddL9l1xEc45cXWaTVLfooXhtQj5HwdxJ90CMjbnAKonwrwIYJ7UkgDM2q2YZdxHbWoAtCeBW4NKuZdxFbEcBxu1ZyG1ziX46d7DHY5Zxl7QdBRjgCHi8/0+bUEs1F8G0b+v7lnEXs60BWAOEJQFc418tHcu4VbYJ8DdD9LKeJkNG6VTOfbnSapHPQzZvJbOMO9f2d5Kk6Age/ppTP3QbFtnquFyWcZe0TYD/HdCW9bD7aGVJAH/kyI13LOPOsf1J8hQdwYOiJQFcfLVYW8ZdzPYRmYsS36FjRgDHLetMvIcPTrrf7aKEf+TpUqpMmWwwI4DjvZBXM5/3TrpPJPw78R1XzdoaBoQT39cG6rV4tH+yrvV3eeMZcIAbByjXPQc4N4ON93eAGwco1z0HODeDjfcnwKwj8xuUcwXuxsNw9yYZYJmTmJ4/RnrAmoQ2mRmi3FvfGSAb5PyhVr9F9w1k1HsHOJqivhVSatGmxHeU35x0v/BL99rZJBPiO2eS0Exswy4nODjR0RXpHv5eZpNUt+iheG1CfLe0TWYFlq6/dK8CGFdEERI2zrOkWdpe4m+pPkXi1gJchIS9MHJL2wtdLtKtSNxRgHGLFGLdnNdV3qEtbc8FW/tYybijACMYAe8hEJhUvzQXQaB7VGRpO+pcRYVicWsA1sShImFrTrRAx9L2AneLdVHFrQFYRmnIM7nSahHfLW2H4r2VrFjcUYD5qjBEFboNi6wK8d3S9q2QDNkpGXcU4MGBHBJ2KIYUmaXtFD9L6xaJWwtwMRL2gixY2l7gbrEuxeI+wqWmie/wz0n3X0q2fCSuj/iOoJx0jySgvcCS9KV7J76f87aun6F+7198Xxes4Wi0f7LCvV3afAYc4OYhynPQAc7LX/O9HeDmIcpz0AHOy1/zvR3g5iHKc9ABzstf870d4OYhynPQAc7LX/O9nfiugAilv65J99rZJBPyOSa/CcHmbCNmzhwtIt2jnxPfkYRoG4r2JoR/sjqwZJPutc/gIiTsaEbDClu1Hc5GolQLcBESdqJvor5V2xJ/1joKcEkSdqqnW7Wdmqc5/SjA6CzU2IfAiYTeKezKgEqWaKu2s5I27qwBeKx/bVtFwr7WOVO+VduqtGkAllEaOqGMsDUS3y3jDuV6kSwKcEkSdqqHW7Wdmqc5/SjAQ+ciJOw5R2aObdX2TEr0h7QAFyNh6127aG7V9iUBuRvaUqUZ+RwBbtI2HlGLvnSPfF1KlSmTDWbkczi8SduoA/DuxVfQxV+6d+I7sre2NtTQnfi+NmBD8Wj/ZIX6uqyDDDjAHYCU46IDnJO9Dvo6wB2AlOPi+DXpI/59Tc91R1bBVOj77WQAmLFGcLjmEQGWL76HdKp8XCVkyGWLMyCEwOAJ/gPMO+U3kXQlGwAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}1 & 0 & 0 & 0\\\\0 & 1 & 0 & 0\\\\0 & 0 & 1 & 0\\\\0 & 0 & 0 & 1\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡1 0 0 0⎤\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢0 1 0 0⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢0 0 1 0⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣0 0 0 1⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 10,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"I_4 = M([[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]])\n",
|
|
|
|
|
"I_4"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 11,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {
|
|
|
|
|
"scrolled": true
|
|
|
|
|
},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"from itertools import product\n",
|
|
|
|
|
"starting_states = [ c0 * c1 * ent * plus_state for ent, c0, c1 in product((I_4, CZ), C_L_0, C_L_1)]\n"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 12,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"True"
|
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 12,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"len(starting_states) == 2 * 24**2"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 13,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"resulting_states = [CZ * s for s in starting_states]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 14,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"starting_states_np = [np.array(s).astype(np.cdouble)[:, 0] for s in starting_states]\n",
|
|
|
|
|
"resulting_states_np = [np.array(s).astype(np.cdouble)[:, 0] for s in resulting_states]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 15,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"from collections import defaultdict"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 16,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"state_indices = defaultdict(list)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"for rs_np in resulting_states_np:\n",
|
|
|
|
|
" for j, ss_np in enumerate(starting_states_np):\n",
|
|
|
|
|
" c = ss_np.transpose().conjugate().dot(rs_np)\n",
|
|
|
|
|
" if(c * c.conjugate() == 1):\n",
|
|
|
|
|
" state_indices[tuple(rs_np)].append(i)\n"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 17,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAB0AAAAPCAYAAAAYjcSfAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAB2klEQVQ4EZ2UgU0CQRBFgVCA2oHQAWAFYgcYO4AOtARDCZRgsAPowEgH0oFCB/jeenO5W+6McZKf3f0zOzM3M3ud0+nUQUauVcBdgEFw7sFcPrj/rl0vdrvdg87ADijulTH6oxtsZixr9y1yxPayquPOsnpm/4LNrl+QX6xiBPbgFTxHQPbKDTAp9dpWZcrhKQiCWZUNWOBjKw9nwmIYQXco71X+IpZ1nOuLALbGREN0vo2ABWn1TLiT+siyxqDW0/yMzSznPOd3tQMqWnv/56AtAe1ZbQg5v4NDk31wUV5rniaTC8oQrDCyh41SlNWAZS8LQ+dih971AXwC/VnN1N/4UpteloO9g2CJppFdvqKzbzU9Z/vmPb92Xr3D2ReSWtTaRwxM5KN6MfbwKak4xwofQaHqM4JuBdLT7LFpEydtENOZGS04n5WeQOlNo/uZ0volv96kJj2cboBEm2iYixOav9WwMWAED666DvzSCWhyfKUl2de+iAS1tbxtjg3a5A86SarCMq+/Z8T6b3IdnH8fDZzus5kIfa7THqSnxJqyrjmAe9QAlBMdTuDi8Tcmqx3iEJZ69n55Ob19jPaUbAnMRLGs9usaXVMJLbf8G2gU7t01+LyFT636Br8XM4/ij6J4AAAAAElFTkSuQmCC\n",
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"text/latex": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"$\\displaystyle 576$"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"576"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 17,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"missing = []\n",
|
|
|
|
|
"for rs, rs_np in zip(resulting_states, resulting_states_np):\n",
|
|
|
|
|
" if(not tuple(rs_np) in state_indices\n",
|
|
|
|
|
" or state_indices[tuple(rs_np)] == []):\n",
|
|
|
|
|
" missing.append(rs)\n",
|
|
|
|
|
"len(missing)"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 18,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAACMAAAAOCAYAAACl66WxAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAB+ElEQVQ4EZWVjW3CMBCFk4oBEN0ANqDtBE03oOoGZQOYgW4AK9ANygYVbFA6AYUN0vdZPnRxHAQnnWy/e7bvx7mUdV0XZVkOi6KYS5FH6Z90LtsOoEu0byzbVLxpFyfFtWcRsYPGkXSh/Xvj4cgap0whSGtpZVhulP1H2tib4xkm7lY6cet+PGMIhqylfSPYKOyI2jodZZtJr3ZG3Hf4mXMI/Av8TpNKelT68NLLRot+LKHHKSvlOUVt2C4sXmXLlf1beMX9OMOle3nG4TlJnYTzJv4qR76AETRvMRW7t+rpUDzOCdHzjhrRKALKs8xt6MKIusvm8AGZaUksAw/bvrDAEQ52koPn19/anAcGEbYseJZlK5TJG2zOo17p0g8D4shnfGt5kiM6l/etzCh6SrDRpY3eIZyv4abyuGstegedp5a1Q8OZeCHvJHWE8vD531qecKP2WXlyb8ewfc/8kyMTzUfeEWE4gTA+aU35vPDIhxHni2y8MU/UnK/WzvMmy8wmdF1ZOJS2fO7CzCV8Oa2GaDzZaIxXdWDxKHOriQqj9FvOLDXB2y8pnnvB47FI/D+yoozgMe/rxQjCSDtO7oQ/GM4oGx2bf95nXMP9lT4LCy2E/wWH5jR4LGKasRCN20MwM+NpzqVLW9sojMtp/2QcpewEHM7/BytymT8mTB2JAAAAAElFTkSuQmCC\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle 24.0$"
|
|
|
|
|
],
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"text/plain": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"24.0"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 18,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"len(missing) / 24 "
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 19,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAADEAAABoCAYAAABYHRgjAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGtUlEQVR4Ae2bXY4cNRDHdxHPKFkkDjBIHGCTnIDJDRKQeN/lBkE5QRRusDkAL5sbBE4A7AGQkgMgEVYcgOH/67gsj6d72naXNzy0JY8/ulxVf1f5q9tzutvtTu4inJ6eXkvOkzlZ0ud0jiZ//kle0aMsAJfi+wYF52KLfJi2tKtqIxAAeFzVqIL40wraJlIBeKaGV9ZY5XvKPw/lTUgvBPLWaFrSczXCHGPxGkstieKLFSIPla/ysurepnV5nueKY/qJdHfCIALE74o/KkKchnci+jmtqMkHK8DjtbVTHco8Nr4qYw3kPlDdjdGlqWgYU3nAPZ+ozTARmCU2qog9tjQvAbjNnhXgqYBC94y/8oDgwbnVlaSih49IdyeLZid6caKXxH/w+5dk0iChrxRT//9ez7HWqBXStlP5ZhBSfiumKHClPL0eQyjTs0ddUXR4AWvHg9i4MbPInSST8fRSCkdXVJnZ6Kh76Dlu9EYxulbKYy6vdtGdlB8GNj7ZNCbUjp6k/aAMqeLBWEiV0nMAxJlPeTqySr7o/UCgnAKzy2ANpUetoOdmARS3SJsqi4jeHQTW+DsoFXs47X3LBzqQ70V7XpqqfQTRPLDFJAYJZh14r/iL4g/xwUhGtPcVD/ZQI6TFVZ7bDpRnEXtXLN2J0A1EsEZcmZ30K2LTBELz+zCaiyRUEOFmFeSRtAlEq7Ao1TnjMrCddapm12SJOSlh2+F+ZpiS2wWEhLHwsa8agkCxmLE9+fJDje9vL3e6lOLbRFV2s+x4WaHdQy8QWOE3d20nGHZxJ7nSq0ze4jNDxm+v2MsSUUhwIZczQ2SaZbqCEICN5DEeOD+np7lMjWXFLu6ESgZAyg/vm4JFblV231t1sUQAwLT6AuUDAMYFO1330MsSrAmc8EhjkBXi2hErHTJdQEjZ+w66FbPo4k7F0p0IVxBOHbmYTZcxgVaakezt318qsvFjU+g+vSKrCwgBYFZ6IaWH46rKw0yltMsZ3H1MSFFepfAOKZ63lWe1psza4R7cQUjDp4pjL4d/Vf02WMUVSA8QW2k4tjLb3onnrsEVRGEvn7kiEDNXEOJnClqvp/qadRjkrsEbRIlyn5cQ1dAAgh76SfHfmoYTtNbbY4/NSqwbHuEPMUHvwZ1g/l3IU9ccwlRK+zGXsTqvBe8ryUFv9zEBTz5xbchkwSxx9BNY1qao2GNMcIfj4Yh0vsvdJNYaIWmrcgchJXnT8V7TLS8HhhCm3m9UuAhVrkmXvZM0pNefS3lzq0cqfy2AN67aB2ZdQASXOfrFyBOMuzt5KlfKawVR2lO96VZL9O7hUv5Y4rNAbGlp249N94UpAIh/QsFSe7Y41Tpxrbibi42C/rR23caEFOeszUWVg9sDeZ0p05p2WeyCMk+lbNUNTAFnp3umdlU73S4gpEzrDUy+uHLe4D5iVTgXNTcEqu4bqbfiJa08L157951Urr6BmfPMy+KJu6p64R1AMTkIuRUCQdevqa7uFHyat3z5WDj6NVXttgKLRzxSW95bVYcqdxL3jeIlZsyj6nn/us3r83Kgixd6Q5lBjVsX3UATXXQn5cvvAIqWHkPRA2GqQ4m9sZArT1mBTuO63aCsUjqFthyiIrCxtmmdaNtAGBMx4IWx6w1M8eNm5jOTMZeKdjEIei1aQ/lZK4iGHo/3A5XHIsOMqDS6EvVzAHiuEEE0rdhiwhtuFiTmdQIuNnmS08AFwLGvqbipvUT4VvnqUDWwrZckBWu43MAUHyzBWxIWyaL1SnTREsqXD2wDYKnaMkABUiTY2nmkkhlBLF0nPtoNTIGIYREI9ShjI34RilzvOFMMQoNzmBLuSj91UPHNzGIQNUzvCqjJaZpirfH/JS22RI3CYSNoawhrBOFC1rz9kPX97QJCKq63MWvt1GtMHD0/1Co5R9/FneT7fKNIA6AW/YMrZZbne1kiytEgZ2/GPmvxP7gi0yzTFUTYvbLDXW9jZh1/UOxiiWCBY+eHA0WWVHQZ2FKI4ytnBNIYNODX25ixN7JMF3fKZHQvriC6d3GhgF4De72NWWiASOY+JrRG8BZivY0Zu7gw424Jyd0qjt1As1Mdz12DKwi5Eqv0XDibI6h97gpCwk1B6/VUH7NOCdC03WzeG8SsQBF0uY1ZIriUxnp7jN6sxNdR1+BqieSVzJjLWF3VN+oStK4ggsD1NmZJz+c07pYIbzrW25h5T8+Vu+xiwwCf/IY3p1Tt8xTEW624efvXUqjpC3/OaElZevFZbTPFAxBMeVMHePfpcEqRmXreXU2G/wBnEHMCwSMhRgAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}\\frac{\\sqrt{2}}{2}\\\\0\\\\\\frac{\\sqrt{2} i}{2}\\\\0\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡ √2 ⎤\n",
|
|
|
|
|
"⎢ ── ⎥\n",
|
|
|
|
|
"⎢ 2 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢ 0 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢√2⋅ⅈ⎥\n",
|
|
|
|
|
"⎢────⎥\n",
|
|
|
|
|
"⎢ 2 ⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣ 0 ⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 19,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"missing[0]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 46,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"starting_matrices = [a * b * c for a, b, c in product(C_L_0, C_L_1, (I_4, CZ))]"
|
2019-11-11 13:36:11 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 47,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"resulting_matrices = [CZ * c for c in starting_matrices]"
|
2019-11-11 13:36:11 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 49,
|
2019-11-11 13:36:11 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
2019-11-11 17:37:44 +00:00
|
|
|
|
{
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAGoAAAAVCAYAAACuVXuDAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADtklEQVRoBe2ZjXEaMRCFjScFMEkHuAMSOsAdJKQDuwOnhqQDu4Q4HWCXQDqADmLTAXmfkDQSJ7jjLGlgxjuzSFrt7Vvt6u+OwWazuQhpMBiMJFuFsvd63QikcnAZuiCFO7XHoSx3XRgN+5INcS431hnbG9lc+CH4RKnjq6SftJr++N6gov6xeCkeBuJGtYPes3Q24oXlpYwsxC+hMfXh7L1ldOfiRpLDZ3LWwQI7p82UrRSOcvCErvrIyZbs1kfwF9RDlgz5oxiHCSb75DDUCZ5v1bO6JAbGFuVPcWRTbVbXY4hj9XhmGspL1YWDb5EPJbAO4aiPmJvYmMSoQSJuDjmifrbFZKLC59r01N86eHTEUfLAEL3CIV6JujAYa/FEteGo/0Z8zxjd1jdT40HCU6GpHHnV0t/dZtkSip5nbEXCWFtWUYa64NiczIjDpX4Iyqnd8kjISo4SsBTtJjCl01f2vdKk7YpDbmYf9PNNbA6vviM79jlNDpa0C/aV6izvv86O6viUInOZCHVTSn1l8ostr8YF4hgccnNNor7UcE4YjkjQbwXbrBYFh4sDt8lryfZOGPWTJHR/iLOT9WMtH4ruLj1wOCunnFEELroaZ49CYFCBICEmSYhtYEhQ20zmgvEg/V88V4BuZbvGOX0sDrkZkqiPYh+4AgHoYpJZzHsTK6ZBkpPEJwXyttGZQSD75naVwdRBEz1xTGxIVDWSo7y08m6wj1jdEdnBsfJKJYnJwatAjS2vD45ZSJxRONgIUBStfA3Ow9Q2izMkw18oaCtJvJlfhUmSzKy6jIHF3kR22VpDMmeilXMDfevZ2BeH3LyQKAJnBh96WajOGZMa8FR40UVCASJQk4Q+yYvOEukyU3tt33oO3AhbbSYJL9dst9ENtC/WsTj4YIlJvLrUD7N4YoWlC/PtLgTRwLmqQj4gkjFxzNcJ1d33PkpkHMY+KZIx43g5PrSlSuVowi7sqRBWA8cDbiufVcypMpuXGrz5nLRbqo/goMgMc9/nkN2Fump31SMJXA5gnqFkRXh8td13RfB2OfVNkissHNkJbXatywb+hPiM3Y9V9SxYbTjOX4s3MsGxjbHrPMdSY2BLfHOiuoy9FpZwmNRmEbH1QXzBLnKrMtbr/HCeretAmbOzBhbnObm5MInSADmc977HoHjKZM+OfzV8rIUlHFYT/7ZvL05uqUvIoTZ37XMq5ffBv2hyjqUWFrkQ+618wCAc2Szyx1x0/XX972WdCCgPfCnh1cC/hP8HWgYo7Ny1PQYAAAAASUVORK5CYII=\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left( 1152, \\ 4, \\ 4\\right)$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"(1152, 4, 4)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 49,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"starting_matrices_np = np.array([np.array(m).astype(np.cdouble) for m in starting_matrices])\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"starting_matrices_np.shape"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 50,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
2019-11-11 13:36:11 +00:00
|
|
|
|
{
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAGoAAAAVCAYAAACuVXuDAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADtklEQVRoBe2ZjXEaMRCFjScFMEkHuAMSOsAdJKQDuwOnhqQDu4Q4HWCXQDqADmLTAXmfkDQSJ7jjLGlgxjuzSFrt7Vvt6u+OwWazuQhpMBiMJFuFsvd63QikcnAZuiCFO7XHoSx3XRgN+5INcS431hnbG9lc+CH4RKnjq6SftJr++N6gov6xeCkeBuJGtYPes3Q24oXlpYwsxC+hMfXh7L1ldOfiRpLDZ3LWwQI7p82UrRSOcvCErvrIyZbs1kfwF9RDlgz5oxiHCSb75DDUCZ5v1bO6JAbGFuVPcWRTbVbXY4hj9XhmGspL1YWDb5EPJbAO4aiPmJvYmMSoQSJuDjmifrbFZKLC59r01N86eHTEUfLAEL3CIV6JujAYa/FEteGo/0Z8zxjd1jdT40HCU6GpHHnV0t/dZtkSip5nbEXCWFtWUYa64NiczIjDpX4Iyqnd8kjISo4SsBTtJjCl01f2vdKk7YpDbmYf9PNNbA6vviM79jlNDpa0C/aV6izvv86O6viUInOZCHVTSn1l8ostr8YF4hgccnNNor7UcE4YjkjQbwXbrBYFh4sDt8lryfZOGPWTJHR/iLOT9WMtH4ruLj1wOCunnFEELroaZ49CYFCBICEmSYhtYEhQ20zmgvEg/V88V4BuZbvGOX0sDrkZkqiPYh+4AgHoYpJZzHsTK6ZBkpPEJwXyttGZQSD75naVwdRBEz1xTGxIVDWSo7y08m6wj1jdEdnBsfJKJYnJwatAjS2vD45ZSJxRONgIUBStfA3Ow9Q2izMkw18oaCtJvJlfhUmSzKy6jIHF3kR22VpDMmeilXMDfevZ2BeH3LyQKAJnBh96WajOGZMa8FR40UVCASJQk4Q+yYvOEukyU3tt33oO3AhbbSYJL9dst9ENtC/WsTj4YIlJvLrUD7N4YoWlC/PtLgTRwLmqQj4gkjFxzNcJ1d33PkpkHMY+KZIx43g5PrSlSuVowi7sqRBWA8cDbiufVcypMpuXGrz5nLRbqo/goMgMc9/nkN2Fump31SMJXA5gnqFkRXh8td13RfB2OfVNkissHNkJbXatywb+hPiM3Y9V9SxYbTjOX4s3MsGxjbHrPMdSY2BLfHOiuoy9FpZwmNRmEbH1QXzBLnKrMtbr/HCeretAmbOzBhbnObm5MInSADmc977HoHjKZM+OfzV8rIUlHFYT/7ZvL05uqUvIoTZ37XMq5ffBv2hyjqUWFrkQ+618wCAc2Szyx1x0/XX972WdCCgPfCnh1cC/hP8HWgYo7Ny1PQYAAAAASUVORK5CYII=\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left( 1152, \\ 4, \\ 4\\right)$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"(1152, 4, 4)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 50,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
2019-11-11 13:36:11 +00:00
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"resulting_matrices_np = np.array([np.array(m).astype(np.cdouble) for m in resulting_matrices])\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"resulting_matrices_np.shape"
|
2019-11-11 13:36:11 +00:00
|
|
|
|
]
|
2019-11-11 17:37:44 +00:00
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 51,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"resulting_matrices_np_invers = np.array([m.transpose().conjugate() for m in resulting_matrices_np])"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 52,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAACcAAAAPCAYAAABnXNZuAAAACXBIWXMAAA7EAAAOxAGVKw4bAAABw0lEQVRIDcWWgVECMRBFOccCbigBOkDtADvQFrQDqUE7wBKkBOjAkQ6gBKADfC8kR+5A0BmQnVmS3fxsfjabHMV6vW6dWoqi6BF3msfFV2K38c9z/8G+5JIC7KEztEy+fe0xHONL1F1/RTXmTlx8HXQYVewYdWOBE/2WOxqhggQYdIecPvQozsBIIqNh/xWtxcSW2CgRifPEOacf7MbgSxysBcoxMchBXHPR5vwYw43urIPPrC/FXNG5lPRZeBlrMecwwSjxdy5JThJzMrTKmWX98jozTtpl508EtE6VLjqESHWD6T+Gkd0fL6UXYnoucpL6SFnxiLBntPf4zNheYVxiYgcBADh/Sg4WesIy8Ve4hLdFfCZmua/Zdxw1w4HTf9acj28nZpFuXfD7lE0g9pxGTk6ORcao7+VPkuqwGgdvfZqtipj2OWrulrgLgzekrQ2B6lJoQ+yBppsTw2fdneWde2chb2dTfNdqlwESXoA78JsLsJ0h4cU5MjdkUYu6OiJsL5BSPR8xO34lJuLD6ObHDPt9fSv4MbXhU4LTI7EmLF7T/ymINsgfcNvngH8iTPaYB8RabSKF47Quzdw+mYK9+QYkdqlSSVVlyAAAAABJRU5ErkJggg==\n",
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"text/latex": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"$\\displaystyle 1152$"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"1152"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 52,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"indices = defaultdict(list)\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"for rm_np, rm_np_i in zip(resulting_matrices_np, resulting_matrices_np_invers):\n",
|
|
|
|
|
" for i, sm_np in enumerate(starting_matrices_np):\n",
|
|
|
|
|
" test = sm_np.dot(rm_np_i)\n",
|
|
|
|
|
" check_indices_zero = [(1,2), (1,3), (1,4), (2,3), (2,4), (3,4)]\n",
|
|
|
|
|
" check_indices_zero += [(k,i) for i,k in check_indices_zero]\n",
|
|
|
|
|
" check_indices_zero = [(i - 1, k - 1) for i, k in check_indices_zero]\n",
|
|
|
|
|
" for check in check_indices_zero:\n",
|
|
|
|
|
" if(not np.allclose(test[check], 0)):\n",
|
|
|
|
|
" continue\n",
|
|
|
|
|
"\n",
|
|
|
|
|
" diag = np.diagonal(test)\n",
|
|
|
|
|
" phases = np.angle(diag)\n",
|
|
|
|
|
" phases = phases - phases[0]\n",
|
|
|
|
|
" if(not np.allclose(phases, np.zeros(4))):\n",
|
|
|
|
|
" continue\n",
|
|
|
|
|
" indices[tuple((tuple(r) for r in rm_np))].append(i)\n",
|
|
|
|
|
" \n",
|
|
|
|
|
"len(indices.keys())"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 54,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAK34AAAAVCAYAAAAlrVKCAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAgAElEQVR4Aey9j9E0t3HtLaocAC1H8EkZ0HYGcgZU3QgsZSCVI3DRGcg3gnvlDKgMfMUMpAxkMQN+p/dtLLHYGeA00MD0vA+2ajn/ehq/Pujuwe4+JL/4yU9+8lu8f4O3vH71ww8/fPdpd/9zK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2ApsBbYCW4GtwFZgK7AV2Apco8AXX3zxW4z8Gx39Vz/Fzj/IiR9++OEXeH+nF5ZtAPRVORjOfYn3z8vz+ziOAjI/eP9a5ioO1SbZCmwFtgLxFNj9Mt6c3IEIebPXR3eYqEWMOx8WCb2H2QpsBT6EAnftqXtN+SHScwd5EwXu2kdE3s1+kyTbmFuBrcBWYCuwFQiowJ3XEQHl3Ehbga3AVmCZAnfu33dmlwm+M/9mX1aie6CtwDIF7lzXy0TaA4VXAHl8679XvTP/Zr+mPO6s+zWK7VGtCiDH9t8nWkXb9luBrcBWYCuwFdgKbAU+kAJ3Xi/emV1S7M78m/0DNYkd6lZgK9ClAPrkbb/n3uxdU+5y09beRcbtZCvwoRXwXqd7+2tNzurxWjyW65vdopafrafu8LV8/ebJ76cq52mzczrNsNraf1J1dc1+7rqv1tNSG97ae/urxbJyrBpHz7XN3qPavsdbgTvnobcW299W4CoFdh1epfznPy5ya/ln8M9dVVZT1u5z12t1fFv31YpfM96d5/nO7DLbd+bf7NfU686brXuPAneu15549z1bga3AqwK7B7zqsY/u/e8A7Plbr8DuIes13yNuBbYCW4HICuC5cMv/7kfU59ld9ZQcvTP758AfuU9stq3AVuCTAlGfPez83Jl/s7Oz7Gt3Z919ldjetgJbga3AVmArsBXYCmwFtgJbga3AvRXA57tb/g5wb9Vt9B9hjj5CjDLrO84vfm7L/m29FdgKbAW2AqUCeJZ89v+e60eIUeZ1x1lm932PP8pc3neGNvlWYCuwFdgKbAW2AluBrcBWYCuQK7A/w+Rq7P2twFZgK3CdAj/88MN/4P0LEPwX3j/7Av/4Bu9vcfKP2D5f0rhx8Bs98SW2cvxNaafXuzcY52+4Wfx/p05kX17/iLG+/7Tr+0+MKX+s8IfWGKydB511LLX/DTRKc+SB8eKjxoRrX8NYNDx7fQ+2vz+72Hu+xiQ+cV3y59/w/ive/4C3HP/BO2/h8+WlXIfzgWtSY/ISJik+qaO/yImZrxpTPi5rl9/Tu18bC9d+n/n9Gfb/FTpN6QEyDsZL85KG/T8YL/Whx48puPA7vfhP2P6PHOc2es1t02JKA7F2yb5nizGoWmLtehiO7sF4zTyBjTyvls1diwnXQ/bLpC/4ms9E2KR6WdLHGCbhZ+2KWA97dbIZ3baYcJ3KT9iFWx8pe1p3SI+org1hn/ImyfrSZ9NJry3Ga+ZyPpbaH+ZDxh4i58FD5QPsqN6d6zC6rzoy6+pqPsAPVRujvOX9Lf5eTdXvYX6VDL3HXuy9MfZwYyx6nmGbcqZZh5ltwnLvNxZ2gYB9sydl3M0YU2C9W29+9Uc/E3q55b4J7FRPHWFO93qyw9fSNSXLrnbNXGD9Je1Gt+x4LH/igX3qTenUZf3Gwg7b5ueoFNDIVpnMn8Vwn/TM0+emXm+uN0bY5V6WH3amPrKCfxZ70lRjOJ2jZNez9WaHv2VrYgM7zQSfqc+EWR/0aor7qrXdky/pHtW++fwRe0ZT9WfuX4nHsvVmT2NncaZTs55RbroLKLiXPKN0LPk8QvGrfXNNb7ET296XJW9kDNhX2a3+ermVhdJdmZp1CDu6p45wp3tZrYx2zTjT+CNbI1OzNuDPtAaKxJ5YEEN6zqZTl/VKw/yEzPkkYNointPnfqb71PUNq2liTtszdvUXql4Tc9qesafracvaJXvLVnVq9pDSZ4sJ16fXqzB588Pfsppl2WFH92/YJt1D1CvLrlosqVdL3rD81voo7S3HqlWzZq3ssE+5k3CufMZSOQ/mcPWaxGP0zGym1qswYSxqTc/yq78lNevNnmJMW/g/XQMlm96tNzv8UbXRy1ve582f/MNvmH7DMoH57v1myXc3TM7AxvRbWZYvIXqlhV/1CNMrLeypNtIW917aK3vYcU+IXmNhh+2yXiNzi/Go9YGFC7ZJ96k1682u/sLUq84PlQ9gD7s+AFvz+QMbKs7Ujzy3GPu0t+FaM5dhIzW0JG/KuM/Ye5nO/JXjehyfjcWyw+6ynJH4G/zMZ/WoedNkL+f/TIvSzuP4bCycp55liQH2qbbTKffvPpLjtB1lx/0hcz7Fl7Znccr1TPep64PEkrY1pmQj2zM7nA9Xrzl3jb3Xrryv9/hM09Jfyw7Xw9Qsyw7my2r2TE+cp9eLmeYh6pVlh91l9Sq5gfEP15UsP5tfpZ3HsRc7/ISpV1Z32F1WrzJ3Z9qneWU0zWxC1CzLDu7LatZD9xSnbFv+ctvR/dZYuE7VIWs3ylveP8qP+y+r2SN2nAv7PXeu/Sg77g9Vr1bdW1rk1733R7VPPPBD1Xay99iOsuP+UPUqmliYYNv8TtND5yMfR9onO4YLNlKz5u940hgj21H2cuyav9J29Lg2Fq418wE2y3IeY9Gf70QXje303+Gy+nPQ2pVfY2zO0Si3jkOxQ1P62Qnb1OOnrufZebaw55rivsPPxLlN776Rvdr/4Mu0futlzu9j+fN7ZL+lqV4/re3SX8+xNzv80bXRw5vf482efMNvqtl0asrvC1b+FheuL1sfeLLD19Ka9WSXBIG/UOsDi56wXVavqhX7jKU0tc5lKuieLTsW7Ch21SPU2sbChDhTnwyxtrGw5/OPOKatbYpxkl7p9OFzRXmaz33WLg02ss3mOrnpYocfqv7TIB5bL/bEwvpL9qPb1ni4TvfwzNfUmk0xZ+OlUy95o+zVNb3cCDu6p6aBRrct9tI/7A/7iMa49O/MWuy4TtWhsjfnp9Ri9LjFn/wb7GRumj01+R3ZGpio50Figd/D/ErXPbYe7PCxdD2fx+3BL/7gJ1S/YTVl7XLNPPZbulv1bPnzYE4+WmMx7LAJnfPgoz5nZFqEWB/IHBnYqRjTvPdse+YZ91SfPbhOr996mNM9M9gz36ZnWbrPsrXwq6bmdQvum/KMncEOn0ueURb2NJ+4p5kPmc3UXjOLP8UqW4wRIW+ofABruH4DJgu7ua7zuWL2LTnDsqdxYd+sjWQ7srVwZUyntWjxN8It91rHgn2IZ6yFHcxSh1QuwzblzOn8yNijL1Z3C3vOhPum9Mk0hic/fC1d03uyZ3pMXxPLWCx74kpb3FfNB70+/Xsblh921PdmWXypbtOpl+8/08mRLcuexoB9kymzCdFvWHZw0z01+ezdMrrDhu4hjL9e1qP72PEMdkvWlRILy6S2dA+E3+o64kjHnnMYh2KCXapVug41hin/TfEUqyc/fFFr/zT26JZhh40pl2Gf5inhufd5cWxgD7Wu9GbX+aFiTBMyuvXUHr5M64go7OCmn2ejzOl+UndKT/gy1XVi6N16sicG+FzSa2Q8b37GX4pzdGsZi9E0s6Gfxb0xeLLD19Kcl5g9+UsN4bv6mbG0tx7PYIfPJTXryQ5fodZl4DE9ezLNp9er5FhLews/bJfWrCd
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[ 3, \\ 5, \\ 11, \\ 14, \\ 17, \\ 18, \\ 40, \\ 44, \\ 49, \\ 62, \\ 66, \\ 88, \\ 92, \\ 97, \\ 110, \\ 114, \\ 125, \\ 136, \\ 140, \\ 158, \\ 162, \\ 184, \\ 188, \\ 206, \\ 210, \\ 213, \\ 215, \\ 217, \\ 232, \\ 236, \\ 239, \\ 241, \\ 251, \\ 254, \\ 258, \\ 280, \\ 284, \\ 302, \\ 306, \\ 309, \\ 311, \\ 313, \\ 328, \\ 332, \\ 335, \\ 336, \\ 338, \\ 340, \\ 342, \\ 344, \\ 346, \\ 348, \\ 350, \\ 352, \\ 354, \\ 356, \\ 358, \\ 360, \\ 362, \\ 364, \\ 365, \\ 366, \\ 368, \\ 369, \\ 370, \\ 372, \\ 374, \\ 376, \\ 378, \\ 380, \\ 382, \\ 385, \\ 398, \\ 402, \\ 424, \\ 428, \\ 432, \\ 434, \\ 436, \\ 438, \\ 440, \\ 442, \\ 444, \\ 446, \\ 448, \\ 450, \\ 452, \\ 454, \\ 456, \\ 458, \\ 460, \\ 462, \\ 464, \\ 466, \\ 468, \\ 470, \\ 472, \\ 474, \\ 476, \\ 478, \\ 489, \\ 493, \\ 494, \\ 498, \\ 515, \\ 519, \\ 520, \\ 524, \\ 537, \\ 541, \\ 542, \\ 546, \\ 563, \\ 567, \\ 568, \\ 572, \\ 585, \\ 589, \\ 590, \\ 594, \\ 611, \\ 615, \\ 616, \\ 620, \\ 638, \\ 642, \\ 653, \\ 661, \\ 664, \\ 667, \\ 668, \\ 677, \\ 686, \\ 687, \\ 690, \\ 699, \\ 712, \\ 716, \\ 734, \\ 738, \\ 760, \\ 764, \\ 782, \\ 783, \\ 786, \\ 808, \\ 812, \\ 830, \\ 834, \\ 837, \\ 839, \\ 841, \\ 856, \\ 860, \\ 863, \\ 878, \\ 882, \\ 891, \\ 904, \\ 908, \\ 926, \\ 930, \\ 933, \\ 935, \\ 937, \\ 952, \\ 956, \\ 959, \\ 960, \\ 962, \\ 964, \\ 966, \\ 968, \\ 970, \\ 972, \\ 974, \\ 976, \\ 978, \\ 980, \\ 982, \\ 984, \\ 986, \\ 988, \\ 990, \\ 992, \\ 994, \\ 996, \\ 998, \\ 1000, \\ 1002, \\ 1004, \\ 1006, \\ 1022, \\ 1026, \\ 1035, \\ 1048, \\ 1051, \\ 1052, \\ 1056, \\ 1058, \\ 1060, \\ 1062, \\ 1064, \\ 1066, \\ 1068, \\ 1070, \\ 1072, \\ 1074, \\ 1076, \\ 1078, \\ 1080, \\ 1082, \\ 1084, \\ 1086, \\ 1088, \\ 1090, \\ 1092, \\ 1094, \\ 1096, \\ 1098, \\ 1100, \\ 1102, \\ 1113, \\ 1117, \\ 1118, \\ 1122, \\ 1139, \\ 1143, \\ 1144, \\ 1148\\right]$"
|
|
|
|
|
],
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"text/plain": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"[3, 5, 11, 14, 17, 18, 40, 44, 49, 62, 66, 88, 92, 97, 110, 114, 125, 136, 140\n",
|
|
|
|
|
", 158, 162, 184, 188, 206, 210, 213, 215, 217, 232, 236, 239, 241, 251, 254, 2\n",
|
|
|
|
|
"58, 280, 284, 302, 306, 309, 311, 313, 328, 332, 335, 336, 338, 340, 342, 344,\n",
|
|
|
|
|
" 346, 348, 350, 352, 354, 356, 358, 360, 362, 364, 365, 366, 368, 369, 370, 37\n",
|
|
|
|
|
"2, 374, 376, 378, 380, 382, 385, 398, 402, 424, 428, 432, 434, 436, 438, 440, \n",
|
|
|
|
|
"442, 444, 446, 448, 450, 452, 454, 456, 458, 460, 462, 464, 466, 468, 470, 472\n",
|
|
|
|
|
", 474, 476, 478, 489, 493, 494, 498, 515, 519, 520, 524, 537, 541, 542, 546, 5\n",
|
|
|
|
|
"63, 567, 568, 572, 585, 589, 590, 594, 611, 615, 616, 620, 638, 642, 653, 661,\n",
|
|
|
|
|
" 664, 667, 668, 677, 686, 687, 690, 699, 712, 716, 734, 738, 760, 764, 782, 78\n",
|
|
|
|
|
"3, 786, 808, 812, 830, 834, 837, 839, 841, 856, 860, 863, 878, 882, 891, 904, \n",
|
|
|
|
|
"908, 926, 930, 933, 935, 937, 952, 956, 959, 960, 962, 964, 966, 968, 970, 972\n",
|
|
|
|
|
", 974, 976, 978, 980, 982, 984, 986, 988, 990, 992, 994, 996, 998, 1000, 1002,\n",
|
|
|
|
|
" 1004, 1006, 1022, 1026, 1035, 1048, 1051, 1052, 1056, 1058, 1060, 1062, 1064,\n",
|
|
|
|
|
" 1066, 1068, 1070, 1072, 1074, 1076, 1078, 1080, 1082, 1084, 1086, 1088, 1090,\n",
|
|
|
|
|
" 1092, 1094, 1096, 1098, 1100, 1102, 1113, 1117, 1118, 1122, 1139, 1143, 1144,\n",
|
|
|
|
|
" 1148]"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 54,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"indices[tuple(((tuple(r) for r in resulting_matrices_np[0])))]"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 56,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAK34AAAAVCAYAAAAlrVKCAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAgAElEQVR4Aey9i7XttpEtKnl0AGp1BM/OwO3OQJ2B/G4EbWdgj46gh24G6gzukzOQMuhrZWBnYPtkoDdrH2AJJPGZBRSK4NngGEskwULVrFkfYHFtSZ9/9tlnf8Dn9/jI8duffvrpx4+X+5+bgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZmAzsBnYDGwGNgObgc3AZuAeBj7//PM/wPLvg/Xf/gIX/yIDP/3006/w+TE82KfFGEDgfn2GhLEv8PnleXzfbwZqDOxcqrGz7jOpdXx+h88X66LcyDYDazGw+51/PHavmsP5zuU5vFpq3TGyZJPTtfsNx9OW2gxsBsYZ2P1mnMNPTcOT130WOyvnGdsVMXn6/wRbO0b+UQLn+32hP+3b4mZgM/AwBnavfFjANlwTBp66L9v1ahL+i5Kn5oM48mTsGvwr+slggsyS30cY7JdC2QObgXfAwJNrY9V+8w7SZjkXn54LT8avwQ7Zpf7OmcXOynkWxoqYPP3ftjYDm4E1GVitz6/J0kZ1ZuCpebPX4nMk9/17YOCp9Sqx2TVrn6EPz4elvptqo8Nyz8pp7Y/IM5ggs39fGCF5z30kA0xtrOjYqvW6Ilcspidz+l6wr1ivLPesHJuvFnIrYrLwa+tYh4EVa3Yddn5G8nSenox/Y78nD5/MuzD2ZPwb+8ecBw/u7z429z/3G8+rT533O3KZjd+nzj3Lg7fc5v1nxp/Mxc9e7CtPBlbuqZ48bFufBgO7B64fxx2j9WPkhXCvP15M99nZ8enj7Wmznhxna+zW+lq5MMMedLr8TdoM7C2+rJ4/GbtwMAP/zpt2ds3gvW3VRsIau7U+Gy+3ls3AXAa8+uRcL7b2zcBHBph83r1+Z0tkYOdCZGKfZzHA9KRZtt+z3l3b7zn62/eVGHhyD2Swr9xrGPwr5crG8iwGnppfK9fsszLgI9qn8/lk/BrskHX5TY3NYRY7K8fatZBbEZOFX1vHOgysVq8aZnZ9aNjasp8KA0+t2V2vn0oGbj+ewsBTe4WG391XNGytLbtjuXZ83hu6d9I/l3pnMyvHmFh+Cv2H8XMWx1vvZmAmA596bn8K/Wdm/FfR/R7i9B58lHzS+AnZx+6VWD9ZuVVqMcXxZOypH/t6M5Bj4Mn9J+dPbmzXcI6VPbYZWIOBT70H7f6zRp5tFJuBzcBmYDOwGdgM+DPwqe/zWEZ/+umn/43PryD/J3y+/Bz/+Aaf7zH4A86vQzaOuPljGPgNzn+Xe8j9GMamnEKgfg87v59iAEqDje9w+a+w8yFnBzLCixx/w0cI+wayf5WBGUcLE57/A3a/wCfyL9dyFH34+Lj/ny1MUXPCVRz6P7PyhMEEmW8jEJy/xOc/SnFO5IYuA65i3jK4hwBkJpcwYdw9lyK8CibJ5//ER+rtX/CR++8Qt0NfwpjZASyxxqPOQ97iufTA2IcEj9xLH7gFE/B8DfvSt0rHB2D759JD7Tjs0TGBrFvNMbYg4xo7BlPkH7LVvItyFmfYkh/5qmud2LGWWwk7fHPvdy0+8bxrfxX0FtcZD94D9mpfhIxrr0r9Dhw1cz7OaXGq1Rf19pxbtvCczmXIuvVk8bWFPfIBuWr/w3PX3p3gqvZK4KLWQ8jRMYq2R8+waYJdcECXS97ADt0DIRtzJvt9DM9d+40Ge+C0Gp8Yf+il5KJ879kav1ZfL26Zp7Wl5TTIT1ljZ2GH3lgfkdrDd5k4OHJmsQe56v5AcLByI5jTuZb4oWvJfsNyCjlqLUv5671meT/rxzzphdk6xLOY79n14Kxr5J7FDzl63ffCb41d4+MI5zLXGnvQudTeBj6q6hDyS+0PWPyQo2vDKG+a60/Ih2YfAXa37yMaW5CtYsfzO9YoS95VteGZNyF3mrXYitEo5jhfkzcM9qDP5fc/FrsGE2Rd+nzgkuoPSj+pOorx7z3PxgT9xT1cL+Y4zxI7dC3ZK1kfhRPILpfzMVbxDIzVfAjP6Xe3Ua/mrOE01VvC3qsv1c1e99qqYF9yX5bykcOOMdd6FTzW3AedLjXLYoecKh8g39wDpbHsubbGzurrwZqbw9qDHMU9K5fDoh2zxA5drjVriT3yBp1L1WuCK34XjEPZ97/AP71eBQDLfQQL+Sp+rb6ot+estdXCfsYA+eoe6CyvubfGDn1UT9JgrMla4oeuJftN9B/4Wjnv9v7DkvfEv2V6ZU8uYM4yvVKDXxvLGK+eM2NLid2t31hi1/jYw/N5jiX2qBs6XepV7DH4g5yqB0Lv9Jq1xs7qi3EaObO2IKfiPWLCvGl7m5AP1LvWIEvlMzBTctFHq3ONKzyLe4Pi79uQobmwwhz1lLAHTM3fDSDn1ucj5ngexS56oOOWnAm2szUWuG/+bsDKRb4szw3uGeyr5k0VO/x2/S5yjtko7yHvlst51s8oF3iY+vtCtBXPJe7j83guyWF8uT4fMcdzBftt9SrYSrgi7njOyWFsyZqNmOM5hz15dkvNljBhXJUPQc8S9cpih9xt9SpxD5xd/m5SgV8Vo5hrFucR7Ji7ZL2yvIfY3VKvwXZ2XxnjCj/id5I4dOtvDBGEBXb4dlvNwvYQ75r8SjmzuB7Bjrm31qtF3gQdt9RsjvseToOe29dYDXbILlWvSuy3ra+lnGfxs3IWvSWnA/YvvVKLCfLL1GuIB/WuFbgpuRxvFmM57qNehlPISM023w1GnZZnA+y31ewoduGRiY8F37Cj4in4Rq89NS5WxA+8bjXLcg85VR0GzukY9cTBGjurrwfreQ5rC3KqPQvkqe9bZzzaexb/WS/mXdbiKJNgL/5+FWVHztbYoU8Vo1WwA7f7dykt95Cv5nPg3mV/YI095kHLxyg3cp6BHTpd9sQMdsjQuQxZt3qVmDH4g1zXug/9xZ46kjMzsANrl489frC8Bz/pXA58L7G3YbFruOjhujQHdqv9O85jOE10Td0fJJjMsEedcg6+Xn5rSWVGrxOuoir1e3fooHtqNGJ1buHH8+W+j0TfSezVv/1IdNF9Kc4ZObewn3VDPrv2hPhQPp519t63sOM59Z4hYHfZU6a+tvBHWUYukVmtV1Z7KnBTMYpcWJwTrqK6Uq+sYo+TWX1RfvTM2mPkILNMvwEW9frD+DjKd5zfsoXnqr1uS1+0a3Fu2WKwQ0YdHwvsoqOFP8jQuQx9so5N3dMLJjmssEOPKr8+Wtf/E3bUcW7xiecu38Gtsffo0zP+8wyNPZZTyLmssTOwCzPQS9f1z0zqrjTYo2bMYfcH03vNDPzQuWLeUD0Q2FXfF2NMtWcN75DVYJ/+Xcoau0aflueSPMupzIdsrNfqdyTITe83AQ+VD9F34Kr2ETx3WWM12LWYWj5GLkbOsEHxzmKHnEufjD5b40/0xvqIQ9nvxfFhz9kaO6uvB+t5Tq8tzMu+MxP9eBY5r/akM5aeexY/5Jr5DBn1d4QezHEOiz2Rj7zGoWwuQ2+1p8bJI2dr7NDnsrcRnxnskFHlAuRd1lcWv1JuuTU24Kc4BfexLjz6jRkm4G72JOHB8mDzlOWU1WfhA2MLMs0+AhlVbTtjp78jwY/pfV58V/DexM7Ex4LvVIcl/lSvXIcYTPsdfAZ26Iz9MrqT3UfEh71nEjvVA6H
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[ 2, \\ 4, \\ 10, \\ 15, \\ 16, \\ 19, \\ 41, \\ 45, \\ 48, \\ 63, \\ 67, \\ 89, \\ 93, \\ 96, \\ 111, \\ 115, \\ 124, \\ 137, \\ 141, \\ 159, \\ 163, \\ 185, \\ 189, \\ 207, \\ 211, \\ 212, \\ 214, \\ 216, \\ 233, \\ 237, \\ 238, \\ 240, \\ 250, \\ 255, \\ 259, \\ 281, \\ 285, \\ 303, \\ 307, \\ 308, \\ 310, \\ 312, \\ 329, \\ 333, \\ 334, \\ 337, \\ 339, \\ 341, \\ 343, \\ 345, \\ 347, \\ 349, \\ 351, \\ 353, \\ 355, \\ 357, \\ 359, \\ 361, \\ 363, \\ 364, \\ 365, \\ 367, \\ 368, \\ 369, \\ 371, \\ 373, \\ 375, \\ 377, \\ 379, \\ 381, \\ 383, \\ 384, \\ 399, \\ 403, \\ 425, \\ 429, \\ 433, \\ 435, \\ 437, \\ 439, \\ 441, \\ 443, \\ 445, \\ 447, \\ 449, \\ 451, \\ 453, \\ 455, \\ 457, \\ 459, \\ 461, \\ 463, \\ 465, \\ 467, \\ 469, \\ 471, \\ 473, \\ 475, \\ 477, \\ 479, \\ 488, \\ 492, \\ 495, \\ 499, \\ 514, \\ 518, \\ 521, \\ 525, \\ 536, \\ 540, \\ 543, \\ 547, \\ 562, \\ 566, \\ 569, \\ 573, \\ 584, \\ 588, \\ 591, \\ 595, \\ 610, \\ 614, \\ 617, \\ 621, \\ 639, \\ 643, \\ 652, \\ 660, \\ 665, \\ 666, \\ 669, \\ 676, \\ 686, \\ 687, \\ 691, \\ 698, \\ 713, \\ 717, \\ 735, \\ 739, \\ 761, \\ 765, \\ 782, \\ 783, \\ 787, \\ 809, \\ 813, \\ 831, \\ 835, \\ 836, \\ 838, \\ 840, \\ 857, \\ 861, \\ 862, \\ 879, \\ 883, \\ 890, \\ 905, \\ 909, \\ 927, \\ 931, \\ 932, \\ 934, \\ 936, \\ 953, \\ 957, \\ 958, \\ 961, \\ 963, \\ 965, \\ 967, \\ 969, \\ 971, \\ 973, \\ 975, \\ 977, \\ 979, \\ 981, \\ 983, \\ 985, \\ 987, \\ 989, \\ 991, \\ 993, \\ 995, \\ 997, \\ 999, \\ 1001, \\ 1003, \\ 1005, \\ 1007, \\ 1023, \\ 1027, \\ 1034, \\ 1049, \\ 1050, \\ 1053, \\ 1057, \\ 1059, \\ 1061, \\ 1063, \\ 1065, \\ 1067, \\ 1069, \\ 1071, \\ 1073, \\ 1075, \\ 1077, \\ 1079, \\ 1081, \\ 1083, \\ 1085, \\ 1087, \\ 1089, \\ 1091, \\ 1093, \\ 1095, \\ 1097, \\ 1099, \\ 1101, \\ 1103, \\ 1112, \\ 1116, \\ 1119, \\ 1123, \\ 1138, \\ 1142, \\ 1145, \\ 1149\\right]$"
|
|
|
|
|
],
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"text/plain": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"[2, 4, 10, 15, 16, 19, 41, 45, 48, 63, 67, 89, 93, 96, 111, 115, 124, 137, 141\n",
|
|
|
|
|
", 159, 163, 185, 189, 207, 211, 212, 214, 216, 233, 237, 238, 240, 250, 255, 2\n",
|
|
|
|
|
"59, 281, 285, 303, 307, 308, 310, 312, 329, 333, 334, 337, 339, 341, 343, 345,\n",
|
|
|
|
|
" 347, 349, 351, 353, 355, 357, 359, 361, 363, 364, 365, 367, 368, 369, 371, 37\n",
|
|
|
|
|
"3, 375, 377, 379, 381, 383, 384, 399, 403, 425, 429, 433, 435, 437, 439, 441, \n",
|
|
|
|
|
"443, 445, 447, 449, 451, 453, 455, 457, 459, 461, 463, 465, 467, 469, 471, 473\n",
|
|
|
|
|
", 475, 477, 479, 488, 492, 495, 499, 514, 518, 521, 525, 536, 540, 543, 547, 5\n",
|
|
|
|
|
"62, 566, 569, 573, 584, 588, 591, 595, 610, 614, 617, 621, 639, 643, 652, 660,\n",
|
|
|
|
|
" 665, 666, 669, 676, 686, 687, 691, 698, 713, 717, 735, 739, 761, 765, 782, 78\n",
|
|
|
|
|
"3, 787, 809, 813, 831, 835, 836, 838, 840, 857, 861, 862, 879, 883, 890, 905, \n",
|
|
|
|
|
"909, 927, 931, 932, 934, 936, 953, 957, 958, 961, 963, 965, 967, 969, 971, 973\n",
|
|
|
|
|
", 975, 977, 979, 981, 983, 985, 987, 989, 991, 993, 995, 997, 999, 1001, 1003,\n",
|
|
|
|
|
" 1005, 1007, 1023, 1027, 1034, 1049, 1050, 1053, 1057, 1059, 1061, 1063, 1065,\n",
|
|
|
|
|
" 1067, 1069, 1071, 1073, 1075, 1077, 1079, 1081, 1083, 1085, 1087, 1089, 1091,\n",
|
|
|
|
|
" 1093, 1095, 1097, 1099, 1101, 1103, 1112, 1116, 1119, 1123, 1138, 1142, 1145,\n",
|
|
|
|
|
" 1149]"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 56,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"indices[tuple(((tuple(r) for r in resulting_matrices_np[1])))]"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"execution_count": 55,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAMQAAABlCAYAAAAS5kwgAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAIa0lEQVR4Ae1d7ZHjRBC1qQvgSMGXALWbAFUsGdxdBrAxUMV/6sjgLgD+QAZABtwFQNVuAPyALQLAvPZKLtuye0ZWS+puPVVpbaulmX6vuzUfGvWut9vtynJbr9cvUd5b7Pco+9ay7KnLyoKFOOo95wXIusHpHy9c8guc+s0FWedwU9ZdI5DACLtlwUIcxy4IPh5wZHN89PkXfH394kDwI77LyYfb4+GP0ncU+AnnfEKlr0vnepdnwUIcHU971zmyWn2NYzufPQyI9yCvVwCcKZiHyIBrBuDjH04VxA1cDu0C4rNTIX+TgSUzwIBYsvWJvcMAA6JDCQ8smQEGxJKtT+wdBhgQHUp4YMkMHM4yDeYBo/UNCrnHLs8iNvj9Hp8PGNnLlG6oLQsW4ujndjLf1D6Ye8Vp137k8ewcDOCm8S2QyGOHNbtMOWxKFEYMMCCMiGQxORgwG0Og2SmuEpQmKQJtWbAQR39vMwuIKM5eQ1EWLMRRY+3jc9hlOuaDvxbOgFkLITyiiZYl3981nMoUrGzf4E719Pw1zt8sWIijn8+ZBgSqfgfnl+cQuw3GkOcQ8q7Fq+cjof5mwUIcPd1OnkPIgHgDZ14N2Zty7toypMzm2E17LMpnFizEUfZpcCTPIeCa25X1GEJahz+wZ9iyYCGOHt5o2mVChJ2+fCHGeMRxeZMu1JYFC3H0czvrFmJfO8YP0hWTt5BCJxoQQFmwEIdYU99GCQgQL2MHeXf1FneocDNMh5RlwUIch1a9/N20yyTVtMQjEOTF7fbu+oTf4d7XzoKFOMQT6zbTFqIhXqZaf8D3G9nxXcYRf9ep4+esLFiIo59PmS7/Bvn/oPpOPia0DiHWMB1SV4MF5whW10nZanAc4vb6vQbHtfbAdfvl39Zdpq9A6F1DqmTui/hAbqc+dP9ccw6QKK1fi7VzE9CunViWwiZT2cM0IKB0mkRlJaeNgjWKniW+S3IrnKZjiJLSlJMB7wwwILxbiPpNygADYlK6WZl3BhgQ3i1E/SZlgAExKd2szDsDDAjvFqJ+kzJgOu2KufkNtJcn0zI/HzpRWckKUbBG0bPEd0luhdP0SXVJacrJgEcGEEz7J9XsMnm0EHWajQEGxGzUs2KPDMgYQlai/oT9vyEKotlxn6gMOsqao9+x91l79KZZFrCnJwJWUda7no7s8SfokhhYcQwhLHBbNAMITI4hFu0BBH+RAetpV+mKpEhUdpGxRtA09+6xRtGzxHdJboXTNCCgdJakWCX+RR4FaxQ9azjXzjHDyURlVyRog2VkEsF9UrYoemLiYrYkeeBIxhBQgYnKwMPVW5QEYFH0vNoQzYUmOE27TIiwNInKStaJgjWKniW+S3IrnKM9mMMgJ02ispIxomCNomeJ75J8CM5RAgIKbaB0ikRlFeSHwLoUmwzFadplEudpFUITFj5RWW0weMe6FJtY4DRtIRqFUiQqqwwG91iXYhNLnF/C+DKF+AXudEOnviRRmZR1tA8td67rgUMeNMqU3MdTHXAsBNYaPXHORZynuL3+rsF5SXdc+z12iLcr6TL9i1229vP513V/UyTFEui448ikwF1DQ2cxIMhTE5k1183+UdKzhHN2APUKDPG9v9pqTMcQID9NorJMWFpjn/vMgtMKh+kY4hzhPEYGIjHAgIhkLeo6OgMMiNEpZgWRGGBARLIWdR2dAQbE6BSzgkgMMCAiWYu6js6A6bRr87QwRaKyTFg0L8qC0woHkwxo3kLZIhhAMDHJwCIsTZC9GeAYojdlvCAzA2ZjCDQ7sqBP3fB4PcR/I/WOBfrJ2qrwCdc84jALiCjOrkZsI/SOBfo9QdXbGizaOXPj9IiDXSbNYyhbHANmLYQw1zSB7pN31Vg5ExYNbxacVjhMAwLEmyWL0ow4kSwTFo2yLDjNcMiLMDIg3qBPN/SNOSnHffKuGpwNJymwaHiz4ByCA9cyURlIKG0mia9KlTiQZ8FpgsO0y4Q7UZpEZZmwaEGXBacVjtFmmTDIka7Ya+yDpwc1g04hy4RF4ysLziE4RgkIKLQB8SkSlWXCUgiGFDYbai/TLpMQ3iqEJix8orJMWGqCIbrNLOxl2kI0CrlP3qU5RyvLhKXFdO4zC04rHKbLv6GUJO86l8Nov4YJ54j8LfZ73JHcji9qsJxzMI/HNM6z4ByCA9ful3+bdpng4GryLlSsJv9y5kxDEl+5gVLivGQzN0DKipjYyzQgSjqD/DCJzCLpqvGeBYeGUWRWOE3HECWlKScD3hlgQHi3EPWblAEGxKR0szLvDDAgvFuI+k3KAANiUrpZmXcGGBDeLUT9JmVg0mlXzIlvgC5EIrNIumoekwWHhlFkVjhNn1SXlKacDHhkAMG0f1LNLpNHC1Gn2RhgQMxGPSv2yEDVGAJNyuDEWCgjTCKzSLpqTuUdh4VfCX5LnFUBgXUigxNjoYz9ilfNiB5kkXTV+PKOw8KvBL8lTnaZNI+ibHEMVLUQVqw0TWSIRGaRdNXskwWHhlFkVjgnDQjobZZMqkSQgTySrhrcLDg0jCIzw2mWqAx9OTXRGZQOk8gskq4a71lwaBhFNgQnrh0tUZkopm0myaS0CgxlkXTVYGfBoWEUmQnOSbtMiOQwicwi6ap5ShYcGkaRWeGcbZYJg6Awicwi6ao5ThYcGkaRDcE5S0BA4Q30DpHILJKumqNkwaFhFNlQnJN2mQ4VRhPnPpFZS24EXTVHyYJDw2jpW+0skwzVT/ef4QzqzFEfOcqXluFX7FJnu0tis5d9ypniXOgURleNjyw4NIwiq8WJ8x6wn/r57reUI8sp2sRh+NrZHnHSb52jVx7AnaqYyOzKos0vi6SrBj4LDg2jyGpx4jyZYj27wdc//A8UdofKoG6y5wAAAABJRU5ErkJggg==\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left[\\begin{matrix}\\frac{1}{2} & \\frac{1}{2} & \\frac{1}{2} & \\frac{1}{2}\\\\\\frac{1}{2} & - \\frac{1}{2} & \\frac{1}{2} & - \\frac{1}{2}\\\\\\frac{1}{2} & \\frac{1}{2} & - \\frac{1}{2} & - \\frac{1}{2}\\\\- \\frac{1}{2} & \\frac{1}{2} & \\frac{1}{2} & - \\frac{1}{2}\\end{matrix}\\right]$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"⎡1/2 1/2 1/2 1/2 ⎤\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2 -1/2 1/2 -1/2⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎢1/2 1/2 -1/2 -1/2⎥\n",
|
|
|
|
|
"⎢ ⎥\n",
|
|
|
|
|
"⎣-1/2 1/2 1/2 -1/2⎦"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 55,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"resulting_matrices[0]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 60,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAGAAAAAVCAYAAAC5d+tKAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADDElEQVRYCe2ZjVHjMBCFzxkKyNx1YDqAowPoAEIHXAdQAyVwLXAdhBZCB6EDfjrIvU+RzFqW73IeyzfKZGc2lp5+vHqydtdOtdlsvlipqqoW9mKxQ3kcBlLczuzU6nCr+onFDuVRGag9x82kzQao4VLoNz39v5pWFdg16YPXla5L6V83iT6MsXONVda8Rdokbp/gQPbD9Va8C5qrtqJsVVgtfYywe2H4rXOLx2W1r+OxcZ8hdc1ZvE1aw0o6Z/2OcFV4Um9iQoQ9ho62Tdg7ajFbVhuuLNcGFG+TuLmRPsBZcEELVX4KjOVcwLuODCfECkdpjiuwIGVhuKcPr0BjS/E2ea4X8DrTDwvqy3og+kUDIDQl8cbQ57pnM1Pjh2D7YhOcL470cyV1wSFmQ0TSlhIXhNX+bBu1mbieLIE33GcfbPJrgfMLXNB3Kf56J/EuBtdzZwcIB/sQQX2nyXYftVySTWbhcO5cOIVLEdfKgPrq6kt/F0BsH2H3UZ1I38qgbPuY5ZJsCuuWzaSia07AV2mfj1fTp+hJw708aZIfn6gLvC6qW2yqckk2RZzgKWo2YCfRQiGZUxKTzzEip/0frqcYmxIkuwefIAxxqWymGSPyOS7Hlnxhzn8J53qmOvm5FQI1b6zgZFKtmGE7DimXZlPiAYXzN9a+lN6qQzIGqA0iW/6dvoyRure51Fi18bLWiQF/GpOaJ4VpjuJt0ho4vctKP3xa4AtoJ+X0TzkbFKepHJ8TjTnWNSkayy4RLy5CB2HsOhvzLPw04P9y3RebtA7iqcs+eRFbi5DOCRBOJgORKe18O2IOCRPbca0TpjZuivaenpQtAYvmju0qySY4qB3pnhCe6M4m5MB0P2LKoA3IYQ9zTmmT7kXcdA/9jDtLcEOt7Mah+X7OtOidUt98JnRmntImEhI4336MExl8iCNjCZkNbVnEx4HXLJMPnHRKmzzHxNztx0+OHyohQC5DPddV9+h89s51r13nndImOJY27rfCyCB+d/ijJfVpOnQ7XAcyIH5JPckMm5fW3+gnifLF8ptzAAAAAElFTkSuQmCC\n",
|
|
|
|
|
"text/latex": [
|
|
|
|
|
"$\\displaystyle \\left( 24, \\ 24, \\ 2\\right)$"
|
|
|
|
|
],
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"(24, 24, 2)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 60,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"from collections import deque\n",
|
|
|
|
|
"lookup_edge_c0_c1_to_number = []\n",
|
|
|
|
|
"r = []\n",
|
|
|
|
|
"c = []\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"i = 0\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"for _ in C_L_0:\n",
|
|
|
|
|
" r = []\n",
|
|
|
|
|
" for _ in C_L_1:\n",
|
|
|
|
|
" c = []\n",
|
|
|
|
|
" for _ in (I_4, CZ):\n",
|
|
|
|
|
" c.append(i)\n",
|
|
|
|
|
" i += 1 \n",
|
|
|
|
|
" r.append(c)\n",
|
|
|
|
|
" lookup_edge_c0_c1_to_number.append(r)\n",
|
|
|
|
|
" \n",
|
|
|
|
|
"lookup_edge_c0_c1_to_number = np.array(lookup_edge_c0_c1_to_number, dtype=np.int)\n",
|
|
|
|
|
"lookup_edge_c0_c1_to_number.shape"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"indices_contracted = {k: contract_indices(v) for k,v in indices.items()}"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 64,
|
2019-11-11 17:37:44 +00:00
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
"text": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"{\n",
|
|
|
|
|
" , { {0, 1}, {2, 3}, {4, 5}, {6, 7}, {8, 9}, {10, 11}, {12, 13}, {14, 15}, {16, 17}, {18, 19}, {20, 21}, {22, 23}, {24, 25}, {26, 27}, {28, 29}, {30, 31}, {32, 33}, {34, 35}, {36, 37}, {38, 39}, {40, 41}, {42, 43}, {44, 45}, {46, 47} }\n",
|
|
|
|
|
" , { {48, 49}, {50, 51}, {52, 53}, {54, 55}, {56, 57}, {58, 59}, {60, 61}, {62, 63}, {64, 65}, {66, 67}, {68, 69}, {70, 71}, {72, 73}, {74, 75}, {76, 77}, {78, 79}, {80, 81}, {82, 83}, {84, 85}, {86, 87}, {88, 89}, {90, 91}, {92, 93}, {94, 95} }\n",
|
|
|
|
|
" , { {96, 97}, {98, 99}, {100, 101}, {102, 103}, {104, 105}, {106, 107}, {108, 109}, {110, 111}, {112, 113}, {114, 115}, {116, 117}, {118, 119}, {120, 121}, {122, 123}, {124, 125}, {126, 127}, {128, 129}, {130, 131}, {132, 133}, {134, 135}, {136, 137}, {138, 139}, {140, 141}, {142, 143} }\n",
|
|
|
|
|
" , { {144, 145}, {146, 147}, {148, 149}, {150, 151}, {152, 153}, {154, 155}, {156, 157}, {158, 159}, {160, 161}, {162, 163}, {164, 165}, {166, 167}, {168, 169}, {170, 171}, {172, 173}, {174, 175}, {176, 177}, {178, 179}, {180, 181}, {182, 183}, {184, 185}, {186, 187}, {188, 189}, {190, 191} }\n",
|
|
|
|
|
" , { {192, 193}, {194, 195}, {196, 197}, {198, 199}, {200, 201}, {202, 203}, {204, 205}, {206, 207}, {208, 209}, {210, 211}, {212, 213}, {214, 215}, {216, 217}, {218, 219}, {220, 221}, {222, 223}, {224, 225}, {226, 227}, {228, 229}, {230, 231}, {232, 233}, {234, 235}, {236, 237}, {238, 239} }\n",
|
|
|
|
|
" , { {240, 241}, {242, 243}, {244, 245}, {246, 247}, {248, 249}, {250, 251}, {252, 253}, {254, 255}, {256, 257}, {258, 259}, {260, 261}, {262, 263}, {264, 265}, {266, 267}, {268, 269}, {270, 271}, {272, 273}, {274, 275}, {276, 277}, {278, 279}, {280, 281}, {282, 283}, {284, 285}, {286, 287} }\n",
|
|
|
|
|
" , { {288, 289}, {290, 291}, {292, 293}, {294, 295}, {296, 297}, {298, 299}, {300, 301}, {302, 303}, {304, 305}, {306, 307}, {308, 309}, {310, 311}, {312, 313}, {314, 315}, {316, 317}, {318, 319}, {320, 321}, {322, 323}, {324, 325}, {326, 327}, {328, 329}, {330, 331}, {332, 333}, {334, 335} }\n",
|
|
|
|
|
" , { {336, 337}, {338, 339}, {340, 341}, {342, 343}, {344, 345}, {346, 347}, {348, 349}, {350, 351}, {352, 353}, {354, 355}, {356, 357}, {358, 359}, {360, 361}, {362, 363}, {364, 365}, {366, 367}, {368, 369}, {370, 371}, {372, 373}, {374, 375}, {376, 377}, {378, 379}, {380, 381}, {382, 383} }\n",
|
|
|
|
|
" , { {384, 385}, {386, 387}, {388, 389}, {390, 391}, {392, 393}, {394, 395}, {396, 397}, {398, 399}, {400, 401}, {402, 403}, {404, 405}, {406, 407}, {408, 409}, {410, 411}, {412, 413}, {414, 415}, {416, 417}, {418, 419}, {420, 421}, {422, 423}, {424, 425}, {426, 427}, {428, 429}, {430, 431} }\n",
|
|
|
|
|
" , { {432, 433}, {434, 435}, {436, 437}, {438, 439}, {440, 441}, {442, 443}, {444, 445}, {446, 447}, {448, 449}, {450, 451}, {452, 453}, {454, 455}, {456, 457}, {458, 459}, {460, 461}, {462, 463}, {464, 465}, {466, 467}, {468, 469}, {470, 471}, {472, 473}, {474, 475}, {476, 477}, {478, 479} }\n",
|
|
|
|
|
" , { {480, 481}, {482, 483}, {484, 485}, {486, 487}, {488, 489}, {490, 491}, {492, 493}, {494, 495}, {496, 497}, {498, 499}, {500, 501}, {502, 503}, {504, 505}, {506, 507}, {508, 509}, {510, 511}, {512, 513}, {514, 515}, {516, 517}, {518, 519}, {520, 521}, {522, 523}, {524, 525}, {526, 527} }\n",
|
|
|
|
|
" , { {528, 529}, {530, 531}, {532, 533}, {534, 535}, {536, 537}, {538, 539}, {540, 541}, {542, 543}, {544, 545}, {546, 547}, {548, 549}, {550, 551}, {552, 553}, {554, 555}, {556, 557}, {558, 559}, {560, 561}, {562, 563}, {564, 565}, {566, 567}, {568, 569}, {570, 571}, {572, 573}, {574, 575} }\n",
|
|
|
|
|
" , { {576, 577}, {578, 579}, {580, 581}, {582, 583}, {584, 585}, {586, 587}, {588, 589}, {590, 591}, {592, 593}, {594, 595}, {596, 597}, {598, 599}, {600, 601}, {602, 603}, {604, 605}, {606, 607}, {608, 609}, {610, 611}, {612, 613}, {614, 615}, {616, 617}, {618, 619}, {620, 621}, {622, 623} }\n",
|
|
|
|
|
" , { {624, 625}, {626, 627}, {628, 629}, {630, 631}, {632, 633}, {634, 635}, {636, 637}, {638, 639}, {640, 641}, {642, 643}, {644, 645}, {646, 647}, {648, 649}, {650, 651}, {652, 653}, {654, 655}, {656, 657}, {658, 659}, {660, 661}, {662, 663}, {664, 665}, {666, 667}, {668, 669}, {670, 671} }\n",
|
|
|
|
|
" , { {672, 673}, {674, 675}, {676, 677}, {678, 679}, {680, 681}, {682, 683}, {684, 685}, {686, 687}, {688, 689}, {690, 691}, {692, 693}, {694, 695}, {696, 697}, {698, 699}, {700, 701}, {702, 703}, {704, 705}, {706, 707}, {708, 709}, {710, 711}, {712, 713}, {714, 715}, {716, 717}, {718, 719} }\n",
|
|
|
|
|
" , { {720, 721}, {722, 723}, {724, 725}, {726, 727}, {728, 729}, {730, 731}, {732, 733}, {734, 735}, {736, 737}, {738, 739}, {740, 741}, {742, 743}, {744, 745}, {746, 747}, {748, 749}, {750, 751}, {752, 753}, {754, 755}, {756, 757}, {758, 759}, {760, 761}, {762, 763}, {764, 765}, {766, 767} }\n",
|
|
|
|
|
" , { {768, 769}, {770, 771}, {772, 773}, {774, 775}, {776, 777}, {778, 779}, {780, 781}, {782, 783}, {784, 785}, {786, 787}, {788, 789}, {790, 791}, {792, 793}, {794, 795}, {796, 797}, {798, 799}, {800, 801}, {802, 803}, {804, 805}, {806, 807}, {808, 809}, {810, 811}, {812, 813}, {814, 815} }\n",
|
|
|
|
|
" , { {816, 817}, {818, 819}, {820, 821}, {822, 823}, {824, 825}, {826, 827}, {828, 829}, {830, 831}, {832, 833}, {834, 835}, {836, 837}, {838, 839}, {840, 841}, {842, 843}, {844, 845}, {846, 847}, {848, 849}, {850, 851}, {852, 853}, {854, 855}, {856, 857}, {858, 859}, {860, 861}, {862, 863} }\n",
|
|
|
|
|
" , { {864, 865}, {866, 867}, {868, 869}, {870, 871}, {872, 873}, {874, 875}, {876, 877}, {878, 879}, {880, 881}, {882, 883}, {884, 885}, {886, 887}, {888, 889}, {890, 891}, {892, 893}, {894, 895}, {896, 897}, {898, 899}, {900, 901}, {902, 903}, {904, 905}, {906, 907}, {908, 909}, {910, 911} }\n",
|
|
|
|
|
" , { {912, 913}, {914, 915}, {916, 917}, {918, 919}, {920, 921}, {922, 923}, {924, 925}, {926, 927}, {928, 929}, {930, 931}, {932, 933}, {934, 935}, {936, 937}, {938, 939}, {940, 941}, {942, 943}, {944, 945}, {946, 947}, {948, 949}, {950, 951}, {952, 953}, {954, 955}, {956, 957}, {958, 959} }\n",
|
|
|
|
|
" , { {960, 961}, {962, 963}, {964, 965}, {966, 967}, {968, 969}, {970, 971}, {972, 973}, {974, 975}, {976, 977}, {978, 979}, {980, 981}, {982, 983}, {984, 985}, {986, 987}, {988, 989}, {990, 991}, {992, 993}, {994, 995}, {996, 997}, {998, 999}, {1000, 1001}, {1002, 1003}, {1004, 1005}, {1006, 1007} }\n",
|
|
|
|
|
" , { {1008, 1009}, {1010, 1011}, {1012, 1013}, {1014, 1015}, {1016, 1017}, {1018, 1019}, {1020, 1021}, {1022, 1023}, {1024, 1025}, {1026, 1027}, {1028, 1029}, {1030, 1031}, {1032, 1033}, {1034, 1035}, {1036, 1037}, {1038, 1039}, {1040, 1041}, {1042, 1043}, {1044, 1045}, {1046, 1047}, {1048, 1049}, {1050, 1051}, {1052, 1053}, {1054, 1055} }\n",
|
|
|
|
|
" , { {1056, 1057}, {1058, 1059}, {1060, 1061}, {1062, 1063}, {1064, 1065}, {1066, 1067}, {1068, 1069}, {1070, 1071}, {1072, 1073}, {1074, 1075}, {1076, 1077}, {1078, 1079}, {1080, 1081}, {1082, 1083}, {1084, 1085}, {1086, 1087}, {1088, 1089}, {1090, 1091}, {1092, 1093}, {1094, 1095}, {1096, 1097}, {1098, 1099}, {1100, 1101}, {1102, 1103} }\n",
|
|
|
|
|
" , { {1104, 1105}, {1106, 1107}, {1108, 1109}, {1110, 1111}, {1112, 1113}, {1114, 1115}, {1116, 1117}, {1118, 1119}, {1120, 1121}, {1122, 1123}, {1124, 1125}, {1126, 1127}, {1128, 1129}, {1130, 1131}, {1132, 1133}, {1134, 1135}, {1136, 1137}, {1138, 1139}, {1140, 1141}, {1142, 1143}, {1144, 1145}, {1146, 1147}, {1148, 1149}, {1150, 1151} }\n",
|
|
|
|
|
"}\n"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2019-11-13 15:08:33 +00:00
|
|
|
|
"print(\"{\")\n",
|
|
|
|
|
"for r in lookup_edge_c0_c1_to_number:\n",
|
|
|
|
|
" print(\" , {\", \", \".join((\"{%d, %d}\" % (c[0], c[1]) for c in r)), \"}\")\n",
|
|
|
|
|
"print(\"}\") "
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 68,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([ 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108,\n",
|
|
|
|
|
" 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121,\n",
|
|
|
|
|
" 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134,\n",
|
|
|
|
|
" 135, 136, 137, 138, 139, 140, 141, 142, 143])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 68,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"I0_indices = lookup_edge_c0_c1_to_number[2, :, :].reshape(2*24)\n",
|
|
|
|
|
"I0_indices"
|
2019-11-11 17:37:44 +00:00
|
|
|
|
]
|
2019-11-13 15:08:33 +00:00
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 69,
|
|
|
|
|
"metadata": {
|
|
|
|
|
"scrolled": true
|
|
|
|
|
},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([ 4, 5, 52, 53, 100, 101, 148, 149, 196, 197, 244,\n",
|
|
|
|
|
" 245, 292, 293, 340, 341, 388, 389, 436, 437, 484, 485,\n",
|
|
|
|
|
" 532, 533, 580, 581, 628, 629, 676, 677, 724, 725, 772,\n",
|
|
|
|
|
" 773, 820, 821, 868, 869, 916, 917, 964, 965, 1012, 1013,\n",
|
|
|
|
|
" 1060, 1061, 1108, 1109])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 69,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"I1_indices = lookup_edge_c0_c1_to_number[:, 2, :].reshape(2*24)\n",
|
|
|
|
|
"I1_indices"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": []
|
2019-11-11 13:36:11 +00:00
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"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.3"
|
|
|
|
|
}
|
|
|
|
|
},
|
|
|
|
|
"nbformat": 4,
|
|
|
|
|
"nbformat_minor": 2
|
|
|
|
|
}
|