用传入数据填充复杂矩阵的最快方法。 Armadillo图书馆

Question

我的情况是，每个1MS都将到达2*200个INT16_T样本的新缓冲区。 （乘以2是因为复杂的数据）

应以最快的方式（小于1 ms）将这些数据添加到Armadillo Complect Double矩阵中。。

有人知道更好的方法吗？

当前的测试代码需要超过1ms（两种情况测试）

nb_of_samples = 200;


    void write(uint32_t rawId, uint32_t slotId, int16_t *samples) 
    {
        size_t currentBlock = slotId * nb_of_samples; // slot identification 
        auto const memPtr = matrix.colptr(currentBlock) + rawId; // pointer to position from where to start write new received samples
        std::for_each(std::execution::par_unseq, m_index.begin(), m_index.end(), [&memPtr, &samples](size_t index)
        {
            new (memPtr + index) std::complex<double>(samples[index], samples[index+1]);
        });
//        second idea, simple for 
//        for(size_t i = 0; i < NB_OF_SAMPLES; ++i)
//        {
//            new (memPtr + i) std::complex<double>(samples[i], samples[i+1]);
//        }
    }

Answer 1

C ++是std :: Complex的混叠和UB规则的特殊例外。这样做是为了兼容广泛使用的C子例程。

https://en.cppreference.com/w/cpp/numeric/complex

因此，您可以使用这样的简单的东西：

void convert(std::complex<double> *pout, const uint16_t* pin)
{
    for (size_t i = 0; i < nb_of_samples*2 ; i++)
        reinterpret_cast<double*>(pout)[i] = pin[i];
}