[petsc-dev] Output MPISBAIJ

2011-12-18 Thread 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

[petsc-dev] Output MPISBAIJ

2011-12-14 Thread Alexander Grayver
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

[petsc-dev] Output MPISBAIJ

2011-12-14 Thread Matthew Knepley
On Wed, Dec 14, 2011 at 10:54 AM, Alexander Grayver agrayver at gfz-potsdam.de 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