Problem with saving the binary matrix via Matview

Matthew Knepley knepley at
Fri Sep 4 14:31:48 CDT 2009

On Fri, Sep 4, 2009 at 2:25 PM, Nguyen, Hung V ERDC-ITL-MS <
Hung.V.Nguyen at> wrote:

> Hello,
> I have a problem with saving a matrix with the binary viewer
> PetscBinaryViewerOpen() for the matrix with nrow=853564, nnzeros=47191472.
> The application ran with 32 pes for more than 3 hours without writing any
> into file.
> However, the same executable ran for small size of application
> (nrows-12856,nnzeros=675744) and were able to write successfully binary
> files.

Much bigger matrices than this are routinely saved. I do not think this has
to do
with MatView() in particular, unless this is a very slow disk. You can check
with the
debugger where the code is currently occupied. This sounds to me more like a
deadlock (not calling the routine with every process).


> Note: using Totalview I found that it hang in the line of
> Matview(matrix->petsc, fd);
> Thank you,
> -hung
> --- code:
>      petsc_analyst_mat(matrix->petsc);
>      sprintf(file[LOAD_MATRIX],"",t_prev);
>      PetscViewerBinaryOpen(PETSC_COMM_WORLD,
>      MatView(matrix->petsc, fd);
>      PetscViewerDestroy(fd);
> ---- Matrix info using salsa/AnaMod module:
> Computed <simple:trace> as <9.362113e+03>
> Computed <simple:trace-abs> as <9.362623e+03>
> Computed <simple:norm1> as <2.800000e+01>
> Computed <simple:normInf> as <2.800000e+01>
> Computed <simple:normF> as <2.404511e+02>
> Computed <simple:diagonal-dominance> as <-3.064463e-02>
> Could not compute <simple:symmetry-snorm>
> Could not compute <simple:symmetry-anorm>
> Could not compute <simple:symmetry-fsnorm>
> Could not compute <simple:symmetry-fanorm>
> Could not compute <structure:n-struct-unsymm>
> Computed <structure:nrows> as <853564>
> Could not compute <structure:symmetry>
> Computed <structure:nnzeros> as <47191472>
> Computed <structure:max-nnzeros-per-row> as <112>
> Computed <structure:min-nnzeros-per-row> as <16>

What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <>

More information about the petsc-users mailing list