自然对数底E值并行计算.ppt

自然对数底E值的并行计算 数学10-02班 朱建伟 摘要 e是数学中最重要的数学符号之一，称为自然常数，是自然对数的底数。它最先由瑞士数学家欧拉在1727 年使用。利用级数的方法求解e的值具有计算精确，速度较快的特点。其基本计算公式为： 本文主要采用多线程技术，对e值的计算进行加速，得出计算时间与加速比，从而判断多线程技术加速效果，进而证明多核技术在数值计算方面的优越性。 并且文中采用不同的并行方式来比较不同结果，查看哪种方式在并行计算中具有优越性，获得最终结果。 关键字 e的计算 多线程 加速比 高性能计算 引言 e值在数学史上具有十分重要的意义，所以对e的计算进行研究也从未间断。已经有多种方式对e值进行快速精确的计算，于是想到依托于现有算法，对其进行改造获得多线程算法，查看其计算速度，获得其加速比。通过对其加速效果的观察，获得多线程计算式中高性能计算方式，可以用于数值计算。 重要常量值 const int numSteps = 2000000;//计算终止值 CRITICAL_SECTION cs;//临界区声明 int k = 1; double e = 1.0; double fact = 1.0; //获得电脑线程数（核数） int numOfProcessors; SYSTEM_INFO SysInfo; GetSystemInfo(SysInfo); numOfProcessors = SysInfo.dwNumberOfProcessors; 串行计算 double temp_e = 1.0; double temp_fact = 1.0; for (int i=1 ; inumSteps ; i++) { temp_fact *= I; temp_e += 1.0 / temp_fact; } 利用临界区进行同步 方式一： EnterCriticalSection(cs); while (k numSteps) { fact *= k; e += 1.0 / fact; k++; } LeaveCriticalSection(cs); 方式二： while (k numSteps) { fact *= k; e += 1.0 / fact; EnterCriticalSection(cs); k++; LeaveCriticalSection(cs); } 利用线程ID号控制 int threadID = *((int*)arg); double fact = 1.0; for (int j=1; j=threadID+1; j++) fact *= j; e += 1.0 / fact; for (int i=threadID; inumSteps; i+=numOfProcessors) { for (int j=2; jnumOfProcessors+2; j++) fact *= (i + j); e += 1.0 / fact; } 利用Open MP并行 double temp_e1 = 0.0; double *temp_fact1 = new double[numOfProcessors]; double temp = 1.0; int j = 0; int end_num = (numOfProcessors-1) *(numSteps/numOfProcessors); //初始化 for (int i=1; i=end_num+1; i++) { if (i == (numSteps/numOfProcessors)*j+1) { temp_fact1[j] = temp; j ++; } temp *= I; } #pragma omp parallel for firstprivate(temp_fact1) reduction(+:temp_e1) for (int i=1 ; i=numSteps ; i++) { temp_fact1[omp_get_thread_num()] *= I; temp_e1 += 1.0 / temp_fact1[omp_get_thread_num()]; } 加速比 序号 计算方式 计算结果 运行时间/s 平均时间/s 加速比 1 串行计算 2.71828 0.671 0.655 0.66300 1.000 0.64 0.686 2 临界区 2.71828 0.39