Determination of the characteristic directions of lossless linear   optical elements
  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. 
Hanno Hammer