
<html>


<head>


<meta http-equiv="Content-Type" content="text/html; charset=utf-8">


</head>


<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">


Hello,


<div><br>


</div>


<div>In order to achieve reasonable performance for Helmholtz with PML, Erlangga in his paper used</div>


<div><br>


</div>


<div>1) Matrix dependent interpolation in the multigrid. The operators are nonlinear, for example an intermediate computation reads something like</div>


<div><span class="Apple-tab-span" style="white-space: pre;"></span>d = max(|a+c|, |b|, …)</div>


<div><br>


</div>


<div>2) Full weighting (This is linear, so I believe I can achieve that with <i>PCMGSetRestriction</i>).</div>


<div><br>


</div>


<div>3) F-cycle with one pre- and postsmoothing with the Jacobi iteration and relaxation factor ω = 0.5.</div>


<div><br>


</div>


<div>I am not sure how to do 1 & 3 in PETSc. Can anyone suggest a way of implementing these?</div>


<div><br>


</div>


<div>Thanks,</div>


<div><br>


</div>


<div>Artur</div>


<div><br>


</div>


<div>PS. for anyone curious, the paper is "Advances in Iterative Methods and Preconditioners for the Helmholtz Equation"</div>


</body>


</html>