Function getIdentityMatrix

tests/utils/qmatrix.cpp:30–39

}

qmatrix getIdentityMatrix(size_t dim) {
    DEMAND( dim >= 1 );

    qmatrix out = getZeroMatrix(dim);

    for (size_t i=0; i<dim; i++)
        out[i][i] = 1;
    
    return out;
}

qmatrix getDiagonalMatrix(qvector v) {
    DEMAND( v.size() >= 1 );