any examples to demonstrate how to Spooles package?

Matthew Knepley knepley at gmail.com
Wed Feb 27 13:12:54 CST 2008


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




More information about the petsc-users mailing list