62
62
ierr = PetscPrintf(PETSC_COMM_SELF,"Reading matrix completes.\n");CHKERRQ(ierr);
63
63
ierr = PetscOptionsGetString(PETSC_NULL,"-fout",fileout,127,PETSC_NULL);CHKERRQ(ierr);
64
ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,fileout,PETSC_BINARY_CREATE,&view);CHKERRQ(ierr);
64
ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,fileout,PETSC_FILE_CREATE,&view);CHKERRQ(ierr);
65
65
ierr = MatView(A,view);CHKERRQ(ierr);
66
66
ierr = VecView(b,view);CHKERRQ(ierr);
67
67
ierr = PetscViewerDestroy(view);CHKERRQ(ierr);