  The invariant formulation of CP violation involves the generation of sets of
invariant constraints for CP conservation, the manipulation of their
expressions and the identification of complete and minimal subsets of such
constraints. In this paper we present a collection of subroutines to deal with
these three tasks in a fast, reliable and systematic way, with examples for the
leptonic sector.
