1c4762a1bSJed Brown 2c4762a1bSJed Brown static char help[] = "Tests the creation of a PC context.\n\n"; 3c4762a1bSJed Brown 4c4762a1bSJed Brown #include <petscpc.h> 5c4762a1bSJed Brown 6c4762a1bSJed Brown int main(int argc,char **args) 7c4762a1bSJed Brown { 8c4762a1bSJed Brown PC pc; 9c4762a1bSJed Brown PetscInt n = 5; 10c4762a1bSJed Brown Mat mat; 11c4762a1bSJed Brown 12*9566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&args,(char*)0,help)); 13*9566063dSJacob Faibussowitsch PetscCall(PCCreate(PETSC_COMM_WORLD,&pc)); 14*9566063dSJacob Faibussowitsch PetscCall(PCSetType(pc,PCNONE)); 15c4762a1bSJed Brown 16c4762a1bSJed Brown /* Vector and matrix must be set before calling PCSetUp */ 17*9566063dSJacob Faibussowitsch PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF,n,n,3,NULL,&mat)); 18*9566063dSJacob Faibussowitsch PetscCall(MatAssemblyBegin(mat,MAT_FINAL_ASSEMBLY)); 19*9566063dSJacob Faibussowitsch PetscCall(MatAssemblyEnd(mat,MAT_FINAL_ASSEMBLY)); 20*9566063dSJacob Faibussowitsch PetscCall(PCSetOperators(pc,mat,mat)); 21*9566063dSJacob Faibussowitsch PetscCall(PCSetUp(pc)); 22*9566063dSJacob Faibussowitsch PetscCall(MatDestroy(&mat)); 23*9566063dSJacob Faibussowitsch PetscCall(PCDestroy(&pc)); 24*9566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 25b122ec5aSJacob Faibussowitsch return 0; 26c4762a1bSJed Brown } 27c4762a1bSJed Brown 28c4762a1bSJed Brown /*TEST 29c4762a1bSJed Brown 30c4762a1bSJed Brown test: 31c4762a1bSJed Brown 32c4762a1bSJed Brown TEST*/ 33