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。）點選“應用”。

開啟選單欄->專案->屬性-> Parallel_Lab1_Martrix-Vecter（視自己新建專案名，最下面一欄）屬性->配置屬性->程式碼生成—>執行庫：選擇多執行緒除錯(/Mtd)。（這一步是為了支援多執行緒。）點選“應用”。

Figure 3 開啟專案的多執行緒支援

開啟選單欄->專案->屬性-> Parallel_Lab1_Martrix-Vecter（視自己新建專案名，最下面一欄）屬性->配置屬性->除錯->命令引數: 在空白欄填入main函式的輸入引數（argv[i]的值）。如需多個引數，用空格分隔每個引數，此處暫編輯輸入“12”。

（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 */

未完待續