contour of potential

Author: Huangzizhou

python/new-formulation-high-order.ipynb

@@ -12,9 +12,9 @@
     "from plotly.subplots import make_subplots\n",
     "import meshio, numpy, copy\n",
     "import meshplot as mp\n",
-    "import sys\n",
-    "import pathlib\n",
-    "sys.path.append(str(pathlib.Path(\"/home/zizhou/smooth-ipc/build\") / \"debug\" / \"python\"))  # noqa\n",
+    "# import sys\n",
+    "# import pathlib\n",
+    "# sys.path.append(str(pathlib.Path(\"/home/zizhou/smooth-ipc/build\") / \"debug\" / \"python\"))  # noqa\n",
     "\n",
     "from ipctk import *"
    ]
@@ -25,8 +25,9 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "V = np.array([[0.5,0],[1.0,0],[1.5,0],[0.9,0.3],[1.01,0.1],[1.1,0.3],[1,-0.5]])\n",
-    "E = np.array([[1,0],[2,1],[3,4],[4,5],[5,3],[0,6],[6,2]], dtype=int)\n",
+    "V = np.array([[0.5,0],[1.0,0],[1.4,0],[0.5,1],[1.01,0.1],[1.5,1],[1,-0.5]])\n",
+    "# E = np.array([[1,0],[2,1],[3,4],[4,5],[5,3],[0,6],[6,2]], dtype=int)\n",
+    "E = np.array([[1,0],[2,1],[0,6],[6,2]], dtype=int)\n",
     "F = np.array([])\n",
     "\n",
     "dhat=0.11\n",
@@ -1985,6 +1986,34 @@
     "fig.show()"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "xs = numpy.linspace(0.3, 1.5, 100)\n",
+    "ys = numpy.linspace(1e-2, 4e-2, 40)\n",
+    "dtypes = numpy.empty((ys.size, xs.size), dtype=float)\n",
+    "for i, x in enumerate(xs):\n",
+    "    for j, y in enumerate(ys):\n",
+    "        V_disp = copy.deepcopy(V)\n",
+    "        V_disp[4, :] = [x, y]\n",
+    "        dtypes[j, i] = potential_single_point(V_disp)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "fig = go.Figure(data=[\n",
+    "    go.Contour(z=dtypes, x=xs, y=ys)\n",
+    "], layout=go.Layout(width=1000, height=500))\n",
+    "fig.show()"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,