<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>   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. </div></div></blockquote><div><br></div><div>Note, I'm not sure why you need the null space of A^T, you want the null space of A.</div><div><br></div><div>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.</div><div><br></div><div>GAMG is use for AMR problems like this a lot in BISICLES.</div><div><br></div><div>You need to use an 'svd' coarse grid solver, or an appropriate iterative solver. LU is the default.</div><div><br></div><div>Mark</div></div></div></div>