diff --git a/src/solvers/solvers.cpp b/src/solvers/solvers.cpp
index 62e267fe21bff0fffa376dd228a0477f19627f68..a17cd70b24e2afebe45c8647a77363c3d37497e2 100644
--- a/src/solvers/solvers.cpp
+++ b/src/solvers/solvers.cpp
@@ -26,6 +26,10 @@ namespace eigenpy
     ConjugateGradientVisitor< ConjugateGradient<MatrixXd,Lower|Upper> >::expose();
     LeastSquaresConjugateGradientVisitor< LeastSquaresConjugateGradient<MatrixXd, LeastSquareDiagonalPreconditionerFix<MatrixXd::Scalar> > >::expose();
     
+    // Conjugate gradient with limited BFGS preconditioner
+    ConjugateGradientVisitor< ConjugateGradient<MatrixXd,Lower|Upper,IdentityPreconditioner > >::expose("IdentityConjugateGradient");
+    ConjugateGradientVisitor< ConjugateGradient<MatrixXd,Lower|Upper,LimitedBFGSPreconditioner<double,Dynamic,Lower|Upper> > >::expose("LimitedBFGSConjugateGradient");
+    
     boost::python::enum_<Eigen::ComputationInfo>("ComputationInfo")
     .value("Success",Eigen::Success)
     .value("NumericalIssue",Eigen::NumericalIssue)