diff --git a/src/solvers/BasicPreconditioners.hpp b/src/solvers/BasicPreconditioners.hpp
index b9a085ef5628fc2273ae96353f23d90a53edd073..00f4199dacff417796b8032a2eda8829eb8deac5 100644
--- a/src/solvers/BasicPreconditioners.hpp
+++ b/src/solvers/BasicPreconditioners.hpp
@@ -47,6 +47,9 @@ namespace eigenpy
       .def("compute",&Preconditioner::template compute<MatrixType>,bp::arg("mat"),
            "Initialize the preconditioner from the matrix value.",
            bp::return_value_policy<bp::reference_existing_object>())
+      .def("factorize",&Preconditioner::template factorize<MatrixType>,bp::arg("mat"),
+           "Initialize the preconditioner from the matrix value, i.e factorize the mat given as input to approximate its inverse.",
+           bp::return_value_policy<bp::reference_existing_object>())
       ;
       
     }