[petsc-dev] Assertion at agg.c:514 tripped

[Mark F. Adams]

I have added this code:

            if(sgid != -1 ) {
              SETERRQ(PETSC_COMM_SELF,PETSC_ERR_PLIB,"Have un-symmetric graph (apparently). Use '-pc_gamg_sym_graph true' to symetrize the graph or '-pc_gamg_threshold 0.0' if the matrix is structurally symmetric.");
            }

I'm  trying to build now.

LuLu:  I've changed the defaults so it looks like you want to add '-pc_gamg_sym_graph true'  to get back to the old semantics.

Mark

On Jul 13, 2012, at 6:20 PM, Jed Brown wrote:

> Mark, in a nonsymmetric (but structurally symmetric) bs=2 problem from Lulu, we are tripping this assertion.
> 
>           else{            /* I'm stealing this local, owned by a ghost */
>             assert(sgid==-1);
> 
> 
> What does this mean? Can't we at least always use SETERRQ for errors?
> 
> 
> Also from that section of code, I notice you are using PETSC_ERR_LIB for internal PETSc logic/corruption errors. They look like they should be PETSC_ERR_PLIB.