Further question about PC with Jaocbi Row Sum

[Shi Jin]

Thank you.
Suppose I have a diagonal matrix, what is the best way to invert it in PETSc?
Do I have to install the external packages superlu_dist or mumps?
I realized that LU or Cholesky decomposition does not work with MPIAIJ matrices.
I also know the best way is probably to directly call Vector operations directly.
However, I want to keep the same KSPSolve structure so that the same code can be used  for non-diagonal MPIAIJ matrices without changing each call to KSPSolve.
Thank you very much.

Shi
> Then you may try direct sparse linear solver,
> sequential run:
> -ksp_type preonly -pc_type cholesky
> parallel run (install external packages superlu_dist or mumps):
> -ksp_type preonly -pc_type lu -mat_type superlu_dist
> or
> -ksp_type preonly -pc_type cholesky -mat_type sbaijmumps
> 
> Hong
> 




__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com