1 2 static char help[] = "Tests PCView() before PCSetup() with -pc_type lu.\n\n"; 3 4 #include <petscmat.h> 5 #include <petscpc.h> 6 7 int main(int argc,char **args) 8 { 9 Mat A; 10 PC pc; 11 12 PetscCall(PetscInitialize(&argc,&args,(char*)0,help)); 13 PetscCall(MatCreate(PETSC_COMM_WORLD,&A)); 14 PetscCall(MatSetSizes(A,1,1,1,1)); 15 PetscCall(MatSetFromOptions(A)); 16 PetscCall(MatSetUp(A)); 17 PetscCall(MatSetValue(A,0,0,1,INSERT_VALUES)); 18 PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY)); 19 PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY)); 20 PetscCall(MatView(A,PETSC_VIEWER_STDOUT_WORLD)); 21 PetscCall(PCCreate(PETSC_COMM_WORLD,&pc)); 22 PetscCall(PCSetOperators(pc,A,A)); 23 PetscCall(PCSetType(pc,PCLU)); 24 PetscCall(PCView(pc,PETSC_VIEWER_STDOUT_WORLD)); 25 PetscCall(PCDestroy(&pc)); 26 PetscCall(MatDestroy(&A)); 27 PetscCall(PetscFinalize()); 28 return 0; 29 } 30 31 /*TEST 32 33 test: 34 35 TEST*/ 36