def METHOD_NAME(a, b, cd, indices, indptr, residual_global):
    """computes residual_global = a*b - cd at given indices in csr_matrix format"""
    N = a.shape[0]
    M = a.shape[1]
    for i in prange(N):
        for j in prange(indptr[i], indptr[i + 1]):
            num = 0.0
            for k in range(M):
                num += a[i, k] * b[k, indices[j]]
            residual_global[j] = num - cd[j]
    return residual_global