[petsc-users] Convergence of transposed linear system.
Jed Brown
jed at jedbrown.org
Wed Nov 26 10:10:51 CST 2014
Gaetan Kenway <gaetank at gmail.com> writes:
> The untransposed system converges about 6 orders of magnitude with
> GRMES(100), ASM (overlap 1) and ILU(1) with RCM reordering. The test is run
> on 128 processors. There are no convergence difficulties.
>
> However, when I try to solve the transpose of the same system, by either
> calling KSPSolveTranspose() or by assembling the transpose of the linear
> system and its preconditioner and calling KSPSolve(), GMRES stagnates after
> a negligible drop in the residual and no further progress is made.
Just a guess here, but the ASM default is "restricted ASM". Can you
compare the nontransposed and transposed convergence with each of
-pc_type <restrict,basic,interpolate>
I.e., 6 runs in total; how does each converge or not?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141126/483a24cb/attachment.pgp>
More information about the petsc-users
mailing list