Function operator *

tests/utils/linalg.cpp:504–514

qvector operator * (const qmatrix& m, const qvector& v) {
    DEMAND( m.size() == v.size() );

    qvector out = getZeroVector(v.size());

    for (size_t r=0; r<v.size(); r++)
        for (size_t c=0; c<v.size(); c++)
            out[r] += m[r][c] * v[c];

    return out;
}