Lisandro,
A followup to our previous discussion. It sounds to me like you
are "actually" solving an n+1 unknown nonlinear problem where the
"special" unknown is kept secret from SNES and managed somehow by the
application code?
You can guess how I feel about this :-). PETSc/SNES is suppo
The problem I encountered is in:
src/mat/impls/aij/seq/aijfact.c
The MatICCFactorSymbolic_SeqAIJ does not check for empty rows when the
number of levels is set to zero. For a larger number of levels it does
get checked.
I changed my code to add the diagonal into the matrix to prevent the