PETSc并行程序设计方法.ppt

#include petscsles.h #undef __FUNCT__ #define __FUNCT__ main int main(int argc,char **args) { Vec x,b,u; /* 近似解，右端向量和分析解 */ Mat A; /* 线性算子/矩阵 */ SLES sles; /* 线性解法器 */ PetscRandom rctx; /* 随机数发生器环境/数据对象 */ PetscReal norm; /* 解误差的范数 */ int i, j, I, J, Istart, Iend, ierr, m = 8, n = 7, its; PetscTruth flg; PetscScalar v, one = 1.0, neg_one = -1.0; KSP ksp; PetscInitialize(argc, args, (char *)0, help); ierr =PetscOptionsGetInt(PETSC_NULL,-m,m,PETSC_NULL); CHKERRQ(ierr); ierr =PetscOptionsGetInt(PETSC_NULL,-n,n,PETSC_NULL); CHKERRQ(ierr); /* - - - - - - - - - - - - - - - - -- - - - - - - - - - -- - - - - - - - - - -- - - - - - - - - - -- - - - - - - - - - - - - - - - - - 填充矩阵算子和右端向量 - - - - - - - - - - - - - - - - - - - - - - - - - - - - -- - - - - - - - - - - - - - - - - - - -- - - - - - - - - - - - - - - - - -- */ /* 创建矩阵对象 */ ierr=MatCreate(PETSC_COMM_WORLD,PETSC_DECIDE,PETSC_DECIDE,m*n,m*n,A); CHKERRQ(ierr); ierr=MatSetFromOptions(A); CHKERRQ(ierr); /*获得局部划分的上下界*/ ierr = MatGetOwnershipRange(A,Istart,Iend); CHKERRQ(ierr); /*给矩阵的每个元素赋值/五点差分格式*/ for (I=Istart; IIend; I++) { MGSetLevels：设置多重网格预条件子的优化级别 MGSetType：设置多重网格预条件子的类型 MGSetCycles：设置W-cycle参数 MGSetNumberSmoothUp：设置向前平滑参数 MGSetNumberSmoothDown：设置向后平滑参数 MGGetCoarseSolve：设置粗网格求解器 MGGetSmoother：设置平滑算子 MGGetSmootherUp：设置向前平滑算子 MGGetSmootherDown：设置向后平滑算子 MGSetInterpolate：设置插值算子 MGSetRestriction：设置约束条件 MGSetResidual：设置多重网格预条件子的残差 MGSetRhs：设置多重网格预条件子的右边向量 MGSetX：设置多重网格预条件子的解向量 MGSetR：设置多重网格预条件子的残差 4.1.4 奇异方程求解 MatNullSpaceCreate：创建一个零子空间对象 MatNullSpaceCreate ( MPI_Comm, PetscTruth hasconstants, int dim, Vec ?basis, MatNullSpace ?nsp ) PCNullSpaceAttach：设置一个零子空间对象 PCNullSpaceAttach ( PC pc, MatNullSpace nsp ) 4.2 非线性方程解法器（SNES） SNES 非线性解法器基于 牛顿迭代法（线性搜索和信赖域方法），并依赖线性解法器 SLES 实现。雅可比矩阵的求解是SNES解法器的重要组成部分。SNES求解以下形式的非线性方程组组 F( x) = 0 (1. 4.2) 其中解算子 F: Rn ? Rn。 4.2.1 基本用法 SNESCreate：创建一个非线性方程