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