[petsc-dev] Output MPISBAIJ

[Barry Smith]

  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