[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: Shi Jin <jinzishuai@xxxxxxxxx>
- Date: Fri, 11 Apr 2008 13:56:56 -0700 (PDT)
- Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=s1024; d=yahoo.com; h=X-YMail-OSG:Received:X-Mailer:Date:From:Subject:To:MIME-Version:Content-Type:Message-ID; b=bcUZEYe/g2ttlzoIOKb3lBT/AZxSSi02WA9jeGSIqEkupRUw7OCfwLPUAcjDMpje+a4Pq+9JXGJYVZUxPjVZ6HGi9sSbHQq6A18Hz3lnjstLlcC2Kd39hNy6Lm5G0TsoYkVQaDmTPudC3r9LpT9kMQ02GMMZcOWrS8FFqZ3+zL0=;
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
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