From 3c66108608461cc69ca20c6ada4a164df7a1c06e Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Tue, 15 Sep 2020 13:15:45 +0200
Subject: [PATCH] Works also for longer time courses.
---
FloryHugg_DiffUnbleached.ipynb | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index 5f5aff1..b24acb8 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -16,7 +16,7 @@
"# domain = ms.Sphere(df.Point(0, 0, 0), 1.0)\n",
"# mesh = ms.generate_mesh(domain, 50)\n",
"mesh = df.UnitIntervalMesh(100000)\n",
- "dt = 0.0000002\n",
+ "dt = 0.000001\n",
"\n",
"F = df.FunctionSpace(mesh, 'CG', 1)"
]
@@ -70,7 +70,7 @@
" t = 0\n",
" # Solve in time\n",
" ti = time.time()\n",
- " for i in range(10):\n",
+ " for i in range(60):\n",
" print(time.time() - ti)\n",
" df.solve(form == 0, c)\n",
" df.assign(c0, c)\n",
@@ -99,7 +99,7 @@
"plt.plot(np.linspace(0, 1, 10000), [c0_1([x]) for x in np.linspace(0, 1, 10000)])\n",
"plt.plot(np.linspace(0, 1, 10000), [c0_9([x]) for x in np.linspace(0, 1, 10000)])\n",
"plt.xlim(0.495, 0.505)\n",
- "plt.ylim(0, 0.3)\n",
+ "# plt.ylim(0, 0.3)\n",
"# 3D:\n",
"# plt.plot(np.linspace(0, 0.5, 1000), [c0([x, 0, 0]) for x in np.linspace(0, 0.5, 1000)])"
]
--
GitLab