どのように私はWindowsでアクセラレーション（リンゴDSPライブラリ）関数を複製できますか？

Question

私はできるだけ簡潔にするつもりです：

私は非常に特定のハードウェアの制約のために、Windowsに移植する必要があるプロジェクトがあります。 AppleのDSPライブラリであるAccelerateを使用してベクトル距離計算を実行するユーティリティクラスがあります。私はこのライブラリを使わずに機能するように書き直す必要がありますが、適切な置き換えを見つけることができませんでした。私の最善の行動は何ですか？

#include <Accelerate/Accelerate.h> 

inline float distBetween(float *x, float *y, unsigned int count) { 
    float *tmp = (float*)malloc(count * sizeof(float)); 
    // float tmp[count]; 
    //t = y - x 
    vDSP_vsub(x, 1, y, 1, tmp, 1, count); 
    //t.squared 
    vDSP_vsq(tmp, 1, tmp, 1, count); 
    //t.sum 
    float sum; 
    vDSP_sve(tmp, 1, &sum, count); 
    delete tmp; 
    return sqrt(sum); 
} 

inline float cosineDistance(float *x, float *y, unsigned int count) { 
    float dotProd, magX, magY; 
    float *tmp = (float*)malloc(count * sizeof(float)); 

    vDSP_dotpr(x, 1, y, 1, &dotProd, count); 

    vDSP_vsq(x, 1, tmp, 1, count); 
    vDSP_sve(tmp, 1, &magX, count); 
    magX = sqrt(magX); 

    vDSP_vsq(y, 1, tmp, 1, count); 
    vDSP_sve(tmp, 1, &magY, count); 
    magY = sqrt(magY); 

    delete tmp; 

    return 1.0 - (dotProd/(magX * magY)); 
} 

Answer 1

ベクトル機能は、通常、特定のアセンブリ言語命令によって実装されます。この実装は非常に遅いです。おそらくSSE命令を使用するライブラリが必要です。

あなたのコードでは、すべての引数stride_x、stride_y、stride_resは1に等しいので、それらを関数の引数から削除することをお勧めします。 Сodeは速くなければなりません。

//t = y - x  
float 
vDSP_vsub(float *x, int stride_x, float *y, int stride_y, float *res, int stride_res, int count) 
{ 
    while(count > 0) 
    { 
     // may be *x - *y ? 
     *res = *y - *x; 
     res += stride_res; 
     x += stride_x; 
     y += stride_y; 
     count--; 
    }  
} 

//t.squared 
float 
vDSP_vsq(float *x, int stride_x, float *res, int stride_res, int count) 
{ 
    while(count > 0) 
    { 
     *res += (*x) * (*x); 
     x += stride_x; 
     res += stride_res; 
     count--; 
    }  
} 

//t.sum 
float 
vDSP_sve(float *x, int stride_x, float *res, int count) 
{ 
    *res = 0.0; 
    while(count > 0) 
    { 
     *res += *x; 
     x += stride_x; 
     count--; 
    }  
} 

float 
vDSP_dotpr(float *x, int stride_x, float *y, int stride_y, float *res, int count) 
{ 
    *res = 0.0; 
    while(count > 0) 
    { 
     *res += (*x) * (*y); 
     x += stride_x; 
     y += stride_y; 
     count--; 
    }  
} 

Answer 2

Intel's IPP librariesをご覧ください。