使用Petsc库构造具有向量的稀疏矩阵

[weixin_38065208]

我正在学习在c/C++中使用petsc以便在其中使用稀疏矩阵分解函数（MatICCFactor）。我做的第一步是将玩具矩阵读入petsc中以创建一个稀疏矩阵并在我做任何事情之前查看矩阵。我是一个用c/C++编写的BTW初学者。 这里是我使用的代码： static char help[] = "Tests Petscmat functions.\n\n"; #include<petscmat.h> #include <petscsys.h> int main(int argc,char **args) { PetscInitialize(&argc,&args,(char*)0, help); Mat *A; PetscInt m = 3; PetscInt n = 3; //to build sparse matrix(CSR) /* * 1 0 0 * 2 0 3 * 4 5 6 */ PetscInt Iv[4]={0,1,3,6}; //i vec size nrow+1 PetscInt Jv[6]={0,0,2,0,1,2}; //j vec size nnz PetscScalar vv[6]={1,2,3,4,5,6}; //v vec size nnz PetscInt nzrv[3]={1,2,3}; // nnz per row PetscViewer viewer; //MatSeqAIJSetPreallocation(A,n,nzrv);//If nnz is given then nz is ignored MatCreateSeqAIJWithArrays(PETSC_COMM_WORLD, m, n, Iv,Jv,vv, A); MatView(A,viewer); PetscViewerDestroy(&viewer); MatDestroy(&A); PetscFinalize(); return 0; } 它编译，但会导致违反分割运行时。我发现MatView()是分段错误的根源。任何人都可以帮助我如何正确使用MatView？

[weixin_38068996]

如果要在控制台中运行代码时查看矩阵，可以将MatView(Mat, PETSC_VIEWER_STDOUT_SELF)用于顺序矩阵，或将MatView(Mat, PETSC_VIEWER_STDOUT_WORLD)用于并行矩阵。 如果你想使用自己的浏览器背景下，程序应该是这样的： PetscViewer viewer; char matrixOutputFile[PETSC_MAX_PATH_LEN];//output file name PetscViewerBinaryOpen(PETSC_COMM_WORLD,matrixOutputFile,FILE_MODE_WRITE,&viewer); PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_INFO_DETAIL); MatView(A,viewer); PetscViewerDestroy(&viewer); 以上代码将矩阵存储到一个二进制文件的PetscViewerBinaryOpen。而且，PetscViewerASCIIOpen()可用于存储ASCII格式的矩阵，PetscViewerDrawOpen()可用于在X窗口显示中显示矩阵，并且PetscViewerSocketOpen()可以将矩阵输出到Socket查看器。