diff --git a/bindings/python/crocoddyl/core/solvers/box-qp.cpp b/bindings/python/crocoddyl/core/solvers/box-qp.cpp
index 024e9500f4..1dfdb47818 100644
--- a/bindings/python/crocoddyl/core/solvers/box-qp.cpp
+++ b/bindings/python/crocoddyl/core/solvers/box-qp.cpp
@@ -67,7 +67,7 @@ void exposeSolverBoxQP() {
           ":param th_grad: gradient tolerance condition (default 1e-9)\n"
           ":param reg: regularization (default 1e-9)"))
       .def("solve", &BoxQP::solve,
-           bp::return_value_policy<bp::return_by_value>(),
+           bp::return_value_policy<bp::copy_const_reference>(),
            bp::args("H", "q", "lb", "ub", "xinit"),
            "Compute the solution of bound-constrained QP based on Newton "
            "projection.\n\n"