Problem with saving the binary matrix via Matview

[Matthew Knepley]

On Fri, Sep 4, 2009 at 2:25 PM, Nguyen, Hung V ERDC-ITL-MS <
Hung.V.Nguyen at usace.army.mil> 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).

  Matt


> 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],"Matrix.at%f",t_prev);
>      PetscViewerBinaryOpen(PETSC_COMM_WORLD,
> file[LOAD_MATRIX],FILE_MODE_WRITE,&fd);
>      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
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090904/c34a7b07/attachment.htm>