SSEを使用する対数、またはFPUに切り替えますか？

Question

統計計算をいくつか行っています。私はそれらを速くする必要があるので、私はSSEを使うためにそのほとんどを書き直しました。私はかなり新しいので、ここで正しいアプローチが何であるか疑問に思っていました：

私の知る限り、SSEにはlog2またはln関数はありません。少なくとも4.1までは最新バージョンです私が使用するハードウェアでサポートされています。

1. エキス4台の山車、そしてenthropyを決定するためにそれらの上にFPUの計算を行う - 私はそれらをまとめる、バックSSEレジスタにそれらの値のいずれかをロードする必要はありません。

   はそれがより良いのです別のフロート

2. は対数関数を実装して何のSSE命令はありませんLOG2

Answer 1

周りにいくつかのSSE log2実装があるようです。 this one。

log2の機能を持つIntel Approximate Maths Libraryもあります。これは古いですが（2000年）、それはSSE2であり、やはり正常に動作するはずです。

---

も参照してください：

• sse_mathfun - SSE vector math library
• avx_mathfun - AVX vector math library
• libmvec - vector math library added in glibc 2.22

Answer 2

を行うSSEための機能を見つけます。しかし、一般的な対数を実行する単一のx86命令もありません。 logやlog10のような対数関数をC標準ライブラリから使用することを考えているなら、libcのようなオープンソースライブラリで使用されている実装を見てみる価値があります。 SSEレジスタのすべての要素で動作する独自の対数近似を容易にロールバックできます。

このような関数は、テイラー級数のような入力引数の特定の領域に対してある精度仕様で有効な多項式近似を使用して実装されることがよくあります。対数プロパティを利用して、汎用入力引数を対数ルーチンの許容入力範囲にラップすることができます。また、あなたが財産を利用することによって対数の底をパラメータ化することができますaはあなたが作成した対数ルーチンのベースである

log_y(x) = log_a(x)/log_a(y) 

。