using PyCall
@pyimport pyphi
@pyimport pyphi.compute as pc

tpm = [
  0 0 0
  0 0 1
  1 0 1
  1 0 0
  1 1 0
  1 1 1
  1 1 1
  1 1 0
]

cm = [
  0 0 1
  1 0 1
  1 1 0
]

labels = ("A", "B", "C")

network = pyphi.Network(tpm, cm=cm, node_labels=labels)

state = (1, 0, 0)
node_indices = (0, 1, 2)
subsystem = pyphi.Subsystem(network, state, node_indices)

pc.phi(subsystem)
pyphi.compute[:phi](subsystem)

pc.sia(subsystem)