[petsc-dev] GASM

[Barry Smith]

  This is why petsc-dev is better than a random list of email addresses.  No one on your random list uses GASM



   Barry

On Mar 27, 2012, at 4:49 PM, Mark F. Adams wrote:

> I'm trying to make general non-overlapping domains in GASM.  I create a vector of IS, each with the global indices of the equations in each block.  
> 
> I use code that looks like:
> 
>        ierr = PCSetType( subpc, PCGASM ); CHKERRQ(ierr);
>        ierr = PCGASMSetLocalSubdomains( subpc, sz, PETSC_NULL, is ); CHKERRQ(ierr);
>        ierr = PetscFree( is ); CHKERRQ(ierr);
>        ierr = PCGASMSetOverlap( subpc, 0 ); CHKERRQ(ierr);
> 
> In parallel, I get errors like this:
> 
> [1]PETSC ERROR: gis_local contains -85 indices outside of gis!
>  <snip>
> [1]PETSC ERROR: PCSetUp_GASM() line 332 in /Users/markadams/Codes/petsc-dev/src/ksp/pc/impls/gasm/gasm.c
> 
> Anyone have any ideas on what I'm doing wrong here?  I've checked the input and I think it is fine,
> 
> Mark