Moved 'height_adjustment' and 'stretch_factor' into 'implicit_burgers_solver().'

Snakefile

@@ -26,7 +26,8 @@ Urgent:
 TODO: Mark Burgers' equation as inviscid -> Done
 TODO: Make sure CFL number is absolute value -> Done
 TODO: Make sure only ghost cells are limited for Dirichlet boundary -> Done
-TODO: Move height adjustment and stretch factor into implicit solver function
+TODO: Move height adjustment and stretch factor into implicit solver
+    function -> Done
 TODO: Check correctness of implicit solver (with burgex.f)
 TODO: Add Burgers class completely
 TODO: Enable choice of equation

scripts/tcd/Equation.py

@@ -483,9 +483,7 @@ class Burgers(Equation):
 
         if self._init_cond.__name__ == 'Sine':
             # u(x,t) = u(x - u0*time, 0)
-            exact = self._init_cond._height_adjustment + \
-                    self._init_cond._stretch_factor * \
-                    self.implicit_burgers_solver(grid, mesh)
+            exact = self.implicit_burgers_solver(grid, mesh)
         elif self._init_cond.__name__ == 'DiscontinuousConstant':
             # u(x,t) = u(x - u0*time, 0)
             exact = self.rarefaction_wave(grid)
@@ -577,7 +575,8 @@ class Burgers(Equation):
 
         burgers_exact = uexact.reshape((1, np.size(grid_values)))
 
-        return burgers_exact
+        return self._init_cond._height_adjustment + \
+            self._init_cond._stretch_factor * burgers_exact
 
     @enforce_boundary()
     def update_right_hand_side(self, projection: ndarray) -> ndarray: