[petsc-users] PCASM subdomains

[Matthew Knepley]

On Thu, Jun 4, 2015 at 3:47 PM, Amneet Bhalla <mail2amneet at gmail.com> wrote:

>
> Hi Folks,
>
> I have a basic question regarding matrix size as created by subdomains in
> PCASM.
>
> In the PCASMSetLocalSubdomains(PC pc,PetscInt n,IS is[],IS is_local[])
> routine, I am defining is[] to be the local+overlapping DOFs and is_local
> [] to be the nonoverlapping local DOFs only. Now
> if I grab the subksp's from PCASM and print out the matrices, would the
> matrix size correspond to size of IS in is[] or in is_local[] dofs? Put in
> another words, do the subdomain matrices are defined for just local DOFs or
> with local+overlapping DOFs?
>

The matrices should definitely be for the overlapping division. If they are
not, something is wrong. If they were
the non-overlapping distribution, we would have block Jacobi.

  Thanks,

     Matt


> Thanks,
> --
> --Amneet
>
>
>
>


-- 
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/20150604/5b967b83/attachment.html>