Function getRandomDiagonalUnitary

tests/deprecated/test_utilities.cpp:667–675

}

QMatrix getRandomDiagonalUnitary(int numQb) {
    DEMAND( numQb >= 1 );

    QMatrix matr = getZeroMatrix(1 << numQb);
    for (size_t i=0; i<matr.size(); i++)
        matr[i][i] = expI(getRandomReal(0,4*M_PI));

    return matr;
}

QMatrix getRandomUnitary(int numQb) {
    DEMAND( numQb >= 1 );