[petsc-users] Neumann BC with non-symmetric matrix

Mark Adams mfadams at lbl.gov
Fri Feb 26 07:12:52 CST 2016


>
>
>    4-4) Along the same lines, I tried a couple of other PCs such as
> {jacobi, sor, gamg, ilu} and none of them were able to converge with bcgs
> as the KSP. However, with gmres, almost all of them converge with the
> exception of gamg.
>

Note, I'm not sure why you need the null space of A^T, you want the null
space of A.

And for singular systems like yours you need to use a pseudo inverse of the
coarse grid because it is singular -- if you represent the null space
exactly.

GAMG is use for AMR problems like this a lot in BISICLES.

You need to use an 'svd' coarse grid solver, or an appropriate iterative
solver. LU is the default.

Mark
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160226/18d16854/attachment.html>


More information about the petsc-users mailing list