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
