🐞 Fix diatomic transform

[Alexsp32]

• Fix flipped indices and other mistakes in internal coordinate transform matrix that was causing transform_to_internal_coordinates to give nonsensical numbers.
• Add a unit test which ensures forward transform ⋅ back transform retain the original matrix.
• Add unit test to ensure transform matrix is unitary