Skip to content
Snippets Groups Projects
Commit 31ab2400 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Revert "Make dirichlet conditions independent of the mesh"

This reverts commit 3d848c50473d40ad89f789b6cb393e33b3a81b6b.

Conflicts:

	src/one-body-sample.cc
	src/one-body-sample.py
parent 5dbb98a3
No related branches found
No related tags found
No related merge requests found
......@@ -135,7 +135,6 @@ int main(int argc, char *argv[]) {
auto const solver_tolerance = parset.get<double>("solver.tolerance");
auto const refinements = parset.get<size_t>("grid.refinements");
size_t const refinement_factor = 1 << refinements;
auto const verbose = parset.get<bool>("verbose");
Solver::VerbosityMode const verbosity =
......@@ -255,10 +254,8 @@ int main(int argc, char *argv[]) {
stiffnessMatrix.mmv(u, rhs);
// Apply Dirichlet condition
for (size_t i = 0; i < finestSize; ++i)
if (ignoreNodes[i].count() == dim) {
if (ignoreNodes[i].count() == dim)
dirichletFunction.evaluate(time, u_diff[i][0]);
u_diff[i][0] /= refinement_factor;
}
for (size_t state_fpi = 0;
state_fpi <
......
......@@ -13,7 +13,7 @@ class neumannCondition:
class dirichletCondition:
def __call__(self, x):
return 3e-3
return 3e-3 / 16
# return 0
fst = 3e-4
snd = 5e-4
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment