7 Subroutine printglobalvec(GlobalVec,filename)
9 #include <finclude/petscsys.h> 10 #include <finclude/petscvec.h> 11 #include <finclude/petscvec.h90> 12 #include <finclude/petscis.h> 16 character(80) :: filename
22 call petscviewerasciiopen(petsc_comm_world,filename,viewer,ierr)
23 call vecview(globalvec,viewer,ierr)
25 call petscviewerdestroy(viewer,ierr)
27 End Subroutine printglobalvec
31 Subroutine printlocalvec(GlobalVec,da,filename)
33 #include <finclude/petscsys.h> 34 #include <finclude/petscvec.h> 35 #include <finclude/petscvec.h90> 36 #include <finclude/petscdmda.h> 37 #include <finclude/petscdmda.h90> 38 #include <finclude/petscis.h> 43 character(80) :: filename
48 call dmgetlocalvector(da,localvec,ierr)
49 call dmglobaltolocalbegin(da,globalvec,insert_values,localvec,ierr)
50 call dmglobaltolocalend(da,globalvec,insert_values,localvec,ierr)
51 call petscviewerasciiopen(petsc_comm_self,filename,viewer,ierr)
52 call vecview(localvec,viewer,ierr)
53 call dmrestorelocalvector(da,localvec,ierr)
56 End Subroutine printlocalvec