PETSc工具箱(簡介、安裝、程序設(shè)計)

基于PETSc的有限元并行程序設(shè)計可移植可擴(kuò)展科學(xué)計算工具箱PETSc（Portable,ExtensibleToolkitforScientificComputation）并行有限元求解程序背景：1、計算規(guī)模日益增大而計算資源不足！2、編寫并行的有限元程序難度很大！1、PETSc工具箱簡介3、基于PETSc的有限元求解程序設(shè)計2、PETSc程序示例主要內(nèi)容1、PETSc工具箱簡介：工具箱結(jié)構(gòu)MatVecBlasLapackMPIKSPSNESTS求解器組件底層用戶開發(fā)平臺數(shù)據(jù)結(jié)構(gòu)組件1、PETSc是系列數(shù)據(jù)結(jié)構(gòu)和功能子程序的集合。2、PETSc提供了高層的程序開發(fā)平臺!分布式的數(shù)據(jù)結(jié)構(gòu)高效的并行求解器矩陣Mat向量Vec非線性求解器SNES時間步進(jìn)求解器TS線性求解器KSP1、PETSc工具箱簡介1、PETSc工具箱簡介：工具箱結(jié)構(gòu)1、PETSc工具箱簡介：下載與編譯主頁：/petsc/petsc-as/Download：源碼（petsc-3.4.4.tar.gz）可調(diào)用的庫文件（包括頭文件）編譯cygwin,winfe,cl...用戶編程PETSc功能

調(diào)用將庫文件（.lib文件）和頭文件添加到用戶的配置中，主要有以下內(nèi)容：1、MPI的頭文件和庫文件MPICH2_DIR/includeMPICH2_DIR/lib（mpi.lib）2、線性代數(shù)運(yùn)算庫（Blas&Lapack）(libf2cblas.lib、libf2clapack.lib)3、PETSc的庫文件和頭文件petsc-3.2-p7\includepetsc-3.2-p7\arch-mswin-c-debug\includepetsc-3.2-p7\arch-mswin-c-debug\lib(libpetsc.lib)注：

由于鏈接庫沖突，通常需要禁用LIBCMT.lib。1、PETSc工具箱簡介：一種Windows下PETSc靜態(tài)庫調(diào)用2、PETSc程序示例：運(yùn)行1、SolveLinearSystem.cpp程序PetscInitialize(&argc,&args,(char*)0,help);功能類似于MPI_Init(&argc,&argv)形成PETSc的數(shù)據(jù)結(jié)構(gòu)（Mat&Vec）調(diào)用PETSc解法器求解（KSP/SNES/TS）PetscFinalize(void);功能類似于MPI_Finalize()并行的數(shù)據(jù)結(jié)構(gòu)（Mat&Vec）1、省去的存儲的細(xì)節(jié)，尤其是在分布式計算中。2、提供了運(yùn)算函數(shù)，如MatMult(A,x,y)。大量的成熟的方程組并行解法，包括直接法，迭代法和預(yù)條件子！2、PETSc程序示例：程序結(jié)構(gòu)1、聲明。2、創(chuàng)建。3、設(shè)置（賦值）。4、使用。5、釋放。矩陣Mat聲明：MatA;創(chuàng)建：MatCreate(A…);設(shè)置(賦值)：MatSet…(A…);聚集：VecAssembly…(A…);使用：如KSP調(diào)用;釋放：MatDestroy(&A);向量Vec聲明：Vecx;創(chuàng)建：VecCreate(x…);設(shè)置(賦值)：VecSet…(x…);聚集：VecAssembly…(x…);使用：如KSP調(diào)用;釋放：VecDestroy(&x);線性解法器KSP聲明：KSPksp;創(chuàng)建：KSPCreate(ksp…);設(shè)置：KSPSet…(ksp…);使用：KSPSolve(ksp…);釋放：KSPDestroy(&ksp);2、PETSc程序示例：組件的使用程序1、”ierr=“和CHKERRQ(ierr)函數(shù);

petsc的錯誤跟蹤，記錄下可能出錯的函數(shù)和錯誤原因。2、MatSetValue&VecSetValue(MatSetValues&VecSetValues)MatSetValue(A,i,j,v,ADD_VALUES);VecSetValue(x,i,v,ADD_VALUES);

為矩陣(向量)賦值。賦值模式：INSERT_VALUE替換

ADD_VALUDE疊加

并行程序中，只要處理器指定其整體行列號即可完成操作，而不必關(guān)心該元素是否在本處理器上。2、PETSc程序示例：幾個PETSc函數(shù)在編譯成可執(zhí)行文件后（**.exe）后，仍能改變程序中的設(shè)置。（程序中的預(yù)置）如：-mat_view將矩陣打印于屏幕(aij格式)

-mat_view_info將矩陣信息打印于屏幕

-ksp_view

顯示迭代法求解器的信息-ksp_typefgmres將迭代方法改為fgmres2、PETSc程序示例：運(yùn)行時選項（Runningoption）基本所有函數(shù)的解釋目錄：petsc-3.3-p5\docs\manualpagesPETSc組件的使用示例目錄：petsc-3.3-p5\src如KSP組件：petsc-3.3-p5\src\ksp\ksp\examples\testsPETSc用戶文檔目錄：petsc-3.3-p5\src\docs\manual.pdf2、PETSc程序示例：幫助文檔線彈性有限元的求解部分歸結(jié)為解線性方程組：[K]{δ}={P}3、基于PETSc的有限元求解程序設(shè)計：一種思路矩陣[A]和向量{p}與{δ}求解方程組Ax=bPETSc的數(shù)據(jù)結(jié)構(gòu)：Mat———矩陣Vec———向量PETSc的求解器組件KSP———解線性方程組1、線彈性有限元程序的一般步驟：前處理建模，形成網(wǎng)格信息包括：單元信息節(jié)點(diǎn)信息材料信息約束信息荷載信息單元分析整體分析組裝總剛矩陣[K]for（ie=0；ie<NE;ie++）{單剛計算；將單剛疊加至總剛}邊界條件修改方程“置大數(shù)法”“劃行劃列”求解支配方程組:[K]{δ}={p}大型線性方程組后處理輸出：位移云圖應(yīng)力云圖...3、基于PETSc的有限元求解程序設(shè)計：有限元流程前處理單元分析整體分析組裝總剛矩陣[K]邊界條件修改方程求解支配方程組:[K]{δ}={p}后處理進(jìn)程1進(jìn)程0進(jìn)程0求解支配方程組:[K]{δ}={p}求解支配方程組:[K]{δ}={p}單元分析整體分析組裝總剛矩陣[K]單元分析整體分析組裝總剛矩陣[K]3、基于PETSc的有限元求解程序設(shè)計：一種并行方法1、一種線彈性有限元程序并行化方法前處理單元分析整體分析組裝總剛矩陣[K]邊界處理求解支配方程組:[K]{u}={p}后處理1、提供成熟高效的方程組求解程序2、存儲總剛的矩陣結(jié)構(gòu)Mat和存儲列陣的向量結(jié)構(gòu)Vec分布式計算環(huán)境下！3、基于PETSc的有限元求解程序設(shè)計：一種并行方法網(wǎng)格信息PETSc求解器求解（KSP）輸出位移形成方程組的系數(shù)矩陣[K]和右端項{P}（PETSc）邊界處理修改方程3、基于PETSc的有限元求解程序設(shè)計：流程類似于數(shù)組空間動態(tài)分配程序一、輸入網(wǎng)格信息（代替前處理）節(jié)點(diǎn)信息NLIST.DAT單元信息ELIST.DAT材料信息NLIST.DAT受約束信息NLIST.DAT節(jié)點(diǎn)力信息NLIST.DAT3、基于PETSc的有限元求解程序設(shè)計：輸入ZFSolver：3、基于PETSc的有限元求解程序設(shè)計：分析主要過程分析：網(wǎng)格信息PETSc求解器求解（KSP）輸出位移形成方程組的系數(shù)矩陣[K]和右端項{P}（PETSc）邊界處理修改方程1、總剛矩陣[K]存儲預(yù)分配2、單剛計算及組集總剛（并行環(huán)境下）3、總剛方程求解（迭代方法）4、邊界修改（置“大數(shù)”法）1、矩陣的種類（Dense&Sparse）2、稀疏矩陣的預(yù)分配空間（類似于數(shù)組存儲空間的動態(tài)分配）。對于大型稀疏矩陣的聚集效率至關(guān)重要。3、提供預(yù)分配參數(shù)MatMPIAIJSetPreallocation(MatB,PetscIntd_nz,constPetscIntd_nnz[],PetscInto_nz,constPetscInto_nnz[])4、做法：對角矩陣和非矩陣對角的劃分存儲預(yù)分配參數(shù)：若每行非零元素個數(shù)基本一樣，給出兩個常數(shù)（d_nz，o_nz）若精確存儲預(yù)分配，給出兩個整形數(shù)組（d_nnz[]，o_nnz[]）一、預(yù)分配矩陣存儲3、基于PETSc的有限元求解程序設(shè)計：過程分析

STLVector容器每個處理器均遍歷網(wǎng)格，按照節(jié)點(diǎn)號的劃分統(tǒng)計出對應(yīng)部分的共單元的節(jié)點(diǎn)，并相應(yīng)的作自由度擴(kuò)展。并行矩陣的預(yù)分配：計算出每行對應(yīng)部分非零元素數(shù)目。在有限元中總剛矩陣中，只有共單元的才會對其有剛度貢獻(xiàn)。一、預(yù)分配矩陣存儲有限元總剛矩陣：高度稀疏、對稱，需要預(yù)分配找出節(jié)點(diǎn)的共單元關(guān)系3、基于PETSc的有限元求解程序設(shè)計：過程分析1、找出所有節(jié)點(diǎn)的共單元關(guān)系

for(intie=0;ie<NE;ie++)//所以單元循環(huán)

{

for（intid=0;id<4;id++）//處理與節(jié)點(diǎn)id共單元ie的節(jié)點(diǎn)

{

for(intjd=0;jd<4;jd++)

{………}

}

}

2、去重，排序、統(tǒng)計

3、自由度擴(kuò)展

4、得到預(yù)分配需要數(shù)組

以平面三個四邊形單元為例。二、單剛計算及組集總剛：每個單元循環(huán){Ele3D8N(EleCoor,iE,iMu,iDense,EK,RE);

//計算單剛Ke；MatSetValues(K,8,id,8,jd,EK,ADD_VALUES);

//Ke疊加至總剛；

}25/共33頁三、總剛方程的求解對于已經(jīng)賦值完畢的矩陣[K]和向量[p]就可以調(diào)用PETSc方程求解器進(jìn)行求解。26/共33頁置“大數(shù)