Now using moore-penrose pseudo-inverse.  I'm about to make a new branch to try to change the implementation of the Cholesky decomposition so I'm not using logged diagonals anymore.