code: def determine_signs(nodes, edges, cutoff=1e-10):
      
      N = len(nodes)
      M = np.zeros((N, N), dtype=float)
      for a, b, w in edges:
          M[a, b] += w
      M = symmetrize(M)
  
      return get_signs(M, cutoff=cutoff, validate=False)
docs: Construct the orientation matrix for the pairs on N molecules.
  
      >>> determine_signs([0, 1, 2], [(0, 1, 1), (0, 2, -1), (1, 2, -1)])
      array([ 1,  1, -1])