From 0515c6762882919a570b64a324bdc93dc3f54231 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Fri, 25 Sep 2020 15:35:31 +0200
Subject: [PATCH] Further simplifying weak form of radial symmetry.
---
FloryHugg_DiffUnbleached.ipynb | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index da0718c..04e3b89 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -51,9 +51,9 @@
"# (1-c_tot)*Ga0*2*X[0]*c.dx(0)*tc+\n",
"# (1-c_tot)*Ga0/c_tot*c*2*X[0]*c_tot.dx(0)*tc) * df.dx\n",
"# # Weak form radial symmetry:\n",
- " form = (df.inner((c-c0)/dt, tc*X[0]**2) +\n",
- " df.inner((1-c_tot)*Ga0*(df.grad(c)-c/c_tot*df.grad(c_tot)),\n",
- " df.grad(tc)*X[0]**2)) * df.dx\n",
+ " form = ((c-c0)/dt*tc*X[0]**2 +\n",
+ " (1-c_tot)*Ga0*(c.dx(0)-c/c_tot*c_tot.dx(0))*\n",
+ " tc.dx(0)*X[0]**2) * df.dx\n",
" \n",
" t = 0\n",
" # Solve in time\n",
--
GitLab