SSE整数除算？

Question

浮動小数点値の除算に_mm_div_psがあります。整数乗算に_mm_mullo_epi16があります。しかし、整数除算（16ビット値）のための何かがありますか？どのように私はそのような分裂を行うことができる？

Answer 1

彼は8ビット、16ビット、32ビット・ワード（ただし64ビット）http://www.agner.org/optimize/#vectorclass

ルックのためにSSE/AVXと整数の除算を行うための高速なアルゴリズムが実装されていAgner霧のvectorclassを参照してください。コードについてはvectori128.hファイル、algoirthmの説明はよく書かれたマニュアルであるVectorClass.pdf

ここには、彼のマニュアルのアルゴリズムを説明する断片があります。

「整数の除算 あり整数ベクトル部門の 便利です何のx86命令セット内の命令とその拡張されず、このような命令は、彼ら が存在する場合、非常に遅くなる。したがって、ベクトルクラスライブラリが使用していますこのアルゴリズムの基本原理は、次の式で表すことができます。 a/b≈a *（2n/b）>> n この計算は、次の手順を実行します。 1.適切なものを見つけるnの値 2. 2n/bを計算する 3.丸め誤差に必要な修正を計算する 4.複数の計算を行うこの問題は、複数の数値を同じ除数で除算すると有利です。 b。被除数aの各値に対してステップ4が繰り返されている間に、ステップ1,2および3は1回だけ実行すればよい。数学的な詳細は、ファイル vectori128.hに記載されています。 （T. GranlundとPLモンゴメリーも参照してください：乗算、SIGPLANの議事録を使用して不変 整数除算を "...

編集：ファイルvectori128.hの終わり近くには、との短い除算を行う方法を示していますスカラー変数 "高速除算に使用するパラメータの計算には、 よりも時間がかかるため、同じ除数オブジェクト を複数回使用すると便利です。たとえば、符号なしの短い整数を10：

short x = 10; 
uint16_t dividends[80], quotients[80];   // numbers to work with 
Divisor_us div10(x);       // make divisor object for dividing by 10 
Vec8us temp;         // temporary vector 
for (int i = 0; i < 80; i += 8) {    // loop for 4 elements per iteration 
    temp.load(dividends+i);     // load 4 elements 
    temp /= div10;        // divide each element by 10 
    temp.store(quotients+i);     // store 4 elements 
} 

「

編集：ショートパンツ

#include <stdio.h> 
#include "vectorclass.h" 

int main() {  
    short numa[] = {10, 20, 30, 40, 50, 60, 70, 80}; 
    short dena[] = {10, 20, 30, 40, 50, 60, 70, 80}; 

    Vec8s num = Vec8s().load(numa); 
    Vec8s den = Vec8s().load(dena); 

    Vec4f num_low = to_float(extend_low(num)); 
    Vec4f num_high = to_float(extend_high(num)); 
    Vec4f den_low = to_float(extend_low(den)); 
    Vec4f den_high = to_float(extend_high(den)); 

    Vec4f qf_low = num_low/den_low; 
    Vec4f qf_high = num_high/den_high; 
    Vec4i q_low = truncate_to_int(qf_low); 
    Vec4i q_high = truncate_to_int(qf_high); 

    Vec8s q = compress(q_low, q_high); 
    for(int i=0; i<8; i++) { 
     printf("%d ", q[i]); 
    } printf("\n"); 
}