diff --git a/classes/sph0/louis++/CMakeLists.txt b/classes/sph0/louis++/CMakeLists.txt
new file mode 100644
index 00000000..d62ebf72
--- /dev/null
+++ b/classes/sph0/louis++/CMakeLists.txt
@@ -0,0 +1,6 @@
+CMAKE_MINIMUM_REQUIRED(VERSION 3.14)
+PROJECT(LOUIS CXX)
+
+FILE(GLOB SRCS *.h *.cpp)
+
+ADD_EXECUTABLE(louis ${SRCS})
diff --git a/classes/sph0/louis++/ParticleManager.cpp b/classes/sph0/louis++/ParticleManager.cpp
new file mode 100644
index 00000000..b61dc2fb
--- /dev/null
+++ b/classes/sph0/louis++/ParticleManager.cpp
@@ -0,0 +1,15 @@
+#include "ParticleManager.h"
+
+ParticleManager::ParticleManager()
+{
+}
+
+void
+ParticleManager::initialisation()
+{
+}
+
+void
+ParticleManager::solver()
+{
+}
diff --git a/classes/sph0/louis++/ParticleManager.h b/classes/sph0/louis++/ParticleManager.h
new file mode 100644
index 00000000..0d8c617a
--- /dev/null
+++ b/classes/sph0/louis++/ParticleManager.h
@@ -0,0 +1,13 @@
+#ifndef PARTICLEMANAGER_H
+#define PARTICLEMANAGER_H
+
+class ParticleManager
+{
+public:
+    ParticleManager();
+
+    void initialisation();
+    void solver();
+};
+
+#endif // PARTICLEMANAGER_H
diff --git a/classes/sph0/louis++/main.cpp b/classes/sph0/louis++/main.cpp
new file mode 100644
index 00000000..b570bc47
--- /dev/null
+++ b/classes/sph0/louis++/main.cpp
@@ -0,0 +1,28 @@
+// !> SPH simulation
+// @n This program is used to solve the Navier-Stokes equations
+//    using the SPH method. A number of files must be given:
+//    paths.txt, *.prm, *.fp and *.mp.
+// @warning The domain must be cubic!
+// @brief   Main program to launch a SPH simulation
+// @author  Louis Goffin
+// @date    2013-05-26
+// @version 1.0.0
+
+#include "ParticleManager.h"
+#include <iostream>
+#include <chrono>
+
+int
+main()
+{
+    std::cout << "============= SPH_simulation (L.Goffin)\n";
+
+    auto t1 = std::chrono::high_resolution_clock::now();
+
+    ParticleManager manager;
+    manager.initialisation();
+    manager.solver();
+
+    auto t2 = std::chrono::high_resolution_clock::now();
+    std::cout << "Elapsed real time = " << std::chrono::duration_cast<std::chrono::duration<double>>(t2 - t1).count() << '\n';
+}