[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Further question about PC with Jaocbi Row Sum
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: Further question about PC with Jaocbi Row Sum
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Fri, 11 Apr 2008 16:04:54 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=5XYKF2gchjYduZC+cxHUqOVGdBZ5ppeDirvzEFDhL+8=; b=aSxBONsXE6rhFjMsAbTXLgURP+ATT3DK7PxpwBb+IYZL7JT2oApYKHODyGZW6px2GjRForWWqhzlxTjnYLL/xE8fknJKQgEsA3h3IFNyD9M16xVlG2G/KqRO74TOTiA6feRnqIdL3oaZh1z/S704OE0d3CYf9HD7WjLe5yvm78I=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=A9axNZUlQy6sPYh5FNWVkZKnlLESWYVDQjk8ANOnXUqJXK3Cfkixkj3fUg8dwgsNrDG1YgHWih8ZfMfRdowoPcGRFcmxNWJ9XNj9g5z9rlJCMhgOmiZeaKx2UAalmV8E+yWDH7dE6jBY/IjfYx1zGWfPln5yPTf4a+klTx1Vj4U=
- In-reply-to: <548693.2032.qm@web36202.mail.mud.yahoo.com>
- References: <548693.2032.qm@web36202.mail.mud.yahoo.com>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
On Fri, Apr 11, 2008 at 3:56 PM, Shi Jin <jinzishuai@xxxxxxxxx> wrote:
> Thank you.
> Suppose I have a diagonal matrix, what is the best way to invert it in PETSc?
If you have a diagonal matrix, you just use -ksp_type preonly -pc_type jacobi
Matt
> 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
>
>
--
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