Stavros Efthymiou

Thanks for opening this issue. Could you please describe in more detail which are the control and target qubits of the CZ gates? In qibo we use the register terminology...

> I think that this method should be easy to implement, and the main problem I see here is computing the conmutator between H_B and H_P, step II subfigure B....

> How is the inner structure of a `TrotterHamiltonian`? I would say that we have a list of many pauli terms and we could do the conmutator for all of...

> I am aware that when we execute the circuit we are collapsing the wave function so once I have executed it I can't do it again, but I was...

The [Unitary](https://qibo.readthedocs.io/en/stable/api-reference/qibo.html#arbitrary-unitary) gate allows you to add any matrix as a gate to a circuit, for example ```Python import numpy as np from qibo import models, gates c = models.Circuit(3)...

It seems that the number of targets is the total number of qubits in the circuit. The `gates.Unitary` approach will work, but you will probably notice some deterioration if nqubits...

Thanks for proposing this, sounds good to me. Here is a function that could serve as a starting point: ```py def light_cone(circuit, observable_qubits): """Reduces circuit to the qubits relevant for...

@bpcarlos thanks for checking. I believe the redundant gates issue is solved in the following: ```py def light_cone(circuit, observable_qubits): """Reduces circuit to the qubits relevant for an observable. Args: circuit...

Thanks for opening this, it is certainly not possible to do directly with Qibo. For now, I will try to find a temporary solution using custom code and will write...

I found a relatively simple way to do this using our existing `TrotterHamiltonian` functionality. It may be considered cheating because the state vector is used in the calculation, however it...