[英特尔多核平台编码优化大赛的优化过程-补充.docxVIP

英特尔多核平台编码优化大赛的优化过程--补充分类：?代码优化2007-01-20 17:20?2521人阅读?评论(2)?收藏?举报 英特尔多核平台编码优化大赛的优化过程--补充??HouSisong@GM 2007.01.20tag: 多核编程,sse2,牛顿迭代,代码优化,优化大赛,invsqrt,开方 主要文章请参看我的《英特尔多核平台编码优化大赛的优化过程》：/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?stdio.h#include?stdlib.h#include?math.h#include?windows.h#includexmmintrin.h//使用SSE1#includeprocess.h#includevector//////////////////////////////////////////////////////////////////////////////////////////////////////////////////作者：侯思松?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?time.h#endif#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?(?(result16)??RAND_MAX?);?}?};?CMyRand?_MyRand;?inline?int?_my_rand?(void){?return?_MyRand._my_rand();?}#else?#define?_my_rand??rand#endifint?rand(?void?);int?computePot