I am using an error mitigation algorithm that uses the readout errors p(0|1) and p(1|0) to correct the errors due to noise in my original circuit. I am trying to measure the readout errors in the same circuit as shown in the image. As shown in the circuit, I am resetting the two qubits twice.
The first part of the circuit corresponds to my original algorithm and in the second and the third parts I am measuring the readout error of the two qubits by observing the counts.
I've tried calculating the readout error rates using a separate circuit but I notice that the readout error rates fluctuate considerably within short periods of time, so they might change between the execution of my original circuit and readout error measurement circuit. Also, it is tricky to map the readout error measurement circuit to the same physical qubits which are being used by my original circuit.
I know that I can use the library function
backend.properties().readout_error(n), but it returns the average of p(0|1) and p(1|0) and it stays constant over time but the actual measured readout error fluctuates over time.
- Is this a reasonable way to calculate the readout error rates of the two qubits?
- Is there a better way to accomplish what I am trying to do?