समस्या निवारण ऑटो vectorize कारण '1200'

Question

MSVC 2013 परम डब्ल्यू/अद्यतन 4

समझने नहीं क्यों मैं इस प्रतीत होता है सरल उदाहरण

जानकारी C5002 पर इस त्रुटि हो रही है: पाश कारण की वजह से '1200' vectorized नहीं

जो

1200 लूप शामिल पाश-ले गए डेटा dependences

मैं नहीं दिख रहा है कि कैसे एल की पुनरावृत्तियों ओओपी एक दूसरे के साथ हस्तक्षेप कर सकता है।

__declspec(align(16)) class PhysicsSystem 
{ 
public: 
    static const int32_t MaxEntities = 65535; 

    __declspec(align(16)) struct VectorizedXYZ 
    { 
     double  mX[ MaxEntities ]; 
     double  mY[ MaxEntities ]; 
     double  mZ[ MaxEntities ]; 

     VectorizedXYZ() 
     { 
      memset(mX, 0, sizeof(mX)); 
      memset(mY, 0, sizeof(mY)); 
      memset(mZ, 0, sizeof(mZ)); 
     } 
    }; 

    void Update(double dt) 
    { 
     for (int32_t i = 0; i < MaxEntities; ++i) <== 1200 
     { 
      mTmp.mX[ i ] = mPos.mX[ i ] + mVel.mX[ i ] * dt; 
      mTmp.mY[ i ] = mPos.mY[ i ] + mVel.mY[ i ] * dt; 
      mTmp.mZ[ i ] = mPos.mZ[ i ] + mVel.mZ[ i ] * dt; 
     } 
    } 

private:  
    VectorizedXYZ mTmp; 
    VectorizedXYZ mPos; 
    VectorizedXYZ mVel; 
}; 

संपादित करें: http://blogs.msdn.com/b/nativeconcurrency/archive/2012/05/08/auto-vectorizer-in-visual-studio-11-rules-for-loop-body.aspx इस से आंकना "उदाहरण 1 - अजीब के ठंड से समानांतर" का एक उदाहरण होने लगते हैं, लेकिन यह काम करता है जैसे कि यह सोचता है कि सरणियों अलियासिंग है, जो मुझे puzzling है से असुरक्षित हैं।

EDIT2: यह अच्छा होगा अगर किसी कारण ऑटो vectorization इस तरह के एक मालूम होता है सरल उदाहरण पर विफल रहता है साझा कर सकता है, लेकिन कुछ समय के लिए यह संवारता के बाद, मैं बजाय का विकल्प चुना शासन काल अपने आप को लेने के लिए

void PhysicsSystem::Update(Real dt) 
{ 
    const __m128d mdt = { dt, dt }; 

    // advance by 2 since we can do 2 at a time at double precision in __m128d 
    for (size_t i = 0; i < MaxEntities; i += 2) 
    { 
     __m128d posX = _mm_load_pd(&mPos.mX[ i ]); 
     __m128d posY = _mm_load_pd(&mPos.mY[ i ]); 
     __m128d posZ = _mm_load_pd(&mPos.mZ[ i ]); 

     __m128d velX = _mm_load_pd(&mVel.mX[ i ]); 
     __m128d velY = _mm_load_pd(&mVel.mY[ i ]); 
     __m128d velZ = _mm_load_pd(&mVel.mZ[ i ]); 

     __m128d velFrameX = _mm_mul_pd(velX, mdt); 
     __m128d velFrameY = _mm_mul_pd(velY, mdt); 
     __m128d velFrameZ = _mm_mul_pd(velZ, mdt); 

     _mm_store_pd(&mPos.mX[ i ], _mm_add_pd(posX, velFrameX)); 
     _mm_store_pd(&mPos.mY[ i ], _mm_add_pd(posX, velFrameY)); 
     _mm_store_pd(&mPos.mZ[ i ], _mm_add_pd(posX, velFrameZ)); 
    } 
} 

Answer 1

सुनिश्चित करें कि आपके संकलक यह समर्थन करता है तो, लेकिन कुछ उचित vectorisation लागू करने के लिए, आप portably ऐसा कर सकते हैं नहीं:

void PhysicsSystem::Update(double dt) { 
    double *tx=mTmp.mX, *ty=mTmp.mY, *tz=mTmp.mZ; 
    double *px=mPos.mX, *py=mPos.mY, *pz=mPos.mZ; 
    double *vx=mVel.mX, *vy=mVel.mY, *vz=mVel.mZ; 
    #pragma omp simd aligned(tx, ty, tz, px, py, pz, vx, vy, vz) 
    for (int i = 0; i < MaxEntities; ++i) { 
     tx[ i ] = px[ i ] + vx[ i ] * dt; 
     ty[ i ] = py[ i ] + vy[ i ] * dt; 
     tz[ i ] = pz[ i ] + vz[ i ] * dt; 
    } 
} 

आप निर्देश के लिए OpenMP समर्थन सक्षम करने की जरूरत है तो ध्यान में रखा जाना करने के लिए।