英特尔多核平台编码优化大赛的优化过程--补充
分类： 代码优化2007-01-20 17:20 2521人阅读 评论(2) 收藏 举报
英特尔多核平台编码优化大赛的优化过程--补充 HouSisong@GMail.com 2007.01.20 tag: 多核编程,sse2,牛顿迭代,代码优化,优化大赛,invsqrt,开方 主要文章请参看我的《英特尔多核平台编码优化大赛的优化过程》：
http://blog.csdn.net/housisong/archive/2007/01/20/1488465.aspx 本文章是其补充；提供一个完整的float实现版本、double到float的手工转换、手工得到invSqrt的粗略起始迭代值 等其它几个不算成功的实现； 我测试和优化过程中用的 CPU：AMD64x2 3600+ (双核CPU) 操作系统：Windows XP 32bit 编译器：Visual Studio 2005 大赛公布的原始代码执行时间 3.97秒
一个float完整实现版本(牺牲了计算精度)，源代码如下：
(如果用汇编实现应该还可以把速度提高一些，或者使用ICC编译器:) /* compute the potential energy of a collection of */ /* particles interacting via pairwise potential */ #include
#include
////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//作者：侯思松 HouSisong@263.net //计算结果的精度为float单精度浮点 版本
#define _IS_FAST
//以牺牲精度的方式加快速度，否则就运算达到float单精度
////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//#define _NEW_TIME_CLOCK
#ifdef _NEW_TIME_CLOCK #define clock_t double double CLOCKS_PER_SEC=0.0; inline double clock() { __int64 result;
if (CLOCKS_PER_SEC==0) {
QueryPerformanceFrequency((LARGE_INTEGER *)&result); CLOCKS_PER_SEC=(double)result; }
QueryPerformanceCounter((LARGE_INTEGER *)&result); return (double)result; } #else
#include
#define _IS_USES_MY_RAND
//单线程执行rand函数，所以使用自定义rand是安全的
const long DefthreadCount=2; //1,2,4,8,16,.. 把计算任务分成多个任务并行执行
float& m128_value(__m128& x,const long index) { return ((float*)(&x))[index]; }
#define NPARTS 1000 #define NITER 201 #define DIMS 3
#ifdef _IS_USES_MY_RAND class CMyRand {
private:
unsigned long _my_holdrand; public:
CMyRand():_my_holdrand(1){} inline int _my_rand (void) {
unsigned long result=_my_holdrand * 214013 + 2531011; _my_holdrand =result;
return ( (result>>16) & RAND_MAX ); } };
CMyRand _MyRand;
inline int _my_rand (void){ return _MyRand._my_rand(); } #else
#define _my_rand rand #endif
int rand( void ); int computePot();
void initPositions(void); void updatePositions(void);
__declspec(align(16)) float r[DIMS][(NPARTS+3)/4*4]; //16byte对齐 double pot;
int main() { int i;
clock_t start, stop;
//char ctmp; std::cin>>ctmp;
initPositions(); updatePositions();
start=clock();
for( i=0; i printf(\ updatePositions(); } pot = 0.0; stop=clock(); printf (\double)(stop-start)/ CLOCKS_PER_SEC); return 0; } void initPositions() { int i, j; for( i=0; i r[i][j]= (float)( 0.5 + _my_rand() *(1.0/RAND_MAX) ); } }
