Implement algorithm from here: https://arxiv.org/abs/1403.2975

The canonical implementation of this algorithm is in Haskell and is GPL-licensed, so cannot be used within Qiskit.

This would be useful e.g. for using the T gate simulator in Aer, studying error correction, etc.