2016-02-15 44 views
-1

我正在寻找一种最佳实践解决方案,以优化嵌套数组的代码。来自一个数组的值作为第二个索引。我认为目前的代码不是CPU缓存友好的代码。我可以使用openmp来优化此代码的性能吗?如果是,如果PICK不是很好,那么最好使用“for”或“sections”10?为嵌套数组值优化代码

#define bino4d(a,b,c,d) (base1[a]+base2[b]+base3[c]+base4[d]) 

static unsigned char _bits_[8] = { 0x80, 0x40, 0x20, 0x10, 0x08, 0x04, 0x02, 0x01 }; 
#define TSTBIT(ch,i) (ch[(i)>>3] & _bits_[(i) & 0x07]) 

lpmax = (PICK*(PICK - 1)*(PICK - 2)*(PICK - 3))/24; 

    for (lp = 0; lp < lpmax; lp++) 
    { 

     j0 = I[LP4_0[lp]]; 
     j1 = I[LP4_1[lp]]; 
     j2 = I[LP4_2[lp]]; 
     j3 = I[LP4_3[lp]]; 
     ul = j0 + j1 + j2 + j3; 

     //j0 = LP4_0[lp]; 
     //j1 = LP4_1[lp]; 
     //j2 = LP4_2[lp]; 
     //j3 = LP4_3[lp]; 
     //ul = bino4d(I[j0], I[j1], I[j2], I[j3]); 
     if (TSTBIT(ToCover, ul)) s4++; 
    } 
+0

什么I的大小,什么是LP4_ * – user3528438

+0

LP4_的范围*范围等于210为使该代码例如PICK = 5 –

回答

1

通过使用一个额外的数组存储的ul的部分结果,您可能使计算更多本地及获得更好的缓存友好性或实现并行。假设我们有一个数组UL[lpmax]=0

for (lp = 0; lp < lpmax; lp++) 
    UL[lp] += I[LP4_0[lp]]; 

for (lp = 0; lp < lpmax; lp++) 
    UL[lp] += I[LP4_1[lp]]; 

for (lp = 0; lp < lpmax; lp++) 
    UL[lp] += I[LP4_2[lp]]; 

for (lp = 0; lp < lpmax; lp++) 
    UL[lp] += I[LP4_3[lp]]; 

for (lp = 0; lp < lpmax; lp++) 
    if (TSTBIT(ToCover, UL[lp])) s4++; 
+0

有趣IDEIA。我会检查它。 –