diff --git a/01-adr-solver/adr-solver.ipynb b/01-adr-solver/adr-solver.ipynb
index 232352af35182c7b994012ca717eb3aa65c2107a..de49e56c2ed5702e40b162708fc68b78563f43f0 100644
--- a/01-adr-solver/adr-solver.ipynb
+++ b/01-adr-solver/adr-solver.ipynb
@@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "markdown",
- "id": "63389be3",
+ "id": "f5c98db2",
"metadata": {},
"source": [
"# Learning Objectives\n",
@@ -18,7 +18,7 @@
},
{
"cell_type": "markdown",
- "id": "e93d2ee3",
+ "id": "543ca6e9",
"metadata": {},
"source": [
"## Introduction\n",
@@ -47,7 +47,7 @@
},
{
"cell_type": "markdown",
- "id": "c4c20a91",
+ "id": "81212ea5",
"metadata": {},
"source": [
"# Nektar++ session files\n",
@@ -79,7 +79,7 @@
},
{
"cell_type": "markdown",
- "id": "4e0115d7",
+ "id": "641e9b33",
"metadata": {},
"source": [
"## Defining a mesh\n",
@@ -109,7 +109,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "fd75912a",
+ "id": "3b35b95a",
"metadata": {},
"outputs": [],
"source": [
@@ -118,7 +118,7 @@
},
{
"cell_type": "markdown",
- "id": "82bc7c81",
+ "id": "dcc84c6a",
"metadata": {},
"source": [
"## Visualising the mesh\n",
@@ -136,7 +136,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "efb78a7a",
+ "id": "9de79b7c",
"metadata": {},
"outputs": [],
"source": [
@@ -157,7 +157,7 @@
},
{
"cell_type": "markdown",
- "id": "28628b0b",
+ "id": "77664457",
"metadata": {},
"source": [
"Some important points to note in the visualisation of high-order meshes and fields:\n",
@@ -168,7 +168,7 @@
},
{
"cell_type": "markdown",
- "id": "c2f1a438",
+ "id": "5caf99c0",
"metadata": {},
"source": [
"# Steady diffusion problem\n",
@@ -193,7 +193,7 @@
},
{
"cell_type": "markdown",
- "id": "635fa17e",
+ "id": "717cdae4",
"metadata": {},
"source": [
"## Setting up the session file\n",
@@ -209,7 +209,7 @@
},
{
"cell_type": "markdown",
- "id": "45210822",
+ "id": "06ad0b05",
"metadata": {},
"source": [
"### Setting the polynomial order\n",
@@ -232,7 +232,7 @@
},
{
"cell_type": "markdown",
- "id": "cad5c53a",
+ "id": "c1e2f69a",
"metadata": {},
"source": [
"### Setting the equations to be solved\n",
@@ -256,7 +256,7 @@
},
{
"cell_type": "markdown",
- "id": "3f198fc7",
+ "id": "bd970cd9",
"metadata": {},
"source": [
"### Setting the boundary conditions\n",
@@ -289,7 +289,7 @@
},
{
"cell_type": "markdown",
- "id": "15cfe985",
+ "id": "5da50b5d",
"metadata": {},
"source": [
"### Setting the forcing function\n",
@@ -307,7 +307,7 @@
},
{
"cell_type": "markdown",
- "id": "d351923f",
+ "id": "4f74f445",
"metadata": {},
"source": [
"## Running the solver\n",
@@ -324,7 +324,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "32e86c90",
+ "id": "477c620f",
"metadata": {},
"outputs": [],
"source": [
@@ -333,7 +333,7 @@
},
{
"cell_type": "markdown",
- "id": "e3ae196c",
+ "id": "b47b2304",
"metadata": {},
"source": [
"## Post-processing\n",
@@ -356,7 +356,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "73015ebf",
+ "id": "34ac9a1b",
"metadata": {},
"outputs": [],
"source": [
@@ -365,7 +365,7 @@
},
{
"cell_type": "markdown",
- "id": "6daa307c",
+ "id": "81d55dc7",
"metadata": {},
"source": [
"Now we can use `PyVista`, combined with `itkwidgets` to visualise this in-browser as follows:"
@@ -374,7 +374,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "308cb160",
+ "id": "ada14670",
"metadata": {},
"outputs": [],
"source": [
@@ -391,7 +391,7 @@
},
{
"cell_type": "markdown",
- "id": "874aef66",
+ "id": "7d0c2e88",
"metadata": {},
"source": [
"## Running in parallel\n",
@@ -402,7 +402,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "46177f60",
+ "id": "cdc11460",
"metadata": {},
"outputs": [],
"source": [
@@ -411,7 +411,98 @@
},
{
"cell_type": "markdown",
- "id": "c885f726",
+ "id": "9e620b9c",
+ "metadata": {},
+ "source": [
+ "## Convergence analysis\n",
+ "\n",
+ "In this example we solved the laplacian equation $$ \\nabla^2 u = f $$ with forcing function $f(x,y) = -8\\pi^2 \\sin(2 \\pi x) \\sin(2\\pi y)$. The exact solution to this problem that satisfies the given boundary conditions is\n",
+ "\n",
+ "$u_e(x,y) = \\sin(2 \\pi x) \\sin(2\\pi y)$\n",
+ "\n",
+ "We can provide this exact solution in the [session file](sDiffSession.xml) and Nekar++ reports $L^2$ and $L^\\infty$ errors with respect to this exact solution at the end of the computation. The exact solution is defined in the [session file](sDiffSession.xml) inside the `<CONDITIONS>...</CONDITIONS>` as follows\n",
+ "\n",
+ "```xml\n",
+ "<FUNCTION NAME=\"ExactSolution\">\n",
+ " <E VAR=\"u\" VALUE=\"sin(2*PI*x)*sin(2*PI*y)\" />\n",
+ "</FUNCTION>\n",
+ "````\n",
+ "\n",
+ "\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "96827143",
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "import subprocess\n",
+ "import re\n",
+ "import numpy as np\n",
+ "import matplotlip.pyplot as plt\n",
+ "\n",
+ "expOrders = [\"exp2.xml\",\"exp3.xml\",\"exp4.xml\",\"exp5.xml\",\"exp6.xml\",\"exp7.xml\",\"exp8.xml\",\"exp9.xml\",\"exp10.xml\",\"exp11.xml\"]\n",
+ "\n",
+ "# first set up the array of Ndof (number of degree of freedom)\n",
+ "Ndof=[]\n",
+ "for e in expOrders:\n",
+ " cmd = \"FieldConvert -m dof AdrMesh.xml sDiffSession-Complete.xml \" + e + \" a.stdout\"\n",
+ " res = subprocess.check_output(cmd, shell=True, text=True)\n",
+ " ndof = ressplit('Total number of DOF:', maxsplit=1)[-1]\\\n",
+ " .split(maxsplit=1)[0]\n",
+ " Ndof.append(ndof)\n",
+ " \n",
+ "Ndof=np.array(Ndof, dtype=int)\n",
+ "\n",
+ "\n",
+ "L2Err=[]\n",
+ "LinfErr=[]\n",
+ "# now solving the equations using several polynomial orders and save L2 and Linf errors\n",
+ "for e in expOrders:\n",
+ " cmd=\"ADRSolver AdrMesh.xml sDiffSession-Complete.xml \" + e \n",
+ " # grab the std out\n",
+ " stdout = subprocess.check_output(cmd, shell=True, text=True)\n",
+ " \n",
+ " # get L2 error\n",
+ " l2 = stdout.split('L 2 error (variable u) :', maxsplit=1)[-1]\\\n",
+ " .split(maxsplit=1)[0]\n",
+ " L2Err.append(l2)\n",
+ " \n",
+ " # get the Linf error\n",
+ " lifn = stdout.split('L inf error (variable u) :', maxsplit=1)[-1]\\\n",
+ " .split(maxsplit=1)[0]\n",
+ " LinfErr.append(l2)\n",
+ "\n",
+ "# convert the results to numeric arrays\n",
+ "L2Err = np.array(L2Err, dtype=np.float)\n",
+ "LinfErr = np.array(LinfErr, dtype=np.float)\n",
+ "\n",
+ "\n",
+ "# Creating figure and axes\n",
+ "fig, (ax1, ax2) = plt.subplots(2, 1, figsize=[7, 11])\n",
+ "\n",
+ "ax1.loglog(Ndof, L2Err, '--r', linewidth=2, label='L2 Error')\n",
+ "ax1.set_title('L2 Error', fontsize=15)\n",
+ "ax1.set_xlabel('Number of Dof', fontsize=13)\n",
+ "ax1.set_ylabel('||E||2', fontsize=13)\n",
+ "ax1.legend()\n",
+ "\n",
+ "ax2.loglog(x, y, '--r', linewidth=2, label='Linf Error')\n",
+ "ax2.set_title('Linf Error', fontsize=15)\n",
+ "ax2.set_xlabel('Number of Dof', fontsize=13)\n",
+ "ax2.set_ylabel('||E||inf', fontsize=13)\n",
+ "ax2.legend()\n",
+ "\n",
+ "plt.tight_layout()\n",
+ "plt.show()"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "1a2ba4a0",
"metadata": {},
"source": [
"## Extra tasks\n",
@@ -433,7 +524,7 @@
},
{
"cell_type": "markdown",
- "id": "0cc9e976",
+ "id": "9b6c6efa",
"metadata": {},
"source": [
"# Unsteady diffusion problem\n",
@@ -457,7 +548,7 @@
},
{
"cell_type": "markdown",
- "id": "38761958",
+ "id": "d146e606",
"metadata": {},
"source": [
"## Changing the session file\n",
@@ -467,7 +558,7 @@
},
{
"cell_type": "markdown",
- "id": "703ad352",
+ "id": "a99318b1",
"metadata": {},
"source": [
"### Setting the solver type and timestepping scheme\n",
@@ -499,7 +590,7 @@
},
{
"cell_type": "markdown",
- "id": "3bc644b4",
+ "id": "fddf599e",
"metadata": {},
"source": [
"### Adding solution parameters\n",
@@ -549,7 +640,7 @@
},
{
"cell_type": "markdown",
- "id": "f8dc31fc",
+ "id": "61de6c04",
"metadata": {},
"source": [
"### Setting the boundary conditions\n",
@@ -562,7 +653,7 @@
},
{
"cell_type": "markdown",
- "id": "1548d346",
+ "id": "83f93165",
"metadata": {},
"source": [
"### Removing the forcing term\n",
@@ -576,7 +667,7 @@
},
{
"cell_type": "markdown",
- "id": "cfa30498",
+ "id": "3e4bb8b4",
"metadata": {},
"source": [
"### Setting initial conditions\n",
@@ -595,7 +686,7 @@
},
{
"cell_type": "markdown",
- "id": "4b90bdf5",
+ "id": "fa7f8437",
"metadata": {},
"source": [
"## Running the solver\n",
@@ -611,7 +702,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "a0b9b829",
+ "id": "d3ee6841",
"metadata": {
"scrolled": true
},
@@ -622,7 +713,7 @@
},
{
"cell_type": "markdown",
- "id": "938b0aff",
+ "id": "4c2c1831",
"metadata": {},
"source": [
"## Post-processing\n",
@@ -635,7 +726,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "d5429236",
+ "id": "ccea59bc",
"metadata": {},
"outputs": [],
"source": [
@@ -645,7 +736,7 @@
},
{
"cell_type": "markdown",
- "id": "a94e02ae",
+ "id": "031156c7",
"metadata": {},
"source": [
"Now that you converted all the results, you can visualise the evolution of `u` over the time by executing the cell below.\n",
@@ -660,7 +751,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "5794c608",
+ "id": "57706b89",
"metadata": {},
"outputs": [],
"source": [
@@ -692,7 +783,7 @@
},
{
"cell_type": "markdown",
- "id": "38abdfe3",
+ "id": "869ea56f",
"metadata": {},
"source": [
"The best to visualize the results is using **Paraview**, **Visit** or Tecplot. Since we converted our results to **vtu** it can be visualized with either **Paraview** or **VisIt**. We recommand you to install [Paraview](https://www.paraview.org/download/), then download the **vtu** files you converted here, Open the files in Paraview together and see the time-evolution of the solution.\n",
@@ -706,7 +797,7 @@
},
{
"cell_type": "markdown",
- "id": "abacc130",
+ "id": "57afcfff",
"metadata": {},
"source": [
"## Extra tasks\n",
@@ -774,7 +865,7 @@
},
{
"cell_type": "markdown",
- "id": "a01814ab",
+ "id": "b1710cc4",
"metadata": {},
"source": [
"# Unsteady advection-diffusion problem\n",
@@ -807,7 +898,7 @@
},
{
"cell_type": "markdown",
- "id": "2913c1fd",
+ "id": "c3d23861",
"metadata": {},
"source": [
"## Generating a new mesh\n",
@@ -822,7 +913,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "703d9d6a",
+ "id": "e65f9edc",
"metadata": {},
"outputs": [],
"source": [
@@ -831,7 +922,7 @@
},
{
"cell_type": "markdown",
- "id": "131ee096",
+ "id": "b79b3e6c",
"metadata": {},
"source": [
"`Gmsh` will produce the mesh and save it with the name [AdvDiff.msh](AdvDiff.msh). The resulting mesh [AdvDiff.msh](AdvDiff.msh) cannot be directly used with the Nektar++ solvers. We need to convert the `msh` format to the appropriate format for Nektar++ using the `NekMesh` utility. Let's convert the mesh and save it to a temporary file `AdvDiff-tmp.xml` as follows:"
@@ -840,7 +931,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "f687433f",
+ "id": "dc1ab3c3",
"metadata": {},
"outputs": [],
"source": [
@@ -849,7 +940,7 @@
},
{
"cell_type": "markdown",
- "id": "a35f56a1",
+ "id": "fdbce79d",
"metadata": {},
"source": [
"If you open the [AdvDiff-tmp.xml](AdvDiff-tmp.xml) file to have a look inside, you will notice that the information regarding definition of verticies and elements are in binary. However, at the end of the file you will see these lines:\n",
@@ -876,7 +967,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "ebb63fd4",
+ "id": "c7e286a7",
"metadata": {},
"outputs": [],
"source": [
@@ -885,7 +976,7 @@
},
{
"cell_type": "markdown",
- "id": "46bdff49",
+ "id": "eab3ade2",
"metadata": {},
"source": [
"Now we have the mesh file which is saved in [AdvDiffPeriodic.xml](AdvDiffPeriodic.xml) and its left and right, and bottom and top boundaries are paired respectively. We will use this mesh for our simulation."
@@ -893,7 +984,7 @@
},
{
"cell_type": "markdown",
- "id": "134071f0",
+ "id": "823fbee1",
"metadata": {},
"source": [
"## Changing the session file\n",
@@ -903,7 +994,7 @@
},
{
"cell_type": "markdown",
- "id": "b113b144",
+ "id": "8a99bedd",
"metadata": {},
"source": [
"### Adjusting the solver settings\n",
@@ -927,7 +1018,7 @@
},
{
"cell_type": "markdown",
- "id": "04ef55dd",
+ "id": "1040e44d",
"metadata": {},
"source": [
"### Solver parameters\n",
@@ -952,7 +1043,7 @@
},
{
"cell_type": "markdown",
- "id": "d7670d04",
+ "id": "9fdb6797",
"metadata": {},
"source": [
"### Setting the boundary conditions\n",
@@ -1011,7 +1102,7 @@
},
{
"cell_type": "markdown",
- "id": "cdd30a67",
+ "id": "4c91ec9a",
"metadata": {},
"source": [
"### Setting the advection velocity\n",
@@ -1029,7 +1120,7 @@
},
{
"cell_type": "markdown",
- "id": "7cdcc5b0",
+ "id": "36d7035f",
"metadata": {},
"source": [
"### Initial conditions\n",
@@ -1044,7 +1135,7 @@
},
{
"cell_type": "markdown",
- "id": "88bbfb4f",
+ "id": "1d6215b5",
"metadata": {},
"source": [
"### Exact solution\n",
@@ -1058,7 +1149,7 @@
},
{
"cell_type": "markdown",
- "id": "e0ff924f",
+ "id": "99b1c95e",
"metadata": {},
"source": [
"## Running the solver\n",
@@ -1068,7 +1159,7 @@
{
"cell_type": "code",
"execution_count": null,
- "id": "28670569",
+ "id": "906d7d44",
"metadata": {},
"outputs": [],
"source": [
@@ -1077,7 +1168,7 @@
},
{
"cell_type": "markdown",
- "id": "b02c57f6",
+ "id": "79907103",
"metadata": {},
"source": [
"## Post-processing\n",
@@ -1087,7 +1178,7 @@
},
{
"cell_type": "markdown",
- "id": "556aadc4",
+ "id": "758e7d58",
"metadata": {},
"source": [
"<div class=\"alert alert-success\">\n",
@@ -1098,7 +1189,7 @@
},
{
"cell_type": "markdown",
- "id": "ca3e378d",
+ "id": "a4a2cb50",
"metadata": {},
"source": [
"## Extra tasks\n",
@@ -1188,7 +1279,7 @@
},
{
"cell_type": "markdown",
- "id": "ba914262",
+ "id": "80907808",
"metadata": {},
"source": [
"### Next Step\n",
@@ -1204,7 +1295,7 @@
},
{
"cell_type": "markdown",
- "id": "825404fc",
+ "id": "e28b75c5",
"metadata": {},
"source": [
"### End of the tutorial. "
diff --git a/01-adr-solver/exp10.xml b/01-adr-solver/exp10.xml
new file mode 100755
index 0000000000000000000000000000000000000000..f0265af063dc269bb9bdf951de8410fcac01ab08
--- /dev/null
+++ b/01-adr-solver/exp10.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="10" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="10" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="10" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp11.xml b/01-adr-solver/exp11.xml
new file mode 100755
index 0000000000000000000000000000000000000000..7f5e95a576ab5ebf856b6c46bb54e9d49ff4a0e8
--- /dev/null
+++ b/01-adr-solver/exp11.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="11" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="11" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="11" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp2.xml b/01-adr-solver/exp2.xml
new file mode 100755
index 0000000000000000000000000000000000000000..6520678efe96f37192ff3bdb9778ce4119a77f36
--- /dev/null
+++ b/01-adr-solver/exp2.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="2" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="2" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="2" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp3.xml b/01-adr-solver/exp3.xml
new file mode 100755
index 0000000000000000000000000000000000000000..15f2b46a67ad1a9fbe038c4fe6d7f7c061413445
--- /dev/null
+++ b/01-adr-solver/exp3.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="3" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="3" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="3" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp4.xml b/01-adr-solver/exp4.xml
new file mode 100755
index 0000000000000000000000000000000000000000..8bf150ec95bf1c59aba9fb375117da5c8bed31ab
--- /dev/null
+++ b/01-adr-solver/exp4.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="4" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="4" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="4" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp5.xml b/01-adr-solver/exp5.xml
new file mode 100755
index 0000000000000000000000000000000000000000..e1e7bbb60420cb740cb09edbba6e1d70d5ec5abc
--- /dev/null
+++ b/01-adr-solver/exp5.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="5" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="5" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="5" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp6.xml b/01-adr-solver/exp6.xml
new file mode 100755
index 0000000000000000000000000000000000000000..47755f79cd6cfa8f496bcba00771556af600f333
--- /dev/null
+++ b/01-adr-solver/exp6.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="6" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="6" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="6" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp7.xml b/01-adr-solver/exp7.xml
new file mode 100755
index 0000000000000000000000000000000000000000..a02fe76f41f423b08d2361ceb5ddbf6d08dca359
--- /dev/null
+++ b/01-adr-solver/exp7.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="7" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="7" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="7" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp8.xml b/01-adr-solver/exp8.xml
new file mode 100755
index 0000000000000000000000000000000000000000..e59e1ff6dfa001456e6a371dffc1b79498bcd0f5
--- /dev/null
+++ b/01-adr-solver/exp8.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="8" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="8" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="8" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>
diff --git a/01-adr-solver/exp9.xml b/01-adr-solver/exp9.xml
new file mode 100755
index 0000000000000000000000000000000000000000..f5f5c68c9544ad257c6a92e3668be9e950cb97f5
--- /dev/null
+++ b/01-adr-solver/exp9.xml
@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<NEKTAR>
+ <EXPANSIONS>
+ <E COMPOSITE="C[1]" NUMMODES="9" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[2]" NUMMODES="9" TYPE="MODIFIED" FIELDS="u" />
+ <E COMPOSITE="C[3]" NUMMODES="9" TYPE="MODIFIED" FIELDS="u" />
+ </EXPANSIONS>
+</NEKTAR>