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>
