static char help[] = "Tests PCView() before PCSetup() with -pc_type lu.\n\n"; #include #include int main(int argc,char **args) { Mat A; PC pc; PetscErrorCode ierr; ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr; CHKERRQ(MatCreate(PETSC_COMM_WORLD,&A)); CHKERRQ(MatSetSizes(A,1,1,1,1)); CHKERRQ(MatSetFromOptions(A)); CHKERRQ(MatSetUp(A)); CHKERRQ(MatSetValue(A,0,0,1,INSERT_VALUES)); CHKERRQ(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY)); CHKERRQ(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY)); CHKERRQ(MatView(A,PETSC_VIEWER_STDOUT_WORLD)); CHKERRQ(PCCreate(PETSC_COMM_WORLD,&pc)); CHKERRQ(PCSetOperators(pc,A,A)); CHKERRQ(PCSetType(pc,PCLU)); CHKERRQ(PCView(pc,PETSC_VIEWER_STDOUT_WORLD)); CHKERRQ(PCDestroy(&pc)); CHKERRQ(MatDestroy(&A)); ierr = PetscFinalize(); return ierr; } /*TEST test: TEST*/