diff --git a/src/scheme/GKSNavier.cpp b/src/scheme/GKSNavier.cpp
index 48b98b66b1375148e80b011389d195d0bb1c5777..a24cb67652d5a3898f0ff41236d9f4637d89b4c7 100644
--- a/src/scheme/GKSNavier.cpp
+++ b/src/scheme/GKSNavier.cpp
@@ -114,21 +114,6 @@ class GKSHandler::GKSNAVIER final : public GKSHandler::IGKSNAVIER
     MatrixP(2, 0) = 0.5 * U_2 - 0.25 * (delta + 1) / lambda;
     MatrixP(2, 1) = -U[0];
 
-    // TinyMatrix<SIZEproblem> MatrixM(zero);
-    // TinyMatrix<SIZEproblem> invMatrixP(zero);
-    // TinyMatrix<SIZEproblem> MatrixMtilde(zero);
-
-    // MatrixMtilde(0, 0) = 1;
-    // MatrixMtilde(1, 1) = 1;
-    // MatrixMtilde(2, 2) = (rho * (delta + 1)) / (8 * lambda * lambda);
-
-    // invMatrixP(0, 0) = sqrt_rho;
-    // invMatrixP(1, 1) = sqrt_rho / (sqrt2 * sqrt_lambda);
-    // invMatrixP(2, 2) = 1;
-    // invMatrixP(1, 0) = sqrt_rho * U[0];
-    // invMatrixP(2, 0) = sqrt_rho * (0.5 * U_2 + 0.25 * (delta + 1) / lambda);
-    // invMatrixP(2, 1) = U[0] * sqrt_rho / (sqrt2 * sqrt_lambda);
-
     invMatrixM = transpose(MatrixP) * invMatrixMtilde * (MatrixP);
 
     return invMatrixM;
@@ -528,6 +513,7 @@ class GKSHandler::GKSNAVIER final : public GKSHandler::IGKSNAVIER
         TinyMatrix<SIZEproblem> InvMInterface =
           this->_computeInvMatrixM(rho_node[node_id], U_node[node_id], lambda_node[node_id], delta);
 
+        // //For the calculation of matrices M using the inverse() routine of pugs for TinyMatrix
         // TinyMatrix<SIZEproblem> InvM_left =
         //   this->_computeInvM_bis(rho_left_node, rhoU_left_node, rhoE_left_node, U_left_node, lambda_left_node,
         //   delta);
@@ -547,19 +533,6 @@ class GKSHandler::GKSNAVIER final : public GKSHandler::IGKSNAVIER
         ar = this->_compute_a_and_b(InvM_right, Gradj_for_a[node_cells[0]]);
         br = this->_compute_a_and_b(InvMInterface, Gradj_for_b[node_cells[0]]);
 
-        // for (size_t l = 0; l < node_cells.size(); l++) {
-        //   if (node_cells.size() == 1) {
-        //     continue;
-        //   }
-        //   if (l == 0) {
-        //     al = this->_compute_a_and_b(InvM_left, Gradj_for_a[node_cells[l]]);
-        //     bl = this->_compute_a_and_b(InvMInterface, Gradj_for_b[node_cells[l]]);
-        //   } else {
-        //     ar = this->_compute_a_and_b(InvM_right, -Gradj_for_a[node_cells[l]]);
-        //     br = this->_compute_a_and_b(InvMInterface, -Gradj_for_b[node_cells[l]]);
-        //   }
-        // }
-
         auto [C0_semi_pos_left, C1_semi_pos_left, C2_semi_pos_left] =
           this->_computeMatricesC0C1C2_semi_moments(rho_left_node, U_left_node, lambda_left_node, u_moments_left,
                                                     xi2_moments_left, 1);