1 2 static char help[] = "Tests MatILUFactorSymbolic() on matrix with missing diagonal.\n\n"; 3 4 #include <petscmat.h> 5 #include <petscpc.h> 6 7 int main(int argc,char **args) 8 { 9 Mat C,A; 10 PetscInt i,j; 11 PetscScalar v; 12 PC pc; 13 Vec xtmp; 14 15 CHKERRQ(PetscInitialize(&argc,&args,(char*)0,help)); 16 CHKERRQ(MatCreate(PETSC_COMM_WORLD,&C)); 17 CHKERRQ(MatSetSizes(C,PETSC_DECIDE,PETSC_DECIDE,3,3)); 18 CHKERRQ(MatSetFromOptions(C)); 19 CHKERRQ(MatSetUp(C)); 20 CHKERRQ(VecCreateSeq(PETSC_COMM_WORLD,3,&xtmp)); 21 i = 0; j = 0; v = 4; 22 CHKERRQ(MatSetValues(C,1,&i,1,&j,&v,INSERT_VALUES)); 23 i = 0; j = 2; v = 1; 24 CHKERRQ(MatSetValues(C,1,&i,1,&j,&v,INSERT_VALUES)); 25 i = 1; j = 0; v = 1; 26 CHKERRQ(MatSetValues(C,1,&i,1,&j,&v,INSERT_VALUES)); 27 i = 1; j = 1; v = 4; 28 CHKERRQ(MatSetValues(C,1,&i,1,&j,&v,INSERT_VALUES)); 29 i = 2; j = 1; v = 1; 30 CHKERRQ(MatSetValues(C,1,&i,1,&j,&v,INSERT_VALUES)); 31 32 CHKERRQ(MatAssemblyBegin(C,MAT_FINAL_ASSEMBLY)); 33 CHKERRQ(MatAssemblyEnd(C,MAT_FINAL_ASSEMBLY)); 34 35 CHKERRQ(MatView(C,PETSC_VIEWER_STDOUT_WORLD)); 36 CHKERRQ(PCCreate(PETSC_COMM_WORLD,&pc)); 37 CHKERRQ(PCSetFromOptions(pc)); 38 CHKERRQ(PCSetOperators(pc,C,C)); 39 CHKERRQ(PCSetUp(pc)); 40 CHKERRQ(PCFactorGetMatrix(pc,&A)); 41 CHKERRQ(MatView(A,PETSC_VIEWER_STDOUT_WORLD)); 42 43 CHKERRQ(PCDestroy(&pc)); 44 CHKERRQ(VecDestroy(&xtmp)); 45 CHKERRQ(MatDestroy(&C)); 46 47 CHKERRQ(PetscFinalize()); 48 return 0; 49 } 50