Description
When the user specifies --circuit_write we write each of the operations/subroutines in their decomposed and Clifford+T forms to QASM files in utils.utils.circuit_estimate()
These files give no indication as to which order they should be applied in order to reconstruct the application circuit. And neither the documentation nor the notebooks address this. I had to look at PyLIQTR.PhaseEstimation.pe.PhaseEstimation.generate_circuit() to get this information to reconstruct the original circuit. It's implicit in the text output that's displayed when running the script, but we never explicitly say this in comments, README, or a Docstring.
Some of this information can be gathered by setting include_nested_resources = True in the gsee_resource_estimation() call.
It could be nice to output the entire application circuit to a single QASM file as well so that someone wanting to do an analysis on the entire circuit doesn't have to reconstruct the circuit. However even if we do the resource estimates of the circuit will not line up with the resource estimates in our json file due to the additions from utils.utils.scale_resources()