  In this paper, we discuss the geometric integration of hamiltonian systems on
Poisson manifolds, in particular, in the case, when the Poisson structure is
induced by a Lie algebra, that is, it is a Lie-Poisson structure.
  A Hamiltonian system on a Poisson manifold $(P, \Pi)$ is a smooth manifold
$P$ equipped with a bivector field $\Pi$ satisfying $[\Pi, \Pi]=0$ (Jacobi
identity), inducing the Poisson bracket on $C^{\infty}(P)$, $\{f, g\}\equiv
\Pi(df, dg)$ where $f, g\in C^{\infty}(P)$. For any $f\in C^{\infty}(P)$ the
Hamiltonian vector field is defined by $X_f(g)=\{g, f\}$. The Hamiltonian
vector fields $X_f$ generate an integrable generalized distribution on $P$ and
the leaves of this foliation are symplectic. The flow of any hamiltonian vector
field preserves the Poisson structure, it fixes each leaf and the hamiltonian
itself is a first integral.
  It is important to characterize numerical methods preserving some of these
fundamental properties of the hamiltonian flow on Poisson manifolds (geometric
integrators). We discuss the difficulties of deriving these Poisson methods
using standard techniques and we present some modern approaches to the problem.
