Use Eigen PermutationMatrix class instead of rolling our own - may save memory and LOC
Introduce a new regularization method based on a 1991 paper by Natarajan