diff --git a/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h b/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h
index 486665226de2b6d97d5f2af407ba68da1bccbb16..6853b0c4b757c162d3fffd2984b0bbf681ec6c97 100644
--- a/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h
+++ b/NumLib/Fem/CoordinatesMapping/ShapeMatrices.h
@@ -14,6 +14,8 @@
 
 #include <ostream>
 
+#include <Eigen/Core>
+
 namespace NumLib
 {
 
@@ -103,6 +105,7 @@ struct ShapeMatrices
      */
     void write (std::ostream& out) const;
 
+    EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
 }; // ShapeMatrices
 
 
diff --git a/Tests/NumLib/SteadyDiffusion2DExample1.h b/Tests/NumLib/SteadyDiffusion2DExample1.h
index 909e515e9f64e681d9f5682fb5a423576e053d4b..f82b59a03837844c935728fb0123442be20ca98c 100644
--- a/Tests/NumLib/SteadyDiffusion2DExample1.h
+++ b/Tests/NumLib/SteadyDiffusion2DExample1.h
@@ -147,4 +147,6 @@ template<typename IndexType>struct SteadyDiffusion2DExample1
 
     LocalMatrixType _localA;
     LocalVectorType _localRhs;
+
+    EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
 };