  We show that the problem of finding the primary and secondary characteristic
directions of a linear lossless optical element can be reformulated in terms of
an eigenvalue problem related to the unimodular factor of the transfer matrix
of the optical device. This formulation makes any actual computation of the
characteristic directions amenable to pre-implemented numerical routines,
thereby facilitating the decomposition of the transfer matrix into equivalent
linear retarders and rotators according to the related Poincare equivalence
theorem. The method is expected to be useful whenever the inverse problem of
reconstruction of the internal state of a transparent medium from optical data
obtained by tomographical methods is an issue.
