My test cases that have a bunch of atoms on a single line keep getting messed up in theobald_rmsd... where the rotation matrix is not calculated correctly. Just calling md.lprmsd works. but trying to get the superposed positions does not...