[petsc-users] AOApplicationToPetscIS

Matthew Knepley knepley at gmail.com
Mon Apr 15 10:55:53 CDT 2013


On Mon, Apr 15, 2013 at 10:52 AM, Hui Zhang <mike.hui.zhang at hotmail.com>wrote:

> 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.
> The task I'm doing is building VecScatter's from subdomains to the global
> domain.  So my program is
>

I do not understand why you would need AOs for this. They are for global
reordering, whereas you seem to only
need a local mapping here.

   Matt


> Step 1.  I can map subdomain petsc ordering to subdomain natural ordering.
>
> Step 2.  I can also map subdomain natural ordering to global domain
> natural ordering.
>
> Step 3.  I have an AO for mapping global domain natural ordering to petsc
> ordering.
>
> Since each subdomain is defined on a sub-communicator of the communicator
> of the global domain.  My question is for
>
> AOApplicationToPetscIS(AO ao,IS is)
>
> can ao and is have different communicators?  Will my program be bad for
> large problems?  How would you do it?




-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130415/e74487c2/attachment.html>


More information about the petsc-users mailing list