From 31ab2400d63a4999e3b282e19e60eb40db5122fd Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Mon, 27 Aug 2012 09:59:17 +0200
Subject: [PATCH] Revert "Make dirichlet conditions independent of the mesh"
This reverts commit 3d848c50473d40ad89f789b6cb393e33b3a81b6b.
Conflicts:
src/one-body-sample.cc
src/one-body-sample.py
---
src/one-body-sample.cc | 5 +----
src/one-body-sample.py | 2 +-
2 files changed, 2 insertions(+), 5 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 2a4932b3..248c1d91 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -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 <
diff --git a/src/one-body-sample.py b/src/one-body-sample.py
index 4f1beca0..ed940280 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -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
--
GitLab