diff --git a/GeoLib/AnalyticalGeometry.h b/GeoLib/AnalyticalGeometry.h
index 34ce39f5e88b485311cc581d78896341cad18d34..3592d5c2f6d074c38fba7a502cf0663524bc99ab 100644
--- a/GeoLib/AnalyticalGeometry.h
+++ b/GeoLib/AnalyticalGeometry.h
@@ -16,9 +16,6 @@
 
 #include <memory>
 
-#include "MathLib/LinAlg/Dense/DenseMatrix.h"
-#include "MathLib/Vector3.h"
-
 #include "Polygon.h"
 
 namespace GeoLib
diff --git a/MeshGeoToolsLib/GeoMapper.cpp b/MeshGeoToolsLib/GeoMapper.cpp
index 97dcd824304a631ac2d14a93c5e264f799c49b63..357365dbf8ade56ae81c84fd655d71346f7d848b 100644
--- a/MeshGeoToolsLib/GeoMapper.cpp
+++ b/MeshGeoToolsLib/GeoMapper.cpp
@@ -89,7 +89,6 @@ void GeoMapper::mapOnMesh(MeshLib::Mesh const*const mesh)
 
     // init grid
     MathLib::Point3d origin(std::array<double,3>{{0,0,0}});
-    MathLib::Vector3 normal(0,0,-1);
     std::vector<MeshLib::Node> flat_nodes;
     flat_nodes.reserve(_surface_mesh->getNumberOfNodes());
     // copy nodes and project the copied nodes to the x-y-plane, i.e. set
diff --git a/MeshLib/ElementCoordinatesMappingLocal.cpp b/MeshLib/ElementCoordinatesMappingLocal.cpp
index 9b3500d9b1236dadaa87c1a03549c6fd60098f77..eb9fcb1e4afa44e738a7dc80ed6c56d3a38f1ff3 100644
--- a/MeshLib/ElementCoordinatesMappingLocal.cpp
+++ b/MeshLib/ElementCoordinatesMappingLocal.cpp
@@ -18,7 +18,6 @@
 #include "MeshLib/Node.h"
 #include "MathLib/MathTools.h"
 #include "MathLib/Point3d.h"
-#include "MathLib/Vector3.h"
 
 namespace detail
 {
diff --git a/MeshLib/Elements/FaceRule.h b/MeshLib/Elements/FaceRule.h
index 70d6e23e670eb13a24a5d304a4e8c9bd1a3b2112..2f4f47de058ed980019126ea6f9a84b4d7e9d43b 100644
--- a/MeshLib/Elements/FaceRule.h
+++ b/MeshLib/Elements/FaceRule.h
@@ -10,7 +10,6 @@
 
 #pragma once
 
-#include "MathLib/Vector3.h"
 #include "Element.h"
 
 namespace MeshLib
diff --git a/MeshLib/Node.h b/MeshLib/Node.h
index 4a3f465321c5cbe4f690f47f497a897c58e826d3..1628c97e98f200b481ce752f09b20c2a8941c97c 100644
--- a/MeshLib/Node.h
+++ b/MeshLib/Node.h
@@ -19,7 +19,6 @@
 #include <vector>
 
 #include "MathLib/Point3dWithID.h"
-#include "MathLib/Vector3.h"
 
 namespace ApplicationUtils
 {
diff --git a/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h b/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
index 47f5b685519297130c4ac7f8cd64a7959e9598ad..facd6e966b9b7dbb0d4260884cd9e3e9945e5718 100644
--- a/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
+++ b/ProcessLib/BoundaryCondition/NormalTractionBoundaryConditionLocalAssembler.h
@@ -11,6 +11,7 @@
 #pragma once
 
 #include "MathLib/LinAlg/Eigen/EigenMapTools.h"
+#include "MathLib/Vector3.h"
 #include "MeshLib/Elements/FaceRule.h"
 #include "NumLib/DOF/DOFTableUtil.h"
 #include "ParameterLib/Parameter.h"
diff --git a/ProcessLib/LIE/Common/Utils.cpp b/ProcessLib/LIE/Common/Utils.cpp
index f18601f1c78b2373deb5b08ebae2501652ca6223..1df5adce7c2a75aed6eb3049f3b690ffa5c2c328 100644
--- a/ProcessLib/LIE/Common/Utils.cpp
+++ b/ProcessLib/LIE/Common/Utils.cpp
@@ -9,6 +9,7 @@
 
 #include "Utils.h"
 
+#include "MathLib/Vector3.h"
 #include "MeshLib/Elements/FaceRule.h"
 
 namespace ProcessLib