Skip to content

Fix ill-posed problem by providing additional data. #144

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 1 commit into from
Jan 8, 2021
Merged

Fix ill-posed problem by providing additional data. #144

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
6 changes: 5 additions & 1 deletion HT/partitioned-heat/fenics-fenics/heat.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@
File, solve, lhs, rhs, grad, inner, dot, dx, ds, interpolate, VectorFunctionSpace, MeshFunction, MPI
from fenicsprecice import Adapter
from errorcomputation import compute_errors
from my_enums import ProblemType, Subcycling
from my_enums import ProblemType, Subcycling, DomainPart
import argparse
import numpy as np
from problem_setup import get_geometry, get_problem_setup
Expand Down Expand Up @@ -211,6 +211,10 @@ def determine_gradient(V_g, u, flux):
precice.store_checkpoint(u_n, t, n)

read_data = precice.read_data()
if problem is ProblemType.DIRICHLET and (domain_part is DomainPart.CIRCULAR or domain_part is DomainPart.RECTANGLE):
# We have to data for an arbitrary point that is not on the circle, to obtain exact solution.
# See https://github.com/precice/fenics-adapter/issues/113 for details.
read_data[(0, 0)] = u_D(0, 0)

# Update the coupling expression with the new read data
precice.update_coupling_expression(coupling_expression, read_data)
Expand Down