<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">

<p></p>

<div>Dear PETSc community,</div>

<div><br>

</div>

<div>I am working on solving a Helmholtz problem with PML. The issue is that I am finding it very hard to deal with the resulting matrix system; I can get the correct solution for coarse meshes, but it takes roughly 2-4 times as long to converge for each successively

 refined mesh. I've noticed that without PML, I do not have problems with convergence speed.</div>

<div><br>

</div>

<div>I am using the GMRES solver with GAMG as the preconditioner (with block-Jacobi preconditioner for the multigrid solves). I have also tried to assemble a separate preconditioning matrix with the complex shift 1+0.5i, that does not seem to improve the results.

 Currently I am running with</div>

<div><br>

</div>

<div>   -ksp_type fgmres \</div>

<div>   -pc_type gamg \</div>

<div>   -mg_levels_pc_type bjacobi \</div>

<div>   -pc_mg_type full \</div>

<div>   -ksp_gmres_restart 150 \</div>

<div><br>

</div>

<div>Can anyone suggest some way of speeding up the convergence? Any help would be appreciated. I am attaching the output from kspview.</div>

<div><br>

</div>

<div>Best,</div>

<div><br>

</div>

<div>Artur</div>

<br>

<p></p>

</div>

</body>

</html>