hhl_2q.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "from logicqubit.logic import *\n",
    "from cmath import *\n",
    "import numpy as np\n",
    "import sympy as sp\n",
    "from scipy.linalg import expm, sinm, cosm, logm, sqrtm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "def qft(qr):\n",
    "    for i in range(len(qr)):\n",
    "        for j in range(i):\n",
    "            qr[i].CU1(qr[j], pi/float(2**(i-j)))\n",
    "        qr[i].H()\n",
    "        \n",
    "def iqft(qr):\n",
    "    for i in range(len(qr)):\n",
    "        for j in range(i):\n",
    "            qr[i].CU1(qr[j], -pi/float(2**(i-j)))\n",
    "        qr[i].H()\n",
    "        \n",
    "def swap(s1, s2):\n",
    "    s2.CX(s1)\n",
    "    s1.CX(s2)\n",
    "    s2.CX(s1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[(1.0000000000000004+4.440892098500627e-16j), 0j], [0j, (1.0000000000000004+4.440892098500627e-16j)]]\n"
     ]
    }
   ],
   "source": [
    "D = np.matrix('1.5707963267948966 0; 0 3.141592653589793')\n",
    "P = np.matrix('1 0; 0 1')/sqrt(2)\n",
    "\n",
    "A = np.matrix('3 1; 1 3')\n",
    "#A = P*D*P.I\n",
    "#A_ad = A.getH()\n",
    "t0 = 2*pi\n",
    "\n",
    "#U = np.matrix(expm(-1j*pi*A))\n",
    "U_A_div2 = np.matrix(expm(1j*A*t0/2))\n",
    "U_A_div4 = np.matrix(expm(1j*A*t0/4))\n",
    "print(U_A_div2.tolist())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 0.375 -0.125]\n"
     ]
    }
   ],
   "source": [
    "b = np.array([1,0])\n",
    "x = np.linalg.solve(A, b)\n",
    "print(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$\\displaystyle \\begin{pmatrix}0.5-2.78426e-17j \\\\ 0.5-2.76685e-17j \\\\ -2.23706e-16+7.5546e-17j \\\\ -1.11856e-17+1.46499e-16j \\\\ -2.45038e-16-4.73355e-17j \\\\ -2.77556e-17+1.76998e-17j \\\\ 0.353553+2.78426e-17j \\\\ -0.353553-1.14097e-17j \\\\ -1.89288e-17+2.78426e-17j \\\\ 9.81308e-18-2.78426e-17j \\\\ -1.26364e-16-3.10767e-17j \\\\ 6.51042e-17+3.10767e-17j \\\\ -1.50275e-16+3.776e-18j \\\\ -1.14967e-17-1.53233e-16j \\\\ 0.353553-2.78426e-17j \\\\ -0.353553-6.70949e-17j \\end{pmatrix}$"
      ],
      "text/plain": [
       "<IPython.core.display.Math object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "{'0000': (0.5-2.7842613392665686e-17j),\n",
       " '0001': (0.4999999999999999-2.7668537838592144e-17j),\n",
       " '0010': (-2.23705968685014e-16+7.554599442488008e-17j),\n",
       " '0011': (-1.1185571585378691e-17+1.4649861050015117e-16j),\n",
       " '0100': (-2.4503807662796223e-16-4.7335537465165435e-17j),\n",
       " '0101': (-2.7755575615628914e-17+1.7699821591488393e-17j),\n",
       " '0110': (0.3535533905932736+2.78426133926657e-17j),\n",
       " '0111': (-0.35355339059327373-1.1409698074428693e-17j),\n",
       " '1000': (-1.8928761761890062e-17+2.7842613392665704e-17j),\n",
       " '1001': (9.813077866773592e-18-2.7842613392665717e-17j),\n",
       " '1010': (-1.2636372189240907e-16-3.107669770100199e-17j),\n",
       " '1011': (6.510420273541852e-17+3.1076697701002e-17j),\n",
       " '1100': (-1.5027461392961003e-16+3.7760034294587894e-18j),\n",
       " '1101': (-1.149673585146546e-17-1.5323324243891858e-16j),\n",
       " '1110': (0.3535533905932737-2.7842613392665674e-17j),\n",
       " '1111': (-0.3535533905932737-6.709492485976009e-17j)}"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "logicQuBit  = LogicQuBit(4)\n",
    "\n",
    "a = Qubit()\n",
    "c1 = Qubit()\n",
    "c2 = Qubit()\n",
    "b = Qubit()\n",
    "r = 2\n",
    "\n",
    "c1.H()\n",
    "c2.H()\n",
    "\n",
    "#b.X()\n",
    "#b.H()\n",
    "\n",
    "# phase estimation\n",
    "b.CU(c2, U_A_div4.tolist())\n",
    "b.CU(c1, U_A_div2.tolist())\n",
    "\n",
    "# iqft\n",
    "swap(c1,c2)\n",
    "c2.H()\n",
    "c2.CS(c1, True)\n",
    "c1.H()\n",
    "swap(c1,c2)\n",
    "\n",
    "a.CRY(c1, pi/2**r)\n",
    "a.CRY(c2, 2*pi/2**r)\n",
    "\n",
    "swap(c1,c2)\n",
    "# qft\n",
    "c1.H()\n",
    "c2.CS(c1)\n",
    "c2.H()\n",
    "swap(c1,c2)\n",
    "\n",
    "\n",
    "# reverse phase estimation\n",
    "b.CU(c1, U_A_div4.conj().T.tolist())\n",
    "b.CU(c2, U_A_div2.conj().T.tolist())\n",
    "\n",
    "c1.H()\n",
    "c2.H()\n",
    "c2.PrintState()\n",
    "c2.getDictPsi()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEDCAYAAAAlRP8qAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAANXklEQVR4nO3df6zdd13H8eeL25UoBgV7CaY/aJUms5gt4rW4hEQXXdIyY8Et2qEsKNAUU4jGJasxmX8guhoCQyhrGmzQv6qJujRSbCL8gQYXe6dlSSddLlXspZDdDbM5nOvuePvHPVsOl3N7vrc79557P3s+kpuc7/f72bnvZN/73Hffe865qSokSevfK8Y9gCRpNAy6JDXCoEtSIwy6JDXCoEtSIwy6JDViw7i+8aZNm2r79u3j+vaStC499NBDj1fV5KBjYwv69u3bmZ6eHte3l6R1KcnXljrmLRdJaoRBl6RGGHRJaoRBl6RGGHRJakSnoCfZk+RCkpkkhwcc/7kkTyY51/u6Z/SjSpKuZujLFpNMAEeBW4BZ4GySU1X1yKKl/1hVv7gCM0qSOuhyhb4bmKmqi1V1BTgJ7FvZsSRJy9XljUWbgUt927PAWwasuynJl4HLwF1VdX7xgiQHgAMA27ZtW/60PdsPf/aa/1m17z/vvXXcIwCep1raSp2jXa7QM2Df4j9z9K/AG6rqRuATwAODnqiqjlfVVFVNTU4OfOeqJOkadQn6LLC1b3sLC1fhL6qqp6rq6d7j08B1STaNbEpJ0lBdgn4W2JlkR5KNwH7gVP+CJK9Pkt7j3b3nfWLUw0qSljb0HnpVzSc5BJwBJoATVXU+ycHe8WPA7cD7k8wDzwD7y78+LUmrqtOnLfZuo5xetO9Y3+NPAp8c7WiSpOXwnaKS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmNMOiS1AiDLkmN6BT0JHuSXEgyk+TwVdb9dJLnk9w+uhElSV0MDXqSCeAosBfYBdyRZNcS644AZ0Y9pCRpuC5X6LuBmaq6WFVXgJPAvgHrPgD8NfDYCOeTJHXUJeibgUt927O9fS9Kshl4B3BsdKNJkpajS9AzYF8t2r4PuLuqnr/qEyUHkkwnmZ6bm+s4oiSpiw0d1swCW/u2twCXF62ZAk4mAdgEvC3JfFU90L+oqo4DxwGmpqYW/0dBkvQSdAn6WWBnkh3A14H9wDv7F1TVjhceJ/kM8HeLYy5JWllDg15V80kOsfDqlQngRFWdT3Kwd9z75pK0BnS5QqeqTgOnF+0bGPKqevdLH0uStFy+U1SSGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRBl2SGmHQJakRnYKeZE+SC0lmkhwecHxfkoeTnEsyneStox9VknQ1G4YtSDIBHAVuAWaBs0lOVdUjfcs+D5yqqkpyA/BXwPUrMbAkabAuV+i7gZmqulhVV4CTwL7+BVX1dFVVb/NVQCFJWlVdgr4ZuNS3Pdvb912SvCPJV4DPAr85mvEkSV11CXoG7PueK/Cq+tuquh54O/ChgU+UHOjdY5+em5tb1qCSpKvrEvRZYGvf9hbg8lKLq+qLwI8l2TTg2PGqmqqqqcnJyWUPK0laWpegnwV2JtmRZCOwHzjVvyDJG5Ok9/jNwEbgiVEPK0la2tBXuVTVfJJDwBlgAjhRVeeTHOwdPwbcBtyZ5DngGeBX+35JKklaBUODDlBVp4HTi/Yd63t8BDgy2tEkScvhO0UlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqREGXZIaYdAlqRGdgp5kT5ILSWaSHB5w/NeSPNz7+lKSG0c/qiTpaoYGPckEcBTYC+wC7kiya9Gy/wB+tqpuAD4EHB/1oJKkq+tyhb4bmKmqi1V1BTgJ7OtfUFVfqqr/7m0+CGwZ7ZiSpGG6BH0zcKlve7a3bynvAT436ECSA0mmk0zPzc11n1KSNFSXoGfAvhq4MLmZhaDfPeh4VR2vqqmqmpqcnOw+pSRpqA0d1swCW/u2twCXFy9KcgPwaWBvVT0xmvEkSV11uUI/C+xMsiPJRmA/cKp/QZJtwN8A76qqR0c/piRpmKFX6FU1n+QQcAaYAE5U1fkkB3vHjwH3AD8MfCoJwHxVTa3c2JKkxbrccqGqTgOnF+071vf4vcB7RzuaJGk5fKeoJDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIwy6JDXCoEtSIzoFPcmeJBeSzCQ5POD49Un+OcmzSe4a/ZiSpGE2DFuQZAI4CtwCzAJnk5yqqkf6ln0L+CDw9pUYUpI0XJcr9N3ATFVdrKorwElgX/+Cqnqsqs4Cz63AjJKkDroEfTNwqW97trdPkrSGdAl6Buyra/lmSQ4kmU4yPTc3dy1PIUlaQpegzwJb+7a3AJev5ZtV1fGqmqqqqcnJyWt5CknSEroE/SywM8mOJBuB/cCplR1LkrRcQ1/lUlXzSQ4BZ4AJ4ERVnU9ysHf8WJLXA9PAq4HvJPltYFdVPbVyo0uS+g0NOkBVnQZOL9p3rO/xN1m4FSNJGhPfKSpJjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjTDoktQIgy5JjegU9CR7klxIMpPk8IDjSfKnveMPJ3nz6EeVJF3N0KAnmQCOAnuBXcAdSXYtWrYX2Nn7OgDcP+I5JUlDdLlC3w3MVNXFqroCnAT2LVqzD/iLWvAg8ENJfmTEs0qSrmJDhzWbgUt927PAWzqs2Qx8o39RkgMsXMEDPJ3kwrKm1VI2AY+Pe4i1IkfGPYEG8Bzt8xLP0TcsdaBL0DNgX13DGqrqOHC8w/fUMiSZrqqpcc8hLcVzdHV0ueUyC2zt294CXL6GNZKkFdQl6GeBnUl2JNkI7AdOLVpzCriz92qXnwGerKpvLH4iSdLKGXrLparmkxwCzgATwImqOp/kYO/4MeA08DZgBvhf4DdWbmQN4G0srXWeo6sgVd9zq1uStA75TlFJaoRBl6RGGHRJaoRBl6RGGHRJaoRBb1SSQe/eldYMz9HRM+gNeOEHo/8HpHw9qtYQz9HVYdAb0PeDsT3J7Un+Psl1Yx1K6uM5ujq6fDiX1qAkeeGHpPeRDL/Mwscy/B/wKuCVwHPjm1Avd56jq8+gr1N9PyjvA3YAPwncBbwbmKuqp8c3neQ5Og7eclnHkvwo8OvAF4D3A18BtgH39Y7771djk2QiyXbgDjxHV4VX6OtMkr0s/KGAp6vq35Pc+sKVTpKtLPwv7LMAVfWd8U2ql6skG6pqvqqeT/I14Jc8R1eHQV9HknwUuAl4ELg5yZ9X1cf6lvwB8FRVzY5lQL3sJbkP+J8kH6+qx3u3XfpvrXiOriCDvk4keRPwC8BNVfXtJD8BnEqyq6re11v2BeDfeutf4dWPVlOSW4F3AV8GNiS5v6r+q+/4dcDngXO9bc/REfPjc9eJJBuATwAfqaqv9vZ9HwufRX+uqn4nyQ8Cz/T+mLe0apJMAG9l4Wp8BrgfuAJ8BHi0/5xM8sqqenYsgzbOX0iscUle27vV8jrgGeDFWyxV9QzwHuAHknx/VT1pzLXakrwW+EPgUeBC7zx8J/AE8MfAriSvSXKvMV9ZBn2Nq6pvAZeAe4CPA08leSTJtt7xi8CNwJvGN6Veznrn6GPAnwCb+/b/LvA54AjwL8B1xnxlectlnUjy88BvAb/PwpszbgMeAH6KhVe8/Mr4ppMgyc3A3cAfVdUX+/afA75aVbf1tl98w5FGy6CvI0neCBwFPgw8D7wGeF1Vnegdn6iq58c4ol7meu+NOA58GvhL4MeB3wPurKryHF1ZBn2d6f3i8zPAPwCf6ns3nj8oWhOSvJqFc/SfgD8Dvt37Y/Mbqmp+rMM1zqCvQ71XvHwUmAc+XFVPjHkk6bv0ztGPsfC5Lfd6jq4Ofym6DvXehfdB4OvAkd6rDKQ1o3eOfgD4Jp6jq8Yr9HUuyQ3ART/oSGuV5+jqMeiS1AhvuUhSIwy6JDXCoEtSIwy6JDXCoEtSI/4fvGmsqkQnMYAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[(0.5+0j), (0.5+0j)]\n"
     ]
    }
   ],
   "source": [
    "res_a = logicQuBit.Measure_One(a)\n",
    "print(res_a)\n",
    "res = logicQuBit.Measure([b])\n",
    "logicQuBit.Plot()\n",
    "print(res)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "def Text2List(table):\n",
    "    list = [int(i, base=2) for i in table]\n",
    "    size = len(table[0])\n",
    "    tmp = sorted(list, key=int, reverse=False)  # values in ascending order\n",
    "    result = [[int(bin(j)[2:].zfill(size)[i]) for i in range(size)] for j in tmp]\n",
    "    return result\n",
    "\n",
    "def BinList(n):\n",
    "    blist = []\n",
    "    for i in iter(range(2 ** n)):\n",
    "        b = bin(i)[2:].zfill(n)  # value in binary, ex: i=1, n=4 -> '0001'\n",
    "        blist.append(b)\n",
    "    return Text2List(blist)\n",
    "\n",
    "#size=3\n",
    "#blist = [value for i,value in enumerate(BinList(size)) if res[i]>0.1] \n",
    "#results = [sum([value/2**(i+1) for i,value in enumerate(j)]) for j in blist]\n",
    "#print(results)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[0.0, 1.5707963267948966, 3.141592653589793, 4.71238898038469]"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# todos possíveis valores\n",
    "size = 2\n",
    "blist = [[i >> bl & 0x1 for bl in range(size)][::-1] for i in range(2**size)]\n",
    "results = [2*pi*sum([value/2**(i+1) for i,value in enumerate(j)]) for j in blist]\n",
    "results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}