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
