SSE to C++ code

SSE to C++ code

本文关键字：code C++ to SSE      更新时间：2023-10-16

我正在尝试将代码从包括SSE指令的c ++源代码重写为仅c ++代码。我知道我会失去表现，但这是一个实验，我正在努力表现。 我想知道是否有一个C++等效于与 、__mm_unpackhi_pd 和 __mm_unpacklo_pd 相同的操作。我对SSE的了解为零。

我正在尝试转换的参考代码片段。任何知识或提示都会有所帮助。谢谢。

for (unsigned chunk = 0; chunk < chunks; chunk++)
{
unsigned start = chunk * chunksize;
unsigned end =
std::min((chunk + 1) * chunksize, (unsigned)2 * w);
__m128d a2b2 =
_mm_load_pd(d_origx +
((2 * init_G_offset + start) & n2_m_1));
unsigned i2_mod_B = 0;
for (unsigned i = start; i < end; i += 2)
{
__m128d ab = a2b2;
a2b2 =
_mm_load_pd(d_origx +
((origx_offset + i) & n2_m_1));
__m128d cd = _mm_load_pd(d_filter + i);
__m128d cc = _mm_unpacklo_pd(cd, cd);
__m128d dd = _mm_unpackhi_pd(cd, cd);
__m128d a0a1 = _mm_unpacklo_pd(ab, a2b2);
__m128d b0b1 = _mm_unpackhi_pd(ab, a2b2);
__m128d ac = _mm_mul_pd(cc, a0a1);
__m128d ad = _mm_mul_pd(dd, a0a1);
__m128d bc = _mm_mul_pd(cc, b0b1);
__m128d bd = _mm_mul_pd(dd, b0b1);
__m128d ac_m_bd = _mm_sub_pd(ac, bd);
__m128d ad_p_bc = _mm_add_pd(ad, bc);
__m128d ab_times_cd = _mm_unpacklo_pd(ac_m_bd, ad_p_bc);
__m128d a2b2_times_cd =
_mm_unpackhi_pd(ac_m_bd, ad_p_bc);
__m128d xy = _mm_load_pd(d_x_sampt + i2_mod_B);
__m128d x2y2 = _mm_load_pd(d_x_sampt + i2_mod_B + 2);
__m128d st = _mm_add_pd(xy, ab_times_cd);
__m128d s2t2 = _mm_add_pd(x2y2, a2b2_times_cd);
_mm_store_pd(d_x_sampt + i2_mod_B, st);
_mm_store_pd(d_x_sampt + i2_mod_B + 2, s2t2);
i2_mod_B += 4;
}
}

下面你找到了这两个函数的描述，我还将每个函数链接到其参考页面。完整的参考资料可在此处获得：https://software.intel.com/sites/landingpage/IntrinsicsGuide/

_mm_unpackhi_p

__m128d _mm_unpackhi_pd (__m128d a, __m128d b)

解包和交错双精度(64 位)浮点 来自 A 和 B 的上半部分的元素，并将结果存储在 DST 中。

_mm_unpacklo_pd

_m128d _mm_unpacklo_pd (__m128d a, __m128d b)

解包和交错双精度(64 位)浮点 来自 A 和 B 的下半部分的元素，并将结果存储在 DST 中。

具体如何实现它取决于您的表示，但基本上您返回一个新值，该值由a的高(或低)一半与b的高(或低)一半连接而成。 例如：

typedef double[2] __m128d;
__m128d _mm_unpackhi_pd(__m128d a, __m128d b) {
__m128d res;
res[0] = a[1];
res[1] = b[1];
return res;
}
__m128d _mm_unpacklo_pd(__m128d a, __m128d b) {
__m128d res;
res[0] = a[0];
res[1] = b[0];
return res;
}

在这个问题上的时机很糟糕... 我在为 SIMDe 实现此功能时发现了这个问题，它只有 17 天的历史。 如果要使用 SIMDe 作为参考，这些功能与许多其他函数一起sse2.h。 SIMDe 中的代码比上面的代码稍微复杂一些，但这主要只是为了匹配其他_mm_unpack*函数的实现。

相关文章：