You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
summary: We solve a simple heat equation. The domain is partitioned and the coupling is established in a Dirichlet-Neumann fashion.
6
6
---
7
7
@@ -25,8 +25,6 @@ This simple case allows us to compare the solution for the partitioned case to a
25
25
26
26
You can either couple a solver with itself or different solvers with each other. In any case you will need to have preCICE and the python bindings installed on your system.
27
27
28
-
* FEniCSx. Install [FEniCSx](https://fenicsproject.org/download/) and the [FEniCSx-adapter](https://github.com/precice/fenicsx-adapter). The code is largely based on this [fenics-tutorial](https://github.com/hplgit/fenics-tutorial/blob/master/pub/python/vol1/ft03_heat.py) from [1] and has been adapted to FEniCSx.
29
-
30
28
* FEniCS. Install [FEniCS](https://fenicsproject.org/download/) and the [FEniCS-adapter](https://github.com/precice/fenics-adapter). The code is largely based on this [fenics-tutorial](https://github.com/hplgit/fenics-tutorial/blob/master/pub/python/vol1/ft03_heat.py) from [1].
@@ -45,18 +43,18 @@ For choosing whether you want to run the Dirichlet-kind and a Neumann-kind parti
45
43
For running the case, open two terminals run:
46
44
47
45
```bash
48
-
cdfenicsx
46
+
cdfenics
49
47
./run.sh -d
50
48
```
51
49
52
50
and
53
51
54
52
```bash
55
-
cdfenicsx
53
+
cdfenics
56
54
./run.sh -n
57
55
```
58
56
59
-
If you want to use FEniCS or Nutils, use `cd fenics`/ `cd nutils` instead of `cd fenicsx`. The FEniCS case also supports parallel runs. Here, you cannot use the `run.sh` script, but must simply execute
57
+
If you want to use Nutils, use `cd nutils` instead of `cd fenics`. The FEniCS case also supports parallel runs. Here, you cannot use the `run.sh` script, but must simply execute
60
58
61
59
```bash
62
60
mpirun -n <N_PROC> heat.py -d
@@ -68,9 +66,7 @@ You can mix the Nutils and FEniCS solver, if you like. Note that the error for a
68
66
69
67
## Visualization
70
68
71
-
Output is written into the folders `fenicsx/out`, `fenics/out` and `nutils`.
72
-
73
-
For FEniCSx you can visualize the content with paraview by opening the `*.xdmf` files. The files `Dirichlet.xdmf` and `Neumann.xdmf` correspond to the numerical solution of the Dirichlet, respectively Neumann, problem, while the files with the prefix `ref` correspond to the analytical reference solution.
69
+
Output is written into the folders `fenics/out` and `nutils`.
74
70
75
71
For FEniCS you can visualize the content with paraview by opening the `*.pvd` files. The files `Dirichlet.pvd` and `Neumann.pvd` correspond to the numerical solution of the Dirichlet, respectively Neumann, problem, while the files with the prefix `ref` correspond to the analytical reference solution, the files with `error` show the error and the files with `ranks` the ranks of the solvers (if executed in parallel).
0 commit comments