Openfoam Pstream類探索

来源:https://www.cnblogs.com/heloveHe/archive/2023/02/23/17148332.html
-Advertisement-
Play Games

對於數值模擬而言,無論是商軟或者開源軟體,並行計算都是非常重要的, 作為一名模擬工程師,如果想把自身數值模擬能力提升一個層次,需要對並行計算有很好的理解與應用 openfoam並行通信主要通過Pstream類完成 Pstream類,類如其名,parallel_stream,並行計算時使用的信息流 O ...


對於數值模擬而言,無論是商軟或者開源軟體,並行計算都是非常重要的,
作為一名模擬工程師,如果想把自身數值模擬能力提升一個層次,需要對並行計算有很好的理解與應用


openfoam並行通信主要通過Pstream類完成

Pstream類,類如其名,parallel_stream,並行計算時使用的信息流
Openfoam對其的介紹是:

Inter-processor communications stream.

處理器間交換信息流

類似的命名方法我們在c++文件讀取時說過,std有fstream類讀取寫入文件/二進位文件,比如說我們要讀取文件,會把讀取內容放入緩存區內進行操作

#include <iostream>
#include <fstream>  // ifstream類需要包含的頭文件。
#include <string>     // getline()函數需要包含的頭文件。
using  namespace std;
 
int main()
{
    string filename = R"(./test.txt)";
 
    //ifstream fin(filename, ios::in);
    ifstream fin;
    fin.open(filename , ios::in);
 
    // 判斷打開文件是否成功。
    // 失敗的原因主要有:1)目錄不存在;2)文件不存在;3)沒有許可權,Linux平臺下很常見。
    if (fin.is_open() == false)
    {
        cout << "打開文件" << filename << "失敗。\n";  return 0;
    }
 
    string buffer;
    while (fin >> buffer)
    {
        cout << buffer << endl;
    }
 
    fin.close();	   // 關閉文件,fin對象失效前會自動調用close()。
 
    cout << "操作文件完成。\n";
}

類似的openfoam也有PstreamBuffers類進行並行通信緩衝
可以這樣使用:

PstreamBuffers pBuffers(Pstream::commsTypes::nonBlocking);

for (label proci = 0; proci < Pstream::nProcs(); proci++)
    { 
        if (proci != Pstream::myProcNo()) 
            {
                someObject vals;
                UOPstream str(proci, pBuffers);
                str << vals; 
            } 
    }

pBuffers.finishedSends(); // no-op for blocking

for (label proci = 0; proci < Pstream::nProcs(); proci++)
    { 
            if (proci != Pstream::myProcNo())
            { 
                UIPstream str(proci, pBuffers);
                someObject vals(str); 
            } 
    }

上面這個程式可以看到,先後使用UOPstream與UIPstream進行緩衝區的文件輸出與讀取,這就很像ofstream類與ifstream類,甚至命名方式上都有幾分相似,我們打開相應的繼承關係圖

image
image

二者分別服務於IPstream類與OPstream類,我們再打開今天文章的主角,Pstream類繼承關係圖
image
發現IPstream類與OPstream類是Pstream類的衍生類,Pstream類是其基礎
打開Pstream類的源碼:

點擊查看代碼
namespace Foam
{

/*---------------------------------------------------------------------------*\
                           Class Pstream Declaration
\*---------------------------------------------------------------------------*/

class Pstream
:
    public UPstream
{

protected:

    // Protected data

        //- Transfer buffer
        DynamicList<char> buf_;

public:

    // Declare name of the class and its debug switch
    ClassName("Pstream");


    // Constructors

        //- Construct given optional buffer size
        Pstream
        (
            const commsTypes commsType,
            const label bufSize = 0
        )
        :
            UPstream(commsType),
            buf_(0)
        {
            if (bufSize)
            {
                buf_.setCapacity(bufSize + 2*sizeof(scalar) + 1);
            }
        }


        // Gather and scatter

            //- Gather data. Apply bop to combine Value
            //  from different processors
            template<class T, class BinaryOp>
            static void gather
            (
                const List<commsStruct>& comms,
                T& Value,
                const BinaryOp& bop,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T, class BinaryOp>
            static void gather
            (
                T& Value,
                const BinaryOp& bop,
                const int tag = Pstream::msgType(),
                const label comm = Pstream::worldComm
            );

            //- Scatter data. Distribute without modification. Reverse of gather
            template<class T>
            static void scatter
            (
                const List<commsStruct>& comms,
                T& Value,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T>
            static void scatter
            (
                T& Value,
                const int tag = Pstream::msgType(),
                const label comm = Pstream::worldComm
            );

        // Combine variants. Inplace combine values from processors.
        // (Uses construct from Istream instead of <<)

            template<class T, class CombineOp>
            static void combineGather
            (
                const List<commsStruct>& comms,
                T& Value,
                const CombineOp& cop,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T, class CombineOp>
            static void combineGather
            (
                T& Value,
                const CombineOp& cop,
                const int tag = Pstream::msgType(),
                const label comm = Pstream::worldComm
            );

            //- Scatter data. Reverse of combineGather
            template<class T>
            static void combineScatter
            (
                const List<commsStruct>& comms,
                T& Value,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T>
            static void combineScatter
            (
                T& Value,
                const int tag = Pstream::msgType(),
                const label comm = Pstream::worldComm
            );

        // Combine variants working on whole List at a time.

            template<class T, class CombineOp>
            static void listCombineGather
            (
                const List<commsStruct>& comms,
                List<T>& Value,
                const CombineOp& cop,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T, class CombineOp>
            static void listCombineGather
            (
                List<T>& Value,
                const CombineOp& cop,
                const int tag = Pstream::msgType(),
                const label comm = Pstream::worldComm
            );

            //- Scatter data. Reverse of combineGather
            template<class T>
            static void listCombineScatter
            (
                const List<commsStruct>& comms,
                List<T>& Value,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T>
            static void listCombineScatter
            (
                List<T>& Value,
                const int tag = Pstream::msgType(),
                const label comm = Pstream::worldComm
            );

        // Combine variants working on whole map at a time. Container needs to
        // have iterators and find() defined.

            template<class Container, class CombineOp>
            static void mapCombineGather
            (
                const List<commsStruct>& comms,
                Container& Values,
                const CombineOp& cop,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class Container, class CombineOp>
            static void mapCombineGather
            (
                Container& Values,
                const CombineOp& cop,
                const int tag = Pstream::msgType(),
                const label comm = UPstream::worldComm
            );

            //- Scatter data. Reverse of combineGather
            template<class Container>
            static void mapCombineScatter
            (
                const List<commsStruct>& comms,
                Container& Values,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class Container>
            static void mapCombineScatter
            (
                Container& Values,
                const int tag = Pstream::msgType(),
                const label comm = UPstream::worldComm
            );



        // Gather/scatter keeping the individual processor data separate.
        // Values is a List of size UPstream::nProcs() where
        // Values[UPstream::myProcNo()] is the data for the current processor.

            //- Gather data but keep individual values separate
            template<class T>
            static void gatherList
            (
                const List<commsStruct>& comms,
                List<T>& Values,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T>
            static void gatherList
            (
                List<T>& Values,
                const int tag = Pstream::msgType(),
                const label comm = UPstream::worldComm
            );

            //- Scatter data. Reverse of gatherList
            template<class T>
            static void scatterList
            (
                const List<commsStruct>& comms,
                List<T>& Values,
                const int tag,
                const label comm
            );

            //- Like above but switches between linear/tree communication
            template<class T>
            static void scatterList
            (
                List<T>& Values,
                const int tag = Pstream::msgType(),
                const label comm = UPstream::worldComm
            );


        // Exchange

            //- Helper: exchange contiguous data. Sends sendData, receives into
            //  recvData. If block=true will wait for all transfers to finish.
            template<class Container, class T>
            static void exchange
            (
                const UList<Container>& sendData,
                const labelUList& recvSizes,
                List<Container>& recvData,
                const int tag = UPstream::msgType(),
                const label comm = UPstream::worldComm,
                const bool block = true
            );

            //- Helper: exchange sizes of sendData. sendData is the data per
            //  processor (in the communicator). Returns sizes of sendData
            //  on the sending processor.
            template<class Container>
            static void exchangeSizes
            (
                const Container& sendData,
                labelList& sizes,
                const label comm = UPstream::worldComm
            );

            //- Exchange contiguous data. Sends sendData, receives into
            //  recvData. Determines sizes to receive.
            //  If block=true will wait for all transfers to finish.
            template<class Container, class T>
            static void exchange
            (
                const UList<Container>& sendData,
                List<Container>& recvData,
                const int tag = UPstream::msgType(),
                const label comm = UPstream::worldComm,
                const bool block = true
            );
};

我們看到Pstream類有一個構造函數,剩下的都是靜態成員函數,而這些成員函數就是並行通訊的工具箱
這裡多問一句,為什麼工具箱的函數都是靜態成員函數


為什麼這裡用靜態成員函數呢

用靜態成員可以變數實現多個對象間的數據共用,比全局變數更安全
這裡我詳細說下,舉個例子

Time mytime1;
mytime1.hour=2;
Time mytime2;
mytime2.hour=4;

這段程式中成員變數是跟著對象走的,他們的對象各自占用不同的記憶體地址,彼此互不影響
那我們想做類內的全局變數滿足相互通信需求,在不同對象mytime1和mytime2中共用一個副本,怎麼辦
這時static關鍵字就派上用場了,增加了static關鍵字或成員函數不隸屬整個對象,而隸屬於整個類
因為這個變數跟著類走,所以調用時用“類名::成員變數名”或“類名::成員變數函數”進行調用(當然也可用“對象名.靜態函數名”),表示明確的隸屬關係,不創建對象也可進行訪問編輯
在Pstream類調用工具箱中函數時,我們常見到這樣的調用方式,而且不創建Pstream對象也可進行調用

// 在head節點收集信息
Pstream::gatherList(nInternalFaces);
Pstream::gatherList(nBoundaries);

因為類的靜態成員脫離了與對象的關係,普通成員變數的記憶體分配是在對象初始化時完成的,對於靜態成員必須在程式的全局區進行清晰的初始化
全局區的初始化過程可由某個.cpp源文件的開頭的靜態成員函數完成,如下所示:

void Time::func(int testValue)
{
	mystatic = testValue ;
}

或者在全局區這樣寫:

int Time::mystatic=10;

這樣能保證這個靜態成員變數能夠被正常使用。
此外靜態成員函數只能調用靜態成員變數,也沒有this指針可以使用
這裡上一張圖可能更方便理解
image

C++程式運行時,靜態變數和全局變數存儲在數據段,所以需要在全局區通過直接分配記憶體或者靜態函數進行分配記憶體
因而靜態成員的生命周期與程式運行周期相同,在程式中只有一份,無論創建對象與否,或者創建多少對象
說到這裡可能大家對Openfoam的並行通信多了一些理解,只要開始了並行計算那麼就可以通過Pstream類內的成員函數進行通信調用,在同樣的數據段副本上進行信息流溝通


接下來依次說下類中各個工具的使用

收發數據
Pstream::gather()與Pstream::scatter()分別有兩個重載,分別是收集以及散佈數據,不如後面Pstream::gatherList()與Pstream::scatterList()常用,這裡不細說了
Pstream::combineGather()、Pstream::combineScatter()重載情況與上同,用於就地集中收集或散佈的數據,不太常用

Pstream::listCombineGather()、Pstream::listCombineScatter()重載情況與上同,用於一次整合list容器中的變數
Pstream::mapCombineGather()、Pstream::mapCombineScatter()重載情況與上同,用於一次整合整個map容器中的變數

Pstream::gatherList()以及Pstream::scatterList()的第二個重載比較常用,

template<class T>
            static void gatherList
            (
                List<T>& Values,
                const int tag = Pstream::msgType(),
                const label comm = UPstream::worldComm
            );
template<class T>
            static void scatterList
            (
                List<T>& Values,
                const int tag = Pstream::msgType(),
                const label comm = UPstream::worldComm
            );

Pstream::gatherList()以及Pstream::scatterList()的輸入第一個參數是Values
這個Values需要自己整合下,Values是UPstream::nProcs()數量大小的List,比如說我要收集內部面可以這樣創建需要收集的List

List<label> nIternalFaces(Pstream::nProcs());
nIternalFaces[Pstream::myProcNo()] = mesh.Cf().size();//比如說看看每個節點分到了多少網格
Pstream::gatherList(nIternalFaces);//在頭結點收集數據

Pstream::scatterList()與之類似
Pstream::gatherList()以及Pstream::scatterList()的輸入第二個參數是Pstream::msgType(),預設為1,可以不輸入

int Foam::UPstream::msgType_(1);

Pstream::gatherList()以及Pstream::scatterList()的輸入第三個參數是Pstream::msgType(),預設為0,可以不輸入

Foam::label Foam::UPstream::worldComm(0);

交換數據
Pstream::exchange()有兩個重載,用於交換連續的數據,一般情況下等待其他所有傳輸完成再傳輸,可通過預設參數block()修改優先權
Pstream::exchangeSizes()用於交換數據的大小


下麵是Pstream類函數相互關係
image


結語

並行開發遠不止收發數據這麼簡單,還有很多類可說的,後續會一一進行介紹,並對openfoam並行計算進行優化

一起探索openfoam也是相當有趣的一件事,非常歡迎私信討論
指正的價值要比打賞更重要,下麵是個人聯繫方式,能結交到志同道合的朋友是我的榮幸
image


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 本文是系列第三篇。系列文章: 現代圖片性能優化及體驗優化指南 - 圖片類型及 Picture 標簽的使用 現代圖片性能優化及體驗優化指南 - 響應式圖片方案 圖片資源,在我們的業務中可謂是占據了非常大頭的一環,尤其是其對帶寬的消耗是十分巨大的。 對圖片的性能優化及體驗優化在今天就顯得尤為重要。本文, ...
  • 一、直接在 el-table-column 外嵌套 el-form 符合表單的校驗習慣,唯一需要註意的地方 el-form 需要綁定 :model="scope.row" // html <div id="app"> <el-table :data="list" border> <el-table- ...
  • 架構的核心是管理複雜度,架構師的核心能力是抽象能力,什麼是抽象能力?抽象能力就是一種化繁為簡的能力。何為化繁為簡?就是把一種複雜的事情變得簡單的能力,比如通過打比喻讓別人很容易聽明白你說的意思就是一種抽象能力。如何鍛煉抽象能力?我覺得有三種方法,第一種是用歸納法找共性,從多個問題中找到共同的問題提煉... ...
  • 1. Scala 1.1. 混合了面向對象和函數式編程的語言 1.2. 直接使用任何一個Java類庫 1.3. 聲明非遞歸的方法時,不需要顯式地返回類型 1.3.1. 會自動地替你推斷生成一個 1.4. 字元串插值 1.4.1. 在字元串的字面量中內嵌變數和表達式 1.5. 以使用val關鍵字替換v ...
  • 說起辦公自動化就離不開對excel表格的處理,現在哪個辦公族不得每天面對著幾份excel過日本子?本次活動同時在稀土掘金、博客園、微信公眾號舉辦,互不幹擾,歡迎各平臺搜索博主 ENG八戒 參與最新活動 ...
  • 自己實現Mybatis底層機制-01 主要實現:封裝SqlSession到執行器+Mapper介面和Mapper.xml+MapperBean+動態代理Mapper的方法 1.Mybatis整體架構分析 對上圖的解讀: 1)mybatis 的核心配置文件 ​ mybatis-config.xml:進 ...
  • VL21 根據狀態轉移表實現時序電路 寫一個簡單的Moore狀態機就可以了,太短就懶得寫三段式了。 `timescale 1ns/1ns module seq_circuit( input A , input clk , input rst_n, output wire Y ); reg [1:0] ...
  • 早期的計算工具 計算尺是基本的手動計算器,易於乘除 1614年,蘇格蘭數學家納皮爾發現利用加減計算乘除的方法,依此發明對數,納皮爾在製作第一張對數表的時候,必需進行大量的乘法運算,而一條物理線的距離或區間可表示真數,於是他設計出計算器納皮爾的骨頭協助計算[2]。到1633年,英國牧師奧特雷德利用對數 ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發(一):使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說,近來想嘗試開發一款安卓APP,考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP,畢竟是使用Visual Studio開發工具,使用起來也比較的順手,結合微軟官方的教程進行了安卓 ...
  • 前言 QuestPDF 是一個開源 .NET 庫,用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫,它徹底改變了我們生成 PDF 文檔的方 ...
  • 項目地址 項目後端地址: https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址: ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址: http://tree ...
  • 話不多說,直接開乾 一.下載 1.官方鏈接下載: https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe,運行開始下載SQL server; 二. ...
  • 前言 隨著物聯網(IoT)技術的迅猛發展,MQTT(消息隊列遙測傳輸)協議憑藉其輕量級和高效性,已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫,為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog支持多種接收器用於日誌存儲,增強器用於添加屬性,LogContext管理動態屬性,支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項,適用於不同需求。 ...
  • 目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據,爬取信息時需要模擬瀏覽器行為,否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下: 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • 1.前言 什麼是熱更新 游戲或者軟體更新時,無需重新下載客戶端進行安裝,而是在應用程式啟動的情況下,在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR,Xlua,ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles,Addressable, ...
  • 本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺,其實像阿裡雲,騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後,它會送 ...
  • 通過以下方式可以高效,並保證數據同步的可靠性 1.API設計 使用RESTful設計,確保API端點明確,並使用適當的HTTP方法(如POST用於創建,PUT用於更新)。 設計清晰的請求和響應模型,以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證,確保接收到的數據符合預期格 ...