From 44eadcae48ef2b666b3f67ec6eeb5aad09d22200 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Thu, 24 Sep 2020 11:09:36 +0200
Subject: [PATCH] Found third value for partitioning/gamma that works, beyond
0.9/0.1, 0.9/0.9.
---
FloryHugg_DiffUnbleached.ipynb | 23 ++++++++++++++++++++---
1 file changed, 20 insertions(+), 3 deletions(-)
diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index b8ed007..da0718c 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -88,12 +88,17 @@
"c0_1 = create_func(F, 'x[0]<0.1 ? 0 :' + p_tot(0.9, 0.9), 1)\n",
"Ga0_1 = create_func(F, p_tot(1, 1/9), 1)\n",
"\n",
+ "c_tot_2 = create_func(F, p_tot(0.9, 0.45), 1)\n",
+ "c0_2 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.45), 1)\n",
+ "Ga0_2 = create_func(F,p_tot(1, 0.08), 1)\n",
+ "\n",
"c_tot_9 = create_func(F, p_tot(0.9, 0.1), 1)\n",
"c0_9 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.1), 1)\n",
"Ga0_9 = create_func(F,p_tot(1, 1), 1)\n",
"\n",
- "c0_1 = calc_sim(c0_1, c_tot_1, Ga0_1, 0.001, 10, 2)\n",
- "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 10, 2)"
+ "c0_1 = calc_sim(c0_1, c_tot_1, Ga0_1, 0.001, 50, 0)\n",
+ "c0_2 = calc_sim(c0_2, c_tot_2, Ga0_2, 0.001, 50, 0)\n",
+ "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 50, 0)"
]
},
{
@@ -104,8 +109,20 @@
"source": [
"x = np.linspace(0, 1, 10000)\n",
"plt.plot(x, eval_func(c0_1, x))\n",
+ "plt.plot(x, eval_func(c0_2, x))\n",
"plt.plot(x, eval_func(c0_9, x))\n",
- "plt.xlim(0, 0.2)"
+ "plt.xlim(0.0, 0.125)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "# Check partitioning\n",
+ "cp = eval_func(c0_9, x)\n",
+ "print('Partitioning: ' + str(np.max(cp[1:1050])/np.min(cp[999:1050])))"
]
},
{
--
GitLab