feellgood · edgar-bonet · May 5, 2026 · May 5, 2026
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -109,50 +109,7 @@ if(${ENABLE_LTO})
     set(CMAKE_INTERPROCEDURAL_OPTIMIZATION TRUE)
 endif()
 
-set(CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} -fopenmp")
-
-SET(HEADERS config.h node.h expression_parser.h mesh.h electrostatSolver.h
-    spinAccumulationSolver.h time_integration.h settings.h tetra.h
-    facette.h linear_algebra.h log-stats.h tags.h chronometer.h element.h)
-
-SET(SOURCES settings.cpp time_integration.cpp solver.cpp electrostatSolver.cpp
-    spinAccumulationSolver.cpp read.cpp save.cpp linear_algebra.cpp recentering.cpp tetra.cpp
-    energy.cpp facette.cpp expression_parser.cpp chronometer.cpp
-    tags.cpp mesh.cpp)
-
-configure_file(config.h.in ./config.h)
-
-add_custom_command(
-    OUTPUT default-settings.o
-    COMMAND ld -r -b binary -z noexecstack default-settings.yml -o default-settings.o
-    DEPENDS default-settings.yml)
-
-#===========================================================================
-# Compiler options
-#===========================================================================
-# -ffast-math -msse4.2 -mfpmath=sse (on X86 -mfpmath=sse set by default) or -mavx
-# -fsanitize=address leak null undefined # both compiler and linker
-#===========================================================================
-
-include_directories (${CMAKE_SOURCE_DIR} ${GMSH_INC})
-
-ADD_EXECUTABLE(feellgood main.cpp ${HEADERS} ${SOURCES} default-settings.o)
-target_compile_options( feellgood PUBLIC -Wall -Wextra -O3 -march=native ${MKL_CFLAGS} )
-target_compile_definitions(feellgood PUBLIC
-        -DNDEBUG -DTBB_SUPPRESS_DEPRECATED_MESSAGES ${EIGEN_MKL_CFLAGS} )
-
-if (MKL_FOUND)
-    #option --no-as-needed for mkl, see
-    # https://www.intel.com/content/www/us/en/developer/tools/oneapi/onemkl-link-line-advisor.html
-    target_link_options(feellgood PUBLIC "-Wl,--no-as-needed" )
-    # not sure we really need it since --no-as-needed is default value, though there might be
-    # a --as-needed before
-    TARGET_LINK_LIBRARIES(feellgood yaml-cpp duktape ann OpenMP::OpenMP_CXX ${TBB_LINK}
-            Eigen3::Eigen ${MKL_TARGET} ${GMSH_LIB})
-else()
-    TARGET_LINK_LIBRARIES(feellgood yaml-cpp duktape ann OpenMP::OpenMP_CXX ${TBB_LINK}
-            Eigen3::Eigen ${GMSH_LIB})
-endif()
+add_subdirectory(src)
 
 #===========================================================================
 # Find where to install the Python modules.
@@ -226,7 +183,7 @@ else()
 endif()
 
 if(${ENABLE_UTESTS})
-    configure_file(./unit-tests/ut_config.h.in ./unit-tests/ut_config.h)
+    configure_file(unit-tests/ut_config.h.in unit-tests/ut_config.h)
     include(CTest)
     add_subdirectory(unit-tests)
 endif()

diff --git a/Doxyfile b/Doxyfile
@@ -765,7 +765,7 @@ WARN_LOGFILE           =
 # spaces. See also FILE_PATTERNS and EXTENSION_MAPPING
 # Note: If this tag is empty the current directory is searched.
 
-INPUT                  = . algebra
+INPUT                  = . src src/algebra
 
 # This tag can be used to specify the character encoding of the source files
 # that doxygen parses. Internally doxygen uses the UTF-8 encoding. Doxygen uses

diff --git a/ci-tests/full_test.py b/ci-tests/full_test.py
@@ -43,15 +43,15 @@ def distance_R3(v1, v2):
 }
 
 # Expected outcome. It depends on whether feeLLGood was compiled in ONE_GAUSS_POINT mode.
-val = subprocess.check_output(["../feellgood","--version"])
+val = subprocess.check_output(["../src/feellgood","--version"])
 if b'ONE_GAUSS_POINT=ON' in val:
     expected = { "m": [0.307432, 0.476202, -0.823843], "E_tot": -4.446980e-19 }
 else:
     expected = { "m": [0.307542, 0.475877, -0.823989], "E_tot": -4.445605e-19 }
 
 # Run the simulation.
 sys.stdout.flush()
-val = subprocess.run(["../feellgood", "--seed", "2", "-"], input=json.dumps(settings), text=True)
+val = subprocess.run(["../src/feellgood", "--seed", "2", "-"], input=json.dumps(settings), text=True)
 
 #(devNote) to avoid -fsanitize=leak ASLR bug, we turn off address randomizer
 #val = subprocess.run(["setarch", "--addr-no-randomize", "../feellgood", "--seed", "2", "-"],

diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
@@ -0,0 +1,44 @@
+SET(HEADERS config.h node.h expression_parser.h mesh.h electrostatSolver.h
+        spinAccumulationSolver.h time_integration.h settings.h tetra.h
+        facette.h linear_algebra.h log-stats.h tags.h chronometer.h element.h)
+
+SET(SOURCES settings.cpp time_integration.cpp solver.cpp electrostatSolver.cpp
+        spinAccumulationSolver.cpp read.cpp save.cpp linear_algebra.cpp recentering.cpp tetra.cpp
+        energy.cpp facette.cpp expression_parser.cpp chronometer.cpp
+        tags.cpp mesh.cpp)
+
+configure_file(config.h.in ./config.h)
+
+set(CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} -fopenmp")
+
+add_custom_command(
+    OUTPUT default-settings.o
+    COMMAND ld -r -b binary -z noexecstack default-settings.yml -o default-settings.o
+    DEPENDS default-settings.yml)
+
+#===========================================================================
+# Compiler options
+#===========================================================================
+# -ffast-math -msse4.2 -mfpmath=sse (on X86 -mfpmath=sse set by default) or -mavx
+# -fsanitize=address leak null undefined # both compiler and linker
+#===========================================================================
+
+include_directories (${CMAKE_SOURCE_DIR} ${GMSH_INC})
+
+ADD_EXECUTABLE(feellgood main.cpp ${HEADERS} ${SOURCES} default-settings.o)
+target_compile_options( feellgood PUBLIC -Wall -Wextra -O3 -march=native ${MKL_CFLAGS} )
+target_compile_definitions(feellgood PUBLIC
+        -DNDEBUG -DTBB_SUPPRESS_DEPRECATED_MESSAGES ${EIGEN_MKL_CFLAGS} )
+
+if (MKL_FOUND)
+    #option --no-as-needed for mkl, see
+    # https://www.intel.com/content/www/us/en/developer/tools/oneapi/onemkl-link-line-advisor.html
+    target_link_options(feellgood PUBLIC "-Wl,--no-as-needed" )
+    # not sure we really need it since --no-as-needed is default value, though there might be
+    # a --as-needed before
+    TARGET_LINK_LIBRARIES(feellgood yaml-cpp duktape ann OpenMP::OpenMP_CXX ${TBB_LINK}
+            Eigen3::Eigen ${MKL_TARGET} ${GMSH_LIB})
+else()
+    TARGET_LINK_LIBRARIES(feellgood yaml-cpp duktape ann OpenMP::OpenMP_CXX ${TBB_LINK}
+            Eigen3::Eigen ${GMSH_LIB})
+endif()
diff --git a/algebra/algebra.h → src/algebra/algebra.h b/algebra/algebra.h → src/algebra/algebra.h
diff --git a/algebra/algebraCore.h → src/algebra/algebraCore.h b/algebra/algebraCore.h → src/algebra/algebraCore.h
diff --git a/algebra/bicg.h → src/algebra/bicg.h b/algebra/bicg.h → src/algebra/bicg.h
diff --git a/algebra/cg.h → src/algebra/cg.h b/algebra/cg.h → src/algebra/cg.h
diff --git a/algebra/iter.h → src/algebra/iter.h b/algebra/iter.h → src/algebra/iter.h
diff --git a/algebra/sparseMat.h → src/algebra/sparseMat.h b/algebra/sparseMat.h → src/algebra/sparseMat.h
@@ -25,7 +25,7 @@ reconstructing the matrix from scratch.
 #include <execution>
 #include <mutex>
 
-#include "config.h"
+#include "../config.h"
 #include "algebraCore.h"
 
 namespace algebra

diff --git a/chronometer.cpp → src/chronometer.cpp b/chronometer.cpp → src/chronometer.cpp
diff --git a/chronometer.h → src/chronometer.h b/chronometer.h → src/chronometer.h
diff --git a/config.h.in → src/config.h.in b/config.h.in → src/config.h.in
diff --git a/default-settings.yml → src/default-settings.yml b/default-settings.yml → src/default-settings.yml
diff --git a/electrostatSolver.cpp → src/electrostatSolver.cpp b/electrostatSolver.cpp → src/electrostatSolver.cpp
diff --git a/electrostatSolver.h → src/electrostatSolver.h b/electrostatSolver.h → src/electrostatSolver.h
diff --git a/element.h → src/element.h b/element.h → src/element.h
diff --git a/energy.cpp → src/energy.cpp b/energy.cpp → src/energy.cpp
diff --git a/expression_parser.cpp → src/expression_parser.cpp b/expression_parser.cpp → src/expression_parser.cpp
diff --git a/expression_parser.h → src/expression_parser.h b/expression_parser.h → src/expression_parser.h
diff --git a/facette.cpp → src/facette.cpp b/facette.cpp → src/facette.cpp
diff --git a/facette.h → src/facette.h b/facette.h → src/facette.h
diff --git a/fem.h → src/fem.h b/fem.h → src/fem.h
diff --git a/fmm_demag.h → src/fmm_demag.h b/fmm_demag.h → src/fmm_demag.h
diff --git a/linear_algebra.cpp → src/linear_algebra.cpp b/linear_algebra.cpp → src/linear_algebra.cpp
diff --git a/linear_algebra.h → src/linear_algebra.h b/linear_algebra.h → src/linear_algebra.h
diff --git a/log-stats.h → src/log-stats.h b/log-stats.h → src/log-stats.h
diff --git a/main.cpp → src/main.cpp b/main.cpp → src/main.cpp
diff --git a/mesh.cpp → src/mesh.cpp b/mesh.cpp → src/mesh.cpp
diff --git a/mesh.h → src/mesh.h b/mesh.h → src/mesh.h
diff --git a/meshUtils.h → src/meshUtils.h b/meshUtils.h → src/meshUtils.h
diff --git a/node.h → src/node.h b/node.h → src/node.h
diff --git a/read.cpp → src/read.cpp b/read.cpp → src/read.cpp
diff --git a/recentering.cpp → src/recentering.cpp b/recentering.cpp → src/recentering.cpp
diff --git a/save.cpp → src/save.cpp b/save.cpp → src/save.cpp
diff --git a/settings.cpp → src/settings.cpp b/settings.cpp → src/settings.cpp
diff --git a/settings.h → src/settings.h b/settings.h → src/settings.h
diff --git a/solver.cpp → src/solver.cpp b/solver.cpp → src/solver.cpp
diff --git a/solver.h → src/solver.h b/solver.h → src/solver.h
diff --git a/spinAccumulationSolver.cpp → src/spinAccumulationSolver.cpp b/spinAccumulationSolver.cpp → src/spinAccumulationSolver.cpp
diff --git a/spinAccumulationSolver.h → src/spinAccumulationSolver.h b/spinAccumulationSolver.h → src/spinAccumulationSolver.h
diff --git a/tags.cpp → src/tags.cpp b/tags.cpp → src/tags.cpp
diff --git a/tags.h → src/tags.h b/tags.h → src/tags.h
diff --git a/tetra.cpp → src/tetra.cpp b/tetra.cpp → src/tetra.cpp
diff --git a/tetra.h → src/tetra.h b/tetra.h → src/tetra.h
diff --git a/time_integration.cpp → src/time_integration.cpp b/time_integration.cpp → src/time_integration.cpp
diff --git a/time_integration.h → src/time_integration.h b/time_integration.h → src/time_integration.h
diff --git a/unit-tests/CMakeLists.txt b/unit-tests/CMakeLists.txt
@@ -16,7 +16,7 @@ if(NOT GMSH_INC)
   message(FATAL_ERROR "Could not find gmsh.h")
 endif()
 
-include_directories (${CMAKE_SOURCE_DIR} ${Boost_INCLUDE_DIRS} ${GMSH_INC})
+include_directories (${CMAKE_SOURCE_DIR}/src ${Boost_INCLUDE_DIRS} ${GMSH_INC})
 
 set( BOOST_TEST_LOG_LEVEL "all" )
 
@@ -32,39 +32,39 @@ add_executable (test_ut_algebra ut_algebra.cpp)
 
 add_executable (test_ut_OMP_solver ut_OMP_solver.cpp)
 
-add_executable(test_ut_dukTapeParser ../expression_parser.cpp ut_dukTapeParser.cpp)
+add_executable(test_ut_dukTapeParser ../src/expression_parser.cpp ut_dukTapeParser.cpp)
 
 add_executable (test_ut_node ut_node.cpp)
 
 add_executable (test_ut_index ut_index.cpp)
 
-SET(SOURCES ../facette.cpp ../tetra.cpp ut_element.cpp)
+SET(SOURCES ../src/facette.cpp ../src/tetra.cpp ut_element.cpp)
 add_executable (test_ut_element ${SOURCES})
 
-SET(SOURCES ../facette.cpp ut_facette.cpp)
+SET(SOURCES ../src/facette.cpp ut_facette.cpp)
 add_executable (test_ut_facette ${SOURCES})
 
-SET(SOURCES ../facette.cpp ut_fac_charges.cpp)
+SET(SOURCES ../src/facette.cpp ut_fac_charges.cpp)
 add_executable (test_ut_fac_charges ${SOURCES})
 
-SET(SOURCES ../tetra.cpp ut_tetra.cpp)
+SET(SOURCES ../src/tetra.cpp ut_tetra.cpp)
 add_executable (test_ut_tetra ${SOURCES})
 
-SET(SOURCES ../tetra.cpp ut_tet_lumping.cpp)
+SET(SOURCES ../src/tetra.cpp ut_tet_lumping.cpp)
 add_executable (test_ut_tet_lumping ${SOURCES})
 
-SET(SOURCES ../tetra.cpp ut_tet_charges.cpp)
+SET(SOURCES ../src/tetra.cpp ut_tet_charges.cpp)
 add_executable (test_ut_tet_charges ${SOURCES})
 
-SET(SOURCES ../tetra.cpp ut_energy.cpp ../facette.cpp)
+SET(SOURCES ../src/tetra.cpp ut_energy.cpp ../src/facette.cpp)
 add_executable (test_ut_energy ${SOURCES})
 
-SET(SOURCES ../tetra.cpp ut_anisotropy.cpp)
+SET(SOURCES ../src/tetra.cpp ut_anisotropy.cpp)
 add_executable (test_ut_anisotropy ${SOURCES})
 
 add_executable (test_ut_tiny ut_tiny.cpp)
 
-set(SOURCES ../tetra.cpp ut_time_int.cpp)
+set(SOURCES ../src/tetra.cpp ut_time_int.cpp)
 add_executable (test_ut_time_int ${SOURCES})
 
 add_executable (test_ut_log-stats ut_log-stats.cpp)

diff --git a/unit-tests/sparse_matrix.h b/unit-tests/sparse_matrix.h
@@ -7,7 +7,7 @@
 
 #include <vector>
 
-#include "../algebra/sparseMat.h"
+#include "algebra/sparseMat.h"
 
 using algebra::SparseMatrix;
 using algebra::MatrixShape;

diff --git a/unit-tests/ut_algebra.cpp b/unit-tests/ut_algebra.cpp
@@ -7,9 +7,9 @@
 
 #include "ut_config.h"  // for tolerance UT_TOL macro
 #include "sparse_matrix.h"
-#include "../algebra/algebra.h"
-#include "../algebra/cg.h"
-#include "../algebra/bicg.h"
+#include "algebra/algebra.h"
+#include "algebra/cg.h"
+#include "algebra/bicg.h"
 
 using namespace algebra;