Pyqrack

QUnit` layer over `QStabilizerHybrid` is theoretically capable of transparently achieving stabilizer simulation efficiency on pure Clifford circuits (as well as some near-Clifford circuits). This capability regressed, however. This release fixes it.

This releases is concurrent with Qrack v7.1.1, which fixes a bug in `QUnit` multi-qubit buffer commutation.

This release fixes an edge case in the Schmidt decomposition "layer" when operating "over" stabilizer layer, which sometimes caused incorrect results. (Specifically, when Schmidt decomposition was previously applied to stabilizer subsystems, it was possible that a two-qubit subsystem would be separated while only updating the probability amplitude caches of one of the two qubits, whereas both need to be updated.)

Work on v0.18.0 prompted attention to the basic implementation of `QStabilizerHybrid` in the underlying C++ Qrack library. This uncovered a conceptual mistake in how we previously attempted measurement through non-Clifford single qubit gate buffers. Now, the minimum necessary edge cases of non-Clifford measurement resort to state vector representation before measurement. This improves accuracy (hopefully to the point of exact universal `QStabilizerHybrid` simulation, up to machine epsilon) at the expense of performance.

Improvements in stabilizer simulation have been made with the "reverse t-injection gadget" in Appendix A of [Fast Estimation of Outcome Probabilities for Quantum Circuits](https://journals.aps.org/prxquantum/abstract/10.1103/PRXQuantum.3.020361). These method improvements can be turned on or off for a `QrackSimulator` with the method `set_t_injection(bool)`. (They are on by default.) While implementing this, other minor improvements were made in `QStabilizerHybrid`.