<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On 22 February 2016 at 23:16, Timothée Nicolas <span dir="ltr"><<a href="mailto:timothee.nicolas@gmail.com" target="_blank">timothee.nicolas@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi all,<br><br></div>It sounds it should be obvious but I can't find it somehow. I would like to use weighted jacobi (more precisely point block jacobi) as a smoother for multigrid, but I don't find the options to set a weight omega, after I set -mg_levels_pc_type pbjacobi  </div></div></div></div></blockquote><div><br></div><div>Why not use  -ksp_type richardson -ksp_richardson_scale XXX ?<br><br></div><div>Cheers,<br></div><div>  Dave<br></div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>So I was thinking I could just artificially scale the second matrix in KSPSetOperators with a weight omega. Since it is the one used in retrieving the diagonal, it should be equivalent to have a weight and should not affect the residuals, right ? However, I notice no change of behaviour whatsoever when I do this.<br><br></div>Is it normal ? Have I grossly misunderstood something ?<br><br></div>Thanks<br><br></div>Timothee<br></div>
</blockquote></div><br></div></div>