Function mv

src/backend/cpu/sparse_blas.cpp:295–316

template<typename T, bool conjugate>
void mv(Param<T> output, CParam<T> values, CParam<int> rowIdx,
        CParam<int> colIdx, CParam<T> right, int 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 < M; ++i) {
        for (int j = rowPtr[i]; j < rowPtr[i + 1]; ++j) {
            // If stride[0] of right is not 1 then rightPtr[colPtr[j]*stride]
            if (conjugate) {
                outPtr[i] += getConjugate(valPtr[j]) * rightPtr[colPtr[j]];
            } else {
                outPtr[i] += valPtr[j] * rightPtr[colPtr[j]];
            }
        }
    }
}

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