  We describe an algorithm for computing the complete set of primitive
orthogonal idempotents in the centralizer ring of the permutation
representation of a wreath product. This set of idempotents determines the
decomposition of the representation into irreducible components. In the
formalism of quantum mechanics, these idempotents are projection operators into
irreducible invariant subspaces of the Hilbert space of a multipartite quantum
system. The C implementation of the algorithm constructs irreducible
decompositions of high-dimensional representations of wreath products. Examples
of computations are given.
