Initial change to FEniCS-Adapter initialization

Author: IshaanDesai

CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py

@@ -116,7 +116,7 @@ def determine_heat_flux(V_g, u, k, flux):
 # Adapter definition and initialization
 precice = Adapter(adapter_config_filename="precice-adapter-config.json")
 
-precice_dt = precice.initialize(coupling_boundary, mesh, V, write_function=f_N_function)
+precice_dt = precice.initialize(coupling_boundary, V, V_g, write_function=f_N_function)
 
 # Create a FEniCS Expression to define and control the coupling boundary values
 coupling_expression = precice.create_coupling_expression()

FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py

@@ -88,7 +88,7 @@ def remaining_boundary(x, on_boundary):
 force_boundary = AutoSubDomain(remaining_boundary)
 
 # Initialize the coupling interface
-precice_dt = precice.initialize(coupling_boundary, mesh, V, dim, fixed_boundary=clamped_boundary_domain)
+precice_dt = precice.initialize(coupling_boundary, V, V, clamped_boundary_domain)
 
 fenics_dt = precice_dt  # if fenics_dt == precice_dt, no subcycling is applied
 # fenics_dt = 0.02  # if fenics_dt < precice_dt, subcycling is applied

FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py

@@ -72,7 +72,7 @@ def neumann_boundary(x, on_boundary):
 force_boundary = AutoSubDomain(neumann_boundary)
 
 # Initialize the coupling interface
-precice_dt = precice.initialize(coupling_boundary, mesh, V, dim, fixed_boundary=clamped_boundary_domain)
+precice_dt = precice.initialize(coupling_boundary, V, V, fixed_boundary=clamped_boundary_domain)
 
 fenics_dt = precice_dt  # if fenics_dt == precice_dt, no subcycling is applied
 # fenics_dt = 0.02  # if fenics_dt < precice_dt, subcycling is applied

HT/partitioned-heat/fenics-fenics/heat.py

@@ -124,10 +124,10 @@ 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, mesh, V, write_function=f_N_function)
+    precice_dt = precice.initialize(coupling_boundary, V, V_g, write_function=f_N_function)
 elif problem is ProblemType.NEUMANN:
     precice = Adapter(adapter_config_filename="precice-adapter-config-N.json")
-    precice_dt = precice.initialize(coupling_boundary, mesh, V_g, write_function=u_D_function)
+    precice_dt = precice.initialize(coupling_boundary, V_g, V, write_function=u_D_function)
 
 boundary_marker = False