Function mtv

src/backend/cpu/sparse_blas.cpp:319–341

template<typename T, bool conjugate>
void mtv(Param<T> output, CParam<T> values, CParam<int> rowIdx,
         CParam<int> colIdx, CParam<T> right, int M) {
    UNUSED(M);

    const T *valPtr   = values.get();
    const int *rowPtr = rowIdx.get();
    const int *colPtr = colIdx.get();
    const T *rightPtr = right.get();
    T *outPtr         = output.get();

    // Output Array Created is a zero value Array
    // Hence, no need to initialize to zero here
    for (int i = 0; i < rowIdx.dims(0) - 1; ++i) {
        for (int j = rowPtr[i]; j < rowPtr[i + 1]; ++j) {
            // If stride[0] of right is not 1 then rightPtr[i*stride]
            if (conjugate) {
                outPtr[colPtr[j]] += getConjugate(valPtr[j]) * rightPtr[i];
            } else {
                outPtr[colPtr[j]] += valPtr[j] * rightPtr[i];
            }
        }
    }
}

template<typename T, bool conjugate>
void mm(Param<T> output, CParam<T> values, CParam<int> rowIdx,