EMsFEA-net/visualization.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "76d6a02c",
   "metadata": {},
   "source": [
    "## test.py可视化"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "aef142c9-a3d0-408c-a87f-695167d49dea",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(61, 181, 2)\n",
      "(60, 180)\n",
      "(432, 25)\n",
      "(432, 4, 2)\n",
      "(432, 72)\n",
      "Total loss: 0.03841373324394226\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAHFCAYAAADcytJ5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA6hUlEQVR4nO3deXxU1f3/8fckhCRkY5MlEBZFQYQgi0hEAQWpVFFUvq6ttmoVRdQf7deq1LpUG5e6tQriUtC2irWIxW9BDSqLomUTZRNcEohC2EnCknXO7w/MMHcyN0wWcu/cvJ6PRx4mN3dmPjf3++15c8655/iMMUYAAAAuFON0AQAAAHYIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKoCHzJw5Uz6fTytWrLA9Jy8vTz6fTzNnzmy8wpqAqr/rn/70J6dLATylmdMFAGhcHTt21KeffqoTTjjB6VIA4KgIKkATEx8fryFDhjhdRo0OHjyoFi1aHJP3PnTokBITE4/JewNoeAz9AE1MuKGf+++/Xz6fT+vWrdOVV16ptLQ0tW/fXtddd50KCwstrzfGaOrUqTr11FOVmJioVq1aafz48fruu+8s5+Xk5Oiiiy5S586dlZCQoB49euimm27Srl27LOdVffaqVas0fvx4tWrV6qi9PR9//LGysrKUkJCgTp066d5779VLL70kn8+nvLy8wHndunXTBRdcoLfeekv9+/dXQkKCHnjgAUnSc889p2HDhqldu3ZKSkpS37599dhjj6m8vNzyWSNGjFCfPn20ZMkSDRkyRImJiYHPrKysDFvfk08+qe7duys5OVlZWVn67LPParweAPboUQEQcOmll+ryyy/X9ddfrzVr1ujuu++WJP31r38NnHPTTTdp5syZuu222/Too49qz549evDBB3XGGWfoiy++UPv27SVJ3377rbKysnTDDTcoLS1NeXl5evLJJ3XmmWdqzZo1iouLs3z2JZdcoiuuuEITJkzQgQMHbGv88ssvde655+qkk07SK6+8ohYtWuj555/X3//+97Dnr1q1Shs2bNDvfvc7de/eXUlJSYH6rrrqKnXv3l3NmzfXF198oYcfflhfffWV5XolqaCgQFdccYXuuusuPfjgg/rPf/6jhx56SHv37tWzzz5rOfe5555Tr1699PTTT0uS7r33Xv30pz9Vbm6u0tLSIrgLACwMAM+YMWOGkWSWL19ue05ubq6RZGbMmBE4dt999xlJ5rHHHrOce8stt5iEhATj9/uNMcZ8+umnRpJ54oknLOfl5+ebxMREc+edd4b9TL/fb8rLy83mzZuNJPPvf/+72mf//ve/j+ga/+d//sckJSWZnTt3Bo5VVlaa3r17G0kmNzc3cLxr164mNjbWbNy4scb3rKysNOXl5ebVV181sbGxZs+ePYHfDR8+vFrNxhjzq1/9ysTExJjNmzcbY478Xfv27WsqKioC5y1btsxIMq+//npE1wfAiqEfAAEXXnih5efMzEyVlJRox44dkqT/+7//k8/n089+9jNVVFQEvjp06KB+/fpp4cKFgdfu2LFDEyZMUEZGhpo1a6a4uDh17dpVkrRhw4Zqn33ppZdGVOOiRYt0zjnnqG3btoFjMTExuuyyy8Ken5mZqZNOOqna8c8//1wXXnih2rRpo9jYWMXFxemaa65RZWWlNm3aZDk3JSWl2t/mqquukt/v1+LFiy3Hzz//fMXGxlo+X5I2b94c0fUBsGLoB0BAmzZtLD/Hx8dLOjwBVZK2b98uY0xgeCfU8ccfL0ny+/0aPXq0tm7dqnvvvVd9+/ZVUlKS/H6/hgwZEni/YB07doyoxt27d4f9fLuawr3vli1bdNZZZ6lnz5565pln1K1bNyUkJGjZsmWaOHFitfrCvXeHDh0C9QQ72t8QQO0QVABErG3btvL5fFqyZEmgAQ5WdWzt2rX64osvNHPmTF177bWB33/zzTe27+3z+SKqoU2bNtq+fXu14wUFBRG/79tvv60DBw7orbfeCvTySNLq1avDvkdNnxcaTAA0LIZ+AETsggsukDFGP/zwgwYNGlTtq2/fvpKOhIPQMDN9+vR61zB8+HB9+OGHlqeH/H6/3nzzzYjfI1x9xhi9+OKLYc8vLi7W3LlzLcdee+01xcTEaNiwYbUpH0At0aMCeNCHH35oeUy3yk9/+tN6ve/QoUN144036pe//KVWrFihYcOGKSkpSdu2bdPHH3+svn376uabb1avXr10wgkn6K677pIxRq1bt9Y777yjnJycen2+JE2ZMkXvvPOORo4cqSlTpigxMVHPP/984EmhmJij//vr3HPPVfPmzXXllVfqzjvvVElJiaZNm6a9e/eGPb9Nmza6+eabtWXLFp100kmaN2+eXnzxRd18883q0qVLva8JgD2CCuBBv/3tb8Mez83Nrfd7T58+XUOGDNH06dM1depU+f1+paena+jQoRo8eLAkKS4uTu+8845uv/123XTTTWrWrJlGjRqlBQsW1Lth79evn3JycvSb3/xG11xzjVq1aqWf//znGj58uH77299G9Ahwr169NHv2bP3ud7/TJZdcojZt2uiqq67S5MmTNWbMmGrnd+jQQc8995x+85vfaM2aNWrdurXuueeewJosAI4dnzHGOF0EANTX6NGjlZeXV+2JnfoaMWKEdu3apbVr1zbo+wKIDD0qAKLO5MmT1b9/f2VkZGjPnj36xz/+oZycHL388stOlwaggRFUAESdyspK/f73v1dBQYF8Pp969+6tv/3tb/rZz37mdGkAGhhDPwAAwLV4PBkAALgWQQUAALgWQQUAALhWVE+m9fv92rp1q1JSUiJefhsAADjLGKPi4mKlp6cfdZHGqA4qW7duVUZGhtNlAACAOsjPz1fnzp1rPCeqg0pKSoqkwxeamprqcDUAACASRUVFysjICLTjNYnqoFI13JOamkpQAQAgykQybYPJtAAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKgAAwLUIKjU4VFYpY4zTZQAA0GQRVGxs2X1Q/f/wvu6bu87pUgAAaLIIKjY2bS9WSblfX3xf6HQpAAA0WQQVG4EBH4Z+AABwDEHFRtXcFGIKAADOIajYqAoodKgAAOAcgoqNqoBi6FMBAMAxBBUAAOBaBBVbP85RoUMFAADHEFRsBIZ+CCoAADiGoGLDhPwXAAA0PoKKjSM9KkQVAACcQlABAACuRVCxwWPJAAA4j6Big8m0AAA4j6Bi48hkWpIKAABOIajYCOz1Q04BAMAxBBUAAOBaBJWjoEMFAADnEFRssI4KAADOI6jYqJpES0wBAMA5BBUbhjX0AQBwHEEFAAC4FkHFRmCOirNlAADQpBFUbARGfphMCwCAYwgqNgILvjlcBwAATRlBxcaRHhVHywAAoEkjqAAAANciqNgJTKalSwUAAKcQVGwEFnwjpwAA4BiCio0jS+g7WwcAAE0ZQcUG+QQAAOcRVAAAgGsRVGywezIAAM4jqNhg92QAAJxHULHBZFoAAJznmqCSnZ0tn8+nO+64w+lSJAWtTEufCgAAjnFFUFm+fLleeOEFZWZmOl0KAABwEceDyv79+3X11VfrxRdfVKtWrZwu5wjDgm8AADjN8aAyceJEnX/++Ro1atRRzy0tLVVRUZHl61gxIf8FAACNr5mTHz5r1iytWrVKy5cvj+j87OxsPfDAA8e4qsOYTAsAgPMc61HJz8/X7bffrr///e9KSEiI6DV33323CgsLA1/5+fnHuEqJPhUAAJzjWI/KypUrtWPHDg0cODBwrLKyUosXL9azzz6r0tJSxcbGWl4THx+v+Pj4RqmPhd4AAHCeY0Fl5MiRWrNmjeXYL3/5S/Xq1Uu//e1vq4WUxhaYo0JeAQDAMY4FlZSUFPXp08dyLCkpSW3atKl23AmBOSrOlgEAQJPm+FM/bnWkR4WoAgCAUxx96ifUwoULnS6hGmIKAADOoUfFBj0pAAA4j6B
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABIsAAAG4CAYAAADFWWfKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAsXklEQVR4nO3de5zVdZ348fdhgMPFYeQSAxMDTiqignjBC6SBImwohtqj1TRDq11dMSW6IJqCho6RGe2SuuaGmKLWI2+tZo0JmA+1wLuWmCsILiLpxgyijVy+vz/68dkduYmema+Mz+fjcR5yvuc738/nDJ/5gi++55xClmVZAAAAAEBEtMl7AgAAAAB8eIhFAAAAACRiEQAAAACJWAQAAABAIhYBAAAAkIhFAAAAACRiEQAAAACJWAQAAABAIhYBAAAAkIhFANCK/P73v48TTjgh+vbtG8ViMSorK2Po0KHx9a9/vcl+u+22W4wdOzanWW7b/Pnzo1AoxPz58/OeSrOaNm1aFAqFXMYeMWJEjBgxIpext2W33XaL008/fbv7fVTWCADkpW3eEwAASuOee+6Jz3zmMzFixIiYMWNG9O7dO1599dVYtGhR3HrrrfH9738/7ynyIXH11VfnPYUtuuOOO6JLly55TwMAPvLEIgBoJWbMmBE1NTXx61//Otq2/d8/4k8++eSYMWNGjjPjw2afffbJewpNvP3229GxY8c44IAD8p4KABBehgYArcYbb7wRPXr0aBKKNmnTZst/5N93331x4IEHRseOHWPAgAHxk5/8ZLN9nn322Rg3blx07do1OnToEPvvv3/MmTMnPZ5lWVRWVsaECRPStg0bNkTXrl2jTZs28dprr6XtV111VbRt2zZWr169w8/v7rvvjqFDh0anTp2ivLw8Ro0aFY888kiTff7yl7/EP//zP0d1dXUUi8X42Mc+Fp/85Cfj/vvvT/s88cQTMXbs2OjZs2cUi8WoqqqKY489Nl555ZWtjj1x4sTo3LlzNDQ0bPbYSSedFJWVlbFu3bqIiLjtttti9OjR0bt37+jYsWPsvffecf7558fatWu3+xwLhUJMmzZts+1bennWypUr48wzz4w+ffpE+/bto6amJi655JJYv379dsd598vQli5dGoVCIa688sq46qqroqamJnbZZZcYOnRoPProo9s9XkTEQw89FEOHDo0OHTrExz/+8bjooovi+uuvj0KhEEuXLm3yXMaOHRu33357HHDAAdGhQ4e45JJLtvo8n3/++fj0pz8dnTp1ih49esRZZ50Va9aseU9zAgDeH1cWAUArMXTo0Lj++uvj3HPPjVNPPTUOPPDAaNeu3Vb3f+qpp+LrX/96nH/++VFZWRnXX399fPnLX4499tgjPvWpT0VExOLFi2PYsGHRs2fP+Nd//dfo3r173HTTTXH66afHa6+9Ft/61reiUCjEUUcd1STILFq0KFavXh0dO3aM3/72t3HKKadERMT9998fBx10UOy666479Nzmzp0bp556aowePTpuueWWaGxsjBkzZsSIESPit7/9bRx++OEREXHaaafF448/Hpdddln0798/Vq9eHY8//ni88cYbERGxdu3aGDVqVNTU1MSPfvSjqKysjJUrV8a8efO2GSC+9KUvxQ9/+MP42c9+Fl/5ylfS9tWrV8ddd90VEyZMSN/rP//5z3HMMcekwPT888/Hd7/73fjDH/4QDzzwwA49761ZuXJlHHLIIdGmTZu4+OKLY/fdd49HHnkkpk+fHkuXLo3Zs2e/r+P+6Ec/igEDBsTMmTMjIuKiiy6KY445JpYsWRIVFRVb/bqnn346Ro0aFf379485c+ZEp06d4tprr42bbrppi/s//vjj8ac//Sm+/e1vR01NTXTu3HmL+7322msxfPjwaNeuXVx99dVRWVkZN998c5xzzjnv6/kBAO9RBgC0Cq+//np2+OGHZxGRRUTWrl27bNiwYVltbW22Zs2aJvv269cv69ChQ/byyy+nbW+//XbWrVu37Mwzz0zbTj755KxYLGbLli1r8vVjxozJOnXqlK1evTrLsiy7/vrrs4hI+02fPj0bMGBA9pnPfCY744wzsizLsnfeeSfr3LlzdsEFF2zzecybNy+LiGzevHlZlmXZhg0bsqqqqmzQoEHZhg0b0n5r1qzJevbsmQ0bNixt22WXXbKJEydu9diLFi3KIiK78847tzmHLTnwwAObjJVlWXb11VdnEZE988wzW/yajRs3ZuvWrcsWLFiQRUT21FNPpcemTp2avfuvYhGRTZ06dbPj9OvXLxs/fny6f+aZZ2a77LJLk9+/LMuyK6+8MouI7Lnnntvmcxk+fHg2fPjwdH/JkiVZRGSDBg3K1q9fn7b/4Q9/yCIiu+WWW7Z5vM997nNZ586ds7/85S9p24YNG7J99tkni4hsyZIlTZ5LWVlZtnjx4u0+z8mTJ2eFQiF78sknm+w3atSoJmsEACgtL0MDgFaie/fu8bvf/S4WLlwYV1xxRYwbNy5eeOGFmDJlSgwaNChef/31Jvvvv//+0bdv33S/Q4cO0b9//3j55ZfTtgceeCBGjhwZ1dXVTb729NNPj7feeiu9DOzoo4+OiEhXF9XV1cWoUaPi6KOPjrq6uoiIeOSRR2Lt2rVp3/dq8eLFsWLFijjttNOavJxul112ic9+9rPx6KOPxltvvRUREYccckjccMMNMX369Hj00UfTS8M22WOPPaJr164xefLkuPbaa+OPf/zje57HGWecEQ8//HAsXrw4bZs9e3YcfPDBMXDgwLTtpZdeilNOOSV69eoVZWVl0a5duxg+fHhERPzpT3/aoee+Nf/5n/8ZRx55ZFRVVcX69evTbcyYMRERsWDBgvd13GOPPTbKysrS/f322y8iosma2JIFCxbEUUcdFT169Ejb2rRpE//4j/+4xf3322+/6N+//3bnM2/evNh3331j8ODBTbZvulINAGgeYhEAtDJDhgyJyZMnx89//vNYsWJFfO1rX4ulS5du9ibX3bt33+xri8VivP322+n+G2+8Eb17995sv6qqqvR4RES/fv1i9913j/vvvz9FpE2x6JVXXonFixfH/fffHx07doxhw4bt0PPZNMbW5rFx48b461//GhF/f7+g8ePHx/XXXx9Dhw6Nbt26xRe/+MVYuXJlRERUVFTEggULYv/9948LLrgg9t1336iqqoqpU6duFpbe7dRTT41isRg33HBDRET88Y9/jIULF8YZZ5yR9nnzzTfjiCOOiN///vcxffr0mD9/fixcuDBuv/32iIgm39sP4rXXXotf/vKX0a5duya3fffdNyJiszD4Xr17TRSLxYjY/rzfeOONqKys3Gz7lrZFbPn3cmvH7dWr12bbt7QNACgd71kEAK1Yu3btYurUqfGDH/wgnn322R3++u7du8err7662fYVK1ZERDS5kmTkyJFx1113xYIFC2Ljxo0xYsSIKC8vj6qqqqirq4v7778/jjjiiBQgdmQOEbHVebRp0ya6du2a5jNz5syYOXNmLFu2LO6+++44//zzY9WqVXHfffdFRMSgQYPi1ltvjSzL4umnn44bbrghLr300ujYsWOcf/75W51H165dY9y4cXHjjTfG9OnTY/bs2dGhQ4f4/Oc/n/Z54IEHYsWKFTF//vx0NVFEvOc39C4Wi9HY2LjZ9k3BbJMePXrEfvvtF5dddtkWj7Mp5rWU7t27N3kj8002Rbp3KxQK7/m4WzrG1o4LAJSGK4sAoJXYUkyJ+N+XPr2fgDBy5MgUQP6vG2+8MTp16hSHHXZY2nb00UfHa6+9FjNnzozDDjssysvL0zHuuOOOWLhw4Q6/BC0iYq+99oqPf/zjMXfu3MiyLG1fu3Zt/OIXv0ifkPZuffv2jXPOOSdGjRoVjz/++GaPFwqFGDx4cPzgBz+IXXfddYv7vNsZZ5wRK1asiHvvvTduuummOOGEE5q8WfemCPLuIPbv//7v7+m57rbbbvH000832fbAAw/Em2++2WTb2LFj49lnn43dd989hgwZstmtpWPR8OHD44EHHmhyRdPGjRvj5z//+Qc
      "text/plain": [
       "<Figure size 1440x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from test import test\n",
    "\n",
    "import numpy as np\n",
    "\n",
    "import torch\n",
    "import torch.nn as nn\n",
    "import torch.nn.functional as F\n",
    "\n",
    "from utils.data_standardizer import standardization\n",
    "from utils.data_loader import data_loader\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# Load datasets\n",
    "# test data select:\n",
    "dataload_mod='mod3' # opt: mod1 mod2 mod3\n",
    "# pretrained model select:\n",
    "pretrained_mod='mod3' # opt: mod1 mod2 mod3\n",
    "\n",
    "dst_path='datasets/top88_'+ dataload_mod + '_xPhys_180_60.npy'\n",
    "U_path='datasets/top88_'+ dataload_mod + '_u_180_60.npy'\n",
    "global_density, global_displace, coarse_density, coarse_displace, fine_displace = data_loader(dst_path, U_path)\n",
    "X = np.hstack((coarse_density[:,:] , coarse_displace[:,:,0] , coarse_displace[:,:,1]))\n",
    "Y = fine_displace[:,:]\n",
    "\n",
    "# Set loss function\n",
    "loss_function = nn.MSELoss()\n",
    "\n",
    "# Predict \n",
    "pred = test('checkpoints/ANN_' + pretrained_mod + '_opt.pt', X)\n",
    "\n",
    "# Calculate loss\n",
    "pred_loss=[]\n",
    "device = f'cuda:{0}' if torch.cuda.is_available() else 'cpu'\n",
    "Y = torch.from_numpy(Y).type(torch.float32).to(device)\n",
    "for i in range(pred.shape[0]):\n",
    "    pred_loss.append(loss_function(pred[i,:],Y[i,:]).item())\n",
    "\n",
    "print('Total loss: '+ str(loss_function(pred,Y).item()))\n",
    "\n",
    "# Plot \n",
    "plt.plot(range(pred.shape[0]),pred_loss)\n",
    "plt.ylabel('Loss')\n",
    "plt.xlabel('Coarse mesh id')\n",
    "plt.title(\"Linear graph\") \n",
    "plt.show()\n",
    "\n",
    "loss_metrix = np.asarray(pred_loss)\n",
    "loss_metrix = loss_metrix.reshape(int(60/5), int(180/5))\n",
    "plt.matshow(loss_metrix)\n",
    "plt.title(\"Show loss value in grid\") \n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "00fa1673",
   "metadata": {},
   "source": [
    "## 拓扑优化 88"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "5afed322",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Minimum compliance problem with OC\n",
      "ndes: 36 x 12\n",
      "volfrac: 0.4, rmin: 5.4, penal: 3.0\n",
      "Filter method: Density based\n",
      "it.: 1 , obj.: 1935.489 Vol.: 0.400, ch.: 0.200\n",
      "it.: 2 , obj.: 1380.004 Vol.: 0.400, ch.: 0.200\n",
      "it.: 3 , obj.: 1134.169 Vol.: 0.400, ch.: 0.200\n",
      "it.: 4 , obj.: 1006.889 Vol.: 0.400, ch.: 0.200\n",
      "it.: 5 , obj.: 919.134 Vol.: 0.400, ch.: 0.176\n",
      "it.: 6 , obj.: 879.596 Vol.: 0.400, ch.: 0.143\n",
      "it.: 7 , obj.: 858.187 Vol.: 0.400, ch.: 0.106\n",
      "it.: 8 , obj.: 848.604 Vol.: 0.400, ch.: 0.082\n",
      "it.: 9 , obj.: 842.698 Vol.: 0.400, ch.: 0.083\n",
      "it.: 10 , obj.: 837.788 Vol.: 0.400, ch.: 0.058\n",
      "it.: 11 , obj.: 833.111 Vol.: 0.400, ch.: 0.065\n",
      "it.: 12 , obj.: 830.436 Vol.: 0.400, ch.: 0.054\n",
      "it.: 13 , obj.: 827.623 Vol.: 0.400, ch.: 0.050\n",
      "it.: 14 , obj.: 825.488 Vol.: 0.400, ch.: 0.054\n",
      "it.: 15 , obj.: 823.347 Vol.: 0.400, ch.: 0.044\n",
      "it.: 16 , obj.: 820.603 Vol.: 0.400, ch.: 0.050\n",
      "it.: 17 , obj.: 819.616 Vol.: 0.400, ch.: 0.055\n",
      "it.: 18 , obj.: 817.396 Vol.: 0.400, ch.: 0.058\n",
      "it.: 19 , obj.: 815.703 Vol.: 0.400, ch.: 0.040\n",
      "it.: 20 , obj.: 814.336 Vol.: 0.400, ch.: 0.030\n",
      "it.: 21 , obj.: 814.397 Vol.: 0.400, ch.: 0.031\n",
      "it.: 22 , obj.: 812.732 Vol.: 0.400, ch.: 0.033\n",
      "it.: 23 , obj.: 811.633 Vol.: 0.400, ch.: 0.034\n",
      "it.: 24 , obj.: 811.338 Vol.: 0.400, ch.: 0.036\n",
      "it.: 25 , obj.: 809.472 Vol.: 0.400, ch.: 0.036\n",
      "it.: 26 , obj.: 808.204 Vol.: 0.400, ch.: 0.036\n",
      "it.: 27 , obj.: 807.856 Vol.: 0.400, ch.: 0.037\n",
      "it.: 28 , obj.: 807.542 Vol.: 0.400, ch.: 0.038\n",
      "it.: 29 , obj.: 807.265 Vol.: 0.400, ch.: 0.040\n",
      "it.: 30 , obj.: 805.490 Vol.: 0.400, ch.: 0.039\n",
      "it.: 31 , obj.: 805.881 Vol.: 0.400, ch.: 0.040\n",
      "it.: 32 , obj.: 805.291 Vol.: 0.400, ch.: 0.023\n",
      "it.: 33 , obj.: 804.307 Vol.: 0.400, ch.: 0.023\n",
      "it.: 34 , obj.: 804.148 Vol.: 0.400, ch.: 0.025\n",
      "it.: 35 , obj.: 802.469 Vol.: 0.400, ch.: 0.023\n",
      "it.: 36 , obj.: 802.718 Vol.: 0.400, ch.: 0.023\n",
      "it.: 37 , obj.: 802.193 Vol.: 0.400, ch.: 0.021\n",
      "it.: 38 , obj.: 801.018 Vol.: 0.400, ch.: 0.022\n",
      "it.: 39 , obj.: 800.864 Vol.: 0.400, ch.: 0.016\n",
      "it.: 40 , obj.: 800.667 Vol.: 0.400, ch.: 0.016\n",
      "it.: 41 , obj.: 800.441 Vol.: 0.400, ch.: 0.017\n",
      "it.: 42 , obj.: 800.184 Vol.: 0.400, ch.: 0.017\n",
      "it.: 43 , obj.: 799.948 Vol.: 0.400, ch.: 0.017\n",
      "it.: 44 , obj.: 799.686 Vol.: 0.400, ch.: 0.017\n",
      "it.: 45 , obj.: 799.443 Vol.: 0.400, ch.: 0.017\n",
      "it.: 46 , obj.: 799.208 Vol.: 0.400, ch.: 0.018\n",
      "it.: 47 , obj.: 797.672 Vol.: 0.400, ch.: 0.016\n",
      "it.: 48 , obj.: 797.913 Vol.: 0.400, ch.: 0.014\n",
      "it.: 49 , obj.: 797.863 Vol.: 0.400, ch.: 0.015\n",
      "it.: 50 , obj.: 797.823 Vol.: 0.400, ch.: 0.015\n",
      "it.: 51 , obj.: 797.734 Vol.: 0.400, ch.: 0.016\n",
      "it.: 52 , obj.: 796.372 Vol.: 0.400, ch.: 0.016\n",
      "it.: 53 , obj.: 796.510 Vol.: 0.400, ch.: 0.016\n",
      "it.: 54 , obj.: 796.227 Vol.: 0.400, ch.: 0.016\n",
      "it.: 55 , obj.: 796.119 Vol.: 0.400, ch.: 0.016\n",
      "it.: 56 , obj.: 795.953 Vol.: 0.400, ch.: 0.017\n",
      "it.: 57 , obj.: 795.817 Vol.: 0.400, ch.: 0.017\n",
      "it.: 58 , obj.: 795.677 Vol.: 0.400, ch.: 0.017\n",
      "it.: 59 , obj.: 795.544 Vol.: 0.400, ch.: 0.017\n",
      "it.: 60 , obj.: 795.414 Vol.: 0.400, ch.: 0.018\n",
      "it.: 61 , obj.: 794.068 Vol.: 0.400, ch.: 0.017\n",
      "it.: 62 , obj.: 794.229 Vol.: 0.400, ch.: 0.018\n",
      "it.: 63 , obj.: 793.962 Vol.: 0.400, ch.: 0.018\n",
      "it.: 64 , obj.: 793.879 Vol.: 0.400, ch.: 0.018\n",
      "it.: 65 , obj.: 793.749 Vol.: 0.400, ch.: 0.018\n",
      "it.: 66 , obj.: 793.755 Vol.: 0.400, ch.: 0.014\n",
      "it.: 67 , obj.: 793.967 Vol.: 0.400, ch.: 0.017\n",
      "it.: 68 , obj.: 792.700 Vol.: 0.400, ch.: 0.016\n",
      "it.: 69 , obj.: 792.880 Vol.: 0.400, ch.: 0.017\n",
      "it.: 70 , obj.: 792.649 Vol.: 0.400, ch.: 0.017\n",
      "it.: 71 , obj.: 792.594 Vol.: 0.400, ch.: 0.013\n",
      "it.: 72 , obj.: 792.555 Vol.: 0.400, ch.: 0.012\n",
      "it.: 73 , obj.: 792.663 Vol.: 0.400, ch.: 0.013\n",
      "it.: 74 , obj.: 792.578 Vol.: 0.400, ch.: 0.013\n",
      "it.: 75 , obj.: 792.538 Vol.: 0.400, ch.: 0.013\n",
      "it.: 76 , obj.: 792.492 Vol.: 0.400, ch.: 0.010\n",
      "[[[ 0.00000000e+00 -7.92492303e+02]\n",
      "  [ 0.00000000e+00 -7.90881464e+02]\n",
      "  [ 0.00000000e+00 -7.89586778e+02]\n",
      "  [ 0.00000000e+00 -7.87987744e+02]\n",
      "  [ 0.00000000e+00 -7.85490451e+02]\n",
      "  [ 0.00000000e+00 -7.80808387e+02]\n",
      "  [ 0.00000000e+00 -7.72183290e+02]\n",
      "  [ 0.00000000e+00 -7.63786319e+02]\n",
      "  [ 0.00000000e+00 -7.60546489e+02]\n",
      "  [ 0.00000000e+00 -7.59546951e+02]\n",
      "  [ 0.00000000e+00 -7.58873647e+02]\n",
      "  [ 0.00000000e+00 -7.57922538e+02]\n",
      "  [ 0.00000000e+00 -7.56516252e+02]]\n",
      "\n",
      " [[-6.57491031e+00 -7.86358023e+02]\n",
      "  [-2.26233032e+00 -7.87875473e+02]\n",
      "  [ 1.07656708e-02 -7.87574712e+02]\n",
      "  [ 1.90747564e+00 -7.86320586e+02]\n",
      "  [ 3.37057042e+00 -7.83993614e+02]\n",
      "  [ 3.94807044e+00 -7.79420594e+02]\n",
      "  [ 2.38042509e+00 -7.70946045e+02]\n",
      "  [-5.64864681e-02 -7.62775011e+02]\n",
      "  [-3.53643946e-01 -7.59621836e+02]\n",
      "  [ 6.55032751e-01 -7.58638816e+02]\n",
      "  [ 2.12493867e+00 -7.57964302e+02]\n",
      "  [ 3.77189157e+00 -7.57007313e+02]\n",
      "  [ 5.54892839e+00 -7.55589191e+02]]\n",
      "\n",
      " [[-1.21895303e+01 -7.80326512e+02]\n",
      "  [-5.71967851e+00 -7.81574272e+02]\n",
      "  [-5.95012104e-01 -7.82093587e+02]\n",
      "  [ 3.44654865e+00 -7.81477367e+02]\n",
      "  [ 6.55988187e+00 -7.79597414e+02]\n",
      "  [ 7.70054669e+00 -7.75460214e+02]\n",
      "  [ 4.47153942e+00 -7.67497735e+02]\n",
      "  [-2.76687070e-01 -7.59807407e+02]\n",
      "  [-7.34615023e-01 -7.56848036e+02]\n",
      "  [ 1.33436380e+00 -7.55896003e+02]\n",
      "  [ 4.29226197e+00 -7.55214995e+02]\n",
      "  [ 7.60027877e+00 -7.54250508e+02]\n",
      "  [ 1.11534040e+01 -7.52830861e+02]]\n",
      "\n",
      " [[-1.74627851e+01 -7.72330634e+02]\n",
      "  [-9.12100085e+00 -7.73610358e+02]\n",
      "  [-1.86461780e+00 -7.74241013e+02]\n",
      "  [ 4.34273120e+00 -7.74070107e+02]\n",
      "  [ 9.17861268e+00 -7.72745345e+02]\n",
      "  [ 1.10098059e+01 -7.69319661e+02]\n",
      "  [ 6.11012281e+00 -7.62098843e+02]\n",
      "  [-7.82677462e-01 -7.54951197e+02]\n",
      "  [-1.13926407e+00 -7.52235248e+02]\n",
      "  [ 2.06667358e+00 -7.51328210e+02]\n",
      "  [ 6.48779119e+00 -7.50640637e+02]\n",
      "  [ 1.14002347e+01 -7.49680564e+02]\n",
      "  [ 1.66596718e+01 -7.48292672e+02]]\n",
      "\n",
      " [[-2.25277724e+01 -7.62643185e+02]\n",
      "  [-1.25009313e+01 -7.63885450e+02]\n",
      "  [-3.45945990e+00 -7.64602804e+02]\n",
      "  [ 4.57485768e+00 -7.64676752e+02]\n",
      "  [ 1.10463862e+01 -7.63841276e+02]\n",
      "  [ 1.36806525e+01 -7.61204653e+02]\n",
      "  [ 7.21676348e+00 -7.54903513e+02]\n",
      "  [-1.55624813e+00 -7.48305804e+02]\n",
      "  [-1.49859047e+00 -7.45861882e+02]\n",
      "  [ 2.88935857e+00 -7.45020345e+02]\n",
      "  [ 8.70323997e+00 -7.44334342e+02]\n",
      "  [ 1.51144203e+01 -7.43391295e+02]\n",
      "  [ 2.19713767e+01 -7.42055620e+02]]\n",
      "\n",
      " [[-2.73927686e+01 -7.51377061e+02]\n",
      "  [-1.58870639e+01 -7.52602601e+02]\n",
      "  [-5.33643883e+00 -7.53385118e+02]\n",
      "  [ 4.22255418e+00 -7.53646888e+02]\n",
      "  [ 1.21090682e+01 -7.53167037e+02]\n",
      "  [ 1.55596731e+01 -7.51132409e+02]\n",
      "  [ 7.66841943e+00 -7.45566534e+02]\n",
      "  [-2.45562233e+00 -7.39823201e+02]\n",
      "  [-1.69486122e+00 -7.37857831e+02]\n",
      "  [ 3.84394712e+00 -7.37117517e+02]\n",
      "  [ 1.09394109e+01 -7.36443639e+02]\n",
      "  [ 1.87121313e+01 -7.35520179e+02]\n",
      "  [ 2.70263175e+01 -7.34235593e+02]]\n",
      "\n",
      " [[-3.20820691e+01 -7.38749696e+02]\n",
      "  [-1.93369247e+01 -7.39952162e+02]\n",
      "  [-7.52670896e+00 -7.40786221e+02]\n",
      "  [ 3.29026058e+00 -7.41178249e+02]\n",
      "  [ 1.23413308e+01 -7.40917468e+02]\n",
      "  [ 1.64690506e+01 -7.39238811e+02]\n",
      "  [ 7.53471100e+00 -7.34444932e+02]\n",
      "  [-3.15174045e+00 -7.29915771e+02]\n",
      "  [-1.60200086e+00 -7.28426931e+02]\n",
      "  [ 4.98955590e+00 -7.27785998e+02]\n",
      "  [ 1.32345446e+01 -7.27125467e+02]\n",
      "  [ 2.22434017e+01 -7.26212845e+02]\n",
      "  [ 3.18833462e+01 -7.24964464e+02]]\n",
      "\n",
      " [[-3.65413224e+01 -7.25036153e+02]\n",
      "  [-2.28285892e+01 -7.26200575e+02]\n",
      "  [-1.00202659e+01 -7.27062560e+02]\n",
      "  [ 1.77937594e+00 -7.27540391e+02]\n",
      "  [ 1.17134368e+01 -7.27432372e+02]\n",
      "  [ 1.63274423e+01 -7.26069583e+02]\n",
      "  [ 7.09016414e+00 -7.22195923e+02]\n",
      "  [-3.41439190e+00 -7.18843028e+02]\n",
      "  [-1.11633735e+00 -7.17752896e+02]\n",
      "  [ 6.38287927e+00 -7.17187414e+02]\n",
      "  [ 1.56452710e+01 -7.16530477e+02]\n",
      "  [ 2.57548140e+01 -7.15620679e+02]\n",
      "  [ 3.65630692e+01 -7.14404383e+02]]\n",
      "\n",
      " [[-4.07449236e+01 -7.10523484e+02]\n",
      "  [-2.63238423e+01 -7.11644909e+02]\n",
      "  [-1.27621582e+01 -7.12519095e+02]\n",
      "  [-2.39014524e-01 -7.13052494e+02]\n",
      "  [ 1.03004388e+01 -7.13055185e+02]\n",
      "  [ 1.52735070e+01 -7.11967289e+02]\n",
      "  [ 6.52241100e+00 -7.08983621e+02]\n",
      "  [-3.15204150e+00 -7.06706560e+02]\n",
      "  [-1.92080570e-01 -7.05972817e+02]\n",
      "  [ 8.05532319e+00 -7.05480792e+02]\n",
      "  [ 1.81911531e+01 -7.04829877e+02]\n",
      "  [ 2.92576029e+01 -7.03923559e+02]\n",
      "  [ 4.10661940e+01 -7.02738624e+02]]\n",
      "\n",
      " [[-4.47184718e+01 -6.95463890e+02]\n",
      "  [-2.97985507e+01 -6.96550091e+02]\n",
      "  [-1.56783728e+01 -6.97430986e+02]\n",
      "  [-2.63399681e+00 -6.98002624e+02]\n",
      "  [ 8.29759936e+00 -6.98086806e+02]\n",
      "  [ 1.35844896e+01 -6.97242333e+02]\n",
      "  [ 5.95650612e+00 -6.95068321e+02]\n",
      "  [-2.41047562e+00 -6.93686769e+02]\n",
      "  [ 1.13771151e+00 -6.93258406e+02]\n",
      "  [ 9.98829256e+00 -6.92843169e+02]\n",
      "  [ 2.08632125e+01 -6.92206432e+02]\n",
      "  [ 3.27520067e+01 -6.91307213e+02]\n",
      "  [ 4.54017436e+01 -6.90152051e+02]]\n",
      "\n",
      " [[-4.85356456e+01 -6.80039296e+02]\n",
      "  [-3.32532899e+01 -6.81106541e+02]\n",
      "  [-1.87012080e+01 -6.81995354e+02]\n",
      "  [-5.26192962e+00 -6.82599888e+02]\n",
      "  [ 5.94913736e+00 -6.82762606e+02]\n",
      "  [ 1.15788852e+01 -6.82165207e+02]\n",
      "  [ 5.45077751e+00 -6.80701732e+02]\n",
      "  [-1.33215033e+00 -6.79992600e+02]\n",
      "  [ 2.76251443e+00 -6.79796915e+02]\n",
      "  [ 1.21169440e+01 -6.79450716e+02]\n",
      "  [ 2.36272669e+01 -6.78831448e+02]\n",
      "  [ 3.62308525e+01 -6.77940464e+02]\n",
      "  [ 4.95916980e+01 -6.76809834e+02]]\n",
      "\n",
      " [[-5.22896259e+01 -6.64349335e+02]\n",
      "  [-3.67021159e+01 -6.65415236e+02]\n",
      "  [-2.17760458e+01 -6.66316358e+02]\n",
      "  [-7.99424547e+00 -6.66957929e+02]\n",
      "  [ 3.47697221e+00 -6.67215729e+02]\n",
      "  [ 9.53062578e+00 -6.66892806e+02]\n",
      "  [ 5.02261649e+00 -6.66053103e+02]\n",
      "  [-1.00726562e-01 -6.65800276e+02]\n",
      "  [ 4.53943748e+00 -6.65747826e+02]\n",
      "  [ 1.43547651e+01 -6.65444409e+02]\n",
      "  [ 2.64421802e+01 -6.64836009e+02]\n",
      "  [ 3.96926660e+01 -6.63949287e+02]\n",
      "  [ 5.36809589e+01 -6.62832076e+02]]\n",
      "\n",
      " [[-5.60693512e+01 -6.48406750e+02]\n",
      "  [-4.01579966e+01 -6.49489796e+02]\n",
      "  [-2.48556014e+01 -6.50411249e+02]\n",
      "  [-1.07234814e+01 -6.51102217e+02]\n",
      "  [ 1.05296111e+00 -6.51484080e+02]\n",
      "  [ 7.62573336e+00 -6.51464308e+02]\n",
      "  [ 4.67820512e+00 -6.51160427e+02]\n",
      "  [ 1.12194771e+00 -6.51192427e+02]\n",
      "  [ 6.34295967e+00 -6.51194281e+02]\n",
      "  [ 1.66238195e+01 -6.50898651e+02]\n",
      "  [ 2.92753258e+01 -6.50288623e+02]\n",
      "  [ 4.31497475e+01 -6.49398429e+02]\n",
      "  [ 5.77370766e+01 -6.48279284e+02]]\n",
      "\n",
      " [[-5.99675632e+01 -6.32139189e+02]\n",
      "  [-4.36469116e+01 -6.33260074e+02]\n",
      "  [-2.79131207e+01 -6.34214664e+02]\n",
      "  [-1.33795941e+01 -6.34970357e+02]\n",
      "  [-1.21960177e+00 -6.35497187e+02]\n",
      "  [ 5.94865962e+00 -6.35776414e+02]\n",
      "  [ 4.41893982e+00 -6.35910923e+02]\n",
      "  [ 2.24151900e+00 -6.36129265e+02]\n",
      "  [ 8.08273712e+00 -6.36132893e+02]\n",
      "  [ 1.88660246e+01 -6.35816685e+02]\n",
      "  [ 3.21061645e+01 -6.35191385e+02]\n",
      "  [ 4.66247131e+01 -6.34287696e+02]\n",
      "  [ 6.18366974e+01 -6.33148967e+02]]\n",
      "\n",
      " [[-6.40648721e+01 -6.15414573e+02]\n",
      "  [-4.72101043e+01 -6.16588549e+02]\n",
      "  [-3.09557561e+01 -6.17587913e+02]\n",
      "  [-1.59363944e+01 -6.18421796e+02]\n",
      "  [-3.30225504e+00 -6.19102560e+02]\n",
      "  [ 4.52279539e+00 -6.19655729e+02]\n",
      "  [ 4.23749937e+00 -6.20142767e+02]\n",
      "  [ 3.22709892e+00 -6.20495676e+02]\n",
      "  [ 9.70768454e+00 -6.20483358e+02]\n",
      "  [ 2.10489381e+01 -6.20134426e+02]\n",
      "  [ 3.49305362e+01 -6.19484624e+02]\n",
      "  [ 5.01511861e+01 -6.18557524e+02]\n",
      "  [ 6.60592733e+01 -6.17382111e+02]]\n",
      "\n",
      " [[-6.84001624e+01 -5.98069916e+02]\n",
      "  [-5.08693587e+01 -5.99301252e+02]\n",
      "  [-3.39956187e+01 -6.00351518e+02]\n",
      "  [-1.84036094e+01 -6.01274044e+02]\n",
      "  [-5.20175074e+00 -6.02114303e+02]\n",
      "  [ 3.33980332e+00 -6.02915669e+02]\n",
      "  [ 4.12743015e+00 -6.03690921e+02]\n",
      "  [ 4.08370427e+00 -6.04151050e+02]\n",
      "  [ 1.11985118e+01 -6.04121627e+02]\n",
      "  [ 2.31600374e+01 -6.03739159e+02]\n",
      "  [ 3.77573257e+01 -6.03061861e+02]\n",
      "  [ 5.37631038e+01 -6.02105300e+02]\n",
      "  [ 7.04649927e+01 -6.00882928e+02]]\n",
      "\n",
      " [[-7.29640389e+01 -5.79934639e+02]\n",
      "  [-5.46263321e+01 -5.81217514e+02]\n",
      "  [-3.70528349e+01 -5.82320068e+02]\n",
      "  [-2.08196609e+01 -5.83337894e+02]\n",
      "  [-6.95927686e+00 -5.84339226e+02]\n",
      "  [ 2.36840099e+00 -5.85369372e+02]\n",
      "  [ 4.09775224e+00 -5.86399448e+02]\n",
      "  [ 4.83512942e+00 -5.86956767e+02]\n",
      "  [ 1.25607683e+01 -5.86911068e+02]\n",
      "  [ 2.51971419e+01 -5.86497090e+02]\n",
      "  [ 4.05923647e+01 -5.85794384e+02]\n",
      "  [ 5.74739612e+01 -5.84810605e+02]\n",
      "  [ 7.50719812e+01 -5.83543274e+02]]\n",
      "\n",
      " [[-7.77067274e+01 -5.60862589e+02]\n",
      "  [-5.84718329e+01 -5.62183819e+02]\n",
      "  [-4.01598960e+01 -5.63333756e+02]\n",
      "  [-2.32478535e+01 -5.64441268e+02]\n",
      "  [-8.64617079e+00 -5.65591982e+02]\n",
      "  [ 1.56117036e+00 -5.66831687e+02]\n",
      "  [ 4.17108852e+00 -5.68107260e+02]\n",
      "  [ 5.52615238e+00 -5.68778052e+02]\n",
      "  [ 1.38224054e+01 -5.68725982e+02]\n",
      "  [ 2.71661189e+01 -5.68284867e+02]\n",
      "  [ 4.34239347e+01 -5.67561327e+02]\n",
      "  [ 6.12620742e+01 -5.66557658e+02]\n",
      "  [ 7.98526166e+01 -5.65256136e+02]]\n",
      "\n",
      " [[-8.25513670e+01 -5.40762957e+02]\n",
      "  [-6.23904894e+01 -5.42104216e+02]\n",
      "  [-4.33508036e+01 -5.43289962e+02]\n",
      "  [-2.57599266e+01 -5.44465059e+02]\n",
      "  [-1.03597046e+01 -5.45729255e+02]\n",
      "  [ 8.50554576e-01 -5.47147744e+02]\n",
      "  [ 4.36843217e+00 -5.48673731e+02]\n",
      "  [ 6.23223659e+00 -5.49497947e+02]\n",
      "  [ 1.50344071e+01 -5.49466782e+02]\n",
      "  [ 2.90862042e+01 -5.49008020e+02]\n",
      "  [ 4.62452883e+01 -5.48266214e+02]\n",
      "  [ 6.50982909e+01 -5.47246960e+02]\n",
      "  [ 8.47609178e+01 -5.45920951e+02]]\n",
      "\n",
      " [[-8.74038777e+01 -5.19621389e+02]\n",
      "  [-6.63586096e+01 -5.20959422e+02]\n",
      "  [-4.66520323e+01 -5.22162186e+02]\n",
      "  [-2.84227794e+01 -5.23374731e+02]\n",
      "  [-1.22134743e+01 -5.24702341e+02]\n",
      "  [ 1.18638148e-01 -5.26237471e+02]\n",
      "  [ 4.70572635e+00 -5.27997921e+02]\n",
      "  [ 7.05452381e+00 -5.29036934e+02]\n",
      "  [ 1.62717556e+01 -5.29068166e+02]\n",
      "  [ 3.10004967e+01 -5.28602145e+02]\n",
      "  [ 4.90734765e+01 -5.27839455e+02]\n",
      "  [ 6.89745055e+01 -5.26801602e+02]\n",
      "  [ 8.97650283e+01 -5.25453613e+02]]\n",
      "\n",
      " [[-9.21601401e+01 -4.97503563e+02]\n",
      "  [-7.03352320e+01 -4.98812961e+02]\n",
      "  [-5.00683977e+01 -5.00011816e+02]\n",
      "  [-3.12895911e+01 -5.01234065e+02]\n",
      "  [-1.43253601e+01 -5.02579800e+02]\n",
      "  [-8.14991836e-01 -5.04157754e+02]\n",
      "  [ 5.16482998e+00 -5.06077290e+02]\n",
      "  [ 8.09940657e+00 -5.07372276e+02]\n",
      "  [ 1.76342370e+01 -5.07500654e+02]\n",
      "  [ 3.29819600e+01 -5.07029806e+02]\n",
      "  [ 5.19558451e+01 -5.06237102e+02]\n",
      "  [ 7.29111566e+01 -5.05171325e+02]\n",
      "  [ 9.48539665e+01 -5.03798571e+02]]\n",
      "\n",
      " [[-9.67257368e+01 -4.74533623e+02]\n",
      "  [-7.42712911e+01 -4.75796506e+02]\n",
      "  [-5.35926560e+01 -4.76979137e+02]\n",
      "  [-3.43978807e+01 -4.78196722e+02]\n",
      "  [-1.67912358e+01 -4.79535484e+02]\n",
      "  [-2.13851388e+00 -4.81102313e+02]\n",
      "  [ 5.61503874e+00 -4.83048304e+02]\n",
      "  [ 9.44938706e+00 -4.84537506e+02]\n",
      "  [ 1.92490917e+01 -4.84763386e+02]\n",
      "  [ 3.51336912e+01 -4.84280024e+02]\n",
      "  [ 5.49630728e+01 -4.83446441e+02]\n",
      "  [ 7.69441288e+01 -4.82343380e+02]\n",
      "  [ 1.00019016e+02 -4.80944411e+02]]\n",
      "\n",
      " [[-1.01063413e+02 -4.50852244e+02]\n",
      "  [-7.81510626e+01 -4.52071715e+02]\n",
      "  [-5.72371931e+01 -4.53234215e+02]\n",
      "  [-3.77825930e+01 -4.54434039e+02]\n",
      "  [-1.96765153e+01 -4.55742659e+02]\n",
      "  [-3.98812843e+00 -4.57251176e+02]\n",
      "  [ 5.80285567e+00 -4.59087353e+02]\n",
      "  [ 1.11150934e+01 -4.60605159e+02]\n",
      "  [ 2.12628290e+01 -4.60880197e+02]\n",
      "  [ 3.75791264e+01 -4.60378924e+02]\n",
      "  [ 5.81711182e+01 -4.59503577e+02]\n",
      "  [ 8.11024266e+01 -4.58362046e+02]\n",
      "  [ 1.05229415e+02 -4.56944642e+02]]\n",
      "\n",
      " [[-1.05246821e+02 -4.26575768e+02]\n",
      "  [-8.20317777e+01 -4.27770927e+02]\n",
      "  [-6.10463605e+01 -4.28910522e+02]\n",
      "  [-4.14787489e+01 -4.30069806e+02]\n",
      "  [-2.30157994e+01 -4.31307052e+02]\n",
      "  [-6.44200798e+00 -4.32681598e+02]\n",
      "  [ 5.45098639e+00 -4.34265199e+02]\n",
      "  [ 1.29909909e+01 -4.35636612e+02]\n",
      "  [ 2.37880132e+01 -4.35938279e+02]\n",
      "  [ 4.04258495e+01 -4.35431775e+02]\n",
      "  [ 6.16318998e+01 -4.34523832e+02]\n",
      "  [ 8.53819782e+01 -4.33350287e+02]\n",
      "  [ 1.10414564e+02 -4.31930492e+02]]\n",
      "\n",
      " [[-1.09373359e+02 -4.01793112e+02]\n",
      "  [-8.59916840e+01 -4.02982765e+02]\n",
      "  [-6.50693369e+01 -4.04096348e+02]\n",
      "  [-4.55040257e+01 -4.05185323e+02]\n",
      "  [-2.67991131e+01 -4.06299779e+02]\n",
      "  [-9.49742390e+00 -4.07466989e+02]\n",
      "  [ 4.39337116e+00 -4.08686993e+02]\n",
      "  [ 1.47969028e+01 -4.09755036e+02]\n",
      "  [ 2.68260576e+01 -4.10086317e+02]\n",
      "  [ 4.37304676e+01 -4.09629619e+02]\n",
      "  [ 6.53632546e+01 -4.08703702e+02]\n",
      "  [ 8.97581897e+01 -4.07504672e+02]\n",
      "  [ 1.15506976e+02 -4.06091487e+02]]\n",
      "\n",
      " [[-1.13525369e+02 -3.76562325e+02]\n",
      "  [-9.01080140e+01 -3.77757809e+02]\n",
      "  [-6.93481883e+01 -3.78838698e+02]\n",
      "  [-4.98445239e+01 -3.79833867e+02]\n",
      "  [-3.09517920e+01 -3.80789354e+02]\n",
      "  [-1.30261170e+01 -3.81706568e+02]\n",
      "  [ 2.69385497e+00 -3.82552244e+02]\n",
      "  [ 1.61217424e+01 -3.83279490e+02]\n",
      "  [ 3.01172550e+01 -3.83577443e+02]\n",
      "  [ 4.74866937e+01 -3.83193960e+02]\n",
      "  [ 6.93718513e+01 -3.82267223e+02]\n",
      "  [ 9.42300979e+01 -3.81046271e+02]\n",
      "  [ 1.20520965e+02 -3.79628801e+02]]\n",
      "\n",
      " [[-1.17771594e+02 -3.50914022e+02]\n",
      "  [-9.44510772e+01 -3.52112736e+02]\n",
      "  [-7.39106490e+01 -3.53132219e+02]\n",
      "  [-5.44657950e+01 -3.54002705e+02]\n",
      "  [-3.53580462e+01 -3.54783006e+02]\n",
      "  [-1.67774956e+01 -3.55462134e+02]\n",
      "  [ 5.82665359e-01 -3.56012156e+02]\n",
      "  [ 1.67038097e+01 -3.56463485e+02]\n",
      "  [ 3.31221339e+01 -3.56668938e+02]\n",
      "  [ 5.15251898e+01 -3.56328704e+02]\n",
      "  [ 7.36460457e+01 -3.55411037e+02]\n",
      "  [ 9.88381417e+01 -3.54151977e+02]\n",
      "  [ 1.25584867e+02 -3.52687584e+02]]\n",
      "\n",
      " [[-1.22190670e+02 -3.24838502e+02]\n",
      "  [-9.90661381e+01 -3.26026238e+02]\n",
      "  [-7.87180372e+01 -3.26958626e+02]\n",
      "  [-5.92507604e+01 -3.27678238e+02]\n",
      "  [-3.98158670e+01 -3.28284588e+02]\n",
      "  [-2.04863062e+01 -3.28791018e+02]\n",
      "  [-1.70617798e+00 -3.29167779e+02]\n",
      "  [ 1.65394840e+01 -3.29432453e+02]\n",
      "  [ 3.53407821e+01 -3.29521960e+02]\n",
      "  [ 5.54239275e+01 -3.29174395e+02]\n",
      "  [ 7.81279309e+01 -3.28223926e+02]\n",
      "  [ 1.03699468e+02 -3.26861086e+02]\n",
      "  [ 1.30971677e+02 -3.25259762e+02]]\n",
      "\n",
      " [[-1.26740313e+02 -2.98330420e+02]\n",
      "  [-1.03913422e+02 -2.99459041e+02]\n",
      "  [-8.36847411e+01 -3.00259941e+02]\n",
      "  [-6.40173599e+01 -3.00807145e+02]\n",
      "  [-4.41379055e+01 -3.01231754e+02]\n",
      "  [-2.40794717e+01 -3.01560805e+02]\n",
      "  [-4.11225819e+00 -3.01773761e+02]\n",
      "  [ 1.57821092e+01 -3.01888115e+02]\n",
      "  [ 3.66413140e+01 -3.01879604e+02]\n",
      "  [ 5.87975276e+01 -3.01497623e+02]\n",
      "  [ 8.27594571e+01 -3.00492522e+02]\n",
      "  [ 1.09050236e+02 -2.98963866e+02]\n",
      "  [ 1.37128434e+02 -2.97107352e+02]]\n",
      "\n",
      " [[-1.31337939e+02 -2.71322523e+02]\n",
      "  [-1.08950530e+02 -2.72285305e+02]\n",
      "  [-8.86939697e+01 -2.72878575e+02]\n",
      "  [-6.86577740e+01 -2.73200018e+02]\n",
      "  [-4.83354579e+01 -2.73377253e+02]\n",
      "  [-2.76300515e+01 -2.73458644e+02]\n",
      "  [-6.68957626e+00 -2.73420581e+02]\n",
      "  [ 1.45443544e+01 -2.73274396e+02]\n",
      "  [ 3.71549078e+01 -2.73068679e+02]\n",
      "  [ 6.14704289e+01 -2.72599556e+02]\n",
      "  [ 8.74427156e+01 -2.71595687e+02]\n",
      "  [ 1.15123350e+02 -2.69966812e+02]\n",
      "  [ 1.44480477e+02 -2.67818140e+02]]\n",
      "\n",
      " [[-1.35944751e+02 -2.43635869e+02]\n",
      "  [-1.14110651e+02 -2.44252219e+02]\n",
      "  [-9.36897254e+01 -2.44515596e+02]\n",
      "  [-7.32411059e+01 -2.44511152e+02]\n",
      "  [-5.25480068e+01 -2.44304218e+02]\n",
      "  [-3.12944439e+01 -2.44001351e+02]\n",
      "  [-9.54762037e+00 -2.43554376e+02]\n",
      "  [ 1.28379615e+01 -2.42949258e+02]\n",
      "  [ 3.69964350e+01 -2.42339955e+02]\n",
      "  [ 6.33814064e+01 -2.41691071e+02]\n",
      "  [ 9.17671684e+01 -2.40753438e+02]\n",
      "  [ 1.21648570e+02 -2.39223256e+02]\n",
      "  [ 1.52867826e+02 -2.37004832e+02]]\n",
      "\n",
      " [[-1.40477281e+02 -2.15027699e+02]\n",
      "  [-1.19314000e+02 -2.14987188e+02]\n",
      "  [-9.87687436e+01 -2.14714291e+02]\n",
      "  [-7.79759537e+01 -2.14228940e+02]\n",
      "  [-5.69240113e+01 -2.13505515e+02]\n",
      "  [-3.52608449e+01 -2.12611109e+02]\n",
      "  [-1.27899116e+01 -2.11551795e+02]\n",
      "  [ 1.06530129e+01 -2.10340801e+02]\n",
      "  [ 3.62490301e+01 -2.09247470e+02]\n",
      "  [ 6.44876096e+01 -2.08378260e+02]\n",
      "  [ 9.51636355e+01 -2.07531031e+02]\n",
      "  [ 1.27573953e+02 -2.06280536e+02]\n",
      "  [ 1.61055685e+02 -2.04352903e+02]]\n",
      "\n",
      " [[-1.44818034e+02 -1.85228363e+02]\n",
      "  [-1.24505606e+02 -1.84090889e+02]\n",
      "  [-1.04066568e+02 -1.82952918e+02]\n",
      "  [-8.30411525e+01 -1.81737560e+02]\n",
      "  [-6.17113739e+01 -1.80276410e+02]\n",
      "  [-3.96706060e+01 -1.78633817e+02]\n",
      "  [-1.64171512e+01 -1.76920073e+02]\n",
      "  [ 8.08630521e+00 -1.75121097e+02]\n",
      "  [ 3.50515362e+01 -1.73633848e+02]\n",
      "  [ 6.48355254e+01 -1.72627590e+02]\n",
      "  [ 9.73868858e+01 -1.71924902e+02]\n",
      "  [ 1.32139957e+02 -1.71046477e+02]\n",
      "  [ 1.67967767e+02 -1.69557762e+02]]\n",
      "\n",
      " [[-1.48738492e+02 -1.54253342e+02]\n",
      "  [-1.29557573e+02 -1.51497817e+02]\n",
      "  [-1.09416730e+02 -1.49230600e+02]\n",
      "  [-8.82555861e+01 -1.47115187e+02]\n",
      "  [-6.67033361e+01 -1.44747265e+02]\n",
      "  [-4.41588386e+01 -1.42269564e+02]\n",
      "  [-2.01186681e+01 -1.39844177e+02]\n",
      "  [ 5.48735518e+00 -1.37515565e+02]\n",
      "  [ 3.36408497e+01 -1.35676049e+02]\n",
      "  [ 6.45325380e+01 -1.34520695e+02]\n",
      "  [ 9.84909371e+01 -1.34038825e+02]\n",
      "  [ 1.35270244e+02 -1.33601548e+02]\n",
      "  [ 1.73381859e+02 -1.32629078e+02]]\n",
      "\n",
      " [[-1.52713191e+02 -1.21495855e+02]\n",
      "  [-1.34335890e+02 -1.17407379e+02]\n",
      "  [-1.14484808e+02 -1.13966040e+02]\n",
      "  [-9.33208214e+01 -1.10843970e+02]\n",
      "  [-7.13392896e+01 -1.07660933e+02]\n",
      "  [-4.80283040e+01 -1.04511097e+02]\n",
      "  [-2.31389388e+01 -1.01491182e+02]\n",
      "  [ 3.50680414e+00 -9.87248217e+01]\n",
      "  [ 3.25136008e+01 -9.64529907e+01]\n",
      "  [ 6.39922742e+01 -9.48307846e+01]\n",
      "  [ 9.85057173e+01 -9.41137862e+01]\n",
      "  [ 1.36887690e+02 -9.41883733e+01]\n",
      "  [ 1.77513749e+02 -9.34978455e+01]]\n",
      "\n",
      " [[-1.57887388e+02 -8.62528954e+01]\n",
      "  [-1.38755782e+02 -8.24988663e+01]\n",
      "  [-1.18813840e+02 -7.84494676e+01]\n",
      "  [-9.74076621e+01 -7.47266759e+01]\n",
      "  [-7.45378295e+01 -7.12168503e+01]\n",
      "  [-5.02699646e+01 -6.77616709e+01]\n",
      "  [-2.45436756e+01 -6.44069960e+01]\n",
      "  [ 2.82224410e+00 -6.11777242e+01]\n",
      "  [ 3.22786555e+01 -5.82047961e+01]\n",
      "  [ 6.39638528e+01 -5.54714552e+01]\n",
      "  [ 9.83845218e+01 -5.30979229e+01]\n",
      "  [ 1.36762723e+02 -5.17637568e+01]\n",
      "  [ 1.81037185e+02 -5.20647500e+01]]\n",
      "\n",
      " [[-1.61848658e+02 -5.27013659e+01]\n",
      "  [-1.42355418e+02 -4.87736600e+01]\n",
      "  [-1.21828277e+02 -4.52016750e+01]\n",
      "  [-9.94443582e+01 -4.22256448e+01]\n",
      "  [-7.54489677e+01 -3.93845421e+01]\n",
      "  [-5.04016634e+01 -3.61803731e+01]\n",
      "  [-2.41648979e+01 -3.25701250e+01]\n",
      "  [ 3.50117751e+00 -2.85215680e+01]\n",
      "  [ 3.31246417e+01 -2.42557769e+01]\n",
      "  [ 6.49279028e+01 -1.96182517e+01]\n",
      "  [ 9.93655574e+01 -1.42942295e+01]\n",
      "  [ 1.37484856e+02 -7.95359408e+00]\n",
      "  [ 1.83656879e+02  0.00000000e+00]]]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi4AAADRCAYAAADrETk6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAbf0lEQVR4nO3df2zU9R3H8dch7VGwrb/oj5MClamTwogDByUCKrMIG9HpHLpEa5hmTHTDStjqj8CWxaIzhDlQ4pRfIxv+UVAMuFEyWuYAI6xOxpCwUGkH7Qhm9sqPXin97I+lF8/rD+7D3ff6uT4fySX07vvu93Of+7T34tvvfd8+Y4wRAACAAwYkewAAAAAXi+ACAACcQXABAADOILgAAABnEFwAAIAzCC4AAMAZBBcAAOAMggsAAHDGwGQPIF46Ojp04sQJZWZmyufzJXs4AADgIhhj1NLSokAgoAEDej+ekjLB5cSJEyooKEj2MAAAgIWGhgYNGzas1+1SJrhkZmZKkkaOHHlRie2L0tLSrPY5cGDs0xfr2C5lX5J02WWXebo/mzrbI2Rez4nta2ezP9sx2q5lr9eJzf7S09Ot9uXlWpbsxzlo0CCruoyMjJhrBg8ebLWv7Oxsq7rO389e1V1++eWe7ct2Lm3rbNeJzbq0/b0Qq2AwqIKCgot+DRIWXF599VX96le/UmNjo4qKirR8+XJNmTKl2+1rampUVlamgwcPKhAIaNGiRZo3b95F76/zzW/AgAExv8HYvjg2dV6++Une/9JO5eDiZZ0rwcXL/XkdXGyfm+04/X6/VZ3NG5lN2LmUOts36SFDhljV2QQXmxrJfoy2dakYXDpd7HtBQk7Ofeutt7RgwQI9++yzqq2t1ZQpUzRz5kzV19d3uX1dXZ1mzZqlKVOmqLa2Vs8884x+/OMfq7KyMhHDAwAAjkpIcFm2bJl+8IMf6NFHH9VNN92k5cuXq6CgQK+99lqX269atUrDhw/X8uXLddNNN+nRRx/V3Llz9fLLLydieAAAwFFxDy5tbW3av3+/SkpKIu4vKSnR7t27u6zZs2dP1PYzZszQvn37dP78+XgPEQAAOCru57icOnVKFy5cUG5ubsT9ubm5ampq6rKmqampy+3b29t16tQp5efnR9WEQiGFQqHw18FgMA6jBwAAfVnCLkD35ZNsjDE9nnjT1fZd3d+poqJC2dnZ4RsfhQYAIPXFPbhcc801uuyyy6KOrpw8eTLqqEqnvLy8LrcfOHCgrr766i5rysvL1dzcHL41NDTE5wkAAIA+K+7BJT09XePHj1dVVVXE/VVVVZo8eXKXNcXFxVHbb9++XRMmTOj2I4l+v19ZWVkRNwAAkNoS8qeisrIyvfHGG1q9erUOHTqkp556SvX19eHrspSXl+vhhx8Obz9v3jwdO3ZMZWVlOnTokFavXq0333xTCxcuTMTwAACAoxJyAbo5c+bos88+0y9+8Qs1NjZqzJgx2rZtm0aMGCFJamxsjLimS2FhobZt26annnpKK1euVCAQ0CuvvKL77rsvEcMDAACO8pnOs2AdFwwGlZ2dreuuuy7mq9N6eUVOry/570IdV86N3764cm40V66ca3tFVJur0tpetdX2kv+2f8q3rbO5fL/tJf+5cu6l63z/bm5uvqjXPGV6FXUKBAIx/6LysldRKr+xSN6GuVSuS/WeVl62y3Al4HodsGy0tbVZ1Z09ezbOI+nZhQsXYq6xvWaY7ZzY7s+2NYFNuwbbNhSxruVYj58k7OPQAAAA8UZwAQAAziC4AAAAZxBcAACAMwguAADAGQQXAADgDIILAABwBsEFAAA4g+ACAACcQXABAADOILgAAABnEFwAAIAzCC4AAMAZKdcdevjw4TG377ZtQ2/TldW2A6wrdTadcX0+n9W+bOtseTlO23153fnathOylx2zbetSeX3Zvm4dHR1WdaFQyKou1q7Bndrb22Ou8brLs80YJbvO17Z1tq/3oEGDYto+1jnkiAsAAHAGwQUAADiD4AIAAJxBcAEAAM4guAAAAGcQXAAAgDMILgAAwBkEFwAA4AyCCwAAcEbcg0tFRYVuueUWZWZmKicnR/fcc48OHz7cY011dbV8Pl/U7ZNPPon38AAAgMPiHlxqamo0f/587d27V1VVVWpvb1dJSYnOnDnTa+3hw4fV2NgYvl1//fXxHh4AAHBY3HsV/fGPf4z4es2aNcrJydH+/fs1derUHmtzcnJ0xRVXxHtIAAAgRSS8yWJzc7Mk6aqrrup125tvvlmtra0aPXq0nnvuOd1+++3dbhsKhSIadwWDQUnSqFGjYm7wZNMs0bbO66aHts/Ntgmbl43pbBuA2TZu83p/XnKhqaPXzRK9ft1s92dT5+W+JO+bM9o0MPS66aHtnHhZ59WcXMxfZL4ooSfnGmNUVlamW2+9VWPGjOl2u/z8fL3++uuqrKzUpk2bdOONN2r69OnatWtXtzUVFRXKzs4O3woKChLxFAAAQB/iMwn8b8X8+fO1detWvf/++xo2bFhMtbNnz5bP59OWLVu6fLyrIy4FBQV6/vnnOeLyBRxxicYRl2gccYnGEZf41dmyfe1sfn/5/X6rfQ0ePNiqLjMzs8/X2T63IUOGxLR9MBjUyJEj1dzcrKysrF63T9gRlyeffFJbtmzRzp07Yw4tkjRp0iQdOXKk28f9fr+ysrIibgAAILXF/RwXY4yefPJJbd68WdXV1SosLLT6PrW1tcrPz4/z6AAAgMviHlzmz5+v3//+93rnnXeUmZmppqYmSVJ2drYyMjIkSeXl5Tp+/LjWr18vSVq+fLlGjhypoqIitbW1acOGDaqsrFRlZWW8hwcAABwW9+Dy2muvSZJuu+22iPvXrFmjRx55RJLU2Nio+vr68GNtbW1auHChjh8/royMDBUVFWnr1q2aNWtWvIcHAAAclpA/FfVm7dq1EV8vWrRIixYtivdQAABAiqFXEQAAcAbBBQAAOIPgAgAAnEFwAQAAziC4AAAAZxBcAACAMxLeHdpr1113Xcz9FWz78tj0D7Ldl23PIdseR7Z9YWx43SXV6y6wNnWu9EXysseR172KbHndC8umznYt2/7suFDndf8mr9elzTppa2uz2tf58+dj2r6lpSWm7TniAgAAnEFwAQAAziC4AAAAZxBcAACAMwguAADAGQQXAADgDIILAABwBsEFAAA4g+ACAACcQXABAADOILgAAABnEFwAAIAzCC4AAMAZKdcdOjc3V0OGDImpxrbrrE2nZy87UV/K/rzsXOp1l2cX6my799IdOn51LnTvta2zXcu23YJDoZCndTbjtP3ZsZ3Lc+fOWdXZsvmdYjv/ra2tMW1/+vTpmLbniAsAAHAGwQUAADgj7sFlyZIl8vl8Ebe8vLwea2pqajR+/HgNGjRI1113nVatWhXvYQEAgBSQkHNcioqKtGPHjvDXPZ1nUVdXp1mzZumxxx7Thg0b9Ne//lWPP/64hg4dqvvuuy8RwwMAAI5KSHAZOHBgr0dZOq1atUrDhw/X8uXLJUk33XST9u3bp5dffpngAgAAIiTkHJcjR44oEAiosLBQDzzwgI4ePdrttnv27FFJSUnEfTNmzNC+fft0/vz5RAwPAAA4Ku7BZeLEiVq/fr3+9Kc/6be//a2ampo0efJkffbZZ11u39TUpNzc3Ij7cnNz1d7erlOnTnW7n1AopGAwGHEDAACpLe7BZebMmbrvvvs0duxYffOb39TWrVslSevWreu25svXS+j8PH1P11GoqKhQdnZ2+FZQUBCH0QMAgL4s4R+HHjJkiMaOHasjR450+XheXp6ampoi7jt58qQGDhyoq6++utvvW15erubm5vCtoaEhruMGAAB9T8KvnBsKhXTo0CFNmTKly8eLi4v17rvvRty3fft2TZgwQWlpad1+X7/fL7/fH9exAgCAvi3uR1wWLlyompoa1dXV6YMPPtB3v/tdBYNBlZaWSvr/kZKHH344vP28efN07NgxlZWV6dChQ1q9erXefPNNLVy4MN5DAwAAjov7EZd///vfevDBB3Xq1CkNHTpUkyZN0t69ezVixAh
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Press any key...\n"
     ]
    }
   ],
   "source": [
    "from utils.topopt_88 import top88\n",
    "\n",
    "# Default input parameters\n",
    "mod_idx='test1'\n",
    "nelx=36\n",
    "nely=12\n",
    "volfrac=0.4\n",
    "rmin=5.4\n",
    "penal=3.0\n",
    "ft=1 # ft==0 -> sens, ft==1 -> dens\n",
    "\n",
    "top88(nelx,nely,volfrac,penal,rmin,ft,mod_idx)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a8f10f7c",
   "metadata": {},
   "source": [
    "## 检查数据编号顺序"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a2946a86-77ac-466f-8869-aa91f5e2602c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(13, 37, 2)\n",
      "(12, 36)\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "mod_idx='test1'\n",
    "nelx=36\n",
    "nely=12\n",
    "global_density = np.load('results/top88_' + mod_idx + '_xPhys_' + str(nelx) + '_' + str(nely) + '.npy')\n",
    "global_displace = np.load('results/top88_' + mod_idx + '_u_' + str(nelx) + '_' + str(nely) + '.npy')\n",
    "global_displace = global_displace.reshape(nelx+1,nely+1,2)\n",
    "global_displace = np.dstack((global_displace[:,:,0].T, global_displace[:,:,1].T))\n",
    "\n",
    "print(global_displace.shape)\n",
    "print(global_density.shape)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "9178f774",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZsAAAGNCAYAAADHF9x+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9eYxkWV7ejX/OuUvsW+5L7dXV1d3VNd3T3TM9VW1gYMzAvAbbGIMX/WxA2MJCsjRiMLZsWRpbmD9ABiQkkO0XeRgjDNh+bbP5hRnegYGpnq27s/Z9z6rcMzJjX+495/dHVETGvTdyj6xt7iO1ZvJWxLk3tvPc7/Y8QmutCREiRIgQIfYQ8klfQIgQIUKEeP4Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiFChAgRYs8Rkk2IECFChNhzhGQTIkSIECH2HCHZhAgRIkSIPUdINiEeO7TWaK2f9GWECBHiMcJ80hcQ4lsLWmuazSbVahXDMDBNE8MwMAwDKcN7nxAhnlcIHd5ihnhMUErRaDRQSlGv1wE6EY6UEtM0Q/IJEeI5RUg2IfYcWmscx8FxHLTWCCFoNBodMtFao5TypNaEEB7yMU0TIcSTegkhQoTYJUKyCbGnUErhOA6u6wItEmmn0tYjj27yaZNTPp8nl8sRi8U6BBSST4gQzw7Cmk2IPUGbMJrNZocw2uTQ3SDQizCEEBiG4Vnr+vXrHD9+vLOWlBLDMLAsq5N2C8knRIinFyHZhOg7utNmgIdodoL2c6WUWJblIbJGo9EhH3/NJySfECGeHoRkE6KvaJPAV7/6Vfbv38/k5GRf1u0mju7Ipx0htZsP6vV6SD4hQjyFCMkmRF+gtcZ1XRzHQSkF9E6R7RTtWk+v40BIPiFCPOUIySbErtEu+LebAKSUSCmfyODmRuRTr9dpNBqdawzJJ0SIx4eQbELsCt2zM1LKzoa9XiSyU+x0vW7yMQyj05ygtaZer3fmfQqFAkNDQ1iWhWmantcSIkSI3SMkmxA7Qjtt1u4282/OG5HDTjfxfpBXd7NCm3yUUpw9e5a3334b27Y7abc28bQHTEPyCRFi5wjJJsS20Stt5t+IN4tEtrtx79VG3yYWoJNWa0c+tVqt85j2YGlIPiFC7Awh2YTYFtrRjD9t5ke/02jQn8hmM/SKfNrRT0g+IULsHCHZhNgS/LMzm22ue1GzeRLYiHzq9Tq1Wq3TEBGST4gQ6yMkmxCboj07025p3opA5rMa2WwG/4Bqm3xc18V13XVbrUPyCfGtjpBsQqyL7kn9zdJm6z2/X3haN+o2+XSLinbPHLX/vdFoEI/HiUQinTbrp/U1hQixFwjJJkRPbKUJYCPsxZzN0xDZbIb1yOeDDz7g2LFj5HK5njWfkHxCPO8IySZEAO1oxnXdXaV/1iOHarXKrVu3SCaTHSXnzfA4NuK9ILNuEmkTS7v+1Va+bpNPt6ho6OUT4nlDSDYhOmi3+y4vL5PL5XZFNOtFNrOzs1y4cIFcLkexWOTq1atEIhFyuVznv0gksu717QUeV0TRJpZekU+bfCA0kgvxfCIkmxDAWtqsUCgwNTXFJz7xiV1vwt3koJTiypUrPHz4kBMnTjAwMAC0WqlXV1fJ5/Pcv3+fS5cuEY/HPeRjWdYzn2LaaMC1F/m0Fa0hJJ8QzwdCsgnhmZ1pp3l2u7l3RzblcpmzZ88CcPr0aWKxWGcjNU2TwcFBBgcHAWg2m6ysrJDP57l9+zYXLlwgmUzSbDYpFosMDQ1hms/v17YX+fjtFPwupiH5hHgW8Pz+akNsivVmZ/qVrtJaMzMzw8WLF5mcnOT48eMdElrvPJZlMTw8zPDwMACNRoN8Ps/169eZnZ1lenqaVCrViXoymYzHaO1pxk4IvJeRXJt82mm30EI7xLOAkGy+ReGfnem+o+4H2WitWVxc5MGDB3zoQx9iZGRkR+vYts3o6CgPHz5kdHSUgYEB8vk8+Xyey5cv02w2SafT5HI5BgYGSKVST+Vdfr8IfD3yuXfvHisrK7z88suhi2mIpxIh2XyLYSO7ZlgbxtxNKq1UKjE7O4sQopM22y3a1xWNRhkfH2d8fBytNdVqtUM+09PTKKXIZrOdyCeZTD7XG63fSM4wjNDFNMRTiZBsvoXgn53pNdvR/nunZDM9Pc3ly5dJJpOkUqm+EE33dfmPxeNx4vE4k5OTaK0pl8sd8rl9+zZCCE+zQTwef6LSN3uF7huH0EguxNOIkGy+RbDV2ZnuwvR24DgOly5dYnFxkddff518Pt9pAugXNrsmIQTJZJJkMsn+/ftRSlEqlcjn8ywuLnLjxg1M0/SQTzQa9RDsXmGvB1J73RyELqYhniaEZPOcw2/XvBUBTaDTmbYVFAoFzp49SyQS4fTp00SjUVZWVp64XI2UknQ6TTqd5uDBgyilKBQKLC8vMzMz45nxAajX632LxHrhcUQ2Wzl/6GIa4kkgJJvnGDuRnNlOZKO15v79+1y9epXDhw9z9OjRPXPq3Oo1bQQpJdlslmw2C7Ravttt1gDvv/++Z8Ynm81i2/ZuL/uxYaceQeu5mHaTT7vZIHQxDbFThGTznGI9u+bNsNWUUrPZ5OLFi+Tzed54443OnEz3Ok86stkMhmF0Znzu3bvHRz7yEWq1WqfeUy6XO5I6bfLZ6YzP40ij7Rbr2Sm0lSWmp6eJx+MMDw+HLqYhto2QbJ4zbGbXvBm602jrYXV1lampKRKJBKdPn+4pL/M0RjYboT2r0mvGpz3nU6vVdjXj86TTaNuFn3yKxWKHhLqN5EIL7RBbQUg2zxGUUjiOs2OlZljbYHpt7Fpr7t69y/Xr1zl69CiHDx9+bE6dT2Lzas/4jI6OAnSinvaMT6PRIJPJdMgnnU6vO+PzJBoE9uIc7ZRbr8gHQhfTEOsjJJvnAG0lgIWFhc7d9m5+3L2IotFocOHCBQqFAm+99VanqL7RGlu57iflj7MTbDbj47quZ8YnlUo91k12r8+llArMZK3nYtomn9DFNEQbIdk842gTTa1W4xvf+Abf9V3ftWvtMCGEJ42Wz+c5e/Ys6XSa06dPb6lo7l9jt3jaNqfNZnzu3r0L0CGf5yWy2ayTsRf5hC6mISAkm2ca/tkZ6M/df1uyRmvN7du3uXnzJseOHePgwYPbbjToJ550ZLMR/DM+WmuKxSL5fJ6lpSWUUrz//vsMDAx0Ip9YLNa39+lxvDfbJbT1FK395NNOu3XruoXk8/whJJtnEL1mZ7r/bbcQomVj/N5771Eul/noRz9KJpPZ9hrPc2SzGYQQnhmfL33pS7z44otUq1Xm5ua4du0atm0HBkx3iscV2exGd24zLx8hBKVSCcuyyGQyoYvpc4aQbJ4xbDQ7068NXmvN2bNnGRwc5PTp01iWte01vtUim83QJp+xsTEOHz7s8fF58OABV65cIRaL7XjGZ7dEsBX4aza7RS/yefDgAbFYDNu2PfI7oYvps4+QbJ4hbDY7I6XcFdlorbl58ybNZpNDhw5x/PjxHW8uz2Jk8zjlagzDYGBgoGMi5zhOZ8D0zp07lEqlbc/4POmazW7RbkwxTRPLsnpGPu0W9dDL59lDSDbPALY6O7Mbe4Barca5c+eo1WpEo1GGh4d33dHWbzzLkc1mME2ToaEhhoaGgFb3X5t8bty4QbVa3XDG51lIo20F3dHTemm30MX02UR
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAGNCAYAAAC4xs6PAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9d5Qk+V3mC38iIk1lVWX56vK2bbXv6e5pN0IagZDYvZg1erWLYJl9kZAAsSCBhAQXgQSCZYUEu7gVL0bau3ClFSCMBIsEWqEZ9fgu7013eZ/l0mdE/N4/qiM6MiuzKl2ZnonnnDlnuirrF5GZEb8nvu55JCGEwIYNGzZs2DgEyId9AjZs2LBh4/ULm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJhw4YNG4cGm4Rs2LBhw8ahwSYhGzZs2LBxaLBJyIYNGzZsHBpsErJx4BBCIIQ47NOwYcPGEYDjsE/AxusLQghisRihUAhFUXA4HCiKgqIoyLL9TGTDxusNkrAfSW0cEHRdJxqNous6kUgEwIyIZFnG4XDYpGTDxusMNgnZ2HcIIVBVFVVVEUIgSRLRaNQkGSEEuq7HpegkSYojJYfDgSRJh/UWbNiwsU+wScjGvkLXdVRVRdM082cTExPMzMxQUlJCeXk55eXleDwek2SspGSQ1tramvk6g5hsUrJh4/GHXROysS8wiCQWi5lEEolE6OnpIRwO097eTjAYZHFxkZGREVwuF+Xl5ZSVlZlkY11rdHSU06dPm2vJsoyiKDidTjN9Z5OSDRuPH2wSspF3WNNvsJ1aW15epre3l+rqai5fvoyu60iShCRJaJrGxsYGa2trzM/PMzw8jNvtjiMl2K4bOZ3OOIKLRqMmKSXWlGxSsmHj6MNOx9nIKwxyeOGFF2hqaqKuro7h4WFmZmY4d+4c9fX1ZoOCQUKJUFXVJKX19XW2trYAqKiooLa2lrKyMtxuN/CosSExfWeTkg0bjwdsErKRFwgh0DQNVVXRdZ2XX36Zmpoa5ufnAbh06RJFRUUAe5JQIlRV5cUXX6SkpIRQKITf76ewsNCsJ5WVleFyuczzMI5hk5ING0cfdjrORs4wZn+M5gNZlonFYoyOjtLc3MypU6dyarc2iKOhoYGKigpisRjr6+usra1x//59AoEARUVFcaTkdDrNcwPMtvBoNGqeo01KNmwcPmwSspETrLM/siyjaRoDAwMEg0Gam5s5c+ZMXo4jSZJJKE6nk+rqaqqrqwGIRqMmKY2PjxMMBikuLo4jJYNsjOhICEEkEjHnlTY3N6mqqsLpdOJwOJBl2SYlGzYOADYJ2cgKRvrN6H6TZZnNzU26u7spKCigrKyM4uLipH+b7eaeKnPscrk4duwYx44dAyASiZikNDo6Sjgcxuv1mk0OZWVlZvRjNDl0d3dz48YNXC6Xmb4zCMkYnLVJyYaN/MMmIRsZIzH9JkkSk5OTjI6O0t7eTnt7O/fu3dtVHy7TDT2T17vdbmpqaqipqQEgHA6bTQ7Dw8NEIhFKSkpMUiopKQEw03NGpBQOh81jGwOzNinZsJFf2CRkIyMY0Y+RfovFYvT29rK1tcW1a9fMdmpr+ixfyHa9goIC6urqqKurAyAUCpmkNDg4aNaJpqamqKqqorS0NC5SMqIlm5Rs2Mg/bBKykRYSZ39kWcbn89HT00NZWRl37twxmwEg/ySUzw3e4/Hg8Xior69HCEEwGOTFF18kFArR39+PqqqUlpaaNSWv15uUlCKRCOFwGFmWdzQ62KRkw0Z6sEnIxp4wZn90XTd/NjY2xoMHDzh9+jRNTU07NtyjFAntBkmSKCwsBOD06dO4XC6CwSBra2usra0xPT2NrutxpFRcXIyiKOY5GfUxTdOIRCJJW8JtUrJhIzlsErKRElZlAiP9Fg6H6enpIRqNcvPmTbxe765/ny8c1AYuSRJFRUUUFRXR2NiIEIJAIGCS0uTkJEIIs55kkJLDsX0rWUlJVVVzFioajVJYWIjb7TbbwW1SsmHDJiEbKZBs9mdpaYm+vj5qamq4evWqufEmgyzLj0UktBckSaK4uJji4mKampoQQuD3+01Sun//PpIkxZFSUVHRDlLq7Ozk5MmTlJeXJ60p2aRk4/UKm4Rs7IAR/WiahizL6LrO0NAQc3NznDt3zizw74XHLRJK53wlScLr9eL1emlubkbXdZOUVldXmZiYQJblOFIy0n3WWSVVVYnFYib5JIqx2l5KNl4vsEnIhgmjLdnn81FeXo4sywQCAbq7u5Flmdu3b5sb6l54nCKhXAhOlmVKSkooKSmhpaUFXdfZ2tpibW2N5eVlxsbGcDgcqKrK6uoqbrfbtKOAR5GSQUrGmrbBn43XC2wSsgE8Sr9tbm7S1dXFm9/8ZmZnZxkcHKS5uZmTJ09mvBE+bpFQPiDLMqWlpZSWltLa2oqu62xsbNDT04PP52N6ejqpbYXV4M/4LnaTGLJJycZrBTYJ2Yib/VEUBV3X6e3tZXV1lStXrlBVVZXxmo9TJLSfkGWZ8vJyHA4Hp06dwuv17mlbUVBQsMN11mpbkeg6a5OSjccZNgm9jpFs9sfv95t6cHfu3DEtE7JdPxl0XWdxcZHCwkK8Xm9aUc7jEgntBqP2U1FRQUVFBRBvWzE7O8vQ0BAFBQVxunfW78BKSkb6zrZCt/E4wyah1ymSzf48ePCA0dFRAK5du5bTRpZqTigYDNLV1UU0GjUL88aGaxTxUx33cYyEDKQ6d4fDQWVlJZWVlcA2KRm6d5OTk/T396e0rTDW1XWdqakp1tfX6ejosF1nbTxWsEnodYZkttvRaJTe3l4CgQCXLl2is7Mz5+NIkhRHcAALCwv09fVRX19Pe3t7XLuztYhvJSXD5ns/hl+PIhwOB1VVVWYKNFPbCiOdarvO2nhcYJPQ6wjJhEdXV1fp6emhoqKC27dvm78zCCpbWElD0zSGh4eZm5vj/Pnz1NTUmCSYrIifrF4SjUbNtOHjimw+z3RtKxwOB7quo2laXOcdPLLbSKXmYJOSjcOETUKvEyTO/gghGBkZYWpqio6ODhoaGuKII9eow1grEAjQ1dUV1+Kdam2jiG+IoFrrJaFQiPHxcRYWFpJGAfnAfkZa+Vo7lW3F1NQUwWCQZ599NqlthfUcbFKycZRgk9BrHIm227IsEwqF6O7uRtM0bt26Fef7Y2w+RqdctpAkiUAgwN27d2lsbOT06dMZd3BZ6yXBYJDS0lI8Hk9cFOD1euNIKZdz3m/sx8Zu2FaEQiECgQDHjx/f1bbCUAgH23XWxtGATUKvYSST3llcXKSvr4+6ujrOnDmzY9O2tgZnC03TWF5eZmtriytXrphP7blCluUdUYAhn2NsuFah0ZKSktdV67IkSWnZVhifUVlZWVLbCsN11kpKRpOD7TprI9+wSeg1ikTbbV3XGRgYYHFxkfPnz1NbW5v074zNJVsS8vv9dHV1oaoqVVVVeSOgZJue2+2mtraW2tpaU+1hbW0Nn8/HzMwMmqbFyeek2w6+H9jvpopU6yfaVhikZLSEp2NbYXy2MzMzFBYWUl1dbbvO2sgbbBJ6jSGZ7bbf76e7uxuHw8Ht27fNjrNksKbjMsXs7CwDAwM0NzfjcDjY2trK+n0kw15OrYkbbqL6NbBDaNR4vwexie7nMdJpJDFsKwoLC2loaDC9lNKxrVAUha2tLZO
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAakAAAGNCAYAAACrEY57AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOz9eZgc13neDf/Oqeq9p3v2fcO+7wAJgHtkinIsMf58JUpMixbzKZQTyaJpUrFjO/ErKzZtS7KkWMrnyI7zUp8kW04iO45eWRJlRaZIkSBBkNiXGWAwmH3tnqX3rqrz/tHonunpFcAABMi6r4uXhO5zTlV199Rdz3Pu536EUkphw4YNGzZs3IaQb/cJ2LBhw4YNG6Vgk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlA0bNmzYuG1hk5QNGzZs2LhtYZOUDRs2bNi4bWGTlI1bDqUUSqm3+zRs2LBxB0B/u0/AxrsLSinS6TTxeBxN09B1HU3T0DQNKe1nJhs2bORDKPuR1sYtgmVZpFIpLMsimUwC5CIqKSW6rtukZcOGjTzYJGXjpkMphWEYGIaBUgohBKlUKkdCSiksy8pLAQoh8khL13WEEG/XJdiwYeNtgk1SNm4qLMvCMAxM0wQy5JNN+ZUineWklSW1cDhMXV0dHo8nR1w2admw8c6HvSdl46YgSzTpdDpHNFlSWS6cKEY0Qgg0Tctbq7+/n02bNuXWklKiaRoOhyOXHrRJy4aNdx5skrKx6lie3gPyCOp6kJ0rpcThcOQRYCqVypHWyj0tm7Rs2LjzYZOUjVVFljyOHDlCV1cXHR0dq7LucsJZHmllI7KsKCOZTNqkZcPGOwg2SdlYFSilME0TwzCwLAsonsq7XmT3soq9DtikZcPGOxQ2Sdm4YWSFEFlxhJQSKeXbUrBbjrSSySSpVCp3jjZp2bBx+8MmKRs3hOW1T1LK3I2+VOQDYJoms7OzBAIBXC5XVccpt16leUCOiLKiDaUUyWQyV6+1sLBAY2MjDocDXdfzrsWGDRtvH2ySsnFdyKb3suq9lTf1UqQSiUQ4fvw4yWSSdDqNz+ejtraWuro66urqcDgcZY95o1gu4siSlmVZnDhxgrvvvhun05lLD2YJK1tYbJOWDRu3HjZJ2bhmFEvvrbyBryQppRSjo6OcO3eO7u5uuru7MU2T+fl5wuEwly9f5vTp0/j9/hxh1dbWout6br2bgSwhAbn0XzbSSiQSuTHZgmKbtGzYuLWwScrGNSEbPa1M763EcpIyDIMzZ84wOzvLnj17aGhoyDlONDU10dTUBEAqlSIcDhMOh+nv7yeRSFBTU0NtbW1eQfDNRLFIKxtt2aRlw8ath01SNqrCytqnSjflLEnNz89z4sQJPB4P99xzDy6Xq2Tazul00tLSQktLCwCJRIJwOMzc3BzJZJKzZ88yMjKSi7QCgUBe0e/NQDnSSiaTJBKJnFDEJi0bNlYfNknZqIhs7VNWWl6t8WsoFKK/v5+1a9eydu3aa75pu91u2traaGtrIxKJ0NbWhpSScDjM2NgYhmEQDAapra2lvr6empqam25Ku7IwOUtapmlimmZJybtNWjZsXB9skrJREsudHSql95YjlUoxNzeHaZrs37+furq6Gz4XIQROp5Pm5mba29tRShGLxXKR1sjICJZlEQwGc5FWTU3NTSeGLGktN8tdXjOWfT+VSuH1enG5XDm5u01aNmxUhk1SNoqiGnFEMYTDYU6cOIEQgu7u7lUhqOXnlIUQAp/Ph8/no7OzE6UU0Wg0t6d15coVgDzloM/ne9tI66233mLDhg3U1dUV3dOyScuGjeKwScpGAbLRk2maVZOTUoqBgQEGBgbYuHEjc3Nzq5p6q3QOQgj8fj9+v5+uri6UUiwuLhIOhwmFQgwMDCClzCMtr9dbkLpbbSwnnywhZff3sk7wWdJabpZr99KyYSMDm6Rs5JCVXYdCIerq6qomqEQiwcmTJ0kkEtx9990EAgEWFxdX/aZ/LesJIQgEAgQCAXp6erAsK0da09PTXLx4EV3Xc4R1retfD7KEVCzSypIW2A0gbdhYDpukbABL6b2FhQWOHz/Oe97znqoIanp6mlOnTtHY2MjevXtzdU3ZNVcLN5oKk1ISDAYJBoP09vZimiYLCwuEw2HGx8cBOHbsGPX19Tnicrvdq3HqQOnPohRpZR3es+duk5aNdytskrKRV/uUTUdVIgXLsujv72doaIitW7fS3t6eN+dmePet5nqapuVFUf/n//wf1q9fTzQaZXR0lPPnz+N2u3Nj6urqcDqdq3b8UihGWivbkqzsWmyTlo13MmySehejVO1TJTKIxWKcOHECy7I4dOgQfr+/5PqlcK1efLdCVFBbW0trayuQKUCem5vLiTDOnDmDz+fLc8MoZ+FUDNdzDcUaQGZJK5seXElauq7bIgwb7xjYJPUuxcrap+VP8OXIY2JigtOnT9Pe3s6mTZtKFtNeryFsOdxKV3Vd12lsbKSxsRGAdDqdI62BgQGi0WhJC6diWK1zL0VaQ0NDzM3NsWXLFrtrsY13FGySepehXFt3WCKXlSk/0zQ5f/484+PjbN++PRdxlIIQIkeAq4GbQXrXAofDkWfhlEwmc6S13MIpS1rBYPCmu2FAYQNITdPsrsU23lGwSepdhJW1T8Vqc7L/Xk5SkUiEEydOIKXk8OHDeL3eisdabVK53W6qLperqIVTOBzm3LlzpFIpAoFAjrSq2ee7ESx/4LAbQNp4J8EmqXcJqq19WqkyW+5cvmHDhqo36K+ZpCLzOL/7PzC278favAuKnN/NjqRuZP3lFk5ZKX+WtMbGxkin01y4cIGmpqacG8Zqih2KkaDdtdjGOwE2Sb3DsbKtezXGsJCxNurr62N2dpbdu3fn0lzV4lpISjv2E9z/6beRoSlcgNXQjHHoPRj3PIy5bR/cYb53Qgg8Hg8ejydn4fTSSy9RX1/P4uIiw8PDWJZFbW1trrj4Ri2cqonU7K7FNu5E2CT1Dsb1WBtln+5ff/11vF4vhw8fvq56oapIKhHH9fzncX7nr/LPYXYKxwt/jbzSjwzPkHr0Q8iWdW/rntSNIJuGa25upqam5qZZOF3v+FJdi5eTVlaEYXcttnGrYZPUOxSl2rqXg1KKoaEhAFpbW9m4ceN134wqkZS8cAr3538DbfRKwXtm1zqEmUI/+yYA7q88x35fgPCDH4B/8VEI1F7XOb2dWOk7WMrCaXZ2NmfhlFUNFrNwKrf+9aJUW5Js+nJkZASv10tTU5PdtdjGLYNNUu8wVGrrXgqpVIrTp0+zsLAAQHd39w3deEqSlGmgf/0/4/rmnyJMI//cNQ1z216088cRKxocOqILNH7/v2NNXMbcspvUP/kwOFfPEeJW3GTLNYgsZuEUCoWKWjjV1dXh8Xjy1rgZwoyVpLW4uJgjr+UNILORlk1aNm4GbJJ6B8GyrLwOttfqXB4IBDh8+DA/+tGPbvjJvBhJiaFL6H/4b5H9ZwrPvbkd5fOhnzlWdL1YUxuay4Hr7FH0s0dxvPj/kHz8aYy733ND53mrcC2f53ILJ6DAwunChQu4XK68SOtmqwdhSeJeLD1ody22cbNgk9Q7AFnniOnp6Vx9zvU4l2ejp9WQj688vvzbbyBf+QFidrpgrLF9P9rAWeTsROE5CoG
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "from mpl_toolkits.mplot3d import Axes3D\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "def surf_plot(dp,ny,nx,UorV='sqrt',non_dp=False):\n",
    "    fig = plt.figure()\n",
    "    ax = fig.add_subplot(projection='3d')\n",
    "\n",
    "    # Make data\n",
    "    # 生成网格数据\n",
    "    X, Y = np.meshgrid(np.arange(nx), np.arange(ny))\n",
    "\n",
    "    if non_dp:\n",
    "        Z = Z_x = Z_y = dp.reshape(ny,nx)\n",
    "    else:\n",
    "        Z = dp.reshape(ny,nx,2)\n",
    "        Z_x = Z[:,:,0].reshape(ny,nx)\n",
    "        Z_y = Z[:,:,1].reshape(ny,nx)\n",
    "\n",
    "\n",
    "    if UorV == 'u':\n",
    "        ax.plot_surface(X, Y, Z_x, rstride = 1, cstride = 1, cmap='rainbow')\n",
    "    elif UorV == 'v':\n",
    "        ax.plot_surface(X, Y, Z_y, rstride = 1, cstride = 1, cmap='rainbow')\n",
    "    else:\n",
    "        ax.plot_surface(X, Y, Z_x*Z_x+Z_y*Z_y, rstride = 1, cstride = 1, cmap='rainbow')\n",
    "\n",
    "    plt.show()\n",
    "\n",
    "surf_plot(global_displace,nely+1,nelx+1,'u')\n",
    "surf_plot(global_displace,nely+1,nelx+1,'v')\n",
    "surf_plot(global_displace,nely+1,nelx+1,'sqrt')\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "language": "python",
   "name": "base"
  },
  "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.10.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}