any examples to demonstrate how to Spooles package?

[Matthew Knepley]

On Wed, Feb 27, 2008 at 1:26 PM, Aron Ahmadia <aja2111 at columbia.edu> wrote:
> Hey Matt,
>
>  You should probably clean up the documentation for MatMatSolve while
>  you're at it, it's indicating that x and b are vectors...  Also,
>  should you reference the factor routine you need to use to get a
>  factored matrix?

The dev has the correct args. I added links to LU and Cholesky.

  Matt

>  ~A
>
>
>
>  On Wed, Feb 27, 2008 at 2:12 PM, 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