C中的矢量化Trig函数？

Question

我正在计算高度并行化的trig函数（在像1024的块中），并且我想利用至少一些现代体系结构所具有的并行性。

当我编译块

for(int i=0; i<SIZE; i++) { 
    arr[i]=sin((float)i/1024); 
} 

GCC不会向量化它，并说

not vectorized: relevant stmt not supported: D.3068_39 = __builtin_sinf (D.3069_38); 

这对我来说很有意义。但是，我想知道是否有一个库来执行并行触发计算。

只有一个简单的泰勒系列上升到第11位，GCC将向量化所有的循环，并且我的速度比无辜罪循环快两倍（具有比特精确的答案，或者9阶系列，对于1600倍数值中的最后两个，只有一个比特关闭，对于> 3倍加速）。我确定以前有人遇到过这样的问题，但是当我google时，我发现没有提及任何库或类似的东西。

A.有没有东西存在？ B.如果不是，建议优化并行触发函数？

编辑：我发现了以下库称为“SLEEF”：http://shibatch.sourceforge.net/它在this论文中描述，并使用SIMD指令来计算几个基本功能。它使用SSE和AVX特定代码，但我认为将其转换为标准C循环并不困难。

Answer 1

我的回答是创建自己的图书馆做的正是这个叫vectrig：https://github.com/jeremysalwen/vectrig

Answer 2

您使用的平台是？这种类型的许多库已经存在：

• 英特尔提供了带有icc的矢量数学库（VML）。
• Apple提供了vForce库作为Accelerate框架的一部分。
• HP为Itanium提供了自己的矢量数学库（也可能是其他架构）。
• Sun提供了libmvec及其编译器工具。
• ...

Answer 3

既然你正在寻找在这里计算的谐波，我有一些code that addressed a similar problem。它的矢量化速度已经超过我发现的任何其他东西。作为一个好处，你可以免费获得余弦。

Answer 4

既然你说像有一些选项你使用GCC它看起来：

• http://gruntthepeon.free.fr/ssemath/
  • 它使用SSE和SSE2指令来实现它。
• http://www.gamasutra.com/view/feature/4248/designing_fast_crossplatform_simd_.php
  • 这具有的替代实施方式。一些评论非常好。

这就是说，我可能会考虑GPGPU解决的办法。也许写在CUDA或OpenCL（如果我没有记错CUDA支持正弦函数）。以下是一些看起来可能会更容易的库。

• https://code.google.com/p/slmath/
• https://code.google.com/p/thrust/

Answer 5

取而代之的是泰勒级数的，我想看看fdlibm使用的算法。他们应该以更少的步骤获得尽可能多的精度。