Function getOuterProduct

tests/utils/linalg.cpp:220–230

qmatrix getOuterProduct(qvector ket, qvector bra) {
    DEMAND( bra.size() == ket.size() );

    qmatrix out = getZeroMatrix(bra.size());

    for (size_t i=0; i<ket.size(); i++)
        for (size_t j=0; j<ket.size(); j++)
            out[i][j] = ket[i] * std::conj(bra[j]);

    return out;
}