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 PetscErrorCode ierr; 12 13 ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr; 14 CHKERRQ(MatCreate(PETSC_COMM_WORLD,&A)); 15 CHKERRQ(MatSetSizes(A,1,1,1,1)); 16 CHKERRQ(MatSetFromOptions(A)); 17 CHKERRQ(MatSetUp(A)); 18 CHKERRQ(MatSetValue(A,0,0,1,INSERT_VALUES)); 19 CHKERRQ(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY)); 20 CHKERRQ(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY)); 21 CHKERRQ(MatView(A,PETSC_VIEWER_STDOUT_WORLD)); 22 CHKERRQ(PCCreate(PETSC_COMM_WORLD,&pc)); 23 CHKERRQ(PCSetOperators(pc,A,A)); 24 CHKERRQ(PCSetType(pc,PCLU)); 25 CHKERRQ(PCView(pc,PETSC_VIEWER_STDOUT_WORLD)); 26 CHKERRQ(PCDestroy(&pc)); 27 CHKERRQ(MatDestroy(&A)); 28 ierr = PetscFinalize(); 29 return ierr; 30 } 31 32 /*TEST 33 34 test: 35 36 TEST*/ 37