Function getSuperOperator

tests/utils/linalg.cpp:484–495

qmatrix getSuperOperator(vector<qmatrix> matrices) {
    DEMAND( matrices.size() > 0 );

    size_t dim = matrices[0].size();

    // out = sum_m conj(m) (x) m
    qmatrix out = getZeroMatrix(dim * dim);
    for (auto& matr : matrices)
        out += getKroneckerProduct(getConjugate(matr), matr);

    return out;
}