pulse.use_ideal = True print "initial state", y0, "| expected order", select_order params.initial_state("quantum", qstate = y0) tic = time.time() if doKitaev: result = Kit.simulateevolution(pulseseq, params, dec) else: data = qc.simulateevolution(pulseseq, params, dec) result = U.displaytracedstates(data.YRPN[-1], pop=True) if not doIdeal: print "runtime: ", time.time()-tic, "sec" else: tic = time.time() if doKitaev: result = Kit.simulateevolution(pulseseq, params, dec) print np.around(result, 3) else: data = qc.simulateevolution(pulseseq, params, dec) result = U.displaytracedstates(data.YRPN[-1], pop=True) print np.around(result, 3) if not doIdeal: print "runtime: ", time.time()-tic, "sec" if saveKitaev: qc.saveRun(pulseseq, params, dec, Kit.data_group, 'OF'+str(select_order)+'-werr.shlv') Kit.getQFTall() Kit.resultSave('OF'+str(select_order)+'-result.npy')
# turn off noise for some pulses if doIdeal: for pulse in pulseseq: pulse.use_ideal = True if doRun: tic = time.time() data = qc.simulateevolution(pulseseq, params, dec) data.RhoPNAll = np.array(data.RhoPNAll) toc = time.time() print("runtime: ", toc - tic, "seconds") if params.savedata: print("saving run data to file") qc.saveRun(pulseseq, params, dec, data, params.savedataname) ### evaluate and plot fidelities if calcfinal: evalobj = evd.EvaluateData() if doRun: evalobj.loadidealdata(idealdata) evalobj.loadsimdata(data) evalobj.calculatePlotFidelities(displ, plottitle) #, grouperrors=100) else: print("no data included, doRun must be set to True") ### include spontaneous decay in measurement if calc_spontdecay_during_measurement: pop_decayed = sim.DecayedPopulations_CCD(data.YRPN[-1], params) print(pop_dcayed)
# turn off noise for some pulses if doIdeal: for pulse in pulseseq: pulse.use_ideal = True if doRun: tic = time.time() data = qc.simulateevolution(pulseseq, params, dec) data.RhoPNAll = np.array(data.RhoPNAll) toc = time.time() print "runtime: ", toc-tic, "seconds" if params.savedata: print "saving run data to file" qc.saveRun(pulseseq, params, dec, data, params.savedataname) ### evaluate and plot fidelities if calcfinal: evalobj = evd.EvaluateData() if doRun: evalobj.loadidealdata(idealdata) evalobj.loadsimdata(data) evalobj.calculatePlotFidelities(displ,plottitle)#, grouperrors=100) else: print "no data included, doRun must be set to True" ### include spontaneous decay in measurement if calc_spontdecay_during_measurement: pop_decayed = sim.DecayedPopulations_CCD(data.YRPN[-1], params) print pop_dcayed