Invariant analysis of CP violation
  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. 
F. del Aguila, J. A. Aguilar-Saavedra, M. Zra{\l}ek