17a3b9f03SLawrence Mitchell 27a3b9f03SLawrence Mitchell static char help[] = "Tests PCView() before PCSetup() with -pc_type lu.\n\n"; 37a3b9f03SLawrence Mitchell 47a3b9f03SLawrence Mitchell #include <petscmat.h> 57a3b9f03SLawrence Mitchell #include <petscpc.h> 67a3b9f03SLawrence Mitchell 77a3b9f03SLawrence Mitchell int main(int argc,char **args) 87a3b9f03SLawrence Mitchell { 97a3b9f03SLawrence Mitchell Mat A; 107a3b9f03SLawrence Mitchell PC pc; 117a3b9f03SLawrence Mitchell 12*9566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&args,(char*)0,help)); 13*9566063dSJacob Faibussowitsch PetscCall(MatCreate(PETSC_COMM_WORLD,&A)); 14*9566063dSJacob Faibussowitsch PetscCall(MatSetSizes(A,1,1,1,1)); 15*9566063dSJacob Faibussowitsch PetscCall(MatSetFromOptions(A)); 16*9566063dSJacob Faibussowitsch PetscCall(MatSetUp(A)); 17*9566063dSJacob Faibussowitsch PetscCall(MatSetValue(A,0,0,1,INSERT_VALUES)); 18*9566063dSJacob Faibussowitsch PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY)); 19*9566063dSJacob Faibussowitsch PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY)); 20*9566063dSJacob Faibussowitsch PetscCall(MatView(A,PETSC_VIEWER_STDOUT_WORLD)); 21*9566063dSJacob Faibussowitsch PetscCall(PCCreate(PETSC_COMM_WORLD,&pc)); 22*9566063dSJacob Faibussowitsch PetscCall(PCSetOperators(pc,A,A)); 23*9566063dSJacob Faibussowitsch PetscCall(PCSetType(pc,PCLU)); 24*9566063dSJacob Faibussowitsch PetscCall(PCView(pc,PETSC_VIEWER_STDOUT_WORLD)); 25*9566063dSJacob Faibussowitsch PetscCall(PCDestroy(&pc)); 26*9566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A)); 27*9566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 28b122ec5aSJacob Faibussowitsch return 0; 297a3b9f03SLawrence Mitchell } 307a3b9f03SLawrence Mitchell 317a3b9f03SLawrence Mitchell /*TEST 327a3b9f03SLawrence Mitchell 337a3b9f03SLawrence Mitchell test: 347a3b9f03SLawrence Mitchell 357a3b9f03SLawrence Mitchell TEST*/ 36