我正在学习在c/C++中使用petsc以便在其中使用稀疏矩阵分解函数(MatICCFactor)。我做的第一步是将玩具矩阵读入petsc中以创建一个稀疏矩阵并在我做任何事情之前查看矩阵。我是一个用c/C++编写的BTW初学者。 这里是我使用的代码:使用Petsc库构造具有向量的稀疏矩阵
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?