-pc_type mg

2009-03-14 Thread (Rebecca) Xuefei YUAN
Dear all, Is there any examples that use mg as a preconditioner? How could I add my own restriction/prolongation matrix? Thanks, Rebecca

Performance of MatMatSolve

2009-03-14 Thread Hong Zhang
David, Yes, MatMatSolve dominates. Can you also send us the output of '-log_summary' from superlu_dist? MUMPS only suppports centralized rhs vector b. Thus, we must scatter petsc distributed b into a seqential rhs vector (stored in root proc) in the petsc interface, which explains why the root

Performance of MatMatSolve

2009-03-14 Thread David Fuentes
Thanks a lot Hong, The switch definitely seemed to balance the load during the SuperLU matmatsolve. Although I'm not completely sure what I'm seeing. Changing the #dof also seemed to affect the load balance of the Mumps MatMatSolve. I need to investigate a bit more. Looking in the profile. The ma

multiple rhs

2009-03-14 Thread Hong Zhang
> Very Many Thanks for your efforts on this Barry. The PLAPACK > website looks like it hasn't been updated since 2007. Maybe PLAPACK > is in need of some maintenance? You said "nonsquare", is plapack working for > you for square matrices ? Yes, it works for square matrices. See ~petsc/src/mat/ex

multiple rhs

2009-03-14 Thread David Fuentes
Very Many Thanks for your efforts on this Barry. The PLAPACK website looks like it hasn't been updated since 2007. Maybe PLAPACK is in need of some maintenance? You said "nonsquare", is plapack working for you for square matrices ? thanks again, df >> >> >> >> [0]PETSC ERROR: -