diff --git a/NumLib/ODESolver/TimeDiscretizedODESystem.cpp b/NumLib/ODESolver/TimeDiscretizedODESystem.cpp
index 87f4be36f09db740a18932ddf4ff0c1cc2bc58f0..f38fa038e19e8995d7e4925129c83ec367855e0c 100644
--- a/NumLib/ODESolver/TimeDiscretizedODESystem.cpp
+++ b/NumLib/ODESolver/TimeDiscretizedODESystem.cpp
@@ -32,6 +32,7 @@ void applyKnownSolutions(std::vector<Solutions> const* const known_solutions,
             MathLib::setVector(x, bc.ids[i], bc.values[i]);
         }
     }
+    MathLib::LinAlg::finalizeAssembly(x);
 }
 }