code: def compare(self, other, r_threshold=1e-3):
          
          return compute_rmsd(self.r, other.r) < r_threshold
docs: Compare two rotations
  
             The RMSD of the rotation matrices is computed. The return value
             is True when the RMSD is below the threshold, i.e. when the two
             rotations are almost identical.