The power method now runs correctly.  The tensor outer product has been made more efficient and more general. Lastly, add an eigendecomposition wrapper around the power method that computes all eigenvectors and eigenvalues.