我想使用 neon 指令 SIMD 和臂组件来计算 4 个点的大小和角度。大多数语言都有一个内置库，在我的例子中是 C++，它计算角度 (atan2)，但仅针对一对浮点变量（x 和 y）。我想利用处理 q 寄存器的 SIMD 指令来计算 4 个值向量的 atan2 。
精度要求不高，速度更重要。
我已经有一些汇编指令可以计算 4 个浮点寄存器的大小，对于我的应用程序来说具有可接受的精度。 q1 包含 4 个“x”值（x1、x2、x3、x4）。 q2 包含 4 个“y”值（y1、y2、y3、y4）。 q7 包含 4 个结果的大小 (x1^2 + y1^2、x2^2 + y2^2、x3^2 + y3^2、x4^2 + y4^2)。
vmul.f32 q7, q1, q1
vmla.f32 q7, q2, q2
vrecpe.f32 q7, q7
vrsqrte.f32 q7, q7
使用 SIMD 指令计算两个向量的近似 atan2 的最快方法是什么？
See 数学霓虹灯对于现有的单值浮点实现。由于它没有（或很少）条件，它应该很好地转换为SIMD执行。
由于 ARM NEON 没有直接计算该值的指令，因此有多种技术可以创建比直接计算更好的近似值泰勒级数。具体来说，min-max方法给出了一个很好的近似多项式候选。min-max指最小化最大误差；与一个切比雪夫近似通常都非常好。
DSP guru有不同方法的具体细节函数逼近。网上还有大量书籍。您可以使用搜索最佳多项式matlab, octave或其他一些工具包。通常，您需要将其与range and 精确。一旦你对单个值有了一个好的算法，就可以将其扩展到SIMD任何形式都应该是微不足道的。
问题计算atan2有参考苹果的atan.c来源。代码中的系数很可能源自我上面给出的内容。这段代码的问题是它不能扩展到SIMD以及atan()近似是分段的，根据范围需要不同的系数。为您SIMD，您将需要在整个范围内使用相同的系数（乘数、除数、方程）。
Abramowitz 和 Stegun：数学函数手册有一个章节是关于循环函数第 4.4.28 节给出了对数公式。这似乎类似于eglibc执行。
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
