[petsc-dev] Output MPISBAIJ

Alexander Grayver agrayver at gfz-potsdam.de
Sun Dec 18 11:18:27 CST 2011

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 


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

More information about the petsc-dev mailing list