Optimization of RMSD calculation

Two types of optimization implemented:

1. Tensor is used instead of Matrix. This allows avoiding the extra
allocations required in Matrix to setup the workspace, since the
needed space is now on the stack rather than on the heap

2. When align==displace, only lower eigenvalue is computed.

Overall, these changes result in a 5-10% speedup of rt46 on my mac