From 68f084c601250957881afcec3b8f67f02571bbe7 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@in.tum.de>
Date: Tue, 12 Jul 2022 16:41:48 +0200
Subject: [PATCH 1/9] Use linear iteration and subcycling.

---
 partitioned-heat-conduction/fenics/heat.py     | 4 ++--
 partitioned-heat-conduction/precice-config.xml | 6 +++---
 2 files changed, 5 insertions(+), 5 deletions(-)

diff --git a/partitioned-heat-conduction/fenics/heat.py b/partitioned-heat-conduction/fenics/heat.py
index e7b7f6529..70c82dcef 100644
--- a/partitioned-heat-conduction/fenics/heat.py
+++ b/partitioned-heat-conduction/fenics/heat.py
@@ -62,7 +62,7 @@ def determine_gradient(V_g, u, flux):
 
 args = parser.parse_args()
 
-fenics_dt = .1  # time step size
+fenics_dt = .01  # time step size
 # Error is bounded by coupling accuracy. In theory we would obtain the analytical solution.
 error_tol = args.error_tol
 
@@ -178,7 +178,7 @@ def determine_gradient(V_g, u, flux):
     if precice.is_action_required(precice.action_write_iteration_checkpoint()):
         precice.store_checkpoint(u_n, t, n)
 
-    read_data = precice.read_data()
+    read_data = precice.read_data(dt(0))
 
     # Update the coupling expression with the new read data
     precice.update_coupling_expression(coupling_expression, read_data)
diff --git a/partitioned-heat-conduction/precice-config.xml b/partitioned-heat-conduction/precice-config.xml
index 280fcb2f2..510916fc9 100644
--- a/partitioned-heat-conduction/precice-config.xml
+++ b/partitioned-heat-conduction/precice-config.xml
@@ -7,7 +7,7 @@
       enabled="true" />
   </log>
 
-  <solver-interface dimensions="2">
+  <solver-interface dimensions="2" experimental="true">
     <data:scalar name="Temperature" />
     <data:scalar name="Heat-Flux" />
 
@@ -25,7 +25,7 @@
       <use-mesh name="Dirichlet-Mesh" provide="yes" />
       <use-mesh name="Neumann-Mesh" from="Neumann" />
       <write-data name="Heat-Flux" mesh="Dirichlet-Mesh" />
-      <read-data name="Temperature" mesh="Dirichlet-Mesh" />
+      <read-data name="Temperature" mesh="Dirichlet-Mesh" waveform-order="1"/>
       <mapping:rbf-thin-plate-splines
         direction="read"
         from="Neumann-Mesh"
@@ -38,7 +38,7 @@
       <use-mesh name="Neumann-Mesh" provide="yes" />
       <use-mesh name="Dirichlet-Mesh" from="Dirichlet" />
       <write-data name="Temperature" mesh="Neumann-Mesh" />
-      <read-data name="Heat-Flux" mesh="Neumann-Mesh" />
+      <read-data name="Heat-Flux" mesh="Neumann-Mesh" waveform-order="0"/>
       <mapping:rbf-thin-plate-splines
         direction="read"
         from="Dirichlet-Mesh"

From 8f02697ff4f8b3831db29f23d7f7a63224fb1911 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Wed, 27 Sep 2023 22:52:14 +0200
Subject: [PATCH 2/9] Update to current state of software stack.

---
 partitioned-heat-conduction/fenics/heat.py    |  13 +-
 .../precice-config.xml                        | 121 +++++++++---------
 2 files changed, 70 insertions(+), 64 deletions(-)

diff --git a/partitioned-heat-conduction/fenics/heat.py b/partitioned-heat-conduction/fenics/heat.py
index 70c82dcef..466447ef9 100644
--- a/partitioned-heat-conduction/fenics/heat.py
+++ b/partitioned-heat-conduction/fenics/heat.py
@@ -101,10 +101,12 @@ def determine_gradient(V_g, u, flux):
 # Initialize the adapter according to the specific participant
 if problem is ProblemType.DIRICHLET:
     precice = Adapter(adapter_config_filename="precice-adapter-config-D.json")
-    precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=f_N_function)
+    precice.initialize(coupling_boundary, read_function_space=V, write_object=f_N_function)
+    precice_dt = precice.get_max_time_step_size()
 elif problem is ProblemType.NEUMANN:
     precice = Adapter(adapter_config_filename="precice-adapter-config-N.json")
-    precice_dt = precice.initialize(coupling_boundary, read_function_space=W, write_object=u_D_function)
+    precice.initialize(coupling_boundary, read_function_space=W, write_object=u_D_function)
+    precice_dt = precice.get_max_time_step_size()
 
 dt = Constant(0)
 dt.assign(np.min([fenics_dt, precice_dt]))
@@ -175,7 +177,7 @@ def determine_gradient(V_g, u, flux):
 while precice.is_coupling_ongoing():
 
     # write checkpoint
-    if precice.is_action_required(precice.action_write_iteration_checkpoint()):
+    if precice.requires_writing_checkpoint():
         precice.store_checkpoint(u_n, t, n)
 
     read_data = precice.read_data(dt(0))
@@ -198,10 +200,11 @@ def determine_gradient(V_g, u, flux):
         # Neumann problem reads flux and writes temperature on boundary to Dirichlet problem
         precice.write_data(u_np1)
 
-    precice_dt = precice.advance(dt(0))
+    precice.advance(dt(0))
+    precice_dt = precice.get_max_time_step_size()
 
     # roll back to checkpoint
-    if precice.is_action_required(precice.action_read_iteration_checkpoint()):
+    if precice.requires_reading_checkpoint():
         u_cp, t_cp, n_cp = precice.retrieve_checkpoint()
         u_n.assign(u_cp)
         t = t_cp
diff --git a/partitioned-heat-conduction/precice-config.xml b/partitioned-heat-conduction/precice-config.xml
index 510916fc9..c6c7bd9aa 100644
--- a/partitioned-heat-conduction/precice-config.xml
+++ b/partitioned-heat-conduction/precice-config.xml
@@ -7,69 +7,72 @@
       enabled="true" />
   </log>
 
-  <solver-interface dimensions="2" experimental="true">
-    <data:scalar name="Temperature" />
-    <data:scalar name="Heat-Flux" />
+  <data:scalar name="Temperature" waveform-degree="1" />
+  <data:scalar name="Heat-Flux" waveform-degree="1" />
 
-    <mesh name="Dirichlet-Mesh">
-      <use-data name="Temperature" />
-      <use-data name="Heat-Flux" />
-    </mesh>
+  <mesh name="Dirichlet-Mesh" dimensions="2">
+    <use-data name="Temperature" />
+    <use-data name="Heat-Flux" />
+  </mesh>
 
-    <mesh name="Neumann-Mesh">
-      <use-data name="Temperature" />
-      <use-data name="Heat-Flux" />
-    </mesh>
+  <mesh name="Neumann-Mesh" dimensions="2">
+    <use-data name="Temperature" />
+    <use-data name="Heat-Flux" />
+  </mesh>
 
-    <participant name="Dirichlet">
-      <use-mesh name="Dirichlet-Mesh" provide="yes" />
-      <use-mesh name="Neumann-Mesh" from="Neumann" />
-      <write-data name="Heat-Flux" mesh="Dirichlet-Mesh" />
-      <read-data name="Temperature" mesh="Dirichlet-Mesh" waveform-order="1"/>
-      <mapping:rbf-thin-plate-splines
-        direction="read"
-        from="Neumann-Mesh"
-        to="Dirichlet-Mesh"
-        constraint="consistent"
-        x-dead="true" />
-    </participant>
+  <participant name="Dirichlet">
+    <provide-mesh name="Dirichlet-Mesh" />
+    <receive-mesh name="Neumann-Mesh" from="Neumann" />
+    <write-data name="Heat-Flux" mesh="Dirichlet-Mesh" />
+    <read-data name="Temperature" mesh="Dirichlet-Mesh" />
+    <mapping:rbf
+      direction="read"
+      from="Neumann-Mesh"
+      to="Dirichlet-Mesh"
+      constraint="consistent"
+      x-dead="true" >
+      <basis-function:thin-plate-splines />
+    </mapping:rbf>
+  </participant>
 
-    <participant name="Neumann">
-      <use-mesh name="Neumann-Mesh" provide="yes" />
-      <use-mesh name="Dirichlet-Mesh" from="Dirichlet" />
-      <write-data name="Temperature" mesh="Neumann-Mesh" />
-      <read-data name="Heat-Flux" mesh="Neumann-Mesh" waveform-order="0"/>
-      <mapping:rbf-thin-plate-splines
-        direction="read"
-        from="Dirichlet-Mesh"
-        to="Neumann-Mesh"
-        constraint="consistent"
-        x-dead="true" />
-    </participant>
+  <participant name="Neumann">
+    <provide-mesh name="Neumann-Mesh" />
+    <receive-mesh name="Dirichlet-Mesh" from="Dirichlet" />
+    <write-data name="Temperature" mesh="Neumann-Mesh" />
+    <read-data name="Heat-Flux" mesh="Neumann-Mesh" />
+    <mapping:rbf
+      direction="read"
+      from="Dirichlet-Mesh"
+      to="Neumann-Mesh"
+      constraint="consistent"
+      x-dead="true" >
+      <basis-function:thin-plate-splines />
+    </mapping:rbf>
+  </participant>
 
-    <m2n:sockets from="Dirichlet" to="Neumann" exchange-directory=".." />
+  <m2n:sockets acceptor="Dirichlet" connector="Neumann" exchange-directory=".." />
 
-    <coupling-scheme:serial-implicit>
-      <participants first="Dirichlet" second="Neumann" />
-      <max-time value="1.0" />
-      <time-window-size value="0.1" />
-      <max-iterations value="100" />
-      <exchange data="Heat-Flux" mesh="Dirichlet-Mesh" from="Dirichlet" to="Neumann" />
-      <exchange
-        data="Temperature"
-        mesh="Neumann-Mesh"
-        from="Neumann"
-        to="Dirichlet"
-        initialize="true" />
-      <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-5" />
-      <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-5" />
-      <acceleration:IQN-ILS>
-        <data name="Temperature" mesh="Neumann-Mesh" />
-        <initial-relaxation value="0.1" />
-        <max-used-iterations value="10" />
-        <time-windows-reused value="5" />
-        <filter type="QR2" limit="1e-3" />
-      </acceleration:IQN-ILS>
-    </coupling-scheme:serial-implicit>
-  </solver-interface>
+  <coupling-scheme:serial-implicit>
+    <participants first="Dirichlet" second="Neumann" />
+    <max-time value="1.0" />
+    <time-window-size value="0.1" />
+    <max-iterations value="100" />
+    <exchange data="Heat-Flux" mesh="Dirichlet-Mesh" from="Dirichlet" to="Neumann" substeps="true" />
+    <exchange
+      data="Temperature"
+      mesh="Neumann-Mesh"
+      from="Neumann"
+      to="Dirichlet"
+      initialize="true"
+      substeps="true" />
+    <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-5" />
+    <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-5" />
+    <acceleration:IQN-ILS>
+      <data name="Temperature" mesh="Neumann-Mesh" />
+      <initial-relaxation value="0.1" />
+      <max-used-iterations value="10" />
+      <time-windows-reused value="5" />
+      <filter type="QR2" limit="1e-3" />
+    </acceleration:IQN-ILS>
+  </coupling-scheme:serial-implicit>
 </precice-configuration>

From 4d552aaf77d002785c31aa4018793af79ffa0b69 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Wed, 27 Sep 2023 23:19:19 +0200
Subject: [PATCH 3/9] Introduce change from #333.

---
 partitioned-heat-conduction/fenics/problem_setup.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/partitioned-heat-conduction/fenics/problem_setup.py b/partitioned-heat-conduction/fenics/problem_setup.py
index b85f30358..ff2996fce 100644
--- a/partitioned-heat-conduction/fenics/problem_setup.py
+++ b/partitioned-heat-conduction/fenics/problem_setup.py
@@ -33,7 +33,7 @@ def inside(self, x, on_boundary):
 
 
 def get_geometry(domain_part):
-    nx = ny = 9
+    nx = ny = 11
 
     if domain_part is DomainPart.LEFT:
         p0 = Point(x_left, y_bottom)

From 39762626625079e60502567c2860551f7575c372 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Thu, 28 Sep 2023 10:16:30 +0200
Subject: [PATCH 4/9] Use constant underrelaxation until
 https://github.com/precice/precice/issues/1689 works.

---
 partitioned-heat-conduction/precice-config.xml | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/partitioned-heat-conduction/precice-config.xml b/partitioned-heat-conduction/precice-config.xml
index 54162e2df..d465dbd3d 100644
--- a/partitioned-heat-conduction/precice-config.xml
+++ b/partitioned-heat-conduction/precice-config.xml
@@ -67,12 +67,15 @@
       substeps="true" />
     <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-5" />
     <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-5" />
-    <acceleration:IQN-ILS>
+    <!-- <acceleration:IQN-ILS>
       <data name="Temperature" mesh="Neumann-Mesh" />
       <initial-relaxation value="0.1" />
       <max-used-iterations value="10" />
       <time-windows-reused value="5" />
       <filter type="QR2" limit="1e-3" />
-    </acceleration:IQN-ILS>
+    </acceleration:IQN-ILS> -->
+    <acceleration:constant>
+      <relaxation value="0.5"/>
+    </acceleration:constant>
   </coupling-scheme:serial-implicit>
 </precice-configuration>

From 330136d66af5e6bf8900ed1dc13adc1e59199e82 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Thu, 28 Sep 2023 19:20:18 +0200
Subject: [PATCH 5/9] Adjust convergence limit for constant underrelaxation.

---
 partitioned-heat-conduction/precice-config.xml | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/partitioned-heat-conduction/precice-config.xml b/partitioned-heat-conduction/precice-config.xml
index d465dbd3d..2d9b9983a 100644
--- a/partitioned-heat-conduction/precice-config.xml
+++ b/partitioned-heat-conduction/precice-config.xml
@@ -65,8 +65,8 @@
       to="Dirichlet"
       initialize="true"
       substeps="true" />
-    <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-5" />
-    <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-5" />
+    <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-10" />
+    <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-10" />
     <!-- <acceleration:IQN-ILS>
       <data name="Temperature" mesh="Neumann-Mesh" />
       <initial-relaxation value="0.1" />

From 6a75a9663ee3c37337051655beb8b4c78d296de0 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Tue, 14 Nov 2023 10:58:20 +0100
Subject: [PATCH 6/9] Also write data at intermediate time steps.

---
 partitioned-heat-conduction/fenics/heat.py | 64 ++++++++++++++++++----
 1 file changed, 54 insertions(+), 10 deletions(-)

diff --git a/partitioned-heat-conduction/fenics/heat.py b/partitioned-heat-conduction/fenics/heat.py
index 1ac0827df..97c379d2c 100644
--- a/partitioned-heat-conduction/fenics/heat.py
+++ b/partitioned-heat-conduction/fenics/heat.py
@@ -157,13 +157,25 @@ def determine_gradient(V_g, u, flux):
 
 # output solution and reference solution at t=0, n=0
 n = 0
-print('output u^%d and u_ref^%d' % (n, n))
-temperature_out << (u_n, t)
-ref_out << u_ref
+print("output u^%d and u_ref^%d" % (n, n))
 ranks << mesh_rank
 
 error_total, error_pointwise = compute_errors(u_n, u_ref, V)
-error_out << error_pointwise
+
+# create buffer for output. We need this buffer, because we only want to write the converged output at the end of the window, but we also want to write the samples that are resulting from substeps inside the window
+u_write = []
+ref_write = []
+error_write = []
+# copy data to buffer and rename
+uu = u_n.copy()
+uu.rename("u","")
+u_write.append((uu, t))
+uu_ref = u_ref.copy()
+uu_ref.rename("u_ref","")
+ref_write.append(uu_ref)
+err = error_pointwise.copy()
+err.rename("err","")
+error_write.append(err)
 
 # set t_1 = t_0 + dt, this gives u_D^1
 # call dt(0) to evaluate FEniCS Constant. Todo: is there a better way?
@@ -180,6 +192,17 @@ def determine_gradient(V_g, u, flux):
     if precice.requires_writing_checkpoint():
         precice.store_checkpoint(u_n, t, n)
 
+        # output solution and reference solution at t_n+1 and substeps (read from buffer)
+        print('output u^%d and u_ref^%d' % (n, n))
+        for sample in u_write:
+            temperature_out << sample
+
+        for sample in ref_write:
+            ref_out << sample
+
+        for sample in error_write:
+            error_out << error_pointwise
+
     dt.assign(np.min([fenics_dt, precice_dt]))
     read_data = precice.read_data(dt)
 
@@ -208,25 +231,46 @@ def determine_gradient(V_g, u, flux):
         u_n.assign(u_cp)
         t = t_cp
         n = n_cp
+        # empty buffer if window has not converged
+        u_write = []
+        ref_write = []
+        error_write = []
     else:  # update solution
         u_n.assign(u_np1)
         t += float(dt)
         n += 1
+        # copy data to buffer and rename
+        uu = u_n.copy()
+        uu.rename("u","u")
+        u_write.append((uu, t))
+        uu_ref = u_ref.copy()
+        uu_ref.rename("u_ref","u_ref")
+        ref_write.append(uu_ref)
+        err = error_pointwise.copy()
+        err.rename("err","err")
+        error_write.append(err)
 
     if precice.is_time_window_complete():
         u_ref = interpolate(u_D, V)
         u_ref.rename("reference", " ")
         error, error_pointwise = compute_errors(u_n, u_ref, V, total_error_tol=error_tol)
-        print('n = %d, t = %.2f: L2 error on domain = %.3g' % (n, t, error))
-        # output solution and reference solution at t_n+1
-        print('output u^%d and u_ref^%d' % (n, n))
-        temperature_out << (u_n, t)
-        ref_out << u_ref
-        error_out << error_pointwise
+        print("n = %d, t = %.2f: L2 error on domain = %.3g" % (n, t, error))
+
 
     # Update Dirichlet BC
     u_D.t = t + float(dt)
     f.t = t + float(dt)
 
+# output solution and reference solution at t_n+1 and substeps (read from buffer)
+print("output u^%d and u_ref^%d" % (n, n))
+for sample in u_write:
+    temperature_out << sample
+
+for sample in ref_write:
+    ref_out << sample
+
+for sample in error_write:
+    error_out << error_pointwise
+
 # Hold plot
 precice.finalize()

From 28c344278d22943354b66bd82bdc561b58129b9d Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Tue, 14 Nov 2023 11:12:04 +0100
Subject: [PATCH 7/9] Fix format.

---
 partitioned-heat-conduction/fenics/heat.py | 17 +++++++++--------
 1 file changed, 9 insertions(+), 8 deletions(-)

diff --git a/partitioned-heat-conduction/fenics/heat.py b/partitioned-heat-conduction/fenics/heat.py
index 14595bcd7..02b27eee8 100644
--- a/partitioned-heat-conduction/fenics/heat.py
+++ b/partitioned-heat-conduction/fenics/heat.py
@@ -168,19 +168,21 @@ def determine_gradient(V_g, u, flux):
 
 error_total, error_pointwise = compute_errors(u_n, u_ref, V)
 
-# create buffer for output. We need this buffer, because we only want to write the converged output at the end of the window, but we also want to write the samples that are resulting from substeps inside the window
+# create buffer for output. We need this buffer, because we only want to
+# write the converged output at the end of the window, but we also want to
+# write the samples that are resulting from substeps inside the window
 u_write = []
 ref_write = []
 error_write = []
 # copy data to buffer and rename
 uu = u_n.copy()
-uu.rename("u","")
+uu.rename("u", "")
 u_write.append((uu, t))
 uu_ref = u_ref.copy()
-uu_ref.rename("u_ref","")
+uu_ref.rename("u_ref", "")
 ref_write.append(uu_ref)
 err = error_pointwise.copy()
-err.rename("err","")
+err.rename("err", "")
 error_write.append(err)
 
 # set t_1 = t_0 + dt, this gives u_D^1
@@ -254,13 +256,13 @@ def determine_gradient(V_g, u, flux):
         n += 1
         # copy data to buffer and rename
         uu = u_n.copy()
-        uu.rename("u","u")
+        uu.rename("u", "")
         u_write.append((uu, t))
         uu_ref = u_ref.copy()
-        uu_ref.rename("u_ref","u_ref")
+        uu_ref.rename("u_ref", "")
         ref_write.append(uu_ref)
         err = error_pointwise.copy()
-        err.rename("err","err")
+        err.rename("err", "")
         error_write.append(err)
 
     if precice.is_time_window_complete():
@@ -269,7 +271,6 @@ def determine_gradient(V_g, u, flux):
         error, error_pointwise = compute_errors(u_n, u_ref, V, total_error_tol=error_tol)
         print("n = %d, t = %.2f: L2 error on domain = %.3g" % (n, t, error))
 
-
     # Update Dirichlet BC
     u_D.t = t + float(dt)
     f.t = t + float(dt)

From 76c43b9e1e452abd56a57b7786b2b6405b0f0ebc Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Mon, 20 Nov 2023 11:22:25 +0100
Subject: [PATCH 8/9] Initialize for both participants.

---
 partitioned-heat-conduction/precice-config.xml | 10 ++--------
 1 file changed, 2 insertions(+), 8 deletions(-)

diff --git a/partitioned-heat-conduction/precice-config.xml b/partitioned-heat-conduction/precice-config.xml
index c1e5e3a28..f1da32d1a 100644
--- a/partitioned-heat-conduction/precice-config.xml
+++ b/partitioned-heat-conduction/precice-config.xml
@@ -43,14 +43,8 @@
     <max-time value="1.0" />
     <time-window-size value="0.1" />
     <max-iterations value="100" />
-    <exchange data="Heat-Flux" mesh="Dirichlet-Mesh" from="Dirichlet" to="Neumann" substeps="true" />
-    <exchange
-      data="Temperature"
-      mesh="Neumann-Mesh"
-      from="Neumann"
-      to="Dirichlet"
-      initialize="true"
-      substeps="true" />
+    <exchange data="Heat-Flux" mesh="Dirichlet-Mesh" from="Dirichlet" to="Neumann" initialize="true" substeps="true" />
+    <exchange data="Temperature" mesh="Neumann-Mesh" from="Neumann" to="Dirichlet" initialize="true" substeps="true" />
     <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-10" />
     <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-10" />
     <!-- <acceleration:IQN-ILS>

From 4bb07e76a31ebf3ce0a1f3be7361863bf721a803 Mon Sep 17 00:00:00 2001
From: Benjamin Rodenberg <benjamin.rodenberg@cit.tum.de>
Date: Mon, 20 Nov 2023 11:32:01 +0100
Subject: [PATCH 9/9] Fix format.

---
 .../precice-config.xml                        | 30 +++++++++++++++----
 1 file changed, 25 insertions(+), 5 deletions(-)

diff --git a/partitioned-heat-conduction/precice-config.xml b/partitioned-heat-conduction/precice-config.xml
index f1da32d1a..baab5d5c3 100644
--- a/partitioned-heat-conduction/precice-config.xml
+++ b/partitioned-heat-conduction/precice-config.xml
@@ -25,7 +25,11 @@
     <receive-mesh name="Neumann-Mesh" from="Neumann" />
     <write-data name="Heat-Flux" mesh="Dirichlet-Mesh" />
     <read-data name="Temperature" mesh="Dirichlet-Mesh" />
-    <mapping:nearest-neighbor direction="read" from="Neumann-Mesh" to="Dirichlet-Mesh" constraint="consistent" />
+    <mapping:nearest-neighbor
+      direction="read"
+      from="Neumann-Mesh"
+      to="Dirichlet-Mesh"
+      constraint="consistent" />
   </participant>
 
   <participant name="Neumann">
@@ -33,7 +37,11 @@
     <receive-mesh name="Dirichlet-Mesh" from="Dirichlet" />
     <write-data name="Temperature" mesh="Neumann-Mesh" />
     <read-data name="Heat-Flux" mesh="Neumann-Mesh" />
-    <mapping:nearest-neighbor direction="read" from="Dirichlet-Mesh" to="Neumann-Mesh" constraint="consistent" />
+    <mapping:nearest-neighbor
+      direction="read"
+      from="Dirichlet-Mesh"
+      to="Neumann-Mesh"
+      constraint="consistent" />
   </participant>
 
   <m2n:sockets acceptor="Dirichlet" connector="Neumann" exchange-directory=".." />
@@ -43,8 +51,20 @@
     <max-time value="1.0" />
     <time-window-size value="0.1" />
     <max-iterations value="100" />
-    <exchange data="Heat-Flux" mesh="Dirichlet-Mesh" from="Dirichlet" to="Neumann" initialize="true" substeps="true" />
-    <exchange data="Temperature" mesh="Neumann-Mesh" from="Neumann" to="Dirichlet" initialize="true" substeps="true" />
+    <exchange
+      data="Heat-Flux"
+      mesh="Dirichlet-Mesh"
+      from="Dirichlet"
+      to="Neumann"
+      initialize="true"
+      substeps="true" />
+    <exchange
+      data="Temperature"
+      mesh="Neumann-Mesh"
+      from="Neumann"
+      to="Dirichlet"
+      initialize="true"
+      substeps="true" />
     <relative-convergence-measure data="Heat-Flux" mesh="Dirichlet-Mesh" limit="1e-10" />
     <relative-convergence-measure data="Temperature" mesh="Neumann-Mesh" limit="1e-10" />
     <!-- <acceleration:IQN-ILS>
@@ -55,7 +75,7 @@
       <filter type="QR2" limit="1e-3" />
     </acceleration:IQN-ILS> -->
     <acceleration:constant>
-      <relaxation value="0.5"/>
+      <relaxation value="0.5" />
     </acceleration:constant>
   </coupling-scheme:serial-implicit>
 </precice-configuration>