Change title of example 8

Author: emmetfrancis

examples/example2/example2.ipynb

@@ -311,7 +311,7 @@
     "avg_X = [X.initial_condition]\n",
     "avg_B = [B.initial_condition]\n",
     "# Set loglevel to warning in order not to pollute notebook output\n",
-    "logger.setLevel(logging.DEBUG)\n",
+    "logger.setLevel(logging.WARNING)\n",
     "\n",
     "while True:\n",
     "    # Solve the system\n",
@@ -333,50 +333,6 @@
     "        break"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "72884c41",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "import petsc4py.PETSc as PETSc\n",
-    "Adof = model_cur.sc[\"A\"].dof_map\n",
-    "Avec = model_cur.sc[\"A\"].sol.vector()[Adof]\n",
-    "Bdof = model_cur.sc[\"B\"].dof_map\n",
-    "Bvec = model_cur.sc[\"B\"].sol.vector()[Bdof]\n",
-    "Xdof = model_cur.sc[\"X\"].dof_map\n",
-    "Xvec = model_cur.sc[\"X\"].sol.vector()[Xdof]\n",
-    "data_dict = dict()\n",
-    "data_dict[\"A\"] = Avec\n",
-    "data_dict[\"B\"] = Bvec\n",
-    "data_dict[\"X\"] = Xvec\n",
-    "\n",
-    "def weak_residual(model_cur, data_dict):\n",
-    "    for s in model_cur.sc:\n",
-    "        if s.name not in data_dict.keys():\n",
-    "            raise ValueError(\"Must provide vector for all species\")\n",
-    "        else:\n",
-    "            cur_data = data_dict[s.name]\n",
-    "        cur_vec = s.u[\"u\"].vector()\n",
-    "        cur_dof = s.dof_map\n",
-    "        cur_vec[cur_dof] = cur_data\n",
-    "        s.u[\"u\"].vector().set_local(cur_vec)\n",
-    "        s.u[\"u\"].vector().apply(\"insert\")\n",
-    "    if len(model_cur.problem.global_sizes) == 1:\n",
-    "        model_cur._ubackend = model_cur.u[\"u\"]._functions[0].vector().vec().copy()\n",
-    "    else:\n",
-    "        model_cur._ubackend = PETSc.Vec().createNest(\n",
-    "            [usub.vector().vec().copy() for usub in model_cur.u[\"u\"]._functions],\n",
-    "            comm=model_cur.mpi_comm_world,\n",
-    "        )\n",
-    "    return model_cur.solver.getFunctionNorm()\n",
-    "\n",
-    "print(f\"Residual is {weak_residual(model_cur, data_dict)}\")\n",
-    "data_dict[\"A\"] += 1.0\n",
-    "print(f\"Perturbed residual is {weak_residual(model_cur, data_dict)}\")"
-   ]
-  },
   {
    "attachments": {},
    "cell_type": "markdown",

examples/example8/example8.ipynb

@@ -6,7 +6,7 @@
    "id": "f65f18d7",
    "metadata": {},
    "source": [
-    "# Example 7: Advection-diffusion in a cylindrical geometry\n",
+    "# Example 8: Advection-diffusion in a cylindrical geometry\n",
     "\n",
     "Here, we consider release of a molecule from a cylindrical wall with Pouiselle flow."
    ]