固有値を使って式をベクトル化する

Question

私のプログラムでは、固有ライブラリを使用していますが、私は2Dベクトルで操作する必要があります。私は、以下の機能を持っている私の内部ループでは：両方x()とy()計算は、同一の経路をたどるよう

static inline double eval(double x, double y, double xi, double yi) 
{ 
    const double invlen2 = 1/(x*x + y*y); 
    const double invlen4 = invlen2*invlen2; 
    const double invlen6 = invlen4*invlen2; 

    const double x2 = x*x, y2 = y*y; 
    const double x3 = x2*x, y3 = y2*y; 
    const double xi2 = xi*xi, yi2 = yi*yi; 

    return x*invlen2 + invlen4*(x2*xi + 2*x*y*yi - xi*y2) 
    + invlen6*(x3*xi2 + 3*x*y2*yi2 + 6*x2*y*xi*yi - 3*x*xi2*y2 - 2*y3*xi*yi - x3*yi2); 
} 

void f(Vector2d& out, const Vector2d& R, const Vector2d& r) 
{ 
    out.x() = eval(R.x(), R.y(), r.x(), r.y()); 
    out.y() = eval(R.y(), R.x(), r.y(), r.x()); 
} 

この表現、厄介なもののは、ベクトル化のための最有力候補のように思えます。私の質問は、手動でアセンブリにドロップする必要なく、Eigenでそれを行う方法です。

Answer 1

この回答はEigenとは関係ありませんが、手作業でアセンブリにドロップダウンすると述べたので、これを追加します。

アセンブリを使用してコードをベクトル化する必要はありません。

http://software.intel.com/sites/products/documentation/studio/composer/en-us/2011/compiler_c/intref_cls/common/intref_overview.htm#intref_overview

言った：固有値はすでにベクトル化のための内部サポートを持っていますが、あなたの例に適用可能であると思われないように見えます組み立てずに手動でベクトル化できるようになるのコンパイラ組み込み関数があります。だから私はなぜそれを手動でやりたいのか分かります。