diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 2a4932b351a36007154e2039b448289c33c8bcbb..248c1d91fb85f20f4a8217fd7c5f23506bce9f50 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 4f1beca0afea4711a4581177942a0a08c99041a1..ed9402805e755c99962fc7341a410b0f2285bc0f 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