Resultados GKB vs. GBB.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "import skprocrustes as skp\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import time\n",
    "import scipy.linalg as sp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "p = skp.ProcrustesProblem((5000,5000,10,10), problemnumber=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = np.diag(2*np.ones((5000,))) + np.diag(-1*np.ones((4999,)), k = 1) + np.diag(-1*np.ones((4999,)), k = -1)\n",
    "B = np.dot(A, np.eye(5000,10))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "p2 = skp.ProcrustesProblem((5000,5000,10,10), matrices = (A, B, np.eye(10,10)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Warning: Maximum number of iterations has been exceeded.\n"
     ]
    }
   ],
   "source": [
    "gkbsolver = skp.GKBSolver(inner_solver=\"spg\", verbose=0)\n",
    "t0 = time.time(); results_gkb = gkbsolver.solve(p2); cpu_gkb = time.time()-t0;\n",
    "gbbsolver = skp.GKBSolver(inner_solver=\"gbb\", halfreorth=True, verbose=0)\n",
    "t0 = time.time(); results_gbb = gbbsolver.solve(p2); cpu_gkbgbb = time.time()-t0;\n",
    "puregbbsolver = skp.GBBSolver(verbose=0)\n",
    "t0 = time.time(); results_puregbb = puregbbsolver.solve(p2); cpu_puregbb = time.time()-t0;\n",
    "spgsolver = skp.SPGSolver(verbose=0)\n",
    "t0 = time.time(); results_spg = spgsolver.solve(p2); cpu_spg = time.time()-t0;"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Resumo"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 1.0, 'Problem 1')"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEICAYAAABF82P+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAFLpJREFUeJzt3X+0ZWV93/H3R0DGCAkiA0EYGCyzaMQfiKPQoomAzRIaA1qIWhKRRZ20C5sYYyOGNJhqTExdwRqtCRUVKxooDTJarCI/1Gi1DoogpcpAxRlxwaBAUBACfvvHeW443nnm3gtz9z33zn2/1pp19n6e5+zzPXvduZ+7f6eqkCRpusdNugBJ0uJkQEiSugwISVKXASFJ6jIgJEldBoQkqcuAkB6FJG9O8uEZ+r+d5EULWZM0FANCO7z2S/v+JD9McnuSDyTZbdJ1bY8kT0/yqSR3JvFiJg3CgNBy8ZKq2g04HHgu8AfTB2Rkqfyf+HvgIuD0SReiHddS+c8gzYuq+i7wSeDpAEmuTvLHSb4A3Ac8NclTkqxP8oMkG5O8ZtpiViS5MMm9Sb6a5Fm9z0ryuCRnJrk5yfeTXJRkz9a3OkklOS3JpiR3JfnXSZ6b5Lokdyd59wzf45tVdR5ww3ysF6nHgNCykmQVcDzwtbHm3wDWAbsDtwIfBTYDTwFOAt6W5Nix8ScA/w3YE/gI8LEku3Q+7reAE4Ffasu6C3jPtDFHAGuAlwPvBM4CXgQcCvxakl96rN9V2l4GhJaLjyW5G/hb4LPA28b6PlhVN1TVQ8DPA88H3lhVP66qa4H3MQqRKddU1cVV9ffAnwMrgCM7n/mbwFlVtbmqHgDeDJyUZOexMW9pn/Np4EfAR6vqjral83ng2fPw3aXHZOfZh0g7hBOr6jPb6Ns0Nv0U4AdVde9Y263A2t74qvpJkqmtjekOBC5J8pOxtoeBfcbmbx+bvr8zv6QPpmtpcwtCgvGzgG4D9kyy+1jbAcB3x+ZXTU20g9r7t/dNtwk4rqr2GPu3om0dSIueASGNqapNwBeBP0myIskzGZ0pdMHYsOckeVnbVfQ64AHgS53F/SXwx0kOBEiyMskJ81FnO+NqBfD4Nr8iya7zsWxpigEhbe2VwGpGWwWXAGdX1eVj/ZcyOqh8F6NjEy9rxyOm+0/AeuDTSe5lFCJHzFONBzLaBTV1FtP9wDfnadkSAPGBQZKkHrcgJEldBoQkqcuAkCR1GRCSpK4lfaHcXnvtVatXr550GZK0pFxzzTV3VtXK2cYt6YBYvXo1GzZsmHQZkrSkJLl1LuPcxSRJ6jIgJEldBoQkqcuAkCR1GRCSpC4DQpLUZUBIkroMCElSlwEhSepa0ldSS0tVMukKJs9H0Sx+bkFIkroMCElSlwEhSeoyICRJXQaEJKnLgJAkdRkQkqQuA0KS1GVASJK6DAhJUpcBIUnqMiAkSV0GhCSpy4CQJHUZEJKkLgNCktRlQEiSugwISVKXASFJ6jIgJEldBoQkqWvQgEjy7STXJ7k2yYbWtmeSy5Pc1F6f1NqT5F1JNia5LsnhQ9YmSZrZQmxBHF1Vh1XV2jZ/JnBFVa0BrmjzAMcBa9q/dcB7F6A2SdI2TGIX0wnA+W36fODEsfYP1ciXgD2S7DuB+iRJDB8QBXw6yTVJ1rW2farqewDtde/Wvh+waey9m1vbT0myLsmGJBu2bNkyYOmStLztPPDyj6qq25LsDVye5P/OMDadttqqoepc4FyAtWvXbtUvSZofg25BVNVt7fUO4BLgecDtU7uO2usdbfhmYNXY2/cHbhuyPknStg0WEEmemGT3qWngl4FvAOuBU9uwU4FL2/R64FXtbKYjgXumdkVJkhbekLuY9gEuSTL1OR+pqv+Z5CvARUlOB74DnNzGXwYcD2wE7gNOG7A2SdIsBguIqroFeFan/fvAsZ32As4Yqh5J0qPjldSSpC4DQpLUZUBIkroMCElSlwEhSeoyICRJXQaEJKnLgJAkdRkQkqQuA0KS1GVASJK6DAhJUpcBIUnqMiAkSV0GhCSpy4CQJHUZEJKkLgNCktRlQEiSugwISVKXASFJ6jIgJEldBoQkqcuAkCR1GRCSpC4DQpLUZUBIkroGD4gkOyX5WpJPtPmDknw5yU1JLkzy+Na+a5vf2PpXD12bJGnbFmIL4reBG8fm3w6cU1VrgLuA01v76cBdVXUwcE4bJ0makEEDIsn+wD8H3tfmAxwDXNyGnA+c2KZPaPO0/mPbeEnSBAy9BfFO4PeAn7T5JwN3V9VDbX4zsF+b3g/YBND672njf0qSdUk2JNmwZcuWIWuXpGVtsIBI8ivAHVV1zXhzZ2jNoe+Rhqpzq2ptVa1duXLlPFQqSerZecBlHwX8apLjgRXAzzLaotgjyc5tK2F/4LY2fjOwCticZGfg54AfDFifJGkGg21BVNWbqmr/qloNvAK4sqpOAa4CTmrDTgUubdPr2zyt/8qq2moLQpK0MCZxHcQbgdcn2cjoGMN5rf084Mmt/fXAmROoTZLUDLmL6R9U1dXA1W36FuB5nTE/Bk5eiHokSbPzSmpJUpcBIUnqMiAkSV0GhCSpy4CQJHUZEJKkLgNCktRlQEiSugwISVLXnAIiyfOTnNamVyY5aNiyJEmTNmtAJDmb0f2T3tSadgE+PGRRkqTJm8sWxEuBXwV+BFBVtwG7D1mUJGny5hIQD7bbbhdAkicOW5IkaTGYS0BclOSvGD3o5zXAZ4D/MmxZkqRJm/V231X1jiT/DPg74BDgD6vq8sErkyRN1JyeB1FVlyf58tT4JHtWlY8DlaQd2KwBkeQ3gf8A3A/8BAij4xFPHbY0SdIkzWUL4g3AoVV159DFSJIWj7kcpL4ZuG/oQiRJi8tctiDeBHyxHYN4YKqxqn5rsKokSRM3l4D4K+BK4HpGxyAkScvAXALioap6/eCVSJIWlbkcg7gqybok+ybZc+rf4JVJkiZqLlsQ/7K9vmmszdNcJWkHN5crqb21tyQtQ9sMiCTHVNWVSV7W66+qvxmuLEnSpM20BfGLjM5eekmnrwADQpJ2YDMFxHUAVXXaAtUiSVpEZjqL6Q+2Z8FJViT530m+nuSGJH/U2g9K8uUkNyW5MMnjW/uubX5j61+9PZ8vSdo+c3om9WP0AHBMVT0LOAx4cZIjgbcD51TVGuAu4PQ2/nTgrqo6GDinjZMkTchMu5j+cZLrOu0BqqqeOdOC21Pofthmd2n/CjiGR06dPR94M/Be4IQ2DXAx8O4kacuRJC2wmQLi/9E/QD1nSXYCrgEOBt7D6MZ/d1fVQ23IZmC/Nr0fsAmgqh5Kcg/wZODOactcB6wDOOCAA7anPEnSDGYKiAer6tbtWXhVPQwclmQP4BLgF3rD2mtm6Btf5rnAuQBr165160KSBjLTMYgvzNeHVNXdwNXAkYyebT0VTPsDt7XpzcAqgNb/c4BPrZOkCdlmQFTVa7dnwUlWti0HkjwBeBFwI3AVcFIbdipwaZte3+Zp/Vd6/EGSJmdOz6R+jPYFzm/HIR4HXFRVn0jyf4C/TvJW4GvAeW38ecB/TbKR0ZbDKwasTZI0i8ECoqquA57dab8FeF6n/cfAyUPVI0l6dGa6F9P0ezAVozOKrq2qewetSpI0cTNtQfROcd0TeGaS06vqyoFqkiQtAtsMiG3dgynJgcBFwBFDFSVJmrxHfauNdm3ELgPUIklaRB51QCQ5hNF9liRJO7CZDlJ/nK2vZN6T0emrvz5kUZKkyZvpIPU7ps0X8H3gpqp6cLiSJEmLwUxXUn8WeBLwXGBFVX2uqm4wHCRpedhmQCT5z8DvMLqj6luS/PsFq0qSNHGzPZP6WVX1cJKfAT4PvGVhypIkTdpMZzE92G7XTVXdR/923JKkHdRcnygX4B+1+Tk9UU6StLTNFBC9h/tIkpaJmQJiF2CfqvqpBwcleQGPPORHkrSDmukYxDuB3l1b7299kqQd2EwBsbo90+GnVNUGYPVgFUmSFoWZAmLFDH1PmO9CJEmLy0wB8ZUkr5nemOR04JrhSpIkLQYzHaR+HXBJklN4JBDWAo8HXjp0YZKkyZrpgUG3A/80ydHA01vz//BJcpK0PMy0BQFAVV0FXLUAtUiSFpFH/cAgSdLyYEBIkroMCElSlwEhSeoyICRJXQaEJKlrsIBIsirJVUluTHJDkt9u7XsmuTzJTe31Sa09Sd6VZGOS65IcPlRtkqTZDbkF8RDwu1X1C8CRwBlJngacCVxRVWuAK9o8wHHAmvZvHfDeAWuTJM1isICoqu9V1Vfb9L3AjcB+wAnA+W3Y+cCJbfoE4EM18iVgjyT7DlWfJGlmC3IMIslq4NnAlxk9hOh7MAoRYO82bD9g09jbNre26ctal2RDkg1btmwZsmxJWtYGD4gkuwH/HXhdVf3dTEM7bbVVQ9W5VbW2qtauXLlyvsqUJE0zaEAk2YVROFxQVX/Tmm+f2nXUXu9o7ZuBVWNv3x8fbSpJEzPkWUwBzgNurKo/H+taD5zapk8FLh1rf1U7m+lI4J6pXVGSpIU3691ct8NRwG8A1ye5trX9PvCnwEXtwUPfAU5ufZcBxwMbgfuA0wasTZI0i8ECoqr+lv5xBYBjO+MLOGOoeiRJj45XUkuSugwISVKXASFJ6jIgJEldBoQkqcuAkCR1GRCSpC4DQpLUZUBIkroMCElSlwEhSeoyICRJXQaEJKnLgJAkdRkQkqQuA0KS1GVASJK6DAhJUpcBIUnqMiAkSV0GhCSpy4CQJHUZEJKkLgNCktRlQEiSugwISVKXASFJ6hosIJK8P8kdSb4x1rZnksuT3NRen9Tak+RdSTYmuS7J4UPVJUmamyG3ID4IvHha25nAFVW1BriizQMcB6xp/9YB7x2wLknSHAwWEFX1OeAH05pPAM5v0+cDJ461f6hGvgTskWTfoWqTJM1uoY9B7FNV3wNor3u39v2ATWPjNre2rSRZl2RDkg1btmwZtFhJWs4Wy0HqdNqqN7Cqzq2qtVW1duXKlQOXJUnL10IHxO1Tu47a6x2tfTOwamzc/sBtC1ybJGnMQgfEeuDUNn0qcOlY+6va2UxHAvdM7YqSJE3GzkMtOMlHgRcCeyXZDJwN/ClwUZLTge8AJ7fhlwHHAxuB+4DThqpLkjQ3gwVEVb1yG13HdsYWcMZQtUiSHr3FcpBakrTIGBCSpC4DQpLUZUBIkroMCElSlwEhSeoa7DRXSRrS1bl60iVM1AvrhYN/hlsQkqQuA0KS1GVASJK6DAhJUpcBIUnqMiAkSV0GhCSpy4CQJHUZEJKkLgNCktRlQEiSugwISVLX8r1ZXzLpCiaratIVSFrk3IKQJHUZEJKkLgNCktS1fI9BaLvkj5b3MZw622M42vG5BSFJ6jIgJEldBoQkqcuAkCR1LaqASPLiJN9MsjHJmZOuR5KWs0UTEEl2At4DHAc8DXhlkqdNtipJWr4WTUAAzwM2VtUtVfUg8NfACROuSZKWrcV0HcR+wKax+c3AEdMHJVkHrGuzP0zyzQWobQh7AXdO7NOX/r2oJrr+8uYlv/5g0utw6a/CCf8f3q53HziXQYspIHpfd6urkarqXODc4csZVpINVbV20nUsVa6/7ec63D7LYf0tpl1Mm4FVY/P7A7dNqBZJWvYWU0B8BViT5KAkjwdeAayfcE2StGwtml1MVfVQktcCnwJ2At5fVTdMuKwhLfndZBPm+tt+rsPts8Ovv5QPjpEkdSymXUySpEXEgJAkdRkQA0iyT5KPJLklyTVJ/leSlyZ5YZJPjI17a5JPJdk1ydXtNiPXJrmxXe+xZE1iHST59STXJbkhydeTvC/JHq1vm8tO8u0k17e+65Ms6Qs0Z1n397TveV2SzyTZu73n1Um2tL4bklyc5Gcm/V0mJclZbT1c19bJEWM/Q19P8oUkh7SxOyd5W5Kb2thrk5w16e8wHwyIeZYkwMeAz1XVU6vqOYzOyNp/2rizgKOAE6vqgdZ8SlUd1trf3s7mWnKGXAdJVie5uvOZLwZ+Bziuqg4FDge+COwzNmymZR/d+k4C3vUYv/rEzWHdf76qDquqZzI6c/CMsbdf2PoOBR4EXr6QtS8WSf4J8CvA4W09vYhHLuI9paqeBZwP/MfW9lbgKcAz2s/QC4BdFrbqYRgQ8+8Y4MGq+suphqq6tar+Ymo+ye8CxwMvqar7O8vYDfgR8PDQxQ5kEuvgLOANVfXd9nkPV9X7q6p3pf1My/5Z4K45fuZiNOu6h38Ikt3pfNckOwNP7PUtE/sCd0790VJVd1bV9GuyPgcc3LayXgP826r6cRt/b1W9eSELHsqiOc11B3Io8NUZ+o8CDgGeU1U/nNZ3QZIHgDXA66pqqQbEJNbBbJ8527Kvar80nwr82hw/czGabT28IMm1wJMZheTvj/W9PMnzGf2C/Bbw8cGqXNw+Dfxhkm8Bn2G0ZfXZaWNeAlwPHAx8p6ruXeAaF4RbEANL8p62z/IrrWkjo9uK/HJn+Cltk/YA4A1J5nS/lMVuPtZBkkvaL7bLgLVj+3pP63zeM1rfzUlePtuym6Or6unAM4B3J9lte7/3YtBZ91O7mFYBHwD+bGz4hW0Xyc8z+uX37xa43EWh/dHyHEb3fNsCXJjk1a37gvZzeBTwhunvTXJa+9nblGTV9P6lxoCYfzcw2v8NQFWdARwLrGxNtzPatXJOkqN7C6iqLYz+CtzqZoVLxLyvg6p6afvldTywof2SO6yqPjD9M6vq+jb2k8ATZlv2tL6bW31L9Vbzs637ceuBX5zeWKOLoz7e61su2i7Kq6vqbOC1wL9oXae0n7sTq2oToz92Dkiye3vfB9rP3j2MLvhd0gyI+XclsCLJvxlr+6mzQarqW8DLgA8nOWz6Atp+zWcDNw9Z6IAmsQ7+BHhHkvED4VuFw2zLbmf1HATcOsfPXWxmXfdjns+21+9MfTu0JIckWTPWdBjb+HmoqvuA8xhtda5o798JWJInmEznMYh5VlWV5ERGfx3/HqNN1B8Bb5w27itt98j6sb+iL0hyP7Ar8MGqumYha58vk1gHVXVZkpXAJ9t/0LuBbzC6dcuUmZZ9VZKHGZ19cmZV3f5ov/diMId1P3UMIoz+yv1XY2+fOgbxOEY3z3z1ghW+uOwG/EVGp0g/xGgrYR1w8TbGnwW8BfhGknuB+xmd5bTkbzbqrTYkSV3uYpIkdRkQkqQuA0KS1GVASJK6DAhJUpcBIUnqMiAkSV3/H7m2KTgQB1goAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "gkb, gbb, puregbb, spg = plt.bar([0,1,2,3], [cpu_gkb, cpu_gkbgbb, cpu_puregbb, cpu_spg])\n",
    "gkb.set_facecolor('r')\n",
    "gbb.set_facecolor('g')\n",
    "puregbb.set_facecolor('b')\n",
    "spg.set_facecolor('m')\n",
    "ax = plt.gca()\n",
    "ax.set_xticks([0,1,2,3])\n",
    "ax.set_xticklabels(['GKB', 'GKB+GBB', 'GBB', 'SPG'])\n",
    "#ax.set_ylim([0, 100])\n",
    "ax.set_ylabel('CPU Time')\n",
    "ax.set_title('Problem 1')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "np.linalg.cond(p.A)"
   ]
  },
  {
   "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.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}