并行计算实验一：VS2017+OpenMP环境配置

1、实验目标

（1） 利用visual studio 2017配置OpenMP环境；

打开VisualStudio 2017菜单栏->文件->新建->项目：新建Visual C++工程 :Windows控制台应用程序

Figure 1 新建Parallel_Lab1_Martrix-Vecter项目
打开菜单栏->项目->Parallel_Lab1_Martrix-Vecter（视自己新建项目名，最下面一栏）属性->配置属性->C/C++->语言->OpenMp 支持：选择是(/openmp)。（这一步是为了支持openmp。）点击“应用”。

Figure 2 打开项目的OpenMP支持

打开菜单栏->项目->属性-> Parallel_Lab1_Martrix-Vecter（视自己新建项目名，最下面一栏）属性->配置属性->代码生成—>运行库：选择多线程调试(/Mtd)。（这一步是为了支持多线程。）点击“应用”。

Figure 3 打开项目的多线程支持

打开菜单栏->项目->属性-> Parallel_Lab1_Martrix-Vecter（视自己新建项目名，最下面一栏）属性->配置属性->调试->命令参数: 在空白栏填入main函数的输入参数（argv[i]的值）。如需多个参数，用空格分隔每个参数，此处暂编辑输入“12”。

Figure 4 设置main函数的输入参数

（2）对以下算法利用OpenMP实现并行化。
i. Martrix-Vecter 乘法;

ii梯形积分法;

Figure 5 梯形积分法示意图
iii奇偶变换排序

2、串行程序代码

Martrix-Vector乘法：

• Function: Mat_vect_mult
• Purpose: Multiply a matrix by a vector
• In args: A: the matrix
• x: the vector being multiplied by A
• m: the number of rows in A and components in y
• n: the number of columns in A components in x
• Out args: y: the product vector Ax

/*------------------------------------------------------------------
void Mat_vect_mult(
                   double  A[]  /* in  */, 
                   double  x[]  /* in  */, 
                   double  y[]  /* out */,
                   int     m    /* in  */, 
                   int     n    /* in  */) {
   int i, j;

   for (i = 0; i < m; i++) {
      y[i] = 0.0;
      for (j = 0; j < n; j++)
         y[i] += A[i*n+j]*x[j];
   }
}  /* Mat_vect_mult */
*/

梯形积分法：

/*------------------------------------------------------------------
 * Function:    Trap
 * Purpose:     Estimate integral from a to b of f using trap rule and
 *              n trapezoids
 * Input args:  a, b, n, h
 * Return val:  Estimate of the integral 
 */
double Trap(double a, double b, int n, double h) {
   double integral;
   int k;

   integral = (f(a) + f(b))/2.0;
   for (k = 1; k <= n-1; k++) {
     integral += f(a+k*h);
   }
   integral = integral*h;

   return integral;
}  /* Trap */

/*------------------------------------------------------------------
 * Function:    f
 * Purpose:     Compute value of function to be integrated
 * Input args:  x
 */
double f(double x) {
   double return_val;

   return_val = x*x;
   return return_val;
}  /* f */

奇偶变换排序：

/*-----------------------------------------------------------------
 * Function:     Odd_even_sort
 * Purpose:      Sort list using odd-even transposition sort
 * In args:      n
 * In/out args:  a
 */
void Odd_even_sort(
      int  a[]  /* in/out */, 
      int  n    /* in     */) {
   int phase, i, temp;

   for (phase = 0; phase < n; phase++) 
      if (phase % 2 == 0) { /* Even phase */
         for (i = 1; i < n; i += 2) 
            if (a[i-1] > a[i]) {
               temp = a[i];
               a[i] = a[i-1];
               a[i-1] = temp;
            }
      } else { /* Odd phase */
         for (i = 1; i < n-1; i += 2)
            if (a[i] > a[i+1]) {
               temp = a[i];
               a[i] = a[i+1];
               a[i+1] = temp;
            }
      }
}  /* Odd_even_sort */

未完待续