<div dir="ltr"><div class="gmail_extra">On Tue, May 14, 2013 at 2:35 PM, Isik Ali Ozcer <span dir="ltr"><<a href="mailto:isik.ozcer@mail.mcgill.ca" target="_blank">isik.ozcer@mail.mcgill.ca</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 All,<br>
<br>
Regarding ILU convergence at subdomain boundaries, I need to add to my colleague's message that we followed a very basic tutorial that shows matrix and rhs assembly, preconditioner and solver setup, and finally a KSPsolve.  We did not prepare any application ordering information or index sets. I am not sure if Petsc in our implementation has all the info it needs to perform proper overlapping operations with ASM. Do we need to worry about setting up AO and IS contexts before attempting KSPsolve with ASM and ILU?<br>
</blockquote><div><br></div><div>No the matrix is enough.</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">
Thanks a lot and have a great day,<br>
<br>
Isik<br>
petsc-newb<br>
<br>
<br>
<br>
-------- Original Message --------<br>
Subject:        Re: [petsc-users] Question on asm<br>
Date:   Tue, 14 May 2013 13:03:38 -0500<br>
From:   Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>><br>
To:     Konstantinos Kontzialis <<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>>,<br>
<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><br>
<br>
<br>
<br>
Konstantinos Kontzialis <<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>> writes:<br>
<br>
> Dear all,<br>
><br>
> I am using asm as a preconditioner for solving the Navier Stokes<br>
> equations. When using ILU as the subdomain preconditioner, I notice<br>
> convergence problems in the partition boundaries.<br>
<br>
I.e., it converges slowly there?  That is natural because the<br>
preconditioner is more accurate in the interior of subdomains.<br>
<br>
> This does not happen with the Jacobi preconditioner.<br>
<br>
Naturally, because Jacobi is identical regardless of the partition.<br>
<br>
> Do I need to provide any connectivity, node mapping, and/or overlap<br>
> information to Petsc before calling KSPsolve?<br>
<br>
No, the matrix is enough.<br>
<br>
</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>