<div dir="ltr">On Mon, Apr 15, 2013 at 10:52 AM, Hui Zhang <span dir="ltr"><<a href="mailto:mike.hui.zhang@hotmail.com" target="_blank">mike.hui.zhang@hotmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I'm implementing a domain decomposition preconditioner.  The dof is ordered by myapp and using AO (and LocalToGlobalMapping for assembly) to map to petsc ordering.<br>

The task I'm doing is building VecScatter's from subdomains to the global domain.  So my program is<br></blockquote><div><br></div><div style>I do not understand why you would need AOs for this. They are for global reordering, whereas you seem to only</div>
<div style>need a local mapping here.</div><div style><br></div><div style>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Step 1.  I can map subdomain petsc ordering to subdomain natural ordering.<br>
<br>
Step 2.  I can also map subdomain natural ordering to global domain natural ordering.<br>
<br>
Step 3.  I have an AO for mapping global domain natural ordering to petsc ordering.<br>
<br>
Since each subdomain is defined on a sub-communicator of the communicator of the global domain.  My question is for<br>
<br>
AOApplicationToPetscIS(AO ao,IS is)<br>
<br>
can ao and is have different communicators?  Will my program be bad for large problems?  How would you do it?</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>