diff --git a/include/eigenpy/decompositions/LDLT.hpp b/include/eigenpy/decompositions/LDLT.hpp
index 347405411f7580a8f99949abb3b286f643f97dd6..2a635a217b028e3d6c10eac97bac814b26d6469a 100644
--- a/include/eigenpy/decompositions/LDLT.hpp
+++ b/include/eigenpy/decompositions/LDLT.hpp
@@ -10,6 +10,7 @@
 
 #include "eigenpy/eigenpy.hpp"
 #include "eigenpy/utils/scalar-name.hpp"
+#include "eigenpy/eigen/EigenBase.hpp"
 
 namespace eigenpy {
 
@@ -36,6 +37,8 @@ struct LDLTSolverVisitor
             bp::args("self", "matrix"),
             "Constructs a LDLT factorization from a given matrix."))
 
+        .def(EigenBaseVisitor<Solver>())
+
         .def("isNegative", &Solver::isNegative, bp::arg("self"),
              "Returns true if the matrix is negative (semidefinite).")
         .def("isPositive", &Solver::isPositive, bp::arg("self"),
diff --git a/include/eigenpy/decompositions/LLT.hpp b/include/eigenpy/decompositions/LLT.hpp
index b0a1e1b57d0165d5de41fb8d97f38a2e3973a5ab..1996d6ccbe3657a44f6a258b77b9adf75eebe08a 100644
--- a/include/eigenpy/decompositions/LLT.hpp
+++ b/include/eigenpy/decompositions/LLT.hpp
@@ -10,6 +10,7 @@
 
 #include "eigenpy/eigenpy.hpp"
 #include "eigenpy/utils/scalar-name.hpp"
+#include "eigenpy/eigen/EigenBase.hpp"
 
 namespace eigenpy {
 
@@ -36,6 +37,8 @@ struct LLTSolverVisitor
             bp::args("self", "matrix"),
             "Constructs a LLT factorization from a given matrix."))
 
+        .def(EigenBaseVisitor<Solver>())
+
         .def("matrixL", &matrixL, bp::arg("self"),
              "Returns the lower triangular matrix L.")
         .def("matrixU", &matrixU, bp::arg("self"),