<div dir="ltr"><div>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.<br>
<br><br>Testing hand-coded Jacobian, if the ratio is<br>O(1.e-8), the hand-coded Jacobian is probably correct.<br>Run with -snes_test_display to show difference<br>of hand-coded and finite difference Jacobian.<br>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[0]PETSC ERROR: No support for this operation for this object type!<br>[0]PETSC ERROR: Mat type nest!<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013 <br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[0]PETSC ERROR: See docs/index.html for manual pages.<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>
[0]PETSC ERROR: Libraries linked from /tools/petsc/petsc-3.4.3/linux-gnu-debug/lib<br>[0]PETSC ERROR: Configure run at Wed Jan 15 15:48:37 2014<br>[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<br>
[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: MatSetValues() line 1085 in /tools/petsc/petsc-3.4.3/src/mat/interface/matrix.c<br>[0]PETSC ERROR: SNESComputeJacobianDefault() line 121 in /tools/petsc/petsc-3.4.3/src/snes/interface/snesj.c<br>
[0]PETSC ERROR: SNESSolve_Test() line 64 in /tools/petsc/petsc-3.4.3/src/snes/impls/test/snestest.c<br>[0]PETSC ERROR: SNESSolve() line 3636 in /tools/petsc/petsc-3.4.3/src/snes/interface/snes.c<br><br clear="all"><br></div>
Thank you in advance,<br><div>Nori<br><br></div></div>