any examples to demonstrate how to Spooles package?

Yujie recrusader at gmail.com
Wed Feb 27 13:40:52 CST 2008


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>


More information about the petsc-users mailing list