optimization - 有效使用英特尔编译器 SVML `__m128 _mm_sincos_ps ()`

Question

我有一个简单的循环:

for (ii = 0; ii < numRows * numCols; ii++)
    {
        mCOmega[ii] = cosf(paramOmega * mI[ii]);
        mSOmega[ii] = sinf(paramOmega * mI[ii]);
    }

}

我想通过使用 __m128 _mm_sincos_ps () 使用英特尔编译器 + SVML 进行优化.

我正在努力寻找从 __m128 _mm_sincos_ps () 的输出加载数据并存储数据的有效方法。 .

使用 __m128 _mm_sincos_ps () 编写此循环的最有效方法是什么？ ？

Answer 1

以下内容在 godbolt 上编译没有任何问题:

#include <stdio.h>
#include <stdlib.h>
#include <mathimf.h>
#include <immintrin.h>

void test (const float * restrict mI, 
           float * restrict mCOmega,
           float * restrict mSOmega,
           float paramOmega,
           int numRows,
           int numCols)
{   
    __m128 x, sin_x, cos_x;
    int ii;
    for (ii = 0; ii < numRows * numCols; ii = ii + 4)
    {   x = _mm_loadu_ps((float*)&mI[ii]);
        x = _mm_mul_ps(x, _mm_set1_ps(paramOmega));
        sin_x = _mm_sincos_ps(&cos_x, x);
        _mm_storeu_ps((float*)&mCOmega[ii], cos_x);
        _mm_storeu_ps((float*)&mSOmega[ii], sin_x);
    }   
}