I am trying to run the QFT algorithm for n=5 (n number of qubits). The number of experiments that it generates is bigger than the one allowed by the IBM devices (i.e. 75). Therefore, I have tried to used the IBMQJobManager but this is not working as my result is not the counts but the density matrix obtain from the tomography process. Is it possible to make it work? I have made the following attempt (qft is a function that I have created that implements the QFT algorithm):
n=5 beta=0.5*np.pi state=initial_state_QFT(n,beta) qr=QuantumRegister(n) cr=ClassicalRegister(n) qc=QuantumCircuit(qr) qc.initialize(state,qr) qft(qc,n) qst = tomo.state_tomography_circuits(qc, qr) job = job_manager.run(qst, backend=device, name='foo',shots=8192) statefit = tomo.StateTomographyFitter(job.results(), qst) p, M, weights = statefit._fitter_data(True, 0.5) M_dg = np.conj(M).T linear_inversion_matrix = np.linalg.inv(M_dg @ M) @ M_dg rho = linear_inversion_matrix @ p rho = np.reshape(rho, (2**n, 2**n))
Thank you in advance for your answer!