I am using the qndiag library to try to find a diagonalisation for 2 given matrices.
The github is here : qndiag libray
The function qndiag is defined like this (n