Hello Barry, Thanks a lot for explanation. I don't really need it. Sometimes during testing I output my system matrix to analyze it outside, but for that purposes I can easily switch to MPIAIJ.
Regards, Alexander On 18.12.2011 18:13, Barry Smith wrote: > Alexander, > > It is not hanging, it is just very very very slow. The code to do the > parallel binary saving of sbaij matrices is total crap. (While for AIJ it is > scalable code). The code will only run at reasonable speed for very very > small matrices. > > Do you really need this functionality? I ask because it is a moderately > big project to fix the code and in many years you are the only one who has > had this problem. > > Barry > > > On Dec 14, 2011, at 10:54 AM, Alexander Grayver wrote: > >> Hello dev-team, >> >> I use the same code to output MPIAIJ and MPISBAIJ matrices, however with >> latter program hangs somewhere in MatView. >> The code is: >> call PetscViewerBinaryOpen(comm,'out',FILE_MODE_WRITE,viewer,ierr); >> CHKERRQ(ierr) >> ... >> call PetscViewerFileSetName(viewer,'A.dat',ierr) >> call MatView(A,viewer,ierr); CHKERRQ(ierr) >> >> When I define A as a MPIAIJ matrix everything works well. >> Any ideas? >> >> Regards, >> Alexander
