[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