This is why I have recompiled PETSc with spooles. spooles can deal with AX=Y(Y is a matrix). However, PETSc only provide the following: 51: PetscErrorCode MatSolve_MPISpooles(Mat A,Vec b,Vec x) I don't set b to a matrix even if I use 178: PetscErrorCode MatFactorNumeric_MPISpooles(Mat A,MatFactorInfo *info,Mat *F) for LU factorization
Could you give me some advice or examples? thanks a lot. Regards, Yujie On Wed, Feb 27, 2008 at 11:32 AM, Matthew Knepley <knepley at gmail.com> wrote: > On Wed, Feb 27, 2008 at 1:21 PM, Yujie <recrusader at gmail.com> wrote: > > Dear Matt: > > > > I checked the codes about MatMatSolve(). However, currently, PETSc > didn't > > realize its parallel version. Is it right? I want to inverse the matrix > > parallelly. could you give me some examples about it? thanks a lot. > > Thats right. The parallel version is not implemented. It looks like this > would > take significant work. > > Matt > > > Regards, > > Yujie > > > > > > > > On 2/27/08, Matthew Knepley <knepley at gmail.com> wrote: > > > On Wed, Feb 27, 2008 at 11:05 AM, Yujie <recrusader at gmail.com> wrote: > > > > Dear Sanjay: > > > > > > > > Thank you for your reply. I don't understand what you said. Now, I > want > > to > > > > use spooles package to inverse a sparse SPD matrix. I have further > > checked > > > > the inferface about spooles in PETSc. I find although spooles can > deal > > with > > > > AX=B (B may be a dense matrix) with parallel LU factorization. > > > > However, PETSc only provide the following: > > > > 51: PetscErrorCode MatSolve_MPISpooles(Mat A,Vec b,Vec x) > > > > I don't set b to a matrix even if I use > > > > 178: PetscErrorCode MatFactorNumeric_MPISpooles(Mat A,MatFactorInfo > > > > *info,Mat *F) for LU factorization. > > > > > > > > Could you have any suggestions about this? thanks a lot. > > > > > > > > > MatMatSolve() > > > > > > > > > http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatMatSolve.html > > > > > > Matt > > > > > > > > > > Regards, > > > > Yujie > > > > > > > > On 2/26/08, Sanjay Govindjee <sanjay at ce.berkeley.edu> wrote: > > > > > from my make file > > > > > > > > > > -@${MPIRUN} -s all -np $(NPROC) $(PROGNAME) -ksp_type preonly > > > > > -ksp_monitor -pc_type cholesky -mat_type mpisbaijspooles > -log_summary > > > > > -on_error_attach_debugger -mat_spooles_symmetryflag 0 > -options_left > > > > > > > > > > > > > > > -sg > > > > > > > > > > > > > > > Yujie wrote: > > > > > > Hi, everyone > > > > > > > > > > > > I have compiled PETSc with spooles. However, I try to find how > to > > use > > > > > > this package in PETSc directory. I can't find any examples for > it. > > > > > > Could you give me some advice? I want to use spooles to inverse > a > > > > > > sparse matrix. thanks a lot. > > > > > > > > > > > > Regards, > > > > > > Yujie > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > -- > > > 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/20080227/0c8d87bf/attachment.htm>
