diff --git a/src/support/mo_math_utilities.cpp b/src/support/mo_math_utilities.cpp
index c2b46dc2f876388b7a15c9a0c27afe982b408fbd..45430b05f343732acaec928ecc5a6d6882a2bf7a 100644
--- a/src/support/mo_math_utilities.cpp
+++ b/src/support/mo_math_utilities.cpp
@@ -8,6 +8,9 @@ void tdma_solver_vec(double *a, double *b, double *c, double *d,
                      int slev, int elev, int startidx, int endidx,
                      double* varout, int opt_acc_queue = -1) {
 
+    // Start timing
+    auto start_time = std::chrono::high_resolution_clock::now();
+
     int acc_queue = (opt_acc_queue == -1) ? 1 : opt_acc_queue; // Use 1 as the default if opt_acc_queue is not provided
 
     // Determine array sizes based on startidx and endidx