code: def hungarian(A, B):
      
  
      # should be kabasch here i think
      distances = cdist(A, B, 'euclidean')
  
      # Perform Hungarian analysis on distance matrix between atoms of 1st
      # structure and trial structure
      indices_a, indices_b = linear_sum_assignment(distances)
  
      return indices_b
docs: Hungarian reordering.
  
      Assume A and B are coordinates for atoms of SAME type only