我建议使用Vc:可移植、零开销的 C++ 类型,用于显式数据并行编程 https://github.com/VcDevel/Vcsimd 库,我听说它的目标是包含到 C++ 标准中。它更容易编写,也更容易阅读。
Example:
#include <iostream>
#include <Vc/Vc>
int main() {
using A = Vc::SimdArray<int, 8>;
A arr1 = A::Random();
A arr2 = A::Random();
std::cout << arr1 << '\n';
std::cout << arr2 << '\n';
std::cout << arr1 / arr2 << '\n';
}
Outputs:
<1513634383 -963914658 1763536262 -1285037745 | -695608406 -35372374 1025922083 444041308>
<824703811 1962744590 1568022524 -293901648 | 549806324 248334095 1663905340 641164273>
[1, 0, 1, 4, -1, 0, 0, 0]
以下功能
using A = Vc::SimdArray<int, 8>;
__attribute__((noinline)) A f(A a0, A a1) {
return a0 / a1;
}
With g++-8.2 -O3 -march=skylake
翻译成以下程序集:
f(Vc_1::SimdArray<int, 8ul, Vc_1::Vector<int, Vc_1::VectorAbi::Avx>, 8ul>, Vc_1::SimdArray<int, 8ul, Vc_1::Vector<int, Vc_1::VectorAbi::Avx>, 8ul>):
vcvtdq2pd ymm3, xmm1
vcvtdq2pd ymm2, xmm0
vextracti128 xmm1, ymm1, 0x1
vextracti128 xmm0, ymm0, 0x1
vcvtdq2pd ymm1, xmm1
vdivpd ymm2, ymm2, ymm3
vcvtdq2pd ymm0, xmm0
vdivpd ymm0, ymm0, ymm1
vcvttpd2dq xmm2, ymm2
vcvttpd2dq xmm0, ymm0
vinserti128 ymm0, ymm2, xmm0, 0x1
ret
请注意,x86 指令集中没有用于整数除法的 simd 指令。