Function isApproxUnitary

tests/utils/linalg.cpp:250–256

bool isApproxUnitary(qmatrix m) {

    // should be identity
    qmatrix md = m * getConjugateTranspose(m);
    qmatrix id = getIdentityMatrix(m.size());
    return doMatricesAgree(md, id);
}


qcomp getTrace(qmatrix m) {