Function mm

src/backend/cpu/sparse_blas.cpp:344–369

template<typename T, bool conjugate>
void mm(Param<T> output, CParam<T> values, CParam<int> rowIdx,
        CParam<int> colIdx, CParam<T> right, int M, int N, int ldb, int ldc) {
    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();

    for (int o = 0; o < N; ++o) {
        for (int i = 0; i < rowIdx.dims(0) - 1; ++i) {
            outPtr[i] = scalar<T>(0);
            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]];
                }
            }
        }
        rightPtr += ldb;
        outPtr += ldc;
    }
}

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