diff --git a/src/scheme/CMakeLists.txt b/src/scheme/CMakeLists.txt
index 019ba8ec45998dd740c2e200c68e30f276e51cbd..78b4603741a462002897b01038ca8f1cefa312b4 100644
--- a/src/scheme/CMakeLists.txt
+++ b/src/scheme/CMakeLists.txt
@@ -13,6 +13,7 @@ add_library(
   HyperelasticSolver.cpp
   GKS.cpp
   GKS2.cpp
+  GKSNavier.cpp
   )
 
 target_link_libraries(
diff --git a/src/scheme/GKS2.cpp b/src/scheme/GKS2.cpp
index 54521fc802fcb385fa987a8a2ce5fa4e4f5842e6..02f61d73647ab10cc1725f898e6a102f86e14f74 100644
--- a/src/scheme/GKS2.cpp
+++ b/src/scheme/GKS2.cpp
@@ -1,5 +1,4 @@
-
-#include <scheme/GKS.hpp>
+#include <scheme/GKS2.hpp>
 
 #include <mesh/Mesh.hpp>
 #include <mesh/MeshData.hpp>
diff --git a/src/scheme/GKSNavier.cpp b/src/scheme/GKSNavier.cpp
index 92eeff2c29fddd57eb8fa0bccd06d49c7a1829e2..9b8d4a615cbfece18aba035855fcabeb503c8bdf 100644
--- a/src/scheme/GKSNavier.cpp
+++ b/src/scheme/GKSNavier.cpp
@@ -1,5 +1,4 @@
-
-#include <scheme/GKS.hpp>
+#include <scheme/GKSNavier.hpp>
 
 #include <mesh/Mesh.hpp>
 #include <mesh/MeshData.hpp>
@@ -243,12 +242,12 @@ class GKSNAVIER
 std::tuple<std::shared_ptr<const DiscreteFunctionVariant>,   // rho
            std::shared_ptr<const DiscreteFunctionVariant>,   // U
            std::shared_ptr<const DiscreteFunctionVariant>>   // E
-gks2(std::shared_ptr<const DiscreteFunctionVariant> rho_v,
-     std::shared_ptr<const DiscreteFunctionVariant> rho_U_v,
-     std::shared_ptr<const DiscreteFunctionVariant> rho_E_v,
-     std::shared_ptr<const DiscreteFunctionVariant> tau,
-     const double delta,
-     const double dt)
+gksNavier(std::shared_ptr<const DiscreteFunctionVariant> rho_v,
+          std::shared_ptr<const DiscreteFunctionVariant> rho_U_v,
+          std::shared_ptr<const DiscreteFunctionVariant> rho_E_v,
+          std::shared_ptr<const DiscreteFunctionVariant> tau,
+          const double delta,
+          const double dt)
 {
   std::shared_ptr mesh_v = getCommonMesh({rho_v, rho_U_v, rho_E_v});
   if (not mesh_v) {