[petsc-users] Max color value exeeds 65535 limit

[Smith, Barry F.]

  Your "sparse" matrix is too dense to realistically use coloring to compute the Jacobian.

  Use instead SNESComputeJacobianDefault



    Barry

> On Aug 14, 2018, at 9:06 AM, Florian Kauer <florian.kauer at koalo.de> wrote:
> 
> Hi,
> I was just testing out if my calculations also work for several 10000 variables and got the following error. Is this a fundamental limit with very large (or complex) systems of equations or does it hint to a bug/wrong configuration in my program? Everything works fine for smaller systems.
> 
> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> [0]PETSC ERROR: Argument out of range
> [0]PETSC ERROR: Max color value exeeds 65535 limit. This number is unrealistic. Perhaps a bug in code?
> Current max: 65535 user rewuested: 65536
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.8.4, unknown 
> [0]PETSC ERROR: ./csma_model on a linux-gnu-c-debug named h3a-srv-crunchie by cto7138 Tue Aug 14 15:59:57 2018
> [0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-mpich --download-fblaslapack
> [0]PETSC ERROR: #1 ISColoringCreate() line 268 in /home/cto7138/petsc/src/vec/is/is/utils/iscoloring.c
> [0]PETSC ERROR: #2 MatColoringPatch_SeqAIJ_Inode() line 2715 in /home/cto7138/petsc/src/mat/impls/aij/seq/inode.c
> [0]PETSC ERROR: #3 MatColoringPatch() line 7566 in /home/cto7138/petsc/src/mat/interface/matrix.c
> [0]PETSC ERROR: #4 MatColoringApply_SL() line 101 in /home/cto7138/petsc/src/mat/color/impls/minpack/color.c
> [0]PETSC ERROR: #5 MatColoringApply() line 360 in /home/cto7138/petsc/src/mat/color/interface/matcoloring.c
> [0]PETSC ERROR: #6 SNESComputeJacobianDefaultColor() line 81 in /home/cto7138/petsc/src/snes/interface/snesj2.c
> [0]PETSC ERROR: #7 SNESComputeJacobian() line 2360 in /home/cto7138/petsc/src/snes/interface/snes.c
> [0]PETSC ERROR: #8 SNESSolve_NEWTONLS() line 222 in /home/cto7138/petsc/src/snes/impls/ls/ls.c
> [0]PETSC ERROR: #9 SNESSolve() line 4131 in /home/cto7138/petsc/src/snes/interface/snes.c
> [0]PETSC ERROR: #10 main() line 401 in /home/cto7138/evaluation/analyticalmultihop/src/csma_model/csma_model.cpp
> [0]PETSC ERROR: PETSc Option Table entries:
> [0]PETSC ERROR: --experiment /home/cto7138/evaluation/results/analytical-csma/inverse-model/21341_95/000/experiment_CSMA.json
> [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint at mcs.anl.gov----------
> 
> Greetings,
> Florian