Skip to content

Commit

Permalink
Removed out-commented code and added clarification comment
Browse files Browse the repository at this point in the history
  • Loading branch information
charnley committed Nov 23, 2024
1 parent 04dc65c commit 9424f25
Showing 1 changed file with 2 additions and 8 deletions.
10 changes: 2 additions & 8 deletions rmsd/calculate_rmsd.py
Original file line number Diff line number Diff line change
Expand Up @@ -2054,24 +2054,18 @@ def main(args: Optional[List[str]] = None) -> str:
if q_reflection is not None:
q_coord_sub = np.dot(q_coord_sub, np.diag(q_reflection))

# q_coord -= q_coord

# Rotate q coordinates
# TODO Should actually follow rotation method
# q_coord = kabsch_rotate(q_coord, p_coord)
# rotation_matrix = kabsch(q_coord, p_coord)
U = kabsch(q_coord_sub, p_coord_sub)

# done and done
if settings.print_only_rmsd_atoms or not use_view:
q_coord_sub = np.dot(q_coord_sub, U)
q_coord_sub += p_cent_sub
return set_coordinates(
q_atoms_sub,
q_coord_sub,
title=f"Rotated {settings.structure_b} to match {settings.structure_a}, with RMSD of {result_rmsd:.8f}",
title=f"Rotated '{settings.structure_b}' to match '{settings.structure_a}', with a RMSD of {result_rmsd:.8f}",
)

# Swap, reflect, rotate and re-center on the full atom and coordinate set
q_coord -= q_cent_sub

if q_swap is not None:
Expand Down

0 comments on commit 9424f25

Please sign in to comment.