專利名稱：：一種基于tbb線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法
技術(shù)領(lǐng)域：
：本發(fā)明涉及一種基于TBB線程構(gòu)建塊(ThreadBuildingBlocking)的連續(xù)系統(tǒng)仿真多核并行方法，屬于計算機仿真與算法優(yōu)化領(lǐng)域。
背景技術(shù)：
：連續(xù)系統(tǒng)仿真廣泛應(yīng)用于航空、航天、氣象和電子等領(lǐng)域。隨著這些學(xué)科技術(shù)的發(fā)展，由此建立的仿真模型也日益復(fù)雜，仿真計算量不斷加大。連續(xù)系統(tǒng)仿真通常使用微分方程組來建立數(shù)學(xué)模型，微分方程組有成熟的數(shù)值解法，如龍格-庫塔法，亞當(dāng)姆斯法等。但傳統(tǒng)的串行方法在解決大規(guī)模問題時計算時間過長，難以滿足仿真的實時性要求，有的復(fù)雜問題使用串行算法甚至無法求解。提高連續(xù)系統(tǒng)仿真速度的一種有效方法是采用并行技術(shù)，人們在連續(xù)系統(tǒng)微分方程組的并行化上做了大量工作，在串行算法的基礎(chǔ)上提出了各種各樣的并行算法。在早期的研究工作中，Nievergelt提出設(shè)置多個初始值進行并行計算，然后對計算結(jié)果進行修正的并行算法，Shampine和Watts提出了隱式板塊方法和預(yù)估-校正執(zhí)行方式，Miranker和Liniger也做了預(yù)估-校正法的并行研究工作，此外還給出了龍格-庫塔法二階、三階的計算格式，并推導(dǎo)出了MIMD型龍格-庫塔法計算公式，具體做法是采用逼近法來計算導(dǎo)數(shù)值以打斷差值之間的相關(guān)鏈，并適當(dāng)選擇參數(shù)構(gòu)造出穩(wěn)定的并行計算格式。但是傳統(tǒng)的并行算法實現(xiàn)方式繁雜，線程管理、數(shù)據(jù)交換等過程需要程序員管理，這使得并行程序的編寫需要耗費大量的時間和精力，這就阻礙了并行技術(shù)在仿真計算中的應(yīng)用。隨著計算機技術(shù)向著多處理器及多核架構(gòu)發(fā)展，單機上的并行數(shù)據(jù)處理成為可能。TBB是最新推出的多核處理器并行計算平臺，基于任務(wù)的算法、容器和同步基元可簡化并行應(yīng)用程序的開發(fā)工作，它提供完善的提取模板、容器和類別支持并行處理器，可以和MPI、OpenMP等結(jié)合使用，是一種簡單高效的并行方法。伹是，TBB和連續(xù)系統(tǒng)仿真結(jié)合的相關(guān)研究工作還處于初始階段。
發(fā)明內(nèi)容本發(fā)明的目的是為了解決上述問題,提出一種基于TBB線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法，提高連續(xù)系統(tǒng)仿真計算的速度，以縮短仿真研究的周期，滿足仿真的實時性要求。連續(xù)系統(tǒng)仿真一般采用微分方程組建立模型，仿真計算的時間主要花費在微分方程組的求解上。本發(fā)明采用基于TBB的預(yù)估-校正并行算法，在該算法中，預(yù)估計算和校正計算是每一步積分計算中的主要內(nèi)容，他們占據(jù)了仿真計算時間中的絕大部分。通過使用TBB并行技術(shù)，將預(yù)估計算和校正計算作為兩個并行的子任務(wù)，為它們各分配一個線程，兩個子任務(wù)在各自的線程上并行執(zhí)行。TBB采用基于任務(wù)概念的封裝，這使得并行方法的實現(xiàn)更為簡易。通過相關(guān)的模版類，將并行計算的部分分配到各個單獨的線程，預(yù)估計算和校正計算兩個子任務(wù)分別在兩個單獨的線程中執(zhí)行，這樣可以充分利用多個處理器內(nèi)核的資源，達到加快計算速度的目的。一種基于TBB線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法，主要包括以下步驟步驟一搭建TBB并行計算平臺；步驟二構(gòu)造并行ParaODES模版類；步驟三模型初始化；步驟四TBB任務(wù)調(diào)度初始化；步驟五仿真時間范圍判斷；步驟六調(diào)用并行模版類ParaODES;步驟七仿真結(jié)果傳遞；步驟八結(jié)束TBB任務(wù)調(diào)度。本發(fā)明的優(yōu)點在于(1)本發(fā)明利用TBB實現(xiàn)連續(xù)系統(tǒng)仿真的并行化處理，提高了仿真運算的速度，從而滿足仿真的實時性要求，縮短了仿真研究周期；(2)本發(fā)明利用TBB進行并行化處理，過程簡單靈活，易于操作，具有很好的擴展性;(3)TBB采用多核架構(gòu)，符合當(dāng)今并行技術(shù)向多核平臺發(fā)展的方向。圖l是本發(fā)明的方法流程圖2是本發(fā)明在兩個處理內(nèi)核P,、P2上并行執(zhí)行示意圖3是傳統(tǒng)的串行方法和本發(fā)明所述的基于TBB的并行方法的運行時間曲線比較圖。具體實施例方式下面將結(jié)合附圖和實例對本發(fā)明作進一步的詳細說明。本發(fā)明的一種基于TBB線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法，流程如圖l所示，包括以下步驟步驟一搭建TBB并行計算平臺；在計算機中置入TBB文件包，在VC目錄中添加TBB的包含文件、庫文件及可執(zhí)行文件，最后添加TBB20—INSTALL_DIR的值為C:\Programfiles\Intel\TBB\2.0到系統(tǒng)環(huán)境變量。步驟二構(gòu)造并行ParaODES模版類；預(yù)估計算和校正計算是可以并行處理的子任務(wù)，根據(jù)并行化處理對象的特點選擇使用parallelJor模版類。parallel_for模版類是TBB中最常用的一個并行模版類，它把一個循環(huán)劃分為多個塊，然后在不同的線程中并行處理，它要求循環(huán)中不存在循環(huán)依賴。根據(jù)parallel—for模版編寫需要的并行類，并將其命名為ParaODES模版類，該名稱是由用戶自己定義的，名稱的異同對結(jié)果沒有影響。ParaODES模版類的調(diào)用最終完成預(yù)估計算和校正計算。具體步驟為1、編寫預(yù)估計算和校正計箅的程序。模版類ParaODES中的預(yù)估計算和校正計算按以下公式進行預(yù)估計算,+1(8尸-5^+4乂:2-D(1)H《)(2)校正計算/=^+^(9,+19",-5乂:2+/:3)(3)r=/(")(4)其中尸=化，力(5)乂〉/K)(6)式中力表示狀態(tài)變量下一步的估計值，乂表示狀態(tài)變量當(dāng)前步的校正值，A表示狀態(tài)變量上一步的校正值，y;p表示狀態(tài)變量導(dǎo)數(shù)當(dāng)前步的估計值，/:表示狀態(tài)變量導(dǎo)數(shù)下一步的估計值，y>/:、>/:分別表示狀態(tài)變量導(dǎo)數(shù)當(dāng)前步、前一步、前兩步和前三步的校正值，A表示當(dāng)前時刻，"表示上一時刻。如圖2所示，將式(1)、(2)和式(3)、(4)分別在兩個處理內(nèi)核P,、&上并行執(zhí)行,在P,中計算(，然后根據(jù)(的值計算厶V，在P2中計算乂，然后根據(jù)乂的值計算7T。PpP2在結(jié)束了一步計算之后進行通訊，而沒有在計算過程中進行通訊，從而提高了計算速度。根據(jù)上述的ParaODES模版類中的預(yù)估計算和校正計算，編寫預(yù)估計算和校正計算程序，程序中預(yù)估計算和校正計算的部分必須是循環(huán)體結(jié)構(gòu)，預(yù)估計算和校正計算本來不是循環(huán)結(jié)構(gòu)，但可以改寫為循環(huán)結(jié)構(gòu)。例如在預(yù)估計算和校正計算程序前各自增加一個i-l或i=2的執(zhí)行條件語句，通過判斷循環(huán)變量i的值分別執(zhí)行即可。預(yù)估計算和校正計算為循環(huán)體結(jié)構(gòu)，是為了方便使用paralleLfor模版類。循環(huán)體結(jié)構(gòu)中的代碼要盡量使用局部變量或模版類ParaODES的成員變量，而不能使用全局變量，否則不僅不能加快運行速度，反而會增加程序運行時間。2、在模版類ParaODES的operator接口中將循環(huán)參數(shù)修改為blocked_range模板類。在模版類ParaODES中，operator接口是并行處理的主要部分，它對并行循環(huán)進行優(yōu)6化，將最外層循環(huán)參數(shù)修改成TBB中定義的blockecLrange模板類，使之能夠支持循環(huán)體內(nèi)任務(wù)的并行劃分，在調(diào)用并行計算的模板類ParaODES時，由模板類參數(shù)指定循環(huán)處理的數(shù)值范圍以及任務(wù)粒度參數(shù)，循環(huán)體內(nèi)實現(xiàn)預(yù)估計算和校正計算過程。并行計算模版類ParaODES構(gòu)造完成。步驟三模型初始化；在仿真開始時，首先對仿真模型進行初始化，即設(shè)置連續(xù)系統(tǒng)微分方程組右函數(shù)、問題規(guī)模n,n為方程組中方程的數(shù)目、仿真總時間Tmax、積分步長h;導(dǎo)數(shù)、狀態(tài)量、中間量的變量聲明以及各變量賦初值。步驟四TBB任務(wù)調(diào)度初始化；TBB由任務(wù)調(diào)度器對象task—scheduler一init實現(xiàn)多任務(wù)的分配和并行計算，支持對多線程的劃分。由于TBB任務(wù)調(diào)度的初始化會帶來相對較大的額外開銷，因此，只需進行一次初始化。在任務(wù)調(diào)度初始化以后，才可能進行連續(xù)系統(tǒng)仿真多核并^^計算。步驟五仿真時間范圍判斷；判斷當(dāng)前時間t是否超出仿真時間范圍Tmax,若KTmax,則執(zhí)行步驟六，否則執(zhí)行步驟八。步驟六調(diào)用并行模版類ParaODES;調(diào)用步驟二中的ParaODES模版類完成預(yù)估-校正的并行計算，通過ParaODES模版類創(chuàng)建一個線程計算,+1的估計值，并根據(jù),+1計算導(dǎo)數(shù)的估計值^;通過ParaODES模版類創(chuàng)建另一個線程并行計算乂的精確值，并根據(jù)乂計算導(dǎo)數(shù)的精確值尸；最終得到仿真結(jié)果乂°在上述過程中，需要用到的導(dǎo)數(shù)值當(dāng)前步導(dǎo)數(shù)估計值y;p、前一步導(dǎo)數(shù)校正值y;:、前兩步導(dǎo)數(shù)校正值/:和前三步導(dǎo)數(shù)校正值0并得到新的導(dǎo)數(shù)估計值y;:和導(dǎo)數(shù)校正值"，計算厶p,和尸是為了能進行下一步仿真而進行的。'在連續(xù)系統(tǒng)仿真中，一般需要進行多步積分計算，通過循環(huán)調(diào)用ParaODES模版類來實現(xiàn)。步驟七仿真結(jié)果傳遞；將步驟六中的當(dāng)前步導(dǎo)數(shù)估計值y;p、前一步導(dǎo)數(shù)校正值y;:、前兩步導(dǎo)數(shù)校正值o前三步導(dǎo)數(shù)校正值《,和狀態(tài)變量當(dāng)前步的校正值,這五個參數(shù)的值傳遞給當(dāng)前變量"后一步，即將/二的值傳遞給《，將《,的值傳遞給《，將r的值傳遞給《，將y;^的值傳遞給〃，將乂的值傳遞給>^,時間t增加一個單位長度。返回步驟五。'步驟八結(jié)束TBB任務(wù)調(diào)度。完成并行計算任務(wù)，結(jié)束TBB任務(wù)調(diào)度。和TBB任務(wù)調(diào)度初始化一樣，TBB任務(wù)調(diào)度的結(jié)束也會帶來較大的額外開銷，所以本方法只進行一次TBB的初始化和結(jié)束，否則額外的時間開銷會降低并行程序的運行速度。實施例并行處理的對象是如下所示的連續(xù)系統(tǒng)一階微分方程組Wo)",o"'=1,2"""。式中>^2，..乂為狀態(tài)變量，乂。為狀態(tài)變量初始值，y,0。)為/。時刻的狀態(tài)變量值，^為初始時間，t為時間，n代表方程組規(guī)模，即方程組所含方程數(shù)目。積分步長h設(shè)為O.Ol,仿真時間范圍從0.04到20。分別采用傳統(tǒng)的串行方法和本發(fā)明的TBB并行方法，在lntel酷睿2雙核1.83GHz處理器上，進行連續(xù)系統(tǒng)仿真，結(jié)果如表1所示表1連續(xù)系統(tǒng)仿真結(jié)果<table>tableseeoriginaldocumentpage8</column></row><table>從表1可以看出分別采用傳統(tǒng)的串行方法和本發(fā)明的TBB并行方法，在方程組規(guī)模分別為200、400、600、800的時候，并行方法加速比分別為1.04、1.22、1.42、1.73，并行方法加速比隨著方程組規(guī)模的增大而增大。本發(fā)明的TBB并行方法微分方程的解算速度比傳統(tǒng)的串行方法速度快。運行時間如圖3所示，橫坐標(biāo)為方程數(shù)目，縱坐標(biāo)為解算時間，從圖中可以看出在方程組數(shù)目相同的情況下，本發(fā)明的TBB并行方法的解算時間要少，當(dāng)方程數(shù)目擴大到800后,本發(fā)明的TBB并行方法的解算時間縮短到原來的57.58%,大大加快了方程的求解速度。通過對CPU利用情況的比較，可以得出并行化處理后，CPU的使用率得到了提高。串行算法的CPU最大使用率為70%,沒有充分利用兩個處理核心的資源，而使用TBB并行化處理后，以方程組規(guī)模為400時為例，CPU的最大使用率達到了100%,這就充分利用了多核CPU的計算能力，從而提高了連續(xù)系統(tǒng)仿真的速度。權(quán)利要求1、一種基于TBB線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法，其特征在于，包括以下步驟步驟一搭建TBB并行計算平臺；在計算機中置入TBB文件包，在VC目錄中添加TBB的包含文件、庫文件及可執(zhí)行文件，最后添加TBB20_INSTALL_DIR的值為C:\Programfiles\Intel\TBB\2.0到系統(tǒng)環(huán)境變量；步驟二構(gòu)造并行ParaODES模版類；預(yù)估計算和校正計算使用parallel_for模版類，根據(jù)parallel_for模版構(gòu)造并行類，并將其命名為ParaODES模版類；具體步驟為一、編寫預(yù)估計算和校正計算的程序；模版類ParaODES中的預(yù)估計算和校正計算按以下公式進行預(yù)估計算<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msubsup><mi>y</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mo>=</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>c</mi></msubsup><mo>+</mo><mfrac><mi>h</mi><mn>3</mn></mfrac><mrow><mo>(</mo><mn>8</mn><msubsup><mi>f</mi><mi>i</mi><mi>p</mi></msubsup><mo>-</mo><mn>5</mn><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>c</mi></msubsup><mo>+</mo><mn>4</mn><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>2</mn></mrow><mi>c</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>3</mn></mrow><mi>c</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0002"num="0002"><math><![CDATA[<mrow><msubsup><mi>f</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>p</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math></maths>校正計算<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msubsup><mi>y</mi><mi>i</mi><mi>c</mi></msubsup><mo>=</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>c</mi></msubsup><mo>+</mo><mfrac><mi>h</mi><mn>24</mn></mfrac><mrow><mo>(</mo><mn>9</mn><msubsup><mi>f</mi><mi>i</mi><mi>p</mi></msubsup><mo>+</mo><mn>19</mn><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>c</mi></msubsup><mo>-</mo><mn>5</mn><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>2</mn></mrow><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>3</mn></mrow><mi>c</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0004"num="0004"><math><![CDATA[<mrow><msubsup><mi>f</mi><mi>i</mi><mi>c</mi></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>,</mo><msubsup><mi>y</mi><mi>i</mi><mi>c</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math></maths>其中<mathsid="math0005"num="0005"><math><![CDATA[<mrow><msubsup><mi>f</mi><mi>i</mi><mi>p</mi></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>,</mo><msubsup><mi>y</mi><mi>i</mi><mi>p</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0006"num="0006"><math><![CDATA[<mrow><msubsup><mi>f</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>c</mi></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mi>c</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math></maths>式中yi+1p表示狀態(tài)變量下一步的估計值，yic表示狀態(tài)變量當(dāng)前步的校正值，yi-1c表示狀態(tài)變量上一步的校正值，fip表示狀態(tài)變量導(dǎo)數(shù)當(dāng)前步的估計值，fi+1p表示狀態(tài)變量導(dǎo)數(shù)下一步的估計值，fic、fi-1c、fi-2c、fi-3c分別表示狀態(tài)變量導(dǎo)數(shù)當(dāng)前步、前一步、前兩步和前三步的校正值，ti表示當(dāng)前時刻，ti-1表示上一時刻；式(1)、(2)和式(3)、(4)分別在兩個處理內(nèi)核P1、P2上并行執(zhí)行，在P1中計算yi+1p，然后根據(jù)yi+1p的值計算fi+1p；在P2中計算yic，然后根據(jù)yic的值計算fic；P1、P2在結(jié)束了一步計算之后進行通訊；二、在模版類ParaODES的operator接口中將循環(huán)參數(shù)修改為blocked_range模板類；將operator接口的最外層循環(huán)參數(shù)修改成TBB中定義的blocked_range模板類，使之能夠支持循環(huán)體內(nèi)任務(wù)的并行劃分；并行計算模版類ParaODES構(gòu)造完成；步驟三模型初始化；在仿真開始時，首先對仿真模型進行初始化，即設(shè)置連續(xù)系統(tǒng)微分方程組右函數(shù)、問題規(guī)模n、仿真總時間Tmax、積分步長h；導(dǎo)數(shù)、狀態(tài)量、中間量的變量聲明以及各變量賦初值，其中n為方程組中方程的數(shù)目；步驟四TBB任務(wù)調(diào)度初始化；步驟五仿真時間范圍判斷；判斷當(dāng)前時間t是否超出仿真時間范圍Tmax，若t＜Tmax，則執(zhí)行步驟六，否則執(zhí)行步驟八；步驟六調(diào)用并行模版類ParaODES；調(diào)用步驟二中的ParaODES模版類完成預(yù)估-校正的并行計算，通過ParaODES模版類創(chuàng)建一個線程計算yi+1p的估計值，并根據(jù)yi+1p計算導(dǎo)數(shù)的估計值fi+1p；通過ParaODES模版類創(chuàng)建另一個線程并行計算yic的精確值，并根據(jù)yic計算導(dǎo)數(shù)的精確值fic；最終得到仿真結(jié)果yic；在上述過程中，需要用到的導(dǎo)數(shù)值當(dāng)前步導(dǎo)數(shù)估計值fip、前一步導(dǎo)數(shù)校正值fi-1c、前兩步導(dǎo)數(shù)校正值fi-2c和前三步導(dǎo)數(shù)校正值fi-3c，并得到新的導(dǎo)數(shù)估計值fi+1p和導(dǎo)數(shù)校正值fic；步驟七仿真結(jié)果傳遞；將步驟六中的當(dāng)前步導(dǎo)數(shù)估計值fip、前一步導(dǎo)數(shù)校正值fi-1p、前兩步導(dǎo)數(shù)校正值fi-2c、前三步導(dǎo)數(shù)校正值fi-3c和狀態(tài)變量當(dāng)前步的校正值yic這五個參數(shù)的值傳遞給當(dāng)前變量的后一步，即將fi-2c的值傳遞給fi-3c，將fi-1c的值傳遞給fi-2c，將fic的值傳遞給fi-1c，將fi+1p的值傳遞給fip，將yic的值傳遞給yi-1c；時間t增加一個單位長度；返回步驟五；步驟八結(jié)束TBB任務(wù)調(diào)度；完成并行計算任務(wù)，結(jié)束TBB任務(wù)調(diào)度。2、根據(jù)權(quán)利要求1所述的一種基于TBB線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法，其特征在于步驟二中所述的編寫預(yù)估計算和校正計算程序，程序中預(yù)估計算和校正計算的部分是循環(huán)體結(jié)構(gòu)，循環(huán)體結(jié)構(gòu)中的代碼使用局部變量或模版類ParaODES的成員變量。3、根據(jù)權(quán)利要求1所述的一種基于TBB線程構(gòu)建塊的連續(xù)系統(tǒng)仿真多核并行方法，其特征在于步驟四所述的TBB任務(wù)調(diào)度初始化具體為TBB由任務(wù)調(diào)度器對象task_scheduler_init實現(xiàn)多任務(wù)的分配和并行計算，支持對多線程的劃分。全文摘要本發(fā)明公開了一種基于TBB的連續(xù)系統(tǒng)仿真多核并行方法，主要包括以下步驟步驟一，搭建TBB并行計算平臺；步驟二，構(gòu)造并行ParaODES模版類；步驟三，模型初始化；步驟四，TBB任務(wù)調(diào)度初始化；步驟五，仿真時間范圍判斷；步驟六，調(diào)用并行模版類ParaODES；步驟七，仿真結(jié)果傳遞；步驟八，結(jié)束TBB任務(wù)調(diào)度。本發(fā)明利用TBB實現(xiàn)連續(xù)系統(tǒng)仿真的并行化處理，提高了仿真運算的速度，從而滿足仿真的實時性要求，縮短了仿真研究周期；利用TBB進行并行化處理，過程簡單靈活，易于操作，具有很好的擴展性；TBB采用多核架構(gòu)，符合當(dāng)今并行技術(shù)向多核平臺發(fā)展的方向。文檔編號G06F9/46GK101639788SQ20091009240公開日2010年2月3日申請日期2009年9月10日優(yōu)先權(quán)日2009年9月10日發(fā)明者唐力勇,妮李申請人:北京航空航天大學(xué)