<div dir="ltr"><div class="gmail_extra">On Tue, May 7, 2013 at 9:20 AM, Frederik Treue <span dir="ltr"><<a href="mailto:frtr@fysik.dtu.dk" target="_blank">frtr@fysik.dtu.dk</a>></span> wrote:<br><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I may be overlooking something very obvious here, but:<br>
<br>
I'm trying to solve a convection-diffusion problem. The method I use<br>
requires me to solve some Helmholtz and Poisson equations, which is the<br>
time consuming part. In order to reduce this time, I try to use<br>
preconditioners. However, no matter which preconditioner I use (I've<br>
tried PCJACOBI,PCBJACOBI,PCPBJACOBI,PCMG with 2 levels, galerkin) it<br>
introduces errors along the edges of the local domains when using<br>
multiple processors. These errors are small, but they don't converge to<br>
0 as a function of resolution. I have checked with PCNONE, which<br>
eliminates the problem, but becomes unbearably slow.<br>
<br>
Is this somehow unavoidable? Or am I making some silly mistake? The code<br>
is somewhat complicated, but if desired, I can try to cook up a<br>
proof-of-(non)concept.</blockquote><div><br></div><div>You have made a mistake somewhere, probably in your parallel function evaluation.</div><div>For Poisson, see SNES ex5.</div><div><br></div><div>   Matt</div><div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
/Frederik Treue<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>