将 SSE 内在函数用于 C 中的二维数组

问题描述

我想用SSE来优化这段代码：

v[0][0] = 2*vi-v[1][0];
for (j=1;j<=jmax;j++) {
    u[0][j] = ui;
    v[0][j] = 2*vi-v[1][j];
}

到目前为止我所做的是：

v[0][0] = 2*vi-v[1][0];
for (j=1;j<=jmax;j++) {
    __m128 uiVec = _mm_set1_ps(ui);
    __m128 viVec = _mm_set1_ps(vi);
    __m128 constant = _mm_set1_ps(2);

    __m128 v1Vec = _mm_load_ps(v[1]+j);

    _mm_store_ps(v[0]+j,uiVec);

    __m128 v0Calc = _mm_mul_ps(constant,viVec);
    __m128 v0Result = _mm_sub_ps(v0Calc,v1Vec);

    _mm_store_ps(v[0]+j,v0Result);
}

但是我遇到了分段错误。当我检查代码时，问题出在 _mm_store_ps 函数中。 我该如何解决？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）