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.
