Hi, I wonder whether "-snes_type test" option can work also with a nested matrix. I got the following error in my program and am not sure if there is a problem in my codes or just PETSc doesn't support it.
Testing hand-coded Jacobian, if the ratio is O(1.e-8), the hand-coded Jacobian is probably correct. Run with -snes_test_display to show difference of hand-coded and finite difference Jacobian. [0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: No support for this operation for this object type! [0]PETSC ERROR: Mat type nest! [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013 [0]PETSC ERROR: See docs/changes/index.html for recent updates. [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. [0]PETSC ERROR: See docs/index.html for manual pages. [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Libraries linked from /tools/petsc/petsc-3.4.3/linux-gnu-debug/lib [0]PETSC ERROR: Configure run at Wed Jan 15 15:48:37 2014 [0]PETSC ERROR: Configure options PETSC_ARCH=linux-gnu-debug --download-f2cblaslapack=1 -with-debugging=1 --download-superlu_dist --download-hypre=1 --download-ml=1 --download-parmetis --download-metis [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: MatSetValues() line 1085 in /tools/petsc/petsc-3.4.3/src/mat/interface/matrix.c [0]PETSC ERROR: SNESComputeJacobianDefault() line 121 in /tools/petsc/petsc-3.4.3/src/snes/interface/snesj.c [0]PETSC ERROR: SNESSolve_Test() line 64 in /tools/petsc/petsc-3.4.3/src/snes/impls/test/snestest.c [0]PETSC ERROR: SNESSolve() line 3636 in /tools/petsc/petsc-3.4.3/src/snes/interface/snes.c Thank you in advance, Nori
