Skip to content

Change FEniCS tutorials to handle new initialization signature of FEniCS-Adapter #119

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 4 commits into from
Nov 28, 2020
Merged

Change FEniCS tutorials to handle new initialization signature of FEniCS-Adapter #119

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
41b1a46
Initial change to FEniCS-Adapter initialization
IshaanDesai Nov 17, 2020
ca1b1dd
Merging write_function_space and write_function into write_object
IshaanDesai Nov 24, 2020
42bd9c5
Explicitly mentioning optional parameter names in adapter initialization
IshaanDesai Nov 25, 2020
860642a
Adding comment to explain why a function space is passed twice to ini…
IshaanDesai Nov 28, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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, read_function_space=V, write_object=f_N_function)

# Create a FEniCS Expression to define and control the coupling boundary values
coupling_expression = precice.create_coupling_expression()
Expand Down
4 changes: 3 additions & 1 deletion FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +88,9 @@ 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)
# Function space V is passed twice as both read and write functions are defined using the same space
precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=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
Expand Down
4 changes: 3 additions & 1 deletion FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,9 @@ 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)
# Function space V is passed twice as both read and write functions are defined using the same space
precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=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
Expand Down
4 changes: 2 additions & 2 deletions HT/partitioned-heat/fenics-fenics/heat.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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, read_function_space=V, write_object=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, read_function_space=V_g, write_object=u_D_function)

boundary_marker = False

Expand Down