隨著多核技術的不斷發展，并行方法已經成為一種處理較大規模問題的手段，同時在許多領域取得了成功地應用。目前，并行算法的實現主要基于兩種標準：

MPI（Message Passing Interface）是一種基于消息傳遞并行編程模型的工業標準，主要用于分布存儲體系結構的現實，是已被證實了的理想的程序設計模型;OpenMP標準，主要用于共享存儲體系結構的并行編程，可實現在SMP 集群系統內多處理器的多線程并行計算。

OpenMP是一個應用程序接口，通過規范一系列的編程制導、運行庫函數和環境變量來說明共享存儲體系結構的并行機制，通常由于其較低的開銷和相對較簡單的編程而受到人們廣泛的關注。

時域有限差分（FDTD）方法是模擬計算電磁場的一種基本算法。自1966年Yee首次提出以來，經過半個多世紀的發展，這一方法得到迅速發展和廣泛應用。但是時域有限差分算法通常由于其串行方式使問題本身復雜化且運算較費時間而采用并行方式計算。因此，在PC機上研究并行FDTD算法問題，具有重要的理論與現實意義，可為大規模工程問題的并行化處理提供一定的方法借鑒與理論依據。本文以一維平面波在自由空間中的傳播為例，討論了采用OpenMP 技術對電磁場FDTD算法程序實現并行化的方法，并將該并行方法在三維瞬態場電偶極子輻射FDTD程序中進行了驗證，也相當于對該并行方法進行了一定的推廣，并通過實驗證明了該并行計算的有效性。

1、電磁場理論簡介

FDTD方法由微分形式的麥克斯韋（Maxwell）旋度方程出發進行離散而得到的一組時域推進公式。一維情況下，設TEM波沿z 軸方向傳播，介質參數和場量均與x，y 無關，即- /-x = 0， -/ -y = 0 ，于是Maxwell方程為：

一維情況E 、H 分量空間節點取樣如圖1所示。

在自由空間中，σ = σm = 0 ，介質為無耗，故而可得場的FDTD迭代方程為：

為了滿足數值計算的Courant穩定性條件，通常選取時間步長為 Δt 空間間隔為 δ ，FDTD截斷邊界條件采用一階近似Mur，設截斷邊界處為Ex 節點，如圖2所示，則在Ex 節點處離散，得：

式中： Ex （k - 1）為截斷邊界內的節點;c 在真空中為光速c0，在介質中則為截斷邊界處波的傳播速度。

激勵源采用高斯脈沖源，其表達式為Ei （t），其中τ 為常數，決定了高斯脈沖的寬度。為了使入射波限制在空間有限區域，根據等效原理，在區域分界面上設置等效面電磁流，并設分界面外的場為零。所以，在總場-散射場區的分界面上（總場邊界）設置入射波電磁場的切向分量便可將入射波只引入到總場區。本實驗通過在一維FDTD 的總場邊界處引入高斯脈沖波，如圖3所示。

2、OpenMP 并行設計

2.1 OpenMP概述

OpenMP是基于共享存儲體系結構的工業標準，它不是一門獨立的語言，而是對基本語言的擴展，如C/C++，Fortran語言。其編程簡單，開銷小，規范并制定了一系列的編譯指導語句、運行庫函數和環境變量。對于傳統的串行代碼，采用OpenMP技術并行化時無需對原程序作大的改動，只需加入一些簡單的編譯指導語句即可。同時，OpenMP 提供了兩種粒度的并行方式：粗粒度并行和細粒度并行。OpenMP的細粒度并行是指利用OpenMP 只求解循環部分計算，又稱為循環級并行。由此可見，細粒度并行是一種最為簡單的并行方法。

2.2 OpenMP并行編程模型

OpenMP采用標準的并行模式--Fork/Join式并行執行模式，如圖4所示，在編譯過程中使用編譯指導語句實現并行化。在程序開始執行時，只有主線程的運行線程存在，在執行的過程中，若遇到OpenMP 的指令要求并行執行時，主線程會派生出子線程來執行并行任務。在并行執行的過程中，由主線程與派生出的子線程組成一個線程組來協同工作。在并行執行結束后，派生出的子線程退出或掛起，不再工作，控制流程回到單獨的主線程中，直到下一個并行區或者程序執行完畢。

2.3 FDTD算法的并行化

在用FDTD算法模擬計算電磁場的過程中，時間步長上的迭代過程是相關聯的、互相影響的，故而不能實現并行化。但是在一次迭代內部，電場與磁場的計算僅需要前一時刻的計算結果，與其他區域的電場或者磁場分量無關，各個計算過程之間沒有影響、相互獨立，可以實現并行化。為此，本文采用OpenMP提供的細粒度并行的方式對該算法實現并行化，即OpenMP+細粒度并行。

3、并行FDTD 性能分析

3.1 仿真算例設計

為驗證FDTD算法并行的性能，以電磁波在自由空間中傳播的一維FDTD 算法為例，采用OpenMP 提供的細粒度技術實現并行化。激勵源采用Gauss脈沖源，其帶寬為10 GHz，τ = 0.2 ns， t0 = 0.8 - τ = 0.16 ns，總場邊界為z0 = 500Δz 處，Δt = τ/ 80 = 0.002 5 ns， Δt = dz/ （2c），吸收邊界采用一階近似Mur，波傳播的空間區域為節點1~1 000.用統計的方法，分別測量了2 000~10 000 之間不同迭代次數的串并行時間。

3.2 并行性能測試

并行程序性能測試可由并行算法的加速比和并行效率來衡量，假設有n 個并行部件，則對加速比、效率的定義如下：

（1）加速比speedup=單一計算機運行整個程序所花費時間同一臺計算機使用n 個并行部件的執行時間;

（2）效率efficiency= speedup n.

程序串并行的運行時間由OpenMP 庫函數提供的OMP_get_wtime（）函數來測量。具體地，在測量串行程序運行時間時，注釋掉了程序中的并行編譯指導語句，并保持時間函數的位置不變。本文采用多次運行程序取穩定值的方法，分別測量并記錄了不同迭代次數下的程序串并行運行時間。

3.3 測試環境

測試環境為Intel（R） Core（TM） 2 Duo CPUT5670@1.8 GHz，內存為2 GB，操作系統為WinXP SP3，開發軟件為Intel Fortran 10.1.014 with vs 2005，測試結果如表1所示。

3.4 并行方法驗證

為了驗證本文所采用的OpenMP 并行算法的可行性，將該并行方法應用于三維瞬態場電偶極子輻射FDTD程序中。設整個計算域空間為真空，垂直點偶極子位于計算域中心，即Ez（0，0，0），FDTD計算空間步長為5 cm，時間步長為83.333 ps，計算域為55×55×55 個元胞，截斷邊界為Mur吸收邊界，輻射源為高斯脈沖，測量并記錄了300~10 000之間不同迭代次數的串并行時間，測量結果見表2.

4、結論

本文從分析OpenMP本身的特點及編程模型入手，結合一維FDTD算法實例，采用OpenMP+細粒度并行的方式實現了并行化，并證明了基于OpenMP 的并行FDTD 算法的有效性，而且并行FDTD 算法在所選測試實例的不同迭代次數上均獲得了超線性的加速比。充分利用了OpenMP共享存儲體系結構的特點，避免了消息傳遞帶來的開銷，取得了較其他并行FDTD算法更快加速比和更高的效率。更值得一提的是，將該并行方法應用在三維電磁場FDTD 程序中也取得了很好的加速比和效率，可以預見，該優化方法在更復雜的算法中應用一定具有更理想的性能提升。但不足之處在于系統的可擴展性差，這是由于OpenMP 本身的特點所限制。

因此，今后的工作將放在對基于SMP 機群的MPI 與OpenMP混合編程模型的研究，從而克服系統擴展性差的缺點，進而提升系統的易用性和可移植性。