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