[petsc-dev] Output MPISBAIJ

Barry Smith bsmith at mcs.anl.gov
Sun Dec 18 11:13:22 CST 2011


     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.


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