def profile_unitary_simulator(self):
"""Profile randomly generated circuits.
Writes profile results to <this_module>.prof as well as recording
to the log file.
number of circuits = 100.
number of operations/circuit in [1, 40]
number of qubits in [1, 5]
"""
nCircuits = 100
maxDepth = 40
maxQubits = 5
pr = cProfile.Profile()
randomCircuits = RandomQasmGenerator(seed=self.seed,
maxDepth=maxDepth,
maxQubits=maxQubits)
randomCircuits.add_circuits(nCircuits, doMeasure=False)
self.qp = randomCircuits.getProgram()
pr.enable()
self.qp.execute(self.qp.get_circuit_names(),
backend='local_unitary_simulator')
pr.disable()
sout = io.StringIO()
ps = pstats.Stats(pr, stream=sout).sort_stats('cumulative')
logging.info('------- start profiling UnitarySimulator -----------')
ps.print_stats()
logging.info(sout.getvalue())
logging.info('------- stop profiling UnitarySimulator -----------')
sout.close()
pr.dump_stats(self.moduleName + '.prof')
def profile_unitary_simulator(self):
"""Profile randomly generated circuits.
Writes profile results to <this_module>.prof as well as recording
to the log file.
number of circuits = 100.
number of operations/circuit in [1, 40]
number of qubits in [1, 5]
"""
nCircuits = 100
maxDepth = 40
maxQubits = 5
pr = cProfile.Profile()
randomCircuits = RandomQasmGenerator(seed=self.seed,
maxDepth=maxDepth,
maxQubits=maxQubits)
randomCircuits.add_circuits(nCircuits, doMeasure=False)
self.qp = randomCircuits.getProgram()
pr.enable()
self.qp.execute(self.qp.get_circuit_names(),
backend='local_unitary_simulator')
pr.disable()
sout = io.StringIO()
ps = pstats.Stats(pr, stream=sout).sort_stats('cumulative')
self.log.info('------- start profiling UnitarySimulator -----------')
ps.print_stats()
self.log.info(sout.getvalue())
self.log.info('------- stop profiling UnitarySimulator -----------')
sout.close()
pr.dump_stats(self.moduleName + '.prof')
def profile_nqubit_speed_constant_depth(self):
"""simulation time vs the number of qubits
where the circuit depth is fixed at 40. Also creates a pdf file
with this module name showing a plot of the results. Compilation
is not included in speed.
"""
import matplotlib.pyplot as plt
from matplotlib.ticker import MaxNLocator
qubitRangeMax = 15
nQubitList = range(1, qubitRangeMax + 1)
maxDepth = 40
minDepth = 40
nCircuits = 10
shots = 1024
seed = 88
maxTime = 30 # seconds; timing stops when simulation time exceeds this number
fmtStr1 = 'profile_nqubit_speed::nqubits:{0}, backend:{1}, elapsed_time:{2:.2f}'
fmtStr2 = 'backend:{0}, circuit:{1}, numOps:{2}, result:{3}'
fmtStr3 = 'minDepth={minDepth}, maxDepth={maxDepth}, num circuits={nCircuits}, shots={shots}'
backendList = ['local_qasm_simulator', 'local_unitary_simulator']
if shutil.which('qasm_simulator'):
backendList.append('local_qasm_cpp_simulator')
else:
self.log.info(
'profile_nqubit_speed::\"qasm_simulator\" executable not in path...skipping'
)
fig = plt.figure(0)
plt.clf()
ax = fig.add_axes((0.1, 0.2, 0.8, 0.6))
for i, backend in enumerate(backendList):
elapsedTime = np.zeros(len(nQubitList))
if backend is 'local_unitary_simulator':
doMeasure = False
else:
doMeasure = True
j, timedOut = 0, False
while j < qubitRangeMax and not timedOut:
nQubits = nQubitList[j]
randomCircuits = RandomQasmGenerator(seed,
minQubits=nQubits,
maxQubits=nQubits,
minDepth=minDepth,
maxDepth=maxDepth)
randomCircuits.add_circuits(nCircuits, doMeasure=doMeasure)
qp = randomCircuits.getProgram()
cnames = qp.get_circuit_names()
qobj = qp.compile(cnames,
backend=backend,
shots=shots,
seed=seed)
start = time.perf_counter()
results = qp.run(qobj)
stop = time.perf_counter()
elapsedTime[j] = stop - start
if elapsedTime[j] > maxTime:
timedOut = True
self.log.info(fmtStr1.format(nQubits, backend, elapsedTime[j]))
if backend is not 'local_unitary_simulator':
for name in cnames:
self.log.info(
fmtStr2.format(backend, name,
len(qp.get_circuit(name)),
results.get_data(name)))
j += 1
ax.xaxis.set_major_locator(MaxNLocator(integer=True))
if backend is 'local_unitary_simulator':
ax.plot(nQubitList[:j],
elapsedTime[:j],
label=backend,
marker='o')
else:
ax.plot(nQubitList[:j],
elapsedTime[:j] / shots,
label=backend,
marker='o')
ax.set_yscale('log', basey=10)
ax.set_xlabel('number of qubits')
ax.set_ylabel('process time/shot')
ax.set_title('profile_nqubit_speed_constant_depth')
fig.text(
0.1, 0.05,
fmtStr3.format(minDepth=minDepth,
maxDepth=maxDepth,
nCircuits=nCircuits,
shots=shots))
ax.legend()
self.pdf.savefig(fig)
def profile_nqubit_speed_constant_depth(self):
"""simulation time vs the number of qubits
where the circuit depth is fixed at 40. Also creates a pdf file
with this module name showing a plot of the results. Compilation
is not included in speed.
"""
import matplotlib.pyplot as plt
from matplotlib.ticker import MaxNLocator
qubitRangeMax = 15
nQubitList = range(1,qubitRangeMax + 1)
maxDepth = 40
minDepth = 40
nCircuits = 10
shots = 1024
seed = 88
maxTime = 30 # seconds; timing stops when simulation time exceeds this number
fmtStr1 = 'profile_nqubit_speed::nqubits:{0}, backend:{1}, elapsed_time:{2:.2f}'
fmtStr2 = 'backend:{0}, circuit:{1}, numOps:{2}, result:{3}'
fmtStr3 = 'minDepth={minDepth}, maxDepth={maxDepth}, num circuits={nCircuits}, shots={shots}'
backendList = ['local_qasm_simulator', 'local_unitary_simulator']
if shutil.which('qasm_simulator'):
backendList.append('local_qasm_cpp_simulator')
else:
self.log.info('profile_nqubit_speed::\"qasm_simulator\" executable not in path...skipping')
fig = plt.figure(0)
plt.clf()
ax = fig.add_axes((0.1, 0.2, 0.8, 0.6))
for i, backend in enumerate(backendList):
elapsedTime = np.zeros(len(nQubitList))
if backend is 'local_unitary_simulator':
doMeasure = False
else:
doMeasure = True
j, timedOut = 0, False
while j < qubitRangeMax and not timedOut:
nQubits = nQubitList[j]
randomCircuits = RandomQasmGenerator(seed,
minQubits=nQubits,
maxQubits=nQubits,
minDepth=minDepth,
maxDepth=maxDepth)
randomCircuits.add_circuits(nCircuits, doMeasure=doMeasure)
qp = randomCircuits.getProgram()
cnames = qp.get_circuit_names()
qobj = qp.compile(cnames, backend=backend, shots=shots, seed=seed)
start = time.perf_counter()
results = qp.run(qobj)
stop = time.perf_counter()
elapsedTime[j] = stop - start
if elapsedTime[j] > maxTime:
timedOut = True
self.log.info(fmtStr1.format(nQubits, backend, elapsedTime[j]))
if backend is not 'local_unitary_simulator':
for name in cnames:
self.log.info(fmtStr2.format(
backend, name, len(qp.get_circuit(name)),
results.get_data(name)))
j += 1
ax.xaxis.set_major_locator(MaxNLocator(integer=True))
if backend is 'local_unitary_simulator':
ax.plot(nQubitList[:j], elapsedTime[:j], label=backend, marker='o')
else:
ax.plot(nQubitList[:j], elapsedTime[:j]/shots, label=backend,
marker='o')
ax.set_yscale('log', basey=10)
ax.set_xlabel('number of qubits')
ax.set_ylabel('process time/shot')
ax.set_title('profile_nqubit_speed_constant_depth')
fig.text(0.1, 0.05,
fmtStr3.format(minDepth=minDepth, maxDepth=maxDepth,
nCircuits=nCircuits, shots=shots))
ax.legend()
self.pdf.savefig(fig)