Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/Qiskit-Extensions/mthree
Matrix-free Measurement Mitigation
https://github.com/Qiskit-Extensions/mthree
Last synced: 17 days ago
JSON representation
Matrix-free Measurement Mitigation
- Host: GitHub
- URL: https://github.com/Qiskit-Extensions/mthree
- Owner: Qiskit-Extensions
- License: apache-2.0
- Created: 2021-08-16T15:47:22.000Z (almost 3 years ago)
- Default Branch: main
- Last Pushed: 2024-04-18T21:16:42.000Z (2 months ago)
- Last Synced: 2024-06-10T16:34:58.060Z (20 days ago)
- Language: Python
- Homepage: https://qiskit.org/ecosystem/mthree/
- Size: 2.92 MB
- Stars: 34
- Watchers: 3
- Forks: 24
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
- License: LICENSE.txt
Lists
- awesome-qiskit - mthree - Matrix-free measurement mitigation (m3). (Extensions)
README
# mthree
[](https://opensource.org/licenses/Apache-2.0)
[](https://badge.fury.io/py/mthree)
[](https://pypi.org/project/mthree/)
Matrix-free Measurement Mitigation (M3).
M3 is a measurement mitigation technique that solves for corrected measurement probabilities using a dimensionality reduction step followed by either direct LU factorization or a preconditioned iterative method that nominally converges in O(1) steps, and can be computed in parallel. For example, M3 can compute corrections on 42 qubit GHZ problems in under two seconds on a quad-core machine (depending on the number of unique bitstrings in the output).
## Documentation
[Online Documentation @ Qiskit.org](https://qiskit.org/ecosystem/mthree/)
## Installation
You can `pip` install M3 in serial mode using PyPi via:
```bash
pip install mthree
```
Alternatively, one can install from source:
```bash
python setup.py install
```
To enable openmp one must have an openmp 3.0+ enabled compiler and install with:
```bash
python setup.py install --openmp
```
### OpenMP on OSX
On OSX one must install LLVM using homebrew (You cannot use GCC):
```bash
brew install llvm
```
after which the following (or the like) must be executed in the terminal:
```bash
export PATH="/usr/local/opt/llvm/bin:$PATH"
```
and
```bash
export LDFLAGS="-L/usr/local/opt/llvm/lib -Wl,-rpath,/usr/local/opt/llvm/lib"
export CPPFLAGS="-I/usr/local/opt/llvm/include"
```
Then installation with openmp can be accomplished using:
```bash
CC=clang CXX=clang python setup.py install --openmp
```
## Usage
### Basic usage
M3 is simple to use:
```python
import mthree
# Specify a mitigator object targeting a given backend
mit = mthree.M3Mitigation(backend)
# Compute the 1Q calibration matrices for the given qubits and given number of shots
# By default it is over all backend qubits at 10000 shots.
mit.cals_from_system(qubits, shots)
# Apply mitigation to a given dict of raw counts over the specified qubits
m3_quasi = mit.apply_correction(raw_counts, qubits)
```
Note that here `qubits` is a list of which qubits are measured to yield the bits in the output.
For example the list `[4,3,1,2,0]` indicates that a measurement on physical qubit 4 was written to
classical bit zero in the output bit-strings, physical qubit 3 maps to classical bit 1, etc.
The fact that the zeroth bit is right-most in the bitstring is handled internally.
### Error bounds
It is possible to compute error bounds in a similarly efficient manner. This is not done
by default, but rather turned on using:
```python
m3_quasi = mit.apply_correction(raw_counts, qubits, return_mitigation_overhead=True)
```
Then the distribution itself can be called to return things like the expectation value and the
standard deviation:
```python
expval, stddev = quasi.expval_and_stddev()
```
### Closest probability distribution
The results of M3 mitigation are quasi-probabilities that nominally contain small negative values.
This is suitable for use in computing corrected expectation values or sampling problems
where one is interested in the highest probability bit-string. However, if one needs
a true probability distribution then it is possible to convert from quasi-probabilites to
the closest true probability distribution in L2-norm using:
```python
closest_probs = m3_quasi.nearest_probability_distribution()
```
## License
[](https://opensource.org/licenses/Apache-2.0)