From 89a9ec55b6ecc4d65d8aa9728db3b7811940fa68 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Wed, 21 Oct 2020 14:27:58 +0200
Subject: [PATCH] Updating Fenics_FloryHugg_Binary examples for Kafa.
---
Fenics_FloryHugg_Binary.ipynb | 28 +++++++++++++++-------------
1 file changed, 15 insertions(+), 13 deletions(-)
diff --git a/Fenics_FloryHugg_Binary.ipynb b/Fenics_FloryHugg_Binary.ipynb
index 11fd2fb..bbba925 100644
--- a/Fenics_FloryHugg_Binary.ipynb
+++ b/Fenics_FloryHugg_Binary.ipynb
@@ -59,7 +59,7 @@
"import random\n",
"from dolfin import *\n",
"\n",
- "fol = '/Users/hubatsch/ownCloud/Dropbox_Lars_Christoph/frap_theory/'"
+ "fol = '/Users/hubatsch/Desktop/frap_theory/'"
]
},
{
@@ -174,7 +174,7 @@
"\n",
"# Step in time\n",
"t = 0.0\n",
- "T = 70*dt\n",
+ "T = 5*dt\n",
"ti = time.time()\n",
"while (t < T):\n",
"# file << (u.split()[0], t)\n",
@@ -208,7 +208,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
- "## Thoughts on stability of algorithm\n",
+ "### Thoughts on stability of algorithm\n",
"- relatively robust to grid size changes\n",
"- initial conditions seem to have strong influence\n",
"- time step also plays a role for Newton solver\n",
@@ -219,7 +219,8 @@
"cell_type": "markdown",
"metadata": {},
"source": [
- "**Example for how to use Spheres class to create n positions of non-overlapping spheres inside a box of dimensions ls.**"
+ "### Example for use of 'Spheres' class \n",
+ "...to create n positions of non-overlapping spheres inside a box of dimensions ls."
]
},
{
@@ -239,19 +240,15 @@
"cell_type": "markdown",
"metadata": {},
"source": [
- "**Example on how to use convert_msh_to_xdmf to convert a mesh created in gmsh via 'gmsh -3 file -format msh2' to xdmf, which can be parallelised, in contrast to .xml meshes.**"
+ "The resulting points in points.txt can now be used in a .geo file, e.g. box_with_spheres. Simply copy them in, then create the mesh (see below)."
]
},
{
- "cell_type": "code",
- "execution_count": null,
+ "cell_type": "markdown",
"metadata": {},
- "outputs": [],
"source": [
- "from fem_utils import convert_msh_to_xdmf\n",
- "path_to_msh = \"/Users/hubatsch/ownCloud/Dropbox_Lars_Christoph/frap_theory/Meshes/box_with_spheres.msh\"\n",
- "path_to_xdmf = \"mesh1.xdmf\"\n",
- "convert_msh_to_xdmf(path_to_msh, path_to_xdmf)"
+ "### Example on how to use convert_msh_to_xdmf\n",
+ "...to convert a mesh created in gmsh via 'gmsh -3 file -format msh2' to xdmf, which can be parallelised, in contrast to .xml meshes."
]
},
{
@@ -259,7 +256,12 @@
"execution_count": null,
"metadata": {},
"outputs": [],
- "source": []
+ "source": [
+ "from fem_utils import convert_msh_to_xdmf\n",
+ "path_to_msh = \"/Users/hubatsch/Desktop/frap_theory/Meshes/box_with_spheres.msh\"\n",
+ "path_to_xdmf = \"Meshes/mesh1.xdmf\"\n",
+ "convert_msh_to_xdmf(path_to_msh, path_to_xdmf)"
+ ]
}
],
"metadata": {
--
GitLab