latin_hypercube_test.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "tribal-custom",
   "metadata": {},
   "outputs": [],
   "source": [
    "import lhsmdu\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "id": "foster-psychiatry",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAS2klEQVR4nO3df4zkdX3H8ef75Ei6Kgv2Tqja3bH+KJ4eGF2pGhuxFytgNoolKbgtEW03tmq1iYbGTa3UbCJpm1wrIpmQxphMJE0F5OrRmqxRGoHCoQcLSClKd0ttA4hdImv0gHf/mNlzb1l2Z28/M9/Z7z4fySaz3/nezCvfnXvvd+c73+8rMhNJ0ta3o+oAkqQyHOiSVBMOdEmqCQe6JNWEA12SauKEqp54165d2Wg0qnp6SdqS7rjjjkczc/dq91U20BuNBocOHarq6SVpS4qIuWe7z7dcJKkmHOiSVBMOdEmqCQe6JNWEA12SasKBLkl90ppt0djfYMdlO2jsb9CabRV9/Mo+tihJ20lrtsXkgUkWjywCMLcwx+SBSQAm9k4UeQ730CWpD6Zmpo4O8yWLRxaZmpkq9hwOdEnqg/mF+Q0tPx4OdEnqg5HhkQ0tPx4OdEnqg+l90wztHDpm2dDOIab3TRd7Dge6JPXBxN4JmuNNRodHCYLR4VGa481iB0QBoqpO0bGxsfTiXJK0MRFxR2aOrXafe+iSVBMOdEmqCQe6JNWEA12SasKBLkk10dVAj4hzIuLfI+KBiPizVe4/JSKui4i7IuK2iHhN+aiSpLWsO9Aj4jnA54FzgT3ARRGxZ8VqnwQOZ+YZwMXA35YOKklaWzd76GcBD2TmDzLz58A1wLtWrLMHmAHIzPuARkScWjSpJGlN3Qz0FwP/tez7hzrLlrsTeA9ARJwFjAIvWflAETEZEYci4tAjjzxyfIklSavqZqDHKstWnl76WeCUiDgMfAT4LvDkM/5RZjMzxzJzbPfu3RvNqnX0+uL5kgZbNwUXDwG/uuz7lwA/XL5CZj4OXAIQEQE82PlSn/Tj4vmSBls3e+i3A6+IiJdGxInAhcANy1eIiJM79wH8AXBTZ8irT/px8XxJg23dPfTMfDIiPgz8C/Ac4O8z856I+GDn/quAVwFfioingHuBD/Qws1bRj4vnSxpsXXWKZuZB4OCKZVctu30L8Iqy0bQRI8MjzC3Mrbpc0vbgmaI10Y+L50sabA70mujHxfMlDTYLLiRpC7HgQpK2AQe6JNWEA12SasKBLkk14UCXpJpwoEtSTTjQJakmHOiSVBMOdEmqiVIl0cMRcSAi7oyIeyLikvJRJUlrKVUS/SHg3sw8Ezgb+Jtl10eXJPVBqZLoBJ7faSt6HvAYq1TQSZJ6p1RJ9BW0Sy5+CMwCH83Mp1c+kCXRktQ7pUqi3wEcBl4EvBa4IiJOesY/siRaknqmm4G+bkk07YLoa7PtAdoF0aeXiShJ6kaRkmhgHtgHEBGnAr8O/KBkUEnS2kqVRH8G+GJEzNJ+i+bSzHy0h7klSSuUKon+IfDbZaNJkjbCM0UlqSYc6JJUEw50SaoJB7ok1YQDXZJqwoEuSTXhQB9wrdkWjf0Ndly2g8b+Bq3ZVtWRJA2orj6Hrmq0ZltMHphk8cgiAHMLc0wemARgYu9EldEkDSD30AfY1MzU0WG+ZPHIIlMzUxUlkjTIHOgDbH5hfkPLJW1vDvQBNjI8sqHlkrY3B/oAm943zdDOoWOWDe0cYnrfdEWJJA2yUiXRn4iIw52vuyPiqYh4Qfm428vE3gma401Gh0cJgtHhUZrjTQ+ISlpVZK4sH1qxQrsk+n7g7bTLLm4HLsrMe59l/XHgTzPzt9Z63LGxsTx06NBxhZak7Soi7sjMsdXuK1USvdxFwJc3HlOStBmlSqIBiIgh4BzgK89yvyXRktQjpUqil4wD387Mx1a705JoSeqdUiXRSy7Et1skqRKlSqKJiGHgrcBXy0aUJHWjVEk0wPnA1zPziZ6llSQ9q3U/ttgrfmxRkjZusx9blCRtAQ50SaoJB7ok1YQDXZJqwoEuSTXhQJekmnCgS1JNONAlqSYc6JJUEw50aZtrzbZo7G+w47IdNPY3aM22qo6k47TutVwk1VdrtsXkgUkWjywCMLcwx+SBSQCrDregIp2inXXO7nSK3hMR3yobU1IvTM1MHR3mSxaPLDI1M1VRIm3GunvonU7Rz7OsUzQibljeKRoRJwNXAudk5nxEvLBHeSUVNL8wv6HlGmylOkXfC1ybmfMAmflw2ZiSemFkeGRDyzXYSnWKvhI4JSK+GRF3RMTFpQJK6p3pfdMM7Rw6ZtnQziGm901XlEib0c1B0W46RU8AXg/sA34JuCUibs3M+495oIhJYBJgZMQ9AKlqSwc+p2ammF+YZ2R4hOl90x4Q3aK6GejddIo+BDzaaSt6IiJuAs4EjhnomdkEmtAuuDje0JLKmdg74QCviVKdol8FfjMiToiIIeA3gO+VjSpJWkuRTtHM/F5E/DNwF/A0cHVm3t3L4JKkY9kpKklbiJ2ikrQNONAlqSYc6JIGihcLO35enEvSwPBiYZvjHrqkgeHFwjbHgS5pYHixsM1xoEsaGF4sbHMc6JIGhhcL2xwHuqSBMbF3guZ4k9HhUYJgdHiU5njTA6Jd8kxRSdpCPFNUkrYBB7ok1USRkuhOQfRCpyT6cER8qnxUSdJa1h3oy0qizwX2ABdFxJ5VVv3XzHxt5+svC+cEPCVYktbSzan/R0uiASJiqST63l4GW8lTgiVpbaVKogHeFBF3RsSNEfHq1R4oIiYj4lBEHHrkkUc2FNRTgiVpbd0M9G5Kor8DjGbmmcDngOtXe6DMbGbmWGaO7d69e0NBPSVYktbWzUBftyQ6Mx/PzJ90bh8EdkbErmIp8ZRgSVpPkZLoiDgtIqJz+6zO4/6oZFBPCZaktRUpiQYuAP4oIp4EfgpcmIVPQV068Dk1M8X8wjwjwyNM75v2gKgkdXjqvyRtIZ76L0nbgANdkmrCgS5JNeFAl6SacKBLUk040CWpJhzoklQTDnRJqgkHuiTVhANdkmrCgS5JNeFAl6SaKFISvWy9N0TEUxFxQbmI0vZid66O17qXz11WEv122mUXt0fEDZl57yrrXU77MruSjoPdudqMbvbQj5ZEZ+bPgaWS6JU+AnwFeLhgPmlbsTtXm1GkJDoiXgycD1y11gNtpiRa2g7sztVmlCqJ3g9cmplPrfVAmymJlrYDu3O1GUVKooEx4JqI+E/adXRXRsS7SwSUthO7c7UZRUqiM/OlmdnIzAbwj8AfZ+b1pcNKdTexd4LmeJPR4VGCYHR4lOZ40wOi6kqpkmhJhUzsnXCA67isO9ABMvMgcHDFslUHeWa+b/OxJEkb5ZmiklQTDnRJqgkHuiTVhANdkmrCgS5JNeFAl6SacKBLUk040CWpJhzoklQTDnRJqgkHuiTVhANdkmqiSEl0RLwrIu6KiMOdRqK3lI8qSVpLqZLoGeCGzMyIOAP4B+D0XgSWJK2uSEl0Zv4kM5dq6Z7LMyvqJEk9VqQkGiAizo+I+4CvAe9f7YEsiVbVWrMtGvsb7LhsB439DVqzraojScWUKokmM6/LzNOBdwOfWe2BLIlWlVqzLSYPTDK3MEeSzC3MMXlg0qGu2ihVEn1UZt4EvCwidm0ym1TU1MwUi0cWj1m2eGSRqZmpihJJZRUpiY6Il0dEdG6/DjgR+FHpsNJmzC/Mb2i5tNWUKon+HeDiiDgC/BT43WUHSaWBMDI8wtzC3KrLpTooUhKdmZcDl5eNJpU1vW+ayQOTx7ztMrRziOl90xWmksrxTFFtGxN7J2iONxkdHiUIRodHaY43mdg7UXU0qYio6p2RsbGxPHToUCXPLUlbVUTckZljq93nHrok1YQDXZJqwoEuSTXhQJekmnCgS1JNONAlqSYc6JJUEw50SaoJB7ok1USpTtGJTqfoXRFxc0ScWT6qJGkt6w70ZZ2i5wJ7gIsiYs+K1R4E3pqZZ9Aut2iWDipJWlupTtGbM/PHnW9vpV2CIUnqo2Kdost8ALhxtTvsFJWk3inWKQoQEW+jPdAvXe1+O0UlqXe6KbjoqlM0Is4ArgbOzUzr5ySpz0p1io4A1wK/n5n3l48pSVpPqU7RTwG/DFzZ6Yp+8tkuwC5J6g0biyRpC7GxSJK2AQe6JNWEA12SasKBLkk14UCXpJpwoEtSTTjQJakmHOiSVBMOdEmqCQe6eq4126Kxv8GOy3bQ2N+gNduqOpJUS91cbVE6bq3ZFpMHJlk8sgjA3MIckwcmAZjYO1FlNKl23ENXT03NTB0d5ksWjywyNTNVUSKpvkqVRJ8eEbdExM8i4uPlY2qrml+Y39ByScevVEn0Y8CfAH9dPKG2tJHhkQ0tl3T8SpVEP5yZtwNHepBRW9j0vmmGdg4ds2xo5xDT+6YrSiTVVy9Kop+VJdHbz8TeCZrjTUaHRwmC0eFRmuNND4hKPdDNp1y6LoleT2Y2gSa0Cy6O5zG09UzsnXCAS33QzR56VyXRkqRqFSmJliRVr0hJdEScBhwCTgKejoiPAXsy8/HeRZckLdfVmaKZeRA4uGLZVctu/y/tt2IkSRXxTFFJqgkHuiTVhANdkmrCgS5JNeFAl6SacKBLUk040CWpJhzoklQTDnRJqgkHuiTVhANdktbRmm3R2N9gx2U7aOxv0JptVR1pVV1dy0WStqvWbIvJA5NHy87nFuaYPDAJMHDX+S9VEh0R8Xed+++KiNeVjypJ/Tc1M3V0mC9ZPLLI1MxURYmeXamS6HOBV3S+JoEvFM4pSZWYX5jf0PIqFSmJ7nz/pWy7FTg5In6lcFZJ6ruR4ZENLa9SqZLoroqkLYmWtNVM75tmaOfQMcuGdg4xvW+6okTPrpuB3k1JdFdF0pnZzMyxzBzbvXt3N/kkqVITeydojjcZHR4lCEaHR2mONwfugCh09ymXbkqiLZKWVFsTeycGcoCvVKok+gbg4s6nXd4ILGTm/xTOKklaQ5GSaNp9o+cBDwCLwCW9iyxJWk2pkugEPlQ2miRpIzz1X5JqwoEuSTUR7XdLKnjiiEeAueP857uARwvGMYMZzFCfDDAYOXqVYTQzV/3cd2UDfTMi4lBmjpnBDGYww6DmqCKDb7lIUk040CWpJrbqQG9WHQAzLDFDmxnaBiEDDEaOvmfYku+hS5KeaavuoUuSVnCgS1JNDOxAH5Tauy5ynB4Rt0TEzyLi4xVlmOhsg7si4uaIOLOCDO/qPP/hzjXv39LvDMvWe0NEPBURF/Q7Q0ScHRELne1wOCI+1e8My3Icjoh7IuJb/c4QEZ9Ytg3u7vw8XtDnDMMRcSAi7uxsh+LXmOoiwykRcV3n/8ZtEfGa0hmOkZkD90X7ImDfB34NOBG4E9izYp3zgBtpX4v9jcC/VZTjhcAbgGng4xVleDNwSuf2uaW3RZcZnscvjsmcAdzX7wzL1vsG7WsPXVDBdjgb+KfSr4MNZjgZuBcYWXqNVvGzWLb+OPCNCrbDJ4HLO7d3A48BJ/Y5w18Bf9G5fTow06vXRmYO7B76oNTerZsjMx/OzNuBI4WfeyMZbs7MH3e+vZX29ej7neEn2XnVAs9llYKTXmfo+AjwFeDhws+/kQy91E2G9wLXZuY8tF+jFWRY7iLgyxVkSOD5ERG0dzgeA57sc4Y9wAxAZt4HNCLi1IIZjjGoA71Y7V0fcvTaRjN8gPZfLn3PEBHnR8R9wNeA9/c7Q0S8GDgfuIre6PZn8abOn/k3RsSrK8jwSuCUiPhmRNwRERdXkAGAiBgCzqH9S7bfGa4AXkW7bGcW+GhmPt3nDHcC7wGIiLOAUcrvcB01qAO9WO1dH3L0WtcZIuJttAf6pVVkyMzrMvN04N3AZyrIsB+4NDOfKvzcG8nwHdrX2jgT+BxwfQUZTgBeD7wTeAfw5xHxyj5nWDIOfDszHyv4/N1meAdwGHgR8Frgiog4qc8ZPkv7l+th2n89fpeyfyUco6vroVdgUGrvBqFar6sMEXEGcDVwbmb+qIoMSzLzpoh4WUTsysxSFyfqJsMYcE37L2x2AedFxJOZeX2/MmTm48tuH4yIKyvYDg8Bj2bmE8ATEXETcCZwfx8zLLmQ8m+3dJvhEuCznbcCH4iIB2m/j31bvzJ0Xg+XQPuDHMCDna/e6OUb9Js42HAC8APgpfziYMOrV6zzTo49KHpbFTmWrftpenNQtJttMUK7LerNFf48Xs4vDoq+Dvjvpe/7/bPorP9Fyh8U7WY7nLZsO5wFzPd7O9B+m2Gms+4QcDfwmn7/LIBh2u9bP7ei1+QXgE93bp/aeU3u6nOGk+kciAX+kPZxv6Lb4pjn6+WDb3JjnUd7j+L7wFRn2QeBD3ZuB/D5zv2zwFhFOU6j/Zv6ceD/OrdP6nOGq4Ef0/7z8jBwqILtcClwT+f5bwHe0u8MK9b9IoUHepfb4cOd7XAn7QPUxX/JdrMdgE/Q/qTL3cDHKsrwPuCa0s+9gZ/Fi4Cvd+bD3cDvVZDhTcB/APcB19L5NFqvvjz1X5JqYlAPikqSNsiBLkk14UCXpJpwoEtSTTjQJakmHOiSVBMOdEmqif8HykpqV6vaNpUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# from sahilm89/lhsmdu on Git\n",
    "\n",
    "k = lhsmdu.sample(5, 10) # Latin Hypercube Sampling with multi-dimensional uniformity, \n",
    "                        # 5 variables with 10 samples each\n",
    "k = np.array(k)\n",
    "fig = plt.figure()\n",
    "ax = fig.gca()\n",
    "ax.set_xticks(numpy.arange(0,1,0.1))\n",
    "ax.set_yticks(numpy.arange(0,1,0.1))\n",
    "plt.scatter(k[0], k[1], color=\"g\", label=\"LHS-MDU\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "negative-interest",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[0.31466083, 0.21596663, 0.17900832, 0.46162907, 0.00448606,\n",
       "        0.6622547 , 0.93089177, 0.53687252, 0.77697314, 0.88109352],\n",
       "       [0.64003566, 0.29471529, 0.94064627, 0.4771836 , 0.14838773,\n",
       "        0.58975323, 0.09077612, 0.80874168, 0.79364933, 0.38255348],\n",
       "       [0.27913107, 0.02009113, 0.68051192, 0.58501327, 0.82891503,\n",
       "        0.99522601, 0.30498601, 0.71499093, 0.15375811, 0.45757925],\n",
       "       [0.16451501, 0.50169691, 0.79603668, 0.40445881, 0.9143062 ,\n",
       "        0.20140145, 0.35668584, 0.69322669, 0.86659472, 0.08226575],\n",
       "       [0.70132812, 0.15420153, 0.3459563 , 0.64989679, 0.50722132,\n",
       "        0.06844354, 0.25032204, 0.87651488, 0.94852906, 0.41493816]])"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "k"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "id": "cleared-pakistan",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUQAAAEzCAYAAABJzXq/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAQXklEQVR4nO3db2jc933A8ffHfwJV2ymhcUvnRFI20qbe3Iz2mpWyP+nMVjvDhEIfJBULC4UjrCl9mDBBt1AE64OBKU1rbiGEwtE8aejikTYMjzaDNGtkSKO4JUVzsaKmEKctHlSw2MlnD/SJK2ty9JN8p/Nd3y8w6H73Pd3ni8I7P93pJ0VmIkmCHYMeQJKuFAZRkopBlKRiECWpGERJKgZRksqGQYyIhyPilYh44RL3R0R8OSIWIuL5iPhQ78eUpP5rcob4CHDwLe4/BNxY/9rA1y5/LEnafhsGMTOfAn75FktuB76eK54Bro6I9/ZqQEnaLr14DXEv8NKq20t1TJKGyq4efI5Y59i61wNGRJuVb6t5+9vf/uGbbrqpB08vSb9x4sSJVzNzz1Ye24sgLgHXr7p9HfDyegszswN0AFqtVs7NzfXg6X/7dOe7zByfYfHsIhPjE8wemGV6//Sgx5KuCBFxequP7cW3zI8Dd9W7zR8Fzmbmz3vwebWO7nyX9rE2p8+eJklOnz1N+1ib7nx30KNJQ6/Jj918A/g+8P6IWIqIz0TEPRFxTy15AjgFLAD/Avxd36YVM8dnWD63fNGx5XPLzByfGdBE0ujY8FvmzLxzg/sT+GzPJtJbWjy7uKnjkprzSpUhMzE+sanjkpoziENm9sAsY7vHLjo2tnuM2QOzA5pIGh0GcchM75+mc7jD5PgkQTA5PknncMd3maUeiEH9CQF/7EZSP0TEicxsbeWxniFKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQMoiQVgyhJxSBKUjGIklQaBTEiDkbEixGxEBH3r3P/eEQci4gfRsTJiLi796NKUn9tGMSI2Ak8CBwC9gF3RsS+Ncs+C/woM28GbgX+OSKu6vGsktRXTc4QbwEWMvNUZr4GPArcvmZNAu+MiADeAfwSON/TSSWpz5oEcS/w0qrbS3Vsta8AHwBeBuaBz2fmG2s/UUS0I2IuIubOnDmzxZElqT+aBDHWOZZrbn8CeA74XeCPgK9ExO/8vwdldjKzlZmtPXv2bHJUSeqvJkFcAq5fdfs6Vs4EV7sbeCxXLAA/BW7qzYiStD2aBPFZ4MaIuKHeKLkDeHzNmkXgAEBEvAd4P3Cql4NKUr/t2mhBZp6PiHuBJ4GdwMOZeTIi7qn7jwJfBB6JiHlWvsW+LzNf7ePcktRzGwYRIDOfAJ5Yc+zoqo9fBv6qt6NJ0vbyShVJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiEDVUuvNdpo5MseOBHUwdmaI73x30SBohuwY9gNRUd75L+1ib5XPLAJw+e5r2sTYA0/unBzmaRoRniBoaM8dnLsTwTcvnlpk5PjOgiTRqDKKGxuLZxU0dlzbLIGpoTIxPbOq4tFkGUUNj9sAsY7vHLjo2tnuM2QOzA5pIo8YgamhM75+mc7jD5PgkQTA5PknncMc3VNQzkZkDeeJWq5Vzc3MDeW5JoysiTmRmayuP9QxRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKTSKIgRcTAiXoyIhYi4/xJrbo2I5yLiZER8r7djSlL/bfgbsyNiJ/Ag8JfAEvBsRDyemT9ateZq4KvAwcxcjIh392leSeqbJmeItwALmXkqM18DHgVuX7Pm08BjmbkIkJmv9HZMSeq/JkHcC7y06vZSHVvtfcA1EfHdiDgREXf1akBJ2i5N/shUrHNs7e8M2wV8GDgAvA34fkQ8k5k/uegTRbSBNsDEhL/lWNKVpckZ4hJw/arb1wEvr7PmO5n568x8FXgKuHntJ8rMTma2MrO1Z8+erc4sSX3RJIjPAjdGxA0RcRVwB/D4mjX/CvxpROyKiDHgj4Ef93ZUSeqvDb9lzszzEXEv8CSwE3g4M09GxD11/9HM/HFEfAd4HngDeCgzX+jn4JLUa/4JAUkjxT8hIEk9YBAlqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESRfpzneZOjLFjgd2MHVkiu58d9AjbZsmvw9R0m+J7nyX9rE2y+eWATh99jTtY20ApvdPD3K0beEZoqQLZo7PXIjhm5bPLTNzfGZAE20vgyjpgsWzi5s6PmoMoqQLJsbX/9Melzo+agyipAtmD8wytnvsomNju8eYPTA7oIm2l0GUdMH0/mk6hztMjk8SBJPjk3QOd34r3lABf2O2pBHjb8yWpB4wiJJUDKIkFYMoScUgSlIxiJJUDKIkFYMoScUgSlIxiJJUDKIkFYMoScUgSlIxiJJUDKIkFYMoScUgSlIxiJJUDKIkFYMoScUgSlIxiJJUDKIkFYMoScUgSlIxiJJUDKIklUZBjIiDEfFiRCxExP1vse4jEfF6RHyqdyNK0vbYMIgRsRN4EDgE7APujIh9l1j3JeDJXg8pSduhyRniLcBCZp7KzNeAR4Hb11n3OeCbwCs9nE+Stk2TIO4FXlp1e6mOXRARe4FPAkd7N5okba8mQYx1juWa20eA+zLz9bf8RBHtiJiLiLkzZ840HFGStseuBmuWgOtX3b4OeHnNmhbwaEQAXAvcFhHnM/NbqxdlZgfoALRarbVRlaSBahLEZ4EbI+IG4GfAHcCnVy/IzBve/DgiHgH+bW0MJelKt2EQM/N8RNzLyrvHO4GHM/NkRNxT9/u6oaSR0OQMkcx8AnhizbF1Q5iZf3v5Y0nS9vNKFUkqIxvE7nyXqSNT7HhgB1NHpujOdwc9kqQrXKNvmYdNd75L+1ib5XPLAJw+e5r2sTYA0/unBzmapCvYSJ4hzhyfuRDDNy2fW2bm+MyAJpI0DEYyiItnFzd1XJJgRIM4MT6xqeOSBCMaxNkDs4ztHrvo2NjuMWYPzA5oIknDYCSDOL1/ms7hDpPjkwTB5PgkncMd31CR9JYiczCXFLdarZybmxvIc0saXRFxIjNbW3nsSJ4hStJWGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpGERJKgZRkopBlKRiECWpNApiRByMiBcjYiEi7l/n/umIeL7+PR0RN/d+VEnqrw2DGBE7gQeBQ8A+4M6I2Ldm2U+BP8/MDwJfBDq9HlSS+q3JGeItwEJmnsrM14BHgdtXL8jMpzPzV3XzGeC63o4pSf3XJIh7gZdW3V6qY5fyGeDb690REe2ImIuIuTNnzjSfUpK2QZMgxjrHct2FER9nJYj3rXd/ZnYys5WZrT179jSfUpK2wa4Ga5aA61fdvg54ee2iiPgg8BBwKDN/0ZvxJGn7NDlDfBa4MSJuiIirgDuAx1cviIgJ4DHgbzLzJ70fU5L6b8MzxMw8HxH3Ak8CO4GHM/NkRNxT9x8FvgC8C/hqRACcz8xW/8aWpN6LzHVfDuy7VquVc3NzA3luSaMrIk5s9YTMK1UkqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESSoGUeqj7nyXqSNT7HhgB1NHpujOdwc9kt5Ckz9DKmkLuvNd2sfaLJ9bBuD02dO0j7UBmN4/PcjRdAmeIUp9MnN85kIM37R8bpmZ4zMDmkgbMYhSnyyeXdzUcQ2eQZT6ZGJ8YlPHNXgGUeqT2QOzjO0eu+jY2O4xZg/MDmgibcQgSn0yvX+azuEOk+OTBMHk+CSdwx3fULmCRWYO5IlbrVbOzc0N5Lklja6IOJGZra081jNESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSkYhAlqRhESSoGUZKKQZSk0iiIEXEwIl6MiIWIuH+d+yMivlz3Px8RH+r9qJLUXxsGMSJ2Ag8Ch4B9wJ0RsW/NskPAjfWvDXytx3NKUt81OUO8BVjIzFOZ+RrwKHD7mjW3A1/PFc8AV0fEe3s8qyT1VZMg7gVeWnV7qY5tdo0kXdF2NVgT6xxb+8ecm6whItqsfEsN8L8R8UKD5x9W1wKvDnqIPnJ/w2uU9wbw/q0+sEkQl4DrV92+Dnh5C2vIzA7QAYiIua3+Melh4P6G2yjvb5T3Biv72+pjm3zL/CxwY0TcEBFXAXcAj69Z8zhwV73b/FHgbGb+fKtDSdIgbHiGmJnnI+Je4ElgJ/BwZp6MiHvq/qPAE8BtwAKwDNzdv5ElqT+afMtMZj7BSvRWHzu66uMEPrvJ5+5scv2wcX/DbZT3N8p7g8vYX6y0TJLkpXuSVPoexFG/7K/B/qZrX89HxNMRcfMg5tyKjfa2at1HIuL1iPjUds53uZrsLyJujYjnIuJkRHxvu2e8HA3+2xyPiGMR8cPa39C89h8RD0fEK5f60b0tdyUz+/aPlTdh/hv4PeAq4IfAvjVrbgO+zcrPMn4U+K9+zjSA/X0MuKY+PjQs+2uyt1Xr/oOV15g/Nei5e/y1uxr4ETBRt9896Ll7vL+/B75UH+8BfglcNejZG+7vz4APAS9c4v4tdaXfZ4ijftnfhvvLzKcz81d18xlWfkZzGDT52gF8Dvgm8Mp2DtcDTfb3aeCxzFwEyMxh2mOT/SXwzogI4B2sBPH89o65NZn5FCvzXsqWutLvII76ZX+bnf0zrPxfaxhsuLeI2At8EjjK8GnytXsfcE1EfDciTkTEXds23eVrsr+vAB9g5SKKeeDzmfnG9ozXd1vqSqMfu7kMPbvs7wrVePaI+DgrQfyTvk7UO032dgS4LzNfXznJGCpN9rcL+DBwAHgb8P2IeCYzf9Lv4Xqgyf4+ATwH/AXw+8C/R8R/Zub/9Hm27bClrvQ7iD277O8K1Wj2iPgg8BBwKDN/sU2zXa4me2sBj1YMrwVui4jzmfmtbZnw8jT9b/PVzPw18OuIeAq4GRiGIDbZ393AP+XKi24LEfFT4CbgB9szYl9trSt9fuFzF3AKuIHfvLD7B2vW/DUXv/j5g0G/YNvj/U2wcgXPxwY9b6/3tmb9IwzXmypNvnYfAI7X2jHgBeAPBz17D/f3NeAf6+P3AD8Drh307JvY4xSXflNlS13p6xlijvhlfw339wXgXcBX60zqfA7BhfUN9za0muwvM38cEd8BngfeAB7KzKH4DU0Nv35fBB6JiHlWwnFfZg7Fb8GJiG8AtwLXRsQS8A/Abri8rniliiQVr1SRpGIQJakYREkqBlGSikGUpGIQJakYREkqBlGSyv8BlxPi1ER9Rh0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# from youtube vid\n",
    "\n",
    "def latin_hypercube_2d_uniform(n):\n",
    "    lower_limits = np.arange(0, n)/n\n",
    "    upper_limits = np.arange(1, n+1)/n\n",
    "    \n",
    "    points = np.random.uniform(low=lower_limits, high=upper_limits, \n",
    "                               size=[2, n]).T\n",
    "    np.random.shuffle(points[:,1])\n",
    "    \n",
    "    return points\n",
    "\n",
    "n = 5\n",
    "p = latin_hypercube_2d_uniform(n)\n",
    "\n",
    "plt.figure(figsize=[5,5])\n",
    "plt.xlim([0,1])\n",
    "plt.ylim([0,1])\n",
    "plt.scatter(p[:,0], p[:,1], c='g')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "pacific-mercy",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[0.01629478, 0.01054196],\n",
       "       [0.23866785, 0.70761893],\n",
       "       [0.59014894, 0.54343629],\n",
       "       [0.68173991, 0.9017412 ],\n",
       "       [0.89333663, 0.20751325]])"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p"
   ]
  }
 ],
 "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.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}