Clifford synthesis breaks when a larger topology is used

[y-richie-y]

from pauliopt.clifford.tableau import CliffordTableau

ct = CliffordTableau(3)

ct.append_h(0)
ct.append_cnot(0, 1)
ct.append_s(1)

from pauliopt.topologies import Topology

topology = Topology.complete(4)

from pauliopt.clifford.tableau_synthesis import synthesize_tableau

qc, perm = synthesize_tableau(ct, topology, include_swaps=False)

    158 def x_out(self, row, col):
    159     """
    160     Get the X operator in row `row` and column `col`.
    161 
   (...)
    164         col (int): Column index.
    165     """
--> 166     return self.tableau[row, col] + 2 * self.tableau[row, col + self.n_qubits]

IndexError: index 6 is out of bounds for axis 1 with size 6