multiple rhs

> > What solver would I use to do a factorization of a dense parallel matrix w/ > plapack? MAT_SOLVER_PLAPACK. See ~petsc-3.0.0/src/mat/examples/tests/ex103.c Hong > > I don't see a > > MPI_SOLVER_PLAPACK > > ? > > > > > > On Thu, 12 Mar 2009, Hong Zhang wrote: > >> Is MatCreateMP

multiple rhs

I'm getting plapack errors in "external library" with MatMatMult_MPIDense_MPIDense with plapack? How is memory handled for a matrix of type MATMPIDENSE? Are all NxN entries allocated and ready for use at time of creation? or do I have to MatInsertValues then Assemble to be ready to use a matri

multiple rhs

Hi Hong, What solver would I use to do a factorization of a dense parallel matrix w/ plapack? I don't see a MPI_SOLVER_PLAPACK ? On Thu, 12 Mar 2009, Hong Zhang wrote: > >>> >>> Is MatCreateMPIDense the recommended matrix type to interface w/ mumps ? >>> Does it use a sparse direct stor

parallel MatMatMult

MatMatMult(X,Y,...,...,Z) where X is MPIDENSE and Y is MPIAIJ seems to work but when X is MPIAIJ and Y is MPIDENSE doesn't ? says its not supported. there seems to be all permuations in the source MatMatMult_MPIAIJ_MPIDense MatMatMult_MPIDense_MPIAIJ MatMatMult_MPIAIJ_MPIAIJ MatMatMult_MPI

Problems with Socket

ranted 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/20090312/3f90ac77/attachment-0001.htm>

multiple rhs

>> >> Is MatCreateMPIDense the recommended matrix type to interface w/ mumps ? >> Does it use a sparse direct storage or allocate the full n x n matrix? > > > No, MUMPS is "sparse direct" so it uses MPIAIJ. For mpi dense matrix, you can use plapack Hong > > >> >> df >> >> On Thu, 12 Mar 2009, Ma

Problems with Socket

e. Do you help me? Thanks. Bye. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090312/cfa7bb95/attachment.htm>

multiple rhs

MatLUFactorSymbolic(Afact,A,isrow,iscol,&info); >>>>> MatLUFactorNumeric(Afact,A,&info); >>>>> MatMatSolve(Afact,B,X); >>>>> >>>>> >>>>> my solve keeps running out of memory >>>>> >>>>> "[0]PETSC ERROR: Memory requested xxx!" >>>>> >>>>> >>>>> is this in bytes? I can't tell if the problem I'm trying to solve >>>>> is too large form my machine or if I just have bug in the call >>>>> sequence. >>>>> >>>>> >>>>> >>>>> >>>>> thank you, >>>>> David Fuentes >>>>> >>>>> >>>>> >>>> >> >> -- >> 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 >> >> -- 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/20090312/8a438ba3/attachment.htm>

multiple rhs

Thanks Matt, Is MatCreateMPIDense the recommended matrix type to interface w/ mumps ? Does it use a sparse direct storage or allocate the full n x n matrix? df On Thu, 12 Mar 2009, Matthew Knepley wrote: > You can try using a sparse direct solver like MUMPS instead of PETSc LU. > > Matt

Problems with Socket

You have to upgrade to 3.0.0 to expect these to work. Barry On Mar 12, 2009, at 8:20 AM, Santolo Felaco wrote: > Hi I use Petsc 2.3.1 (I don't use the last version because I work on > a software developed with this version). > I tried to used the Petsc socket, but the client and server

Problems with Socket

> Do you help me? Thanks. >>> Bye. >>> >> -- >> 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 >> > > -- 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/20090312/a08e9943/attachment.htm>

multiple rhs

; >>> "[0]PETSC ERROR: Memory requested xxx!" >>> >>> >>> is this in bytes? I can't tell if the problem I'm trying to solve >>> is too large form my machine or if I just have bug in the call sequence. >>> >>> >>> >>> >>> thank you, >>> David Fuentes >>> >>> >> -- 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/20090312/cbf2390e/attachment.htm>

multiple rhs

Thanks Hong, The complete error message is attached. I think I just had too big of a matrix. The matrix i'm trying to factor is 327680 x 327680 [0]PETSC ERROR: - Error Message [0]PETSC ERROR: Out of memory. This could be due to allocating

multiple rhs

David, I do not see any problem with the calling sequence. The memory is determined in MatLUFactorSymbolic(). Does your code crashes within MatLUFactorSymbolic()? Please send us complete error message. Hong On Wed, 11 Mar 2009, David Fuentes wrote: > > Hello, > > I have a sparse matrix, A, wi

Problems with Socket

rbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090312/d47bf614/attachment.htm>